[SciPy-user] scipy.sparse.linalg.cg not thread safe?
Tue Jan 27 13:11:49 CST 2009
On Mon, Jan 26, 2009 at 6:29 PM, Robert Kern <email@example.com> wrote:
> It's probably not threadsafe.
I don't know Fortran, so I can't say:
Anway, here's a pure-SciPy CG implementation (also BSD-licensed):
It should be a drop-in replacement for sparse.linalg.cg() and have
comparable speed. The only dependency is to the norm() function in
PyAMG, but you can swipe that easily too.
In time we should replace all of the Fortran implementations of the
iterative methods with pure-Python code. This would be a nice target
for SciPy 0.8.
Nathan Bell firstname.lastname@example.org
More information about the SciPy-user