[Scipy-tickets] [SciPy] #1189: Define operations on LinearOperator instances

SciPy Trac scipy-tickets@scipy....
Tue Jun 1 10:02:57 CDT 2010


#1189: Define operations on LinearOperator instances
------------------------------------------------+---------------------------
 Reporter:  barbey                              |       Owner:  wnbell
     Type:  enhancement                         |      Status:  new   
 Priority:  normal                              |   Milestone:  0.8.0 
Component:  scipy.sparse.linalg                 |     Version:  0.7.0 
 Keywords:  LinearOperator, sparse, operations  |  
------------------------------------------------+---------------------------
 The LinearOperator class is very convenient to handle large linear
 operator which can be spare or defined analytically but cannot be stored
 as matrices. This is a very common situation in large linear inverse
 problem. As examples you can think of tomography (radon transform),
 Fourier transform (with non uniform step for instance). I personally work
 on sky map estimation using data from the Herschel Space Telescope. This
 is a huge linear inverse problem which can be solved using
 scipy.sparse.linalg.cgs type of methods.

 Unfortunately as of now the LinearOperator class lack some functionalities
 that I think would be very convenient for my purpose. I would like to have
 the possibility to add or multiply 2 linear models. Or multiply a linear
 model by a scalar for instance. For this purpose, I modified __mul__ and
 created __add__ methods of the LinearOperator class and their derivative.
 I join the code in attachment for review.

 Please be patient since this is my first attempt of contribution to scipy
 (or any open source project).

-- 
Ticket URL: <http://projects.scipy.org/scipy/ticket/1189>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list