[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