[SciPy-dev] Tri-diagonal LAPACK Routines - Shall I interface them?
Wed Dec 2 02:43:12 CST 2009
2009/12/2 David Warde-Farley <email@example.com>:
> On 1-Dec-09, at 8:34 PM, Simon Clift wrote:
>> Is there a broader reason for not completing the interface?
> Nope. I think they were just wrapped on an as-needed basis. This came
> up a little while ago w.r.t. general banded systems as well:
>> Would anyone object if I went ahead and patched scipy/linalg/
>> generic_flapack.pyf to allow
>> the call?
> Doubtful. :) Although if you do, it might be a good idea to add a
> function to scipy.linalg with a comprehensible name (like
> solve_tridiagonal to go along with solve_banded [which actually does
> wrap dgbsv, my post there is wrong in that regard]).
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
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.
> SciPy-Dev mailing list
More information about the SciPy-Dev