[SciPy-user] Not a number in linalg.svdvals

[Pietro Berkes]

I think we're getting nearer... now try substituting your prints with

print a, a1, overwrite_a, a1 is not a

The only relevant difference I see between your and mine installation is the 
Python version: 2.1 vs. 2.2, although I don't see anything 2.2-specific in 
the code. Is there anybody using scipy with Python 2.1 who is able reproduce 
the bug?

Pietro.


> Compute singular values using svd
>
> [[-0.01064368  0.93221664 -0.35956157  0.03967781]
>  [ 0.02760821  0.35939368  0.93265212  0.01529682]
>  [ 0.35956157 -0.03967781 -0.01064368  0.93221664]
>  [-0.93265212 -0.01529682  0.02760821  0.35939368]] [ 1.29446282
> 1.00042554  0.99970394  0.72808433] [[-0.40511624  0.91346254
> -0.01377785  0.0357378 ]
>  [ 0.02830277  0.01255214  0.93261333  0.35954662]
>  [ 0.01552803 -0.03501285 -0.35945512  0.932376  ]
>  [-0.91369506 -0.40521936  0.02888879  0.01113738]] 0
> 140
>
> Singular values using svd
>
> [ 1.29446282  1.00042554  0.99970394  0.72808433]
>
> Compute singular values using svdvals
>
> [ [ 0.]] [              nan               nan
> nan               nan] [ [ 0.]] 0
> 268
>
> Singular values using svdvals
>
> [              nan               nan               nan
> nan]
>
> Nils
>