A user-friendly hybrid sparse matrix class in C++
From MaRDI portal
Publication:1662286
DOI10.1007/978-3-319-96418-8_50zbMATH Open1397.65343DBLPconf/icms/SandersonC18arXiv1805.03380OpenAlexW2802799492WikidataQ61013238 ScholiaQ61013238MaRDI QIDQ1662286FDOQ1662286
Authors: Conrad Sanderson, Ryan R. Curtin 
Publication date: 17 August 2018
Abstract: When implementing functionality which requires sparse matrices, there are numerous storage formats to choose from, each with advantages and disadvantages. To achieve good performance, several formats may need to be used in one program, requiring explicit selection and conversion between the formats. This can be both tedious and error-prone, especially for non-expert users. Motivated by this issue, we present a user-friendly sparse matrix class for the C++ language, with a high-level application programming interface deliberately similar to the widely used MATLAB language. The class internally uses two main approaches to achieve efficient execution: (i) a hybrid storage framework, which automatically and seamlessly switches between three underlying storage formats (compressed sparse column, coordinate list, Red-Black tree) depending on which format is best suited for specific operations, and (ii) template-based meta-programming to automatically detect and optimise execution of common expression patterns. To facilitate relatively quick conversion of research code into production environments, the class and its associated functions provide a suite of essential sparse linear algebra functionality (eg., arithmetic operations, submatrix manipulation) as well as high-level functions for sparse eigendecompositions and linear equation solvers. The latter are achieved by providing easy-to-use abstractions of the low-level ARPACK and SuperLU libraries. The source code is open and provided under the permissive Apache 2.0 license, allowing unencumbered use in commercial products.
Full work available at URL: https://arxiv.org/abs/1805.03380
Recommendations
- Dense and sparse matrix classes using the \(C++\) Standard Template Library
- Sparse Matrices in MATLAB: Design and Implementation
- SparseX: a library for high-performance sparse matrix-vector multiplication on multicore platforms
- Computational Science – ICCS 2005
- Object-oriented techniques for sparse matrix computations in Fortran 2003
Graphs and linear algebra (matrices, eigenvalues, etc.) (05C50) Packaged methods for numerical algorithms (65Y15) Numerical linear algebra (65Fxx)
Cited In (15)
- Efficient alternating least squares algorithms for low multilinear rank approximation of tensors
- Methods for suspensions of passive and active filaments
- Multilevel convergence analysis of multigrid-reduction-in-time
- KAPPA: kinetic approach to physical processes in atmospheres library in C++
- Direct statistical inference for finite Markov jump processes via the matrix exponential
- Systematic large flavor fTWA approach to interaction quenches in the Hubbard model
- Implicit-explicit multirate infinitesimal stage-restart methods
- Energy-conserving hyper-reduction and temporal localization for reduced order models of the incompressible Navier-Stokes equations
- Time- and ensemble-average statistical mechanics of the Gaussian network model
- Dense and sparse matrix classes using the \(C++\) Standard Template Library
- FELINE: finite element solver for hydrodynamic lubrication problems using the inexact Newton method
- Stokes-Cahn-Hilliard formulation in sliding bi-periodic frames for the simulation of two-phase flows
- FELICITY: a Matlab/C++ toolbox for developing finite element methods and simulation modeling
- Quantum illumination with noisy probes: conditional advantages of non-Gaussianity
- Compressed sensing reconstruction using expectation propagation
Uses Software
This page was built for publication: A user-friendly hybrid sparse matrix class in C++
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1662286)
