[SciPy-dev] problem with linalg.cholesky?

Travis Oliphant oliphant.travis at ieee.org
Fri Dec 2 15:43:07 CST 2005


Andrew Jaffe wrote:

>I still do get the problem with the latest version (rev 1550). I thought 
>I was using the optimized lapack -- it's automatic on OS X, isn't it? 
>Does the fact that the problem is occuring during a lapack_lite call 
>imply something else?
>  
>
O.K.  That's good to know.

>If there's no problem there, I guess the questions are: does 
>_castCopyAndTranspose result in a non-contiguous array on all platforms? 
>And does lapack_lite.dpotrf() balk at that, and why?
>  
>
No, it doesn't, that's what has me stumped.  It should be returning a 
contiguous array (because of the .astype method).

So, figuring out why the .astype() method is not returning a contiguous 
array is the key to understanding this, I think.

If you could determine what the parameters being passed to the .astype() 
method are and try to reproduce the problem of not getting a contiguous 
parameter on output, that would be very helpful.

Sorry for your trouble.

-Travis




More information about the Scipy-dev mailing list