[SciPy-user] scipy.sparse.linalg.cg not thread safe?

Nathan Bell wnbell@gmail....
Tue Jan 27 13:11:49 CST 2009


On Mon, Jan 26, 2009 at 6:29 PM, Robert Kern <robert.kern@gmail.com> wrote:
>
> It's probably not threadsafe.
>

I don't know Fortran, so I can't say:
http://projects.scipy.org/scipy/scipy/browser/trunk/scipy/sparse/linalg/isolve/iterative/CGREVCOM.f.src


Anway, here's a pure-SciPy CG implementation (also BSD-licensed):
http://code.google.com/p/pyamg/source/browse/trunk/pyamg/krylov/cg.py

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 wnbell@gmail.com
http://graphics.cs.uiuc.edu/~wnbell/


More information about the SciPy-user mailing list