Abstract:There are many classes of mathematical problems which give rise to matrices, where a large number of the elements are zero. In this case it makes sense to have a special matrix type to handle this class of problems where only the non-zero elements of the matrix are stored. Not only does this reduce the amount of memory to store the matrix, but it also means that operations on this type of matrix can take advantage of the a-priori knowledge of the positions of the non-zero elements to accelerate their calculations. A matrix type that stores only the non-zero elements is generally called sparse. Until recently Octave has lacked a full implementation of sparse matrices. This article address the implementation of sparse matrices within Octave, including their storage, creation, fundamental algorithms used, their implementations and the basic operations and functions implemented for sparse matrices. Mathematical issues such as the return types of sparse operations, matrix fill-in and reordering for sparse matrix factorization is discussed in the context of a real example. Benchmarking of Octave's implementation of sparse operations compared to their equivalent in Matlab are given and their implications discussed. Results are presented for multiplication and linear algebra operations for various matrix orders and densities. Furthermore, the use of Octave's sparse matrix implementation is demonstrated using a real example of a finite element model (FEM) problem. Finally, the method of using sparse matrices with Octave's oct-files is discussed. The means of creating, using and returning sparse matrices within oct-files is discussed as well as the differences between Octave's Sparse and Array classes.

Awkward Arrays in Python, C++, and Numba

The Awkward World of Python and C++

Awkward Just-In-Time (JIT) Compilation: A Developer's Experience

Awkward to RDataFrame and back

Bridging Worlds: Achieving Language Interoperability between Julia and Python in Scientific Computing

Array Programming with NumPy

Describe Data to get Science-Data-Ready Tooling: Awkward as a Target for Kaitai Struct YAML

An array-oriented Python interface for FastJet

Finch: Sparse and Structured Array Programming with Control Flow

PARRAY

mdspan in C++: A Case Study in the Integration of Performance Portable Features into International Language Standards

Efficiently manipulating Pauli strings with PauliArray

Array Requirements for Scientific Applications and an Implementation for Microsoft SQL Server

Enabling MPI communication within Numba/LLVM JIT-compiled Python code using numba-mpi v1.0

BioNumPy: array programming for biology

Associative Arrays: Unified Mathematics for Spreadsheets, Databases, Matrices, and Graphs

Verifying spatial properties of array computations

Sparse Matrix Implementation in Octave

PyArmadillo: a streamlined linear algebra library for Python

Difference Covering Arrays and Pseudo-Orthogonal Latin Squares

DArray: A High Performance RDMA-Based Distributed Array.