[SciPy-user] determinants

dmitrey openopt@ukr....
Thu Sep 13 12:10:46 CDT 2007


from scipy import linalg
from scipy import rand
for N in [10, 50, 100, 200, 400, 500, 1000]:
    print linalg.det(rand(N, N))

So typical A[i,j] values are 0.0 ... 1.0 (that is much less than your 
-13...54) and
typical output is
0.00356531521304
210176.6131
-6.3723083425e+24
1.45560352703e+80
8.58872496027e+217
-2.76708200113e+296
-inf

Same in MATLAB/Octave
So I don't see any bugs here.
Regards, D.


Alastair Basden wrote:
> Hi,
> I wonder whether there is a bug in scipy.linalg.det?  If I call this with 
> my matrix, it returns inf.  The matrix is 2048x2048 in size, with no inf 
> elements (max 54, min -13), float64 type.
> I can't work out why it would think the determinant is infinity!
>
> Thanks...
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>
>
>
>   



More information about the SciPy-user mailing list