Re: [sage-support] Re: how to change show() defaults for 2D plots

2010-07-12 Thread Johannes Huisman

Mike Hansen wrote:

On Sun, Jul 11, 2010 at 5:48 AM, Johannes johannes.huis...@gmail.com wrote:

Sorry for posting twice. I just noticed that
my first message got through all right, though I did
not receive it myself! Is this a common behaviour
of google groups?


Yes, this is a common behavior of Google Groups and Gmail.

You can set defaults in the plot.options dictionary.

sage: plot.options['frame'] = False


Thanks, but I did try that. It doesn't make
any difference, in any case in sage version 4.4.3
that I'm running. In fact, the sage dictionary
plot.options does not contain the keyword 'frame'.
So, setting it to False is, of course possible,
but seemingly without any effect.

Johan

--
http://pageperso.univ-brest.fr/~huisman

--
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


Re: [sage-support] Re: how to change show() defaults for 2D plots

2010-07-12 Thread William Stein
On Mon, Jul 12, 2010 at 8:57 AM, Johannes Huisman
johannes.huis...@gmail.com wrote:
 Mike Hansen wrote:

 On Sun, Jul 11, 2010 at 5:48 AM, Johannes johannes.huis...@gmail.com
 wrote:

 Sorry for posting twice. I just noticed that
 my first message got through all right, though I did
 not receive it myself! Is this a common behaviour
 of google groups?

 Yes, this is a common behavior of Google Groups and Gmail.

 You can set defaults in the plot.options dictionary.

 sage: plot.options['frame'] = False

 Thanks, but I did try that. It doesn't make
 any difference, in any case in sage version 4.4.3

The above changes the options for the plot command.
The poster is asking about changing the default options for
the show command.I don't think Sage supports this yet.

We implemented 3d plotting after 2d plotting, hence 3d plotting got some
features that 2d plotting obviously should have...

William

 that I'm running. In fact, the sage dictionary
 plot.options does not contain the keyword 'frame'.
 So, setting it to False is, of course possible,
 but seemingly without any effect.

 Johan

 --
 http://pageperso.univ-brest.fr/~huisman

 --
 To post to this group, send email to sage-support@googlegroups.com
 To unsubscribe from this group, send email to
 sage-support+unsubscr...@googlegroups.com
 For more options, visit this group at
 http://groups.google.com/group/sage-support
 URL: http://www.sagemath.org




-- 
William Stein
Professor of Mathematics
University of Washington
http://wstein.org

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] Re: Pattern matching of a symbolic function acting on a symbolic variable

2010-07-12 Thread David Sanders


On Jul 11, 12:22 pm, William Stein wst...@gmail.com wrote:
 On Sat, Jul 10, 2010 at 10:04 PM, David Sanders dpsand...@gmail.com wrote:

  On Jul 10, 10:01 pm, David Sanders dpsand...@gmail.com wrote:
  Hi,

  Following up from a couple of my previous posts, I am now wondering
  how to do symbolic pattern matching for an expression of the following
  form:

  Apologies for replying to my own post, but I was wondering if there's
  any documentation I should be reading about this stuff,
  rather than posting every naive question here?
  (Although maybe it's actually useful to post naive questions here!)

 Maybe

    http://sagemath.org/doc/reference/calculus.html

 I'm curious if you had any trouble finding this.  It's very easy --
 just click Documentation, then Reference from the sage homepage.

Thanks, I had seen this, I think via a Google search,
but I did not explore it enough. I guess I was confused by the title
'calculus' instead of 'symbolic manipulation' -- I assumed it was just
about differentiation etc.

By the way, a minor but important point: the typographical conventions
on this page (and many others) are not consistent: some of the titles
have every word with an initial capital, whereas some have only the
first word and proper names capitalized (my preference). Some have
full stops (points) at the end, and others don't (my preference).

Personally I have never liked monolithic pages with all possible
information, such as
http://sagemath.org/doc/reference/sage/symbolic/expression.html
It is very difficult to navigate and find what you are looking for.
Perhaps the examples could somehow be hidden until you hit a relevant
link, at which point they reveal themselves.

I suppose the page is probably automatically generated, but it would
be more helpful to split it up into bitesize pieces.
For example, there is a lot of noise generated by obvious
functions like sin, cos etc.
Finding the relevant bit about how to do pattern matching is not easy,
for example. (Especially if you don't know that what you need is
pattern matching!)

In general, I find that the documentation is difficult to navigate for
a Sage beginner. It is not clear what level each document is at.
Even the order of the documents in the Sage standard documentation
hides A tour of Sage, for example, amid more advanced documents.
Of course, I realise that writing documentation is difficult!

Perhaps it would be useful to have a kind of beginner's page, which
lists a suggested order in which to read the different types of
documentation.

David.




 William

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] Re: Pattern matching of a symbolic function acting on a symbolic variable

