[Numeric] Announce: comfort-array, lapack

Henning Thielemann lemming at henning-thielemann.de
Sat May 25 14:42:51 UTC 2019


I like two announce two of my packages:

1. comfort-array
http://hackage.haskell.org/package/comfort-array

It provides Boxed and Storable arrays with very liberal shape definitions. 
You may use ranges of indices like in 'array' or zero-based indexing like 
in 'repa', but you can also use arbitrary index Sets, enumerations, 
concatenation of arrays and more. E.g. an array with (Shape.Enumeration 
Ordering) has three elements with indices LT, EQ, GT.

2. lapack
http://hackage.haskell.org/package/lapack

A high-level interface to the numerical Linear Algebra package LAPACK. It 
is based on comfort-array, which means that a vector can have a rich 
structure. E.g. matrix vector multiplication has the signature:
    (#*|) :: Matrix.General height width a -> Vector width a -> Vector height a

Thus, a 'Matrix.General (Shape.Enumeration Ordering) (h,w) a' would map a 
rectangular array with dimensions h and w to a three-element vector with 
indices LT, EQ, GT.

The LAPACK interface provides solvers for simultaneous linear equations, 
linear least-squares problems, eigenvalue and singular value problems for 
square, triangular, symmetric, Hermitian, banded and banded Hermitian 
matrices, as LAPACK supports them.


For a motivating example refer to:
    http://code.henning-thielemann.de/bob2019/main.pdf


More information about the Numeric mailing list