[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