2010-07-12 Thread David Sanders


On Jul 11, 8:13 pm, ma...@mendelu.cz ma...@mendelu.cz wrote:
 On 11 čnc, 12:22, William Stein wst...@gmail.com wrote:

  On Sat, Jul 10, 2010 at 10:04 PM, David Sanders dpsand...@gmail.com wrote:

  Maybe

     http://sagemath.org/doc/reference/calculus.html

 Perhaps 
 alsohttp://www.ginac.de/tutorial/Pattern-matching-and-advanced-substituti...


OK, that is also useful, thanks.
In particular, I notice that there is the concept of indexed objects
in GiNaC.
Is this accessible from Sage?


 Or pass to Maxima and use pattern matching from Maxima, which is well
 documented in documantation to Maxima.

I am now confused about which system is used for the symbolics in
Sage? Is it GiNaC, or Maxima, or a mixture?  How can I find out which
system is being used for which operation?
I would very much prefer not to have to learn Maxima if I can help it,
since the whole point is that Sage is supposed to provide the nice,
coherent interface which makes this unnecessary!

I also note that after reading the documentation, I am still left
without an answer to my original question, which is how to do pattern
matching in Sage (or if it's even possible) for something of the
form
f(i) !

David.


 Robert

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] How to plot regions defined by inequalities in 3D?

2010-07-12 Thread David Sanders
Hi,

I have been playing around with the implicit_plot3d command, and it's
very nice.

Is there something similar to plot regions defined by inequalities in
3D, along the lines of the Mathematica
RegionPlot3D command?  I see that there is an old discussion from 2
years ago about this.

It seems to me (in my ignorance) that the kind of algorithm required
to do this should not be too different from that for implicit_plot3d
with the region option. (Though implicit_plot3d draws surfaces,
whereas region_plot3d would draw volumes. But apparently the marching
cubes algorithm is used for both?)

Thanks and best wishes,
David.

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] Re: Notebooks forgets assignments

2010-07-12 Thread Joaquim Puig
Hi,

Thanks for the answers and sorry for the late reply. The parameters
for running the notebook server are:

/usr/local/sage-4.3.2/sage -python notebook.py

on notebook.py:

from sage.all import *

notebook(directory='/home/nbuser/nbfiles', port=8000, accounts=True,
address='', ulimit='-u 100 -v 300 -t 3600', open_viewer=False,
timeout=120, secure=True, server_pool=['nbu...@localhost'])

I guess that the problem is with the timeout parameter. Which would be
a reasonable value?

Thanks again,

Joaquim Puig


On Jul 10, 7:07 am, William Stein wst...@gmail.com wrote:
 On Sat, Jul 10, 2010 at 2:48 AM, Mike Hansen mhan...@gmail.com wrote:
  2010/7/9 Jorge E. ´Sanchez Sanchez hnr...@hotmail.com:
  Which could be the problem? Thanks in advance,

  It sounds like the issue could be that the notebook closes the
  worksheet session after that amount of inactivity.  Then, when you
  execute the next command it starts up a new session in which a
  hasn't been defined.

 That's what it sounds like.  How *exactly* do you run the notebook server?

 William



  --Mike

  --
  To post to this group, send email to sage-support@googlegroups.com
  To unsubscribe from this group, send email to 
  sage-support+unsubscr...@googlegroups.com
  For more options, visit this group 
  athttp://groups.google.com/group/sage-support
  URL:http://www.sagemath.org

 --
 William Stein
 Professor of Mathematics
 University of Washingtonhttp://wstein.org

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


Re: [sage-support] Re: Notebooks forgets assignments

2010-07-12 Thread William Stein
On Mon, Jul 12, 2010 at 11:29 AM, Joaquim Puig joaquim.p...@gmail.com wrote:
 Hi,

 Thanks for the answers and sorry for the late reply. The parameters
 for running the notebook server are:

 /usr/local/sage-4.3.2/sage -python notebook.py

 on notebook.py:

 from sage.all import *

 notebook(directory='/home/nbuser/nbfiles', port=8000, accounts=True,
 address='', ulimit='-u 100 -v 300 -t 3600', open_viewer=False,
 timeout=120, secure=True, server_pool=['nbu...@localhost'])

 I guess that the problem is with the timeout parameter. Which would be
 a reasonable value?

