[SciPy-user] Problems with linalg.cholesky?
Andrew Jaffe
a.h.jaffe at gmail.com
Fri Nov 18 07:20:19 EST 2005
Hi all,
In the newest incarnation of scipy_core, I am having trouble with the
cholesky(a) routine:
/Library/Frameworks/Python.framework/Versions/2.4/lib/python2.4/site-packages/scipy/linalg/basic_lite.py
in cholesky_decomposition(a)
115 else:
116 lapack_routine = lapack_lite.dpotrf
--> 117 results = lapack_routine('L', n, a, m, 0)
118 if results['info'] > 0:
119 raise LinAlgError, 'Matrix is not positive definite -
Cholesky decomposition cannot be computed'
LapackError: Parameter a is not contiguous in lapack_lite.dpotrf
But this isn't true; I get this error even when I pass trivial and
contiguous matrices such as the output of identity(). Other linalg
routines (included complicated ones like singular_value_decomp) seem to
work fine.
Any ideas?
Andrew
More information about the SciPy-User
mailing list