sparse-linear-algebra: Numerical computation in native Haskell

[ gpl, library, numeric ] [ Propose Tags ]

Currently it provides iterative linear solvers, matrix decompositions, eigenvalue computations and related utilities. Please see for details

[Skip to Readme]
Versions [faq],,,,,,,,,,,,,,,, 0.2.9,,,,,,,,,, 0.3, 0.3.1 (info)
Dependencies base (>=4.7 && <5), containers, hspec, mtl (>=2.2.1), mwc-random, primitive (>=, QuickCheck [details]
License GPL-3.0-only
Copyright 2016 Marco Zocca
Author Marco Zocca
Maintainer zocca.marco gmail
Category Numeric
Home page
Source repo head: git clone
Uploaded by ocramz at Tue Oct 18 03:58:21 UTC 2016
Distributions LTSHaskell:0.3.1, NixOS:0.3.1, Stackage:0.3.1
Executables sparse-linear-algebra
Downloads 6672 total (373 in the last 30 days)
Rating 2.0 (votes: 1) [estimated by rule of succession]
Your Rating
  • λ
  • λ
  • λ
Status Hackage Matrix CI
Docs available [build log]
Last success reported on 2016-10-18 [all 1 reports]




Maintainer's Corner

For package maintainers and hackage trustees

Readme for sparse-linear-algebra-

[back to package description]


Numerical computation in native Haskell

TravisCI : Build Status

This library provides common numerical analysis functionality, without requiring any external bindings. It is not optimized for performance (yet), but it serves as an experimental platform for scientific computation in a purely functional setting.

Algorithms :

  • Iterative linear solvers

    • Conjugate Gradient Squared (CGS)

    • BiConjugate Gradient Stabilized (BiCGSTAB) (non-Hermitian systems)

  • Matrix decompositions

    • QR factorization
  • Eigenvalue algorithms

    • QR algorithm

    • Rayleigh quotient iteration

  • Utilities : Vector and matrix norms, matrix condition number, Givens rotation, Householder reflection

  • Predicates : Matrix orthogonality test (A^T A ~= I)

This is also an experiment in principled scientific programming :

  • set the stage by declaring typeclasses and some useful generic operations (normed linear vector spaces, i.e. finite-dimensional spaces equipped with an inner product that induces a distance function),

  • define appropriate data structures, and how they relate to those properties (sparse vectors and matrices, defined internally via Data.IntMap, are made instances of the VectorSpace and Additive classes respectively). This allows to decouple the algorithms from the actual implementation of the backend,

  • implement the algorithms, following 1:1 the textbook [1]




Inspired by


[1] : Y. Saad, Iterative Methods for Sparse Linear Systems, 2nd ed., 2000