[sage-support] Re: adding more points to an existing polyhedron

2010-07-13 Thread tvn
On Jul 13, 2:16 pm, TianWei wrote: > > Hi,  if I have an existing polyhedron p and want to add a vertex v > > (that was not in p) to p, is it possible to modify it (e.g., adding > > new vertices v to p) instead of creating a new polyhedron p'  that > > contains all the previous vertices of p and

[sage-support] Re: adding more points to an existing polyhedron

2010-07-13 Thread TianWei
> Hi, if I have an existing polyhedron p and want to add a vertex v > (that was not in p) to p, is it possible to modify it (e.g., adding > new vertices v to p) instead of creating a new polyhedron p' that > contains all the previous vertices of p and v ? It might be worth pointing out the doc

Re: [sage-support] An argument for more direct axiom integration

2010-07-13 Thread William Stein
On Tue, Jul 13, 2010 at 7:24 PM, Dr. David Kirkby wrote: > On 07/13/10 08:05 AM, William Stein wrote: >> >> On Tue, Jul 13, 2010 at 7:46 AM, Robert Bradshaw >>  wrote: >>> >>> On Mon, Jul 5, 2010 at 4:39 AM, Dr. David Kirkby >>>  wrote: On 07/ 5/10 06:18 AM, William Stein wrote: >

[sage-support] Re: Superimpose various region_plot's

2010-07-13 Thread David Sanders
On Jul 13, 4:47 pm, Jason Grout wrote: > On 7/13/10 4:46 AM, David Sanders wrote: > > > > > Hi, > > > I need to superimpose several region_plot's. These have regions > > colored in different colors which may or may not overlap. However, if > > I do something like > > > var('p q') > > plot1 = reg

Re: [sage-support] An argument for more direct axiom integration

2010-07-13 Thread Dr. David Kirkby
On 07/13/10 08:05 AM, William Stein wrote: On Tue, Jul 13, 2010 at 7:46 AM, Robert Bradshaw wrote: On Mon, Jul 5, 2010 at 4:39 AM, Dr. David Kirkby wrote: On 07/ 5/10 06:18 AM, William Stein wrote: Great idea - you could add an algorithm="axiom" option to sage's integrate command. Perso

[sage-support] adding more points to an existing polyhedron

2010-07-13 Thread tvn
Hi, if I have an existing polyhedron p and want to add a vertex v (that was not in p) to p, is it possible to modify it (e.g., adding new vertices v to p) instead of creating a new polyhedron p' that contains all the previous vertices of p and v ? Thanks, -- To post to this group, send email

Re: [sage-support] Re: sage.numerical.optimize.minimize() returns without returning a value

2010-07-13 Thread Bouazizi Faiçal
cvvcv 2010/7/13 8fjm39j > Following up with another data point: > > Sage 4.3.4 under some version of Redhat Linux on an intel computer > with 64GB RAM. > > The case that failed on the 4GB macbookpro succeeds here: minimize() > returns a value. > > A larger case fails in the same manner, minimize

[sage-support] Re: sage.numerical.optimize.minimize() returns without returning a value

2010-07-13 Thread Harald Schilly
On 13 Jul., 17:32, 8fjm39j wrote: > One minor question:  in your definition of the gradient function, >     sage: gradfun = lambda x:np.array(map(lambda f:f(*x), > eq.gradient())) > should that be 'eq.gradient()' or should it be 'req.gradient()' ? > (recall that req = eq.change_ring(RDF)) I only

[sage-support] Re: sage.numerical.optimize.minimize() returns without returning a value

2010-07-13 Thread 8fjm39j
Thanks for your suggestion. I've now running a worksheet using your method to minimize a large polynomial. So far, minimize() has not yet returned. One minor question: in your definition of the gradient function, sage: gradfun = lambda x:np.array(map(lambda f:f(*x), eq.gradient())) should

[sage-support] Re: Superimpose various region_plot's

2010-07-13 Thread Jason Grout
On 7/13/10 4:46 AM, David Sanders wrote: Hi, I need to superimpose several region_plot's. These have regions colored in different colors which may or may not overlap. However, if I do something like var('p q') plot1 = region_plot([p+q<1, p+q>-1], (p,-2,2), (q,-2,2), incol='red') plot2 = region_