You set it to *2 minutes*.  That's probably not so reasonable.   Maybe
timeout = 3600?


 Thanks again,

 Joaquim Puig


 On Jul 10, 7:07 am, William Stein wst...@gmail.com wrote:
 On Sat, Jul 10, 2010 at 2:48 AM, Mike Hansen mhan...@gmail.com wrote:
  2010/7/9 Jorge E. ´Sanchez Sanchez hnr...@hotmail.com:
  Which could be the problem? Thanks in advance,

  It sounds like the issue could be that the notebook closes the
  worksheet session after that amount of inactivity.  Then, when you
  execute the next command it starts up a new session in which a
  hasn't been defined.

 That's what it sounds like.  How *exactly* do you run the notebook server?

 William



  --Mike

  --
  To post to this group, send email to sage-support@googlegroups.com
  To unsubscribe from this group, send email to 
  sage-support+unsubscr...@googlegroups.com
  For more options, visit this group 
  athttp://groups.google.com/group/sage-support
  URL:http://www.sagemath.org

 --
 William Stein
 Professor of Mathematics
 University of Washingtonhttp://wstein.org

 --
 To post to this group, send email to sage-support@googlegroups.com
 To unsubscribe from this group, send email to 
 sage-support+unsubscr...@googlegroups.com
 For more options, visit this group at 
 http://groups.google.com/group/sage-support
 URL: http://www.sagemath.org




-- 
William Stein
Professor of Mathematics
University of Washington
http://wstein.org

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


Re: [sage-support] Re: Pattern matching of a symbolic function acting on a symbolic variable

2010-07-12 Thread Burcin Erocal
Hi David,

On Mon, 12 Jul 2010 02:24:41 -0700 (PDT)
David Sanders dpsand...@gmail.com wrote:

 On Jul 11, 8:13 pm, ma...@mendelu.cz ma...@mendelu.cz wrote:
  On 11 čnc, 12:22, William Stein wst...@gmail.com wrote:
 
   On Sat, Jul 10, 2010 at 10:04 PM, David Sanders
   dpsand...@gmail.com wrote:
 
   Maybe
 
      http://sagemath.org/doc/reference/calculus.html
 
  Perhaps
  alsohttp://www.ginac.de/tutorial/Pattern-matching-and-advanced-substituti...
 
 
 OK, that is also useful, thanks.
 In particular, I notice that there is the concept of indexed objects
 in GiNaC.
 Is this accessible from Sage?

There is an experimental patch. I'm really busy these days, but this is
close to the top of my list. :)

http://groups.google.com/group/sage-devel/t/69ab50fe11672111

  Or pass to Maxima and use pattern matching from Maxima, which is
  well documented in documantation to Maxima.
 
 I am now confused about which system is used for the symbolics in
 Sage? Is it GiNaC, or Maxima, or a mixture?  How can I find out which
 system is being used for which operation?

It's a mixture, though we are trying to move as much as possible to
GiNaC and native Sage. AFAIK, the only way to tell what is being used
is to read the code. As a general rule, basic arithmetic and pattern
matching is done with GiNaC, more advanced functionality, limits,
simplification, factorization, etc. calls maxima.

 I would very much prefer not to have to learn Maxima if I can help it,
 since the whole point is that Sage is supposed to provide the nice,
 coherent interface which makes this unnecessary!

I can totally understand that. We definitely need to improve the
interface to cover this functionality. Thanks for pointing it out.

 I also note that after reading the documentation, I am still left
 without an answer to my original question, which is how to do pattern
 matching in Sage (or if it's even possible) for something of the
 form
 f(i) !

I don't think this is supported by GiNaC expressions at the moment. If
using wildcards for functions is available in GiNaC, we can wrap it
easily. Implementing it would take more time though. Can you ask the
GiNaC list if this is possible purely using GiNaC (from C++)?


Cheers,
Burcin

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] Re: Notebooks forgets assignments

2010-07-12 Thread Joaquim Puig
Thanks! I will try it this way.

Joaquim Puig

