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

Ergun Bicici ebicici@ku.edu...
Tue Jan 27 15:36:16 CST 2009


Sounds good. +1 :)

Ergun Bicici
Koc University


On Tue, Jan 27, 2009 at 11:10 PM, Robert Kern <robert.kern@gmail.com> wrote:

> 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
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://projects.scipy.org/pipermail/scipy-user/attachments/20090127/67218628/attachment.html 


More information about the SciPy-user mailing list