[sage-support] Re: Using ipython (magic) shell commands in notebook interface

2010-07-13 Thread orca
Hi Tim Thanks for the prompt reply! I have tried and, of course, confirmed your suggestions. However, I am a little concerned about the best way to streamline my experience with Sage. In fact, I find the embedding of the figures (all kinds of plots) and even Latex in the notebook quite amazing an

[sage-support] Re: sage.numerical.optimize.minimize() returns without returning a value

2010-07-13 Thread 8fjm39j
Following up with another data point: Sage 4.3.4 under some version of Redhat Linux on an intel computer with 64GB RAM. The case that failed on the 4GB macbookpro succeeds here: minimize() returns a value. A larger case fails in the same manner, minimize() returning without a value: S is a

[sage-support] Re: sage.numerical.optimize.minimize() returns without returning a value

2010-07-13 Thread Harald Schilly
On Jul 13, 3:07 pm, 8fjm39j wrote: > Any help would be much appreciated. I'm not sure if problems of this size work. Also, you should add the gradient to the minimize method. Here is a snippet that might help you. You do not need the SR as far as i can see. sage: RQ = PolynomialRing(QQ, 30, '

[sage-support] sage.numerical.optimize.minimize() returns without returning a value

2010-07-13 Thread 8fjm39j
I am trying to minimize numerically (over the reals) a polynomial S with rational coefficients. The function sage.numerical.optimize.minimize() returns without returning a value. Here's the relevant fragment of Sage code (executed using the menu item 'evaluate all') = BEGIN CODE FRAGMENT =

[sage-support] Superimpose various region_plot's

2010-07-13 Thread David Sanders
Hi, I need to superimpose several region_plot's. These have regions colored in different colors which may or may not overlap. However, if I do something like var('p q') plot1 = region_plot([p+q<1, p+q>-1], (p,-2,2), (q,-2,2), incol='red') plot2 = region_plot([p-q<1, p-q>-1], (p,-2,2), (q,-2,2), i

Re: [sage-support] Using ipython (magic) shell commands in notebook interface

2010-07-13 Thread Tim Joseph Dumol
Hi, You can use: {{{ %sh cd ~/foobar ls echo "Blah" }}}, but any directory changes are limited to the cell only (due to architectural reasons). Otherwise, os.listdir('.') and and os.chdir('~/foo') are your best choices. On Tue, Jul 13, 2010 at 5:53 PM, orca wrote: > Hi there, > > I'm a newbie

[sage-support] Using ipython (magic) shell commands in notebook interface

2010-07-13 Thread orca
Hi there, I'm a newbie Sage user and my question is related to using its notebook interface in a more convenient ipython-like way in my single machine (be it a laptop or desktop), if possible. What I mean, explicitly, for example, is being able to type: cd ~orca/research and being able to chang

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

2010-07-13 Thread Burcin Erocal
On Mon, 12 Jul 2010 08:46:16 -0700 (PDT) David Sanders wrote: > > > 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 > > >

Re: [sage-support] An argument for more direct axiom integration

2010-07-13 Thread William Stein
On Tue, Jul 13, 2010 at 7:46 AM, Robert Bradshaw wrote: > On Mon, Jul 5, 2010 at 4:39 AM, Dr. David Kirkby > wrote: >> On 07/ 5/10 06:18 AM, William Stein wrote: >>> >>> Great idea - you could add an algorithm="axiom" option to sage's >>> integrate command. >> >> Personally, and I am going to dar

Re: [sage-support] Re: Which python for sage?

2010-07-13 Thread William Stein
On Tue, Jul 13, 2010 at 8:11 AM, Robert Bradshaw wrote: > On Sun, Jul 4, 2010 at 9:47 AM, Greg Kuperberg > wrote: >> On Jul 3, 8:33 am, William Stein wrote: >>> Sage can't switch to Python 3 until every single Python package in >>> Sage is ported to Python 3. >>> This is far from done.  It's po

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

2010-07-13 Thread rickhg12hs
Found the problem was with my system's tar command. After specifying tar="internal" (meaning that it uses R's builtin tar) everything worked fine. Thanks for your consideration though. On Jul 12, 11:14 pm, kcrisman wrote: > On Jul 12, 10:57 pm, rickhg12hs wrote: > > > I've tried to use R's unt