On Jul 12, 11:42 am, William Stein wst...@gmail.com wrote:
 On Mon, Jul 12, 2010 at 11:29 AM, Joaquim Puig joaquim.p...@gmail.com wrote:
  Hi,

  Thanks for the answers and sorry for the late reply. The parameters
  for running the notebook server are:

  /usr/local/sage-4.3.2/sage -python notebook.py

  on notebook.py:

  from sage.all import *

  notebook(directory='/home/nbuser/nbfiles', port=8000, accounts=True,
  address='', ulimit='-u 100 -v 300 -t 3600', open_viewer=False,
  timeout=120, secure=True, server_pool=['nbu...@localhost'])

  I guess that the problem is with the timeout parameter. Which would be
  a reasonable value?

 You set it to *2 minutes*.  That's probably not so reasonable.   Maybe
 timeout = 3600?





  Thanks again,

  Joaquim Puig

  On Jul 10, 7:07 am, William Stein wst...@gmail.com wrote:
  On Sat, Jul 10, 2010 at 2:48 AM, Mike Hansen mhan...@gmail.com wrote:
   2010/7/9 Jorge E. ´Sanchez Sanchez hnr...@hotmail.com:
   Which could be the problem? Thanks in advance,

   It sounds like the issue could be that the notebook closes the
   worksheet session after that amount of inactivity.  Then, when you
   execute the next command it starts up a new session in which a
   hasn't been defined.

  That's what it sounds like.  How *exactly* do you run the notebook server?

  William

   --Mike

   --
   To post to this group, send email to sage-support@googlegroups.com
   To unsubscribe from this group, send email to 
   sage-support+unsubscr...@googlegroups.com
   For more options, visit this group 
   athttp://groups.google.com/group/sage-support
   URL:http://www.sagemath.org

  --
  William Stein
  Professor of Mathematics
  University of Washingtonhttp://wstein.org

  --
  To post to this group, send email to sage-support@googlegroups.com
  To unsubscribe from this group, send email to 
  sage-support+unsubscr...@googlegroups.com
  For more options, visit this group 
  athttp://groups.google.com/group/sage-support
  URL:http://www.sagemath.org

 --
 William Stein
 Professor of Mathematics
 University of Washingtonhttp://wstein.org

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] Re: Pattern matching of a symbolic function acting on a symbolic variable

2010-07-12 Thread David Sanders


On Jul 12, 12:19 pm, Burcin Erocal bur...@erocal.org wrote:
 Hi David,

 On Mon, 12 Jul 2010 02:24:41 -0700 (PDT)



 David Sanders dpsand...@gmail.com wrote:
  On Jul 11, 8:13 pm, ma...@mendelu.cz ma...@mendelu.cz wrote:
   On 11 čnc, 12:22, William Stein wst...@gmail.com wrote:

On Sat, Jul 10, 2010 at 10:04 PM, David Sanders
dpsand...@gmail.com wrote:

Maybe

   http://sagemath.org/doc/reference/calculus.html

   Perhaps
   alsohttp://www.ginac.de/tutorial/Pattern-matching-and-advanced-substituti...

  OK, that is also useful, thanks.
  In particular, I notice that there is the concept of indexed objects
  in GiNaC.
  Is this accessible from Sage?

 There is an experimental patch. I'm really busy these days, but this is
 close to the top of my list. :)


OK great, thanks!

 http://groups.google.com/group/sage-devel/t/69ab50fe11672111

   Or pass to Maxima and use pattern matching from Maxima, which is
   well documented in documantation to Maxima.

  I am now confused about which system is used for the symbolics in
  Sage? Is it GiNaC, or Maxima, or a mixture?  How can I find out which
  system is being used for which operation?

 It's a mixture, though we are trying to move as much as possible to
 GiNaC and native Sage. AFAIK, the only way to tell what is being used
 is to read the code. As a general rule, basic arithmetic and pattern
 matching is done with GiNaC, more advanced functionality, limits,
 simplification, factorization, etc. calls maxima.


