Re: An error of matrix inversion using NumPy

2007-04-04 Thread Travis Oliphant
lancered wrote: > > >So, can you tell me what goes wrong? Is this a bug in > Numpy.linalg? How to deal with this situation? If you need, I can > post the matrix I used below, but it is so long,so not at the moment. > As you discovered, it is very likely your problem is a very high co

Re: An error of matrix inversion using NumPy

2007-04-04 Thread Lou Pecora
In article <[EMAIL PROTECTED]>, "lancered" <[EMAIL PROTECTED]> wrote: > Here is the eigenvalues of KK I obtained: > > >>> linalg.eigvals(KK) > array([ 1.11748411e+05, 3.67154458e+04, 3.41580846e+04, > 2.75272440e+04, 2.09790868e+04, 1.86242332e+04, > 8.68628325e+03,

Re: An error of matrix inversion using NumPy

2007-04-04 Thread Robin Becker
lancered wrote: > Here is the eigenvalues of KK I obtained: > > >>> linalg.eigvals(KK) > array([ 1.11748411e+05, 3.67154458e+04, 3.41580846e+04, > 2.75272440e+04, 2.09790868e+04, 1.86242332e+04, > 8.68628325e+03, 6.66127732e+03, 6.15547187e+03, > 4.6862619

Re: An error of matrix inversion using NumPy

2007-04-04 Thread lancered
Here is the eigenvalues of KK I obtained: >>> linalg.eigvals(KK) array([ 1.11748411e+05, 3.67154458e+04, 3.41580846e+04, 2.75272440e+04, 2.09790868e+04, 1.86242332e+04, 8.68628325e+03, 6.66127732e+03, 6.15547187e+03, 4.68626197e+03, 3.17838339e+03, 2.8

Re: An error of matrix inversion using NumPy

2007-04-04 Thread Robin Becker
Robin Becker wrote: > lancered wrote: h. If > your matrix is symmetric then you should certainly be using a qr decomposition I meant to say :) -- Robin Becker -- http://mail.python.org/mailman/listinfo/python-list

Re: An error of matrix inversion using NumPy

2007-04-04 Thread Robin Becker
lancered wrote: > Hi dear all, .. > matrices are correct. > >So, can you tell me what goes wrong? Is this a bug in > Numpy.linalg? How to deal with this situation? If you need, I can > post the matrix I used below, but it is so long,so not at the moment. ... presumably the

Re: An error of matrix inversion using NumPy

2007-04-04 Thread Lou Pecora
In article <[EMAIL PROTECTED]>, "lancered" <[EMAIL PROTECTED]> wrote: > Hi dear all, > I am using Python2.4.2+NumPy1.0.1 to deal with a parameter > estimation problem with the least square methods. During the > calculations, I use NumPy package to deal with matrix operations, > mostly ma

Re: An error of matrix inversion using NumPy

2007-04-04 Thread Dayong Wang
Here below is the matrix KK I used: [[ 1939.33617572 -146.94170404 0. 0. 0. 0. 0. 0. 0. -1172.61032101 0. 0. -193.69962687 -426.08452381 0. 0. 0. 0.

Re: An error of matrix inversion using NumPy

2007-04-04 Thread BJörn Lindqvist
On 4 Apr 2007 06:15:18 -0700, lancered <[EMAIL PROTECTED]> wrote: > During the calculation, I noticed an apparent error of > inverion of a 19x19 matrix. Denote this matrix as KK, U=KK^ -1, I > found the product of U and KK is not equivalent to unit matrix! This > apparently violate the defi

An error of matrix inversion using NumPy

2007-04-04 Thread lancered
Hi dear all, I am using Python2.4.2+NumPy1.0.1 to deal with a parameter estimation problem with the least square methods. During the calculations, I use NumPy package to deal with matrix operations, mostly matrix inversion and trasposition. The dimentions of the matrices used are about 29x

An error of matrix inversion using NumPy

2007-04-04 Thread lancered
Hi dear all, I am using Python2.4.2+NumPy1.0.1 to deal with a parameter estimation problem with the least square methods. During the calculations, I use NumPy package to deal with matrix operations, mostly matrix inversion and trasposition. The dimentions of the matrices used are about 29x