lazylinop 1.6.0a2

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

Lazylinop philosophy – what is a LazyLinOp object?

Starting from a numpy array, a scipy matrix, a Faust object, or potentially many other compatible linear operators with efficient implementatons, a .LazyLinOp object follows the lazy evaluation paradigm. In short, one can aggregate low-level .LazyLinOp objects into higher-level ones using many classical operations (addition, concatenation, adjoint, real part, etc.), without actually building arrays. The actual effect of these operations is delayed until the resulting linear operator is actually applied to a vector (or to a collection of vectors, seen as a matrix).

The main interest of this paradigm is to enable the construction of processing pipelines that exploit as building blocks efficient implementations of “low-level” linear operators.

LazyLinOperators are complementary to other “lazy” objects such as LazyTensors in Kheops, or the ones of lazyarray, WeldNumpy an 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 is more likely to propose Python functions in the scipy style (like :py`.kron`).

Warning

lazylinop is currently in a beta stage of development. It is thus intensively evolving and we invite you to consult this documentation regularly.

Getting started with lazylinop

Quick installation

Using Anaconda

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

or

Using PIP

pip install lazylinop

For a straightforward installation of lazylinop go to this guide. Then for a quick introduction to lazylinop go to this notebook. The raw notebook is available here.