[Numpy-discussion] linalg.qr versus scilab qr

2007-01-08 Thread R. David
Hello, In scilab, the QR decomposition can produce "sorted" Q and R matrices, where diagonal elements of R are decreasing. Thus, a permutation vector E is returned, so that X*E = Q*R. Is there a way to compute such a permutation QR in numpy ? It seems that numpy.linalg.qr does not sort the diago

Re: [Numpy-discussion] lapack_lite dgesv

2006-12-11 Thread R. David
Hello, > Try replacing 'int' with intc (or numpy.intc if you are not using > 'import *'). The following 'works' for me in the sense that it doesn't > throw any errors (although I imagine the results are nonsense): Thanks, it works now !! Sorry for non including the whole code, I just not wante

Re: [Numpy-discussion] lapack_lite dgesv

2006-12-11 Thread R. David
Hello Tim, > > The problem is probably your definition of ipiv. "(DIM)" is just a > parenthesized scalar, what you probably want is "(DIM,)", which is a > one-tuple. Personally, I'd recommend using list notation ("[nbrows, > nbcols]", "[DIM]") rather than tuple notation since it's both easier t

[Numpy-discussion] lapack_lite dgesv

2006-12-11 Thread R. David
a clue for this ? Regards, Romaric -- ---------- R. David - [EMAIL PROTECTED] Tel. : 03 90 24 45 48 (Fax 45 47) -- ___ Numpy-discussion mailing list Numpy-discussion@scipy.or