A Haskell package for multidimensional arrays, simple tensor computations and multilinear algebra.
Array dimensions have an "identity" which is preserved in data manipulation. Indices are explicitly selected by name in expressions, and Einstein's summation convention for repeated indices is automatically applied.
The library has a purely functional interface: arrays are immutable, and operations work on whole structures which can be assembled and decomposed using simple primitives. Arguments are automatically made conformable by replicating them along extra dimensions appearing in an operation.
There is preliminary support for geometric algebra, multidimensional linear systems of equations, and tensor decompositions.
Application to Multiview Geometry:
Source repository: https://github.com/AlbertoRuiz/hTensor
We need GHC and hmatrix. See this page for installation help.
$ sudo apt-get install libgsl0-dev liblapack-dev
$ cabal update
$ cabal install hTensor
$ ghci
> import Numeric.LinearAlgebra.Exterior
> printA "%4.0f" $ leviCivita 4 !"pqrs" * cov (leviCivita 4)!"qrsu"
p^4 x u_4
u
-6 0 0 0
p 0 -6 0 0
0 0 -6 0
0 0 0 -6