lazylinop 1.8.0

A package dedicated to lazy linear operators based on diverse backends/libraries.

Lazylinop philosophy – what is a LazyLinOp object?

A LazyLinOp is a high-level linear operator based on an arbitrary underlying implementation, such as:

• custom Python functions,

• a NumPy array,

• a SciPy matrix,

• a Faust object,

• any Python linear operator.

Thanks to the Lazylinop API, this operator can be easily manipulated, transformed or aggregated with other linear operators to form more complex LazyLinOp objects. Thus, many operations are available such as the addition, concatenation, adjoint etc.

These operations are all ruled by the lazy paradigm: their evaluation is delayed until the resulting LazyLinOp is actually applied to a vector (or to a collection of vectors, seen as a matrix).

Main benefits of Lazylinop

Lazylinop aims at providing:

• a unified NumPy/SciPy-like API for any underlying implementation (including transparent use of efficient third-party libraries for advanced linear operators such as wavelets, etc.),

• an easy way to combine existing operators into more complex operators with all the LazyLinOp operations at your fingertips (and in particular the adjoint),

• potential performance gains & memory sparing offered by the lazy paradigm,

• many pre-built operators and related APIs, delivered in specialized Lazylinop sub-packages (polynomials of operators, signal-processing oriented operators, etc.).

• Compatibility with SciPy algorithms implemented for LinearOperator-s: svds, eigsh, etc.

Lazylinop is also an alternative or a complement to other related libraries that harness the power of the lazy paradigm. See the section Related lazy libraries for more details.

Beta stage admonition warning

Lazylinop is currently in a beta stage of development. It is thus intensively evolving and we invite you to consult this documentation regularly. We welcome all beta testers and any feedback.

Getting started with lazylinop

Quick install using Anaconda

conda config --add channels conda-forge
conda config --add channels pyfaust
conda config --add channels lazylinop
conda install lazylinop

Quick install using PIP

pip install lazylinop

Python beginners willing to safely install Lazylinop with a virtual environment can consult this guide.

Quick introduction

For a quick introduction to the use of Lazylinop go to this notebook. The raw notebook is available here.

Related lazy libraries

LazyLinOp-s are complementary to other “lazy” objects such as LazyTensors in KeOps, or the ones of lazyarray, WeldNumpy and PyLops libraries, which, to the best of our knowledge, primarily provide compact descriptions of arrays which entries can be evaluated efficiently on the fly.

It’s worth noting that the way PyLops exhibits SciPy LinearOperator-s in its API is different. Indeed, PyLops provides Python classes (like Kronecker) while Lazylinop rather provides Python functions in the SciPy style (as kron ).