[SciPy-User] Ax = b for symmetric positive definite matrix
Dag Sverre Seljebotn
Wed Dec 8 12:03:39 CST 2010
On 12/08/2010 04:16 PM, Mathieu Blondel wrote:
> On Wed, Dec 8, 2010 at 10:47 PM, Dag Sverre Seljebotn
> <email@example.com> wrote:
>> It is not supported in SciPy yet, but if you want to play with calling
>> LAPACK directly, this is supported through sppsv/dppsv/cppsv/zppsv. For
>> instance you could modify the scipy/linalg/*.pyf files to include the
>> function, rebuild SciPy. Or, search for "tokyo cython lapack" and play
>> with that.
> Ideally, I would like to include the relevant Fortran or C files
> directly in my project so that it is easy to compile for other users.
> In this regard, maybe CLAPACK + Cython seem like a better choice than
> LAPACK + f2py?
To have any performance at all you *need* to link to whatever LAPACK the
user is using (ATLAS, Intel MKL, etc.). The most pain-free way to do
that currently is to patch/improve SciPy I believe.
More information about the SciPy-User