[sage-devel] Re: diagonal matrix constructor

2008-09-04 Thread Mike Hansen
> 3 1 3 > 100 0 99 > > There doesn't seem to be any limit on the number of arguments. That's because you didn't try hard enough :-) sage: def f(*args): : return sum(args) : sage: f(0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 2

[sage-devel] Re: diagonal matrix constructor

2008-09-04 Thread Ondrej Certik
On Thu, Sep 4, 2008 at 11:19 PM, Harald Schilly <[EMAIL PROTECTED]> wrote: > > > > On Sep 4, 10:54 pm, Jason Grout <[EMAIL PROTECTED]> wrote: >> Inhttp://trac.sagemath.org/sage_trac/ticket/3704, the diagonal matrix >> constructor is modified (and simplified/enhanced) to let someone do: >> >> sage:

[sage-devel] syntax could be more consistent

2008-09-04 Thread Jason Merrill
There are many operations in Sage that take a function/expression as the first argument, and one or more variables with ranges as subsequent arguments. The prototype is plot(x^2*sin(x),-4*pi,4*pi). It would be nice if the syntax were consistent across all such functions. There are already ticket

[sage-devel] Re: Sage 3.1.2.alpha4 released

2008-09-04 Thread mabshoff
On Sep 4, 5:16 pm, "William Stein" <[EMAIL PROTECTED]> wrote: > On Tue, Sep 2, 2008 at 6:06 AM, mabshoff <[EMAIL PROTECTED]> wrote: > > > Hello folks, > > > this is alpha4 and more progress toward the final 3.1.2. As some > > people have noticed this release is getting larger than the original >

[sage-devel] Re: Sage 3.1.2.alpha4 released

2008-09-04 Thread William Stein
On Tue, Sep 2, 2008 at 6:06 AM, mabshoff <[EMAIL PROTECTED]> wrote: > > Hello folks, > > this is alpha4 and more progress toward the final 3.1.2. As some > people have noticed this release is getting larger than the original > 3.1, but I guess such is life. THe m4ri issues have been sorted out, >

[sage-devel] Re: diagonal matrix constructor

2008-09-04 Thread Harald Schilly
On Sep 4, 11:45 pm, Jason Grout <[EMAIL PROTECTED]> wrote: > Note that there was a patch up on trac to make vector also behave this > way, so that vector(1,2,3) gives a vector of three elements. yeahr, well, i think it was somewhere else too, but i don't remember. the point is, i think, that this

[sage-devel] Re: diagonal matrix constructor

2008-09-04 Thread Jason Grout
Carl Witty wrote: > On Sep 4, 1:54 pm, Jason Grout <[EMAIL PROTECTED]> wrote: >> Inhttp://trac.sagemath.org/sage_trac/ticket/3704, the diagonal matrix >> constructor is modified (and simplified/enhanced) to let someone do: >> >> sage: diagonal_matrix(1,2,3) >> [1 0 0] >> [0 2 0] >> [0 0 3] >> >> t

[sage-devel] Re: diagonal matrix constructor

2008-09-04 Thread Carl Witty
On Sep 4, 1:54 pm, Jason Grout <[EMAIL PROTECTED]> wrote: > Inhttp://trac.sagemath.org/sage_trac/ticket/3704, the diagonal matrix > constructor is modified (and simplified/enhanced) to let someone do: > > sage: diagonal_matrix(1,2,3) > [1 0 0] > [0 2 0] > [0 0 3] > > to get a diagonal matrix with

[sage-devel] Re: diagonal matrix constructor

2008-09-04 Thread Jason Grout
Harald Schilly wrote: > > > On Sep 4, 10:54 pm, Jason Grout <[EMAIL PROTECTED]> wrote: >> Inhttp://trac.sagemath.org/sage_trac/ticket/3704, the diagonal matrix >> constructor is modified (and simplified/enhanced) to let someone do: >> >> sage: diagonal_matrix(1,2,3) > > hi, i didn't read the pa

[sage-devel] Re: diagonal matrix constructor

2008-09-04 Thread Harald Schilly
On Sep 4, 10:54 pm, Jason Grout <[EMAIL PROTECTED]> wrote: > Inhttp://trac.sagemath.org/sage_trac/ticket/3704, the diagonal matrix > constructor is modified (and simplified/enhanced) to let someone do: > > sage: diagonal_matrix(1,2,3) hi, i didn't read the patch and don't want to comment much,

[sage-devel] diagonal matrix constructor

2008-09-04 Thread Jason Grout
In http://trac.sagemath.org/sage_trac/ticket/3704, the diagonal matrix constructor is modified (and simplified/enhanced) to let someone do: sage: diagonal_matrix(1,2,3) [1 0 0] [0 2 0] [0 0 3] to get a diagonal matrix with diagonal entries 1, 2, 3. However, this conflicts with the previous sy

[sage-devel] Re: polyhedra with float coordinates

2008-09-04 Thread Marshall Hampton
OK, I think I have fixed this.A patch is up at http://trac.sagemath.org/sage_trac/ticket/4060 awaiting review. I am cross-posting this to sage-devel to hopefully attract a reviewer. In the meantime you could patch your own source if you'd like. If you aren't comfortable doing that, you coul

[sage-devel] Re: Bug in Matrix multiplication, Eigenvalue calculation?

2008-09-04 Thread chris75de
Ok, thanks. On 4 Sep., 20:53, "John Cremona" <[EMAIL PROTECTED]> wrote: > Your construction gives a rank 1 matrix, so there's a double > eigenvalue of zero.  The two eigenvectors which are found for that > eigenvalue need not be orthogonal, they are just some basis for the > eigenspace.  The entr