I see, thanks for the clarification. The structure is gradly becoming
clearer to me.
I have been looking at GiNaC a bit, and it seems to be very clean.

  I would very much prefer not to have to learn Maxima if I can help it,
  since the whole point is that Sage is supposed to provide the nice,
  coherent interface which makes this unnecessary!

 I can totally understand that. We definitely need to improve the
 interface to cover this functionality. Thanks for pointing it out.

  I also note that after reading the documentation, I am still left
  without an answer to my original question, which is how to do pattern
  matching in Sage (or if it's even possible) for something of the
  form
  f(i) !

 I don't think this is supported by GiNaC expressions at the moment. If
 using wildcards for functions is available in GiNaC, we can wrap it
 easily. Implementing it would take more time though. Can you ask the
 GiNaC list if this is possible purely using GiNaC (from C++)?

Actually, from browsing the documentation, I can't even find symbolic
functions in GiNaC like
f = function('f')
in Sage.
Does this exist in GiNaC?  Does the above Sage statement reference
something in GiNaC? If so, what?

(I tried to answer this question myself using
function?
function??

in Sage, but this just said that it was a built-in function, and gave
me neither the name of a file, nor the source code.
Is there another way  of getting information about this?)

Thanks,
David.


 Cheers,
 Burcin

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] Non-Linear Programming in Sage

2010-07-12 Thread Alexander Strachan
After reading through the documentation, it doesn't look like sage can
solve a non-linear program.  Is this correct?

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] Re: How to plot regions defined by inequalities in 3D?

2010-07-12 Thread Carl Witty
On Jul 12, 2:28 am, David Sanders dpsand...@gmail.com wrote:
 Hi,

 I have been playing around with the implicit_plot3d command, and it's
 very nice.

 Is there something similar to plot regions defined by inequalities in
 3D, along the lines of the Mathematica
 RegionPlot3D command?  I see that there is an old discussion from 2
 years ago about this.

 It seems to me (in my ignorance) that the kind of algorithm required
 to do this should not be too different from that for implicit_plot3d
 with the region option. (Though implicit_plot3d draws surfaces,
 whereas region_plot3d would draw volumes. But apparently the marching
 cubes algorithm is used for both?)

I haven't looked at this stuff in more than a year, but I think this
is all accurate:

Our plotting framework doesn't really understand volumes, only
surfaces.  So it would be a major overhaul to produce a plot that
showed (via some sort of volumetric shading, say) the difference
between the inside and the outside of your region.

However, if you want to produce a plot of the surface of your region,
that's pretty easy.  If your region is defined by a single inequality
F(x,y,z)0, then you can just implicit_plot3d F(x,y,z).  If your
region is defined as a boolean combination of inequalities, then
arrange all the inequalities to be of the form F(x,y,z)  0, then drop
all the  0, replace and with max_symbolic, replace or with
min_symbolic, and replace not F(x,y,z) with -F(x,y,z).  Also,
when you plot, because of http://trac.sagemath.org/sage_trac/ticket/9483
you need to add smooth=False.

Here's a complete example.  This forms the intersection between a cube
and the union of two cylinders.

sage: var('x,y,z')
(x, y, z)
sage: implicit_plot3d(max_symbolic(min_symbolic(x*x+y*y-1, x*x+z*z-2),
x-1.8, y-1.8, z-1.8, -x-1.8, -y-1.8, -z-1.8), (x, -2, 2), (y, -2, 2),
(z, -2, 2), smooth=False)

It would be great to put all of this into a region_plot3d command, but
as far as I know, Sage does not yet support symbolic conjunctions and
disjunctions (ands and ors); so it would be difficult to tell
region_plot3d about any region more complicated than a single
inequality.

Carl

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


RE: [sage-support] Non-Linear Programming in Sage

2010-07-12 Thread Jorge E . ´Sanchez Sanchez

Sage comes with the library cvxopt preinstalled, so you can solve non-linear 
problems.

Regards
Jorge

 Date: Mon, 12 Jul 2010 10:09:40 -0700
 Subject: [sage-support] Non-Linear Programming in Sage
 From: alexander.r.strac...@gmail.com
 To: sage-support@googlegroups.com
 
 After reading through the documentation, it doesn't look like sage can
 solve a non-linear program.  Is this correct?
 
 -- 
 To post to this group, send email to sage-support@googlegroups.com
 To unsubscribe from this group, send email to 
 sage-support+unsubscr...@googlegroups.com
 For more options, visit this group at 
 http://groups.google.com/group/sage-support
 URL: http://www.sagemath.org
  
_
¿Qué signigica Messenger para ti?
www.vivirmessenger.com

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] Re: Non-Linear Programming in Sage

2010-07-12 Thread Harald Schilly


On Jul 12, 7:09 pm, Alexander Strachan
alexander.r.strac...@gmail.com wrote:
 After reading through the documentation, it doesn't look like sage can
 solve a non-linear program.  Is this correct?

