[SciPy-dev] Tri-diagonal LAPACK Routines - Shall I interface them?

David Goldsmith d.l.goldsmith@gmail....
Sun Dec 6 20:31:54 CST 2009


On Sun, Dec 6, 2009 at 6:13 PM, Nathan Bell <wnbell@gmail.com> wrote:

> On Wed, Dec 2, 2009 at 3:43 AM, Benny Malengier
> <benny.malengier@gmail.com> wrote:
> >
> > Interesting, this is exactly the function I needed for my problem, but
> > I was looking in scipy.sparse.linalg, so did not notice banded matrix
> > solver was present in scipy.linalg.
> >
> > In my logic, the "matrix diagonal orded form" of
> >
> http://docs.scipy.org/doc/scipy/reference/generated/scipy.linalg.solve_banded.html#scipy.linalg.solve_banded
> > would be a type of sparse matrix one can manipulate. This would allow
> > things like changing matrix diagonal orded form sparse matrix to a csr
> > matrix, adding some extra elements off the diagonals, and then calling
> > a more generic solver.
> >
>
> Can't you do that already with scipy.sparse.dia_matrix?  If I'm not
> mistaken, dia_matrix is (slightly) more general than the banded format
> but similarly efficient.
>
> In an ideal world scipy.sparse.spsolve() would detect the case that A
> was a dia_matrix (with small bandwidth) and invoke the LAPACK method
> in scipy.linalg instead of using the general sparse LU solver.
>

Perhaps you might file this as an "enhancement" "issue"?

DG


>
> --
> Nathan Bell wnbell@gmail.com
> http://www.wnbell.com/
> _______________________________________________
> SciPy-Dev mailing list
> SciPy-Dev@scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-dev
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/scipy-dev/attachments/20091206/d56e0235/attachment-0001.html 


More information about the SciPy-Dev mailing list