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

Robert Kern robert.kern@gmail....
Tue Jan 27 15:10:15 CST 2009


On Tue, Jan 27, 2009 at 13:11, Nathan Bell <wnbell@gmail.com> wrote:
> 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.

+1

-- 
Robert Kern

"I have come to believe that the whole world is an enigma, a harmless
enigma that is made terrible by our own mad attempt to interpret it as
though it had an underlying truth."
  -- Umberto Eco


More information about the SciPy-user mailing list