Re: [sage-support] identity_matrix(103).det() fails?

2010-12-28 Thread Volker Braun
The fact that this dies somewhere deep in the atlas internals probably means that your CPU doesn't support one of the SSE instructions that atlas was compiled with. I'd bet that if you compile Sage from sources on your machine then it'll work. -- To post to this group, send email to sage-suppo

Re: [sage-support] cant explain this behavior - simple aritmethics / extended euclidean algorithm

2010-12-28 Thread Johannes
Ok, from this point of view see the point. But I dont know if it's documented this way somewhere. greatz Am 28.12.2010 23:00, schrieb Iwao Kimura: > Hi Johannes, > > I'm not sure but if we think gcd(2, 0) is the generator of the ideal > generated by 2 and 0 in the rationals, we see gcd(2,0) is 1

Re: [sage-support] cant explain this behavior - simple aritmethics / extended euclidean algorithm

2010-12-28 Thread Iwao Kimura
Hi Johannes, I'm not sure but if we think gcd(2, 0) is the generator of the ideal generated by 2 and 0 in the rationals, we see gcd(2,0) is 1 because the ideal is not zero (and QQ is a field :) Best regards. On Wed, Dec 29, 2010 at 6:33 AM, Johannes wrote: > I got it solved. It's a type error:

Re: [sage-support] identity_matrix(103).det() fails?

2010-12-28 Thread Iwao Kimura
Hi Volker, Thank you for your message. I followed your suggestion and run Sage with gdb. I agree with you that linbox and blas (ATLAS) are suspected. Best regards -- | Sage Version 4.6, Release Date: 2010-10-30

Re: [sage-support] cant explain this behavior - simple aritmethics / extended euclidean algorithm

2010-12-28 Thread Johannes
I got it solved. It's a type error: the matrix sing is a matrix over QQ. If its a matrix over ZZ it works fine. but why is gcd(2,0) = 1 for 2,0 as Rationals? in my eyes its not well defined or should be the maximum of those values. greatz Am 28.12.2010 22:08, schrieb Johannes: > Hi list, > I've

[sage-support] cant explain this behavior - simple aritmethics / extended euclidean algorithm

2010-12-28 Thread Johannes
Hi list, I've a very confusing problem with some simple algorithm: following setting: a matrix sing: sage: sing [ 3 0 1] [-2 -1 -2] [ 0 1 0] sage: l = sings.column(2);l (1, -2, 0) #my algorithem, code see below: sage: extgcd(l) #some printoutput for debugging: l: (1, -2, 0) tlist: (-2, 0) #th

[sage-support] Re: Differential forms

2010-12-28 Thread jvkersch
Hi Dox, I'll try to come up with something in the next month or so. Also, a good way to help out with the development is to review the patch when it comes out. All Sage code has to go through a review process before it comes out, so this is a pretty important part of development. You can just p

[sage-support] Unimodular Matrix Transformation

2010-12-28 Thread Santanu Sarkar
I have a matrix M=(A : B), where A is a (36, 24) matrix and B (36, 9) is a matrix. I want reduce B by LLL algorithm. And also want to change A (elementary row operation) according to the changes of B. As an example during LLL if 10th row and 12th row of B are interchanged, we also have to interch

Re: [sage-support] Re: Hermite Normal Form

2010-12-28 Thread Santanu Sarkar
Thank you very much. On 28 December 2010 23:56, luisfe wrote: > On Dec 28, 6:23 pm, Santanu Sarkar > wrote: > > Size of my matrix is (90, 36) with entries are around 2^1000. What is the > > fastest > > method to compute Hermite Normal Form? > > In that case, the fastest may be the default one

[sage-support] Re: Hermite Normal Form

2010-12-28 Thread luisfe
On Dec 28, 6:23 pm, Santanu Sarkar wrote: > Size of my matrix is (90, 36) with entries are around 2^1000. What is the > fastest > method to compute Hermite Normal Form? In that case, the fastest may be the default one you are already using. Note that computing the Hermite form is fast, the hard

Re: [sage-support] Re: Hermite Normal Form

2010-12-28 Thread Santanu Sarkar
Size of my matrix is (90, 36) with entries are around 2^1000. What is the fastest method to compute Hermite Normal Form? In my matrix number of rows greater than number of columns. That is A= random_matrix(ZZ, 90, 36). Then how can I calculate transformation matrix of LLL? On 28 December 20

[sage-support] Re: Hermite Normal Form

2010-12-28 Thread luisfe
On Dec 28, 5:27 pm, Santanu Sarkar wrote: > Is there any faster method to compute Hermite Normal Form > of a matrix A and corresponding transformation matrix? I use > A.hermite_form(transformation=true). However it is very slow. > > Also is there any transformation matrix corresponding to the L

[sage-support] Hermite Normal Form

2010-12-28 Thread Santanu Sarkar
Is there any faster method to compute Hermite Normal Form of a matrix A and corresponding transformation matrix? I use A.hermite_form(transformation=true). However it is very slow. Also is there any transformation matrix corresponding to the LLL algorithm. -- To post to this group, send email

Re: [sage-support] identity_matrix(103).det() fails?

2010-12-28 Thread Volker Braun
You are not running out of memory in that computation if you have 3GB of ram. It seems like something in linbox (which is used for matrices of that size) dies during the computation. I also recently ran into some issues with linbox and blas, and was wondering if they are related. Can you run sag

[sage-support] Re: windows installation 3

2010-12-28 Thread emil
On 27 Dez., 20:37, Cyrille Piatecki wrote: > Quoting emil : > > Dear Emil, > first of all, I leave for the end of the week tomorrow so if you   > answer this message, don't be astonished if I do not react. > > Yesterday, it was impossible to download from the swiss domain so I went on >  http://b

Re: [sage-support] identity_matrix(103).det() fails?

2010-12-28 Thread Iwao Kimura
Hi Dan, Thank you for your interest. My Ubuntu box has 3GB memory. I tried several times (via Web interface, via terminal, exit Sage and invoke new Sage session...), and I always can reproduce same run-time error. I can install more 3GB memory after holiday and I'll try again. Best regards. On

[sage-support] 3D scatterplots in sage

2010-12-28 Thread bgbg
1 down vote favorite 2 Is it possible to create 3D scatterplots in sage? By scatterplot I mean graph like this: http://stackoverflow.com/questions/4439894/3d-scatterplots-in-sage -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to

Re: [sage-support] identity_matrix(103).det() fails?

2010-12-28 Thread Dan Drake
On Tue, 28 Dec 2010 at 03:41PM +0900, Iwao Kimura wrote: > On Save v4.6, (Ubuntu 10.04, 64bit), I found that > {{{ > p = 103; identity_matrix(p).det() > Traceback (most recent call last): Hrm, that works fine for me -- using Ubuntu 10.10, 64bit, and on Ubuntu 10.04, 32 bit. I increased p up to abo