[Scipy-tickets] [SciPy] #1766: Unify scipy.linalg.solve_banded and scipy.sparse.dia_matrix
SciPy Trac
scipy-tickets@scipy....
Mon Nov 5 05:26:20 CST 2012
#1766: Unify scipy.linalg.solve_banded and scipy.sparse.dia_matrix
----------------------------+-----------------------------------------------
Reporter: johntyree | Owner: somebody
Type: enhancement | Status: new
Priority: normal | Milestone: Unscheduled
Component: Other | Version: 0.11.0
Keywords: sparse, linalg |
----------------------------+-----------------------------------------------
Creating dia_matrix objects and then deconstructing them manually to pass
to solve_banded is unintuitive and error prone.
Is it possible to (at a minimum), wrap solve_banded and make it a method
of dia_matrix? Otherwise you have to convert the "offset" attribute to
solve_banded's (l,u) tuple of "upper and lower diagonals". The "data"
attribute of dia_matrix objects is already in the correct form.
Even something as simple as:
{{{
L = dia_matrix([ones(10)], 1, 10, 10)
R = ones(10)
L.solve(R)
}}}
would make sense.
In fact, it seems like this would be a reasonable method to have on any
Matrix object.
--
Ticket URL: <http://projects.scipy.org/scipy/ticket/1766>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.
More information about the Scipy-tickets
mailing list