PyLops – A Linear-Operator Python Library for large scale optimization

07/29/2019
by   Matteo Ravasi, et al.
0

Linear operators and optimisation are at the core of many algorithms used in signal and image processing, remote sensing, and inverse problems. For small to medium-scale problems, existing software packages (e.g., MATLAB, Python numpy and scipy) allow for explicitly building dense (or sparse) matrices and performing algebraic operations (e.g., computation of matrix-vector products and manipulation of matrices) with syntax that closely represents their corresponding analytical forms. However, many real application, large-scale operators do not lend themselves to explicit matrix representations, usually forcing practitioners to forego of the convenient linear-algebra syntax available for their explicit-matrix counterparts. PyLops is an open-source Python library providing a flexible and scalable framework for the creation and combination of so-called linear operators, class-based entities that represent matrices and inherit their associated syntax convenience, but do not rely on the creation of explicit matrices. We show that PyLops operators can dramatically reduce the memory load and CPU computations compared to explicit-matrix calculations, while still allowing users to seamlessly use their existing knowledge of compact matrix-based syntax that scales to any problem size because no explicit matrices are required.

READ FULL TEXT

page 2

page 3

page 4

page 5

page 6

page 7

page 9

page 11

research
04/22/2021

PyArmadillo: a streamlined linear algebra library for Python

PyArmadillo is a linear algebra library for the Python language, with th...
research
08/31/2017

Algorithmic patterns for H-matrices on many-core processors

In this work, we consider the reformulation of hierarchical (H) matrix a...
research
10/24/2016

Large Scale Parallel Computations in R through Elemental

Even though in recent years the scale of statistical analysis problems h...
research
09/12/2021

H2Opus: A distributed-memory multi-GPU software package for non-local operators

Hierarchical ℋ^2-matrices are asymptotically optimal representations for...
research
11/21/2019

Decomposition of arrow type positive semidefinite matrices with application to topology optimization

Decomposition of large matrix inequalities for matrices with chordal spa...
research
12/31/2018

Complexity of Linear Operators

Let A ∈{0,1}^n × n be a matrix with z zeroes and u ones and x be an n-di...
research
05/23/2022

Nancy: an efficient parallel Network Calculus library

This paper describes Nancy, a Network Calculus (NC) library that allows ...

Please sign up or login with your details

Forgot password? Click here to reset