[SciPy-User] solving large sparse linear system with Laplacian matrix

David Warde-Farley dwf@cs.toronto....
Fri Oct 30 05:38:26 CDT 2009


On 30-Oct-09, at 5:56 AM, Pauli Virtanen wrote:

> Fri, 30 Oct 2009 05:28:13 -0400, David Warde-Farley wrote:
> [clip]
>> I'm not sure they exist in scipy.sparse.linalg, but looking into your
>> problem a little further it appears that several options exist in  
>> LAPACK
>> (but are not wrapped in SciPy, yet). Conveniently, they have built-in
>> support for multiple right-hand-sides.
>>
>> Solve a general banded system of linear equations:
>> 	http://www.netlib.org/lapack/double/dgbsv.f
>
> This is available as scipy.linalg.solve_banded, but I guess the rest  
> of
> the gang is missing. (They'd be great to have *wink wink nudge  
> nugde* :)

Oh, duh. I was looking in scipy.linalg.fblas when I should've been  
looking in scipy.linalg.flapack.

David


More information about the SciPy-User mailing list