This is a high-level interface to LAPACK. It provides solvers for simultaneous linear equations, linear least-squares problems, eigenvalue and singular value problems for matrices with certain kinds of structures. Features: Based on comfort-array: Allows to precisely express one-column or one-row matrices, as well as dense, square, triangular, banded, symmetric, Hermitian, banded Hermitian, blocked and LU or QR decomposed matrices. Support all data types that are supported by LAPACK, i.e. Float, Double, Complex Float, Complex Double No need for c2hs, hsc, Template Haskell or C helper functions Dependency only on BLAS and LAPACK, no GSL Works with matrices and vectors with zero dimensions. This one is tricky and still leads to surprises since different LAPACK implementations consider different situations as corner cases. No automatic (and dangerous) implicit expansion of singleton vectors or matrices. Instead there are special operators for scaling of vectors and matrices. Separate formatting operator (##): Works better for tuples of matrices and vectors than show. Show is used for code one-liners that can be copied back into Haskell modules. Support for nice formatting in HyperHaskell. See also: hmatrix. For examples see the Example modules and the packages sum-pyramid, wuerfelschlange, magico, linear-circuit, resistor-cube.