There is minimize and minimize_constraint wrapping different solvers.
http://www.sagemath.org/doc/reference/sage/numerical/optimize.html

H

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] Re: How to plot regions defined by inequalities in 3D?

2010-07-12 Thread David Sanders


On Jul 12, 7:29 pm, Carl Witty carl.wi...@gmail.com wrote:
 On Jul 12, 2:28 am, David Sanders dpsand...@gmail.com wrote:

  Hi,

  I have been playing around with the implicit_plot3d command, and it's
  very nice.

  Is there something similar to plot regions defined by inequalities in
  3D, along the lines of the Mathematica
  RegionPlot3D command?  I see that there is an old discussion from 2
  years ago about this.

  It seems to me (in my ignorance) that the kind of algorithm required
  to do this should not be too different from that for implicit_plot3d
  with the region option. (Though implicit_plot3d draws surfaces,
  whereas region_plot3d would draw volumes. But apparently the marching
  cubes algorithm is used for both?)

 I haven't looked at this stuff in more than a year, but I think this
 is all accurate:

 Our plotting framework doesn't really understand volumes, only
 surfaces.  So it would be a major overhaul to produce a plot that
 showed (via some sort of volumetric shading, say) the difference
 between the inside and the outside of your region.

Yes, I agree, although that overhaul could well be worth it!


 However, if you want to produce a plot of the surface of your region,
 that's pretty easy.  If your region is defined by a single inequality
 F(x,y,z)0, then you can just implicit_plot3d F(x,y,z).  If your
 region is defined as a boolean combination of inequalities, then
 arrange all the inequalities to be of the form F(x,y,z)  0, then drop
 all the  0, replace and with max_symbolic, replace or with
 min_symbolic, and replace not F(x,y,z) with -F(x,y,z).  Also,
 when you plot, because ofhttp://trac.sagemath.org/sage_trac/ticket/9483
 you need to add smooth=False.

 Here's a complete example.  This forms the intersection between a cube
 and the union of two cylinders.

 sage: var('x,y,z')
 (x, y, z)
 sage: implicit_plot3d(max_symbolic(min_symbolic(x*x+y*y-1, x*x+z*z-2),
 x-1.8, y-1.8, z-1.8, -x-1.8, -y-1.8, -z-1.8), (x, -2, 2), (y, -2, 2),
 (z, -2, 2), smooth=False)

 It would be great to put all of this into a region_plot3d command, but
 as far as I know, Sage does not yet support symbolic conjunctions and
 disjunctions (ands and ors); so it would be difficult to tell
 region_plot3d about any region more complicated than a single
 inequality.

OK, thanks, I had come to a similar conclusion, though the trick with
max_ and min_symbolic is neat.
For the moment I think this will do most of the things that I need,
but the fact that there are no symbolic and's and or's
makes it quite (=very) messy if there are multiple intersecting
volumes which restrict each other.

Thanks,
David.


 Carl

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] How to use R's untar in Sage?

2010-07-12 Thread rickhg12hs
I've tried to use R's untar command a number of different ways from
Sage's commandline interface and I am having trouble listing the
contents of a LZMA compressed TAR file. Actually, I don't seem to get
any desirable output from r.untar so I am thinking I'm not even close
to using the correct syntax.

If I have a LZMA compressed TAR file, how can I list/read it with Sage
and/or r.untar in Sage?

Cheers!

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] Re: How to use R's untar in Sage?

2010-07-12 Thread kcrisman


On Jul 12, 10:57 pm, rickhg12hs rickhg1...@gmail.com wrote:
 I've tried to use R's untar command a number of different ways from
 Sage's commandline interface and I am having trouble listing the
 contents of a LZMA compressed TAR file. Actually, I don't seem to get
 any desirable output from r.untar so I am thinking I'm not even close
 to using the correct syntax.

Syntax with R commands is not always obvious for optional keywords.
Can you give us the *exact* R command (with options) you are trying to
use?  Note the documentation:  What options are supported will depend
on the ‘tar’ used.   So it is possible that the R untar in question
on Sage behaves differently from a 'native' implementation.  But most
likely it is the use of options that is the problem.  Unfortunately I
can't remember offhand how I've used them in the past :( but perhaps
someone else will.

 If I have a LZMA compressed TAR file, how can I list/read it with Sage
 and/or r.untar in Sage?

If you are just trying to use R from within Sage, without integrating
it with the rest of Sage, you can also use r_console(), or use %r in
the notebook, and just do things that way.

- kcrisman

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org