[sage-devel] Re: Bug in Matrix multiplication, Eigenvalue calculation?

2008-09-04 Thread John Cremona
Your construction gives a rank 1 matrix, so there's a double eigenvalue of zero. The two eigenvectors which are found for that eigenvalue need not be orthogonal, they are just some basis for the eigenspace. The entries you thought should be zero are the inner product of those eigenvectors with e

[sage-devel] Re: Bug in Matrix multiplication, Eigenvalue calculation?

2008-09-04 Thread David Joyner
On Thu, Sep 4, 2008 at 2:16 PM, chris75de <[EMAIL PROTECTED]> wrote: > > Hi, > > the following code > > b=matrix(3,1,[random() for i in [1..3]]) > B=b*b.transpose() > lb,Sb=B.right_eigenvectors() > print lb > print > print Sb.transpose()*B*Sb > print > print Sb.transpose()*Sb > > give's in my sage

[sage-devel] Bug in Matrix multiplication, Eigenvalue calculation?

2008-09-04 Thread chris75de
Hi, the following code b=matrix(3,1,[random() for i in [1..3]]) B=b*b.transpose() lb,Sb=B.right_eigenvectors() print lb print print Sb.transpose()*B*Sb print print Sb.transpose()*Sb give's in my sage version 3.0.6 the following result: [1.844328223, 4.46105473802e-17, -3.42287439238e-16] [

[sage-devel] Re: Hashing and equality of number field ideals

2008-09-04 Thread Jason Grout
John Cremona wrote: > 2008/9/4 William Stein <[EMAIL PROTECTED]>: >> On Thu, Sep 4, 2008 at 1:06 AM, John Cremona <[EMAIL PROTECTED]> wrote: >>> Perhaps a trivial contribution, but I would like to give primes ideals >>> "human-readable" names such as P2a,P2b,.. for primes above 2, etc; so >>> the

[sage-devel] Re: Hashing and equality of number field ideals

2008-09-04 Thread John Cremona
2008/9/4 William Stein <[EMAIL PROTECTED]>: > > On Thu, Sep 4, 2008 at 1:06 AM, John Cremona <[EMAIL PROTECTED]> wrote: >> >> Perhaps a trivial contribution, but I would like to give primes ideals >> "human-readable" names such as P2a,P2b,.. for primes above 2, etc; so >> then when I print a fact

[sage-devel] Re: Sciface Software GmbH (makers of MuPAD) bought by Mathworks (i.e. the MATLAB people)

2008-09-04 Thread William Stein
On Thu, Sep 4, 2008 at 3:32 AM, Jason Grout <[EMAIL PROTECTED]> wrote: > > mabshoff wrote: >> In German from heise.de: >> >> http://www.heise.de/newsticker/MuPAD-verschwindet-aber-nur-als-Produktname--/meldung/115390 >> >> the gist: >> >> * MuPAD will only be sold until September 28th 2008 >> *

[sage-devel] Re: Hashing and equality of number field ideals

2008-09-04 Thread William Stein
On Thu, Sep 4, 2008 at 1:06 AM, John Cremona <[EMAIL PROTECTED]> wrote: > > Perhaps a trivial contribution, but I would like to give primes ideals > "human-readable" names such as P2a,P2b,.. for primes above 2, etc; so > then when I print a factorizaed idea I see something like > P2a^3*P2b*P17. >

[sage-devel] Re: Sage 3.1.2.alpha4 released

2008-09-04 Thread Jason Merrill
On Sep 3, 8:35 pm, Andrzej Giniewicz <[EMAIL PROTECTED]> wrote: > Just noticed one thing, during close (by ctrl+c) of sage that run > notebook (with sage --notebook) I got unhandled error (it's 100% > reproducible for me): > > 2008-09-04 02:30:56+0200 [-] Saving notebook... > 2008-09-04 02:30:56+0

[sage-devel] Re: Sciface Software GmbH (makers of MuPAD) bought by Mathworks (i.e. the MATLAB people)

2008-09-04 Thread Harald Schilly
On Sep 4, 2:59 am, mabshoff <[EMAIL PROTECTED]> wrote: > http://www.heise.de/newsticker/MuPAD-verschwindet-aber-nur-als-Produk... > in google english: http://translate.google.com/translate?hl=en&u=http%3A%2F%2Fwww.heise.de%2Fnewsticker%2FMuPAD-verschwindet-aber-nur-als-Produktname--%2Fmeldung%2F1

[sage-devel] Re: Sciface Software GmbH (makers of MuPAD) bought by Mathworks (i.e. the MATLAB people)

2008-09-04 Thread Jason Grout
mabshoff wrote: > In German from heise.de: > > http://www.heise.de/newsticker/MuPAD-verschwindet-aber-nur-als-Produktname--/meldung/115390 > > the gist: > > * MuPAD will only be sold until September 28th 2008 > * MuPAD (Pro) licenses will remain valid > * MuPAD will be available as part of t

[sage-devel] Re: Hashing and equality of number field ideals

2008-09-04 Thread John Cremona
Perhaps a trivial contribution, but I would like to give primes ideals "human-readable" names such as P2a,P2b,.. for primes above 2, etc; so then when I print a factorizaed idea I see something like P2a^3*P2b*P17. I suppose this idea would apply to other complicated structures too. John 2008/9

[sage-devel] Re: how to make tut and ref

2008-09-04 Thread John Cremona
Thanks, I'll take a look! John 2008/9/3 Alex Ghitza <[EMAIL PROTECTED]>: > > John Cremona wrote: >> On another point, how exactly did you merge a whole sequence of >> patches into one? I tried that using a very error-prone manual way >> (create newclone, copy all files which I know have changed