Есть ли в стандартной библиотеке модуль (пакет?) для работы с матрицами?

matrix :: BLAS3 e => (Int, Int) -> [((Int, Int), e)] -> Matrix (n, p) e

numRows :: MatrixShaped a => a mn e -> Int

-- Матрицы, как списки списков:
 
-- Элемент по индексу
 
getByIndex :: Num a => [[a]] -> Int -> Int -> a
getByIndex m i j = if (i > (length m)) then error "!" else if (j > (length row)) then error "!" else row !! j where row = m !! i
 
-- дать строку
 
getRow :: Num a => [[a]] -> Int -> [a]
getRow m i = if (i > (length m)) then error "!" else m !! i
 
-- Дать столбец
 
getCol :: Num a => [[a]] -> Int -> [a]
getCol m j = map (\ x -> x !! j) m

import Numeric.LinearAlgebra
 
testMatrix :: Matrix Double
testMatrix = (3 >< 4) [1..]
 
test :: Show a => String -> a -> IO ()
test str val = do
  putStr $ str ++ ": "
  print val
 
getRow :: Element a => Int -> Matrix a -> Vector a
getRow rowNum = flatten . extractRows [rowNum]
 
getColumn :: Element a => Int -> Matrix a -> Vector a
getColumn colNum mat = flatten $ subMatrix (0, colNum) (rows mat, 1) mat
 
main :: IO ()
main = do
  test "Initial matrix" $ testMatrix
  test "Element at index (2, 1)" $ atIndex testMatrix (2, 1)
  test "Element at index (2, 1) (different syntax)" $ testMatrix @@> (2, 1)
  test "Number of rows" $ rows testMatrix
  test "Number of columns" $ cols testMatrix
  test "Submatrix at index (1, 1) of size (2, 3)" $
    subMatrix (1, 1) (2, 3) testMatrix
  test "Extract row 1 as vector" $ getRow 1 testMatrix
  test "Extract column 2 as vector" $ getColumn 2 testMatrix

GHCi, version 7.6.3: [url]http://www.haskell.org/ghc/[/url]  :? for help
Loading package ghc-prim ... linking ... done.
Loading package integer-gmp ... linking ... done.
Loading package base ... linking ... done.
Prelude> :m + Data.Packed.Matrix
Prelude Data.Packed.Matrix> (2><3)[1..6]
Loading package array-0.4.0.1 ... linking ... done.
Loading package deepseq-1.3.0.1 ... linking ... done.
Loading package primitive-0.5.0.1 ... linking ... done.
Loading package vector-0.10.0.1 ... linking ... done.
Loading package containers-0.5.0.0 ... linking ... done.
Loading package bytestring-0.10.0.2 ... linking ... done.
Loading package binary-0.5.1.1 ... linking ... done.
Loading package Win32-2.3.0.0 ... linking ... done.
Loading package filepath-1.3.0.1 ... linking ... done.
Loading package old-locale-1.0.0.5 ... linking ... done.
Loading package time-1.4.0.1 ... linking ... done.
Loading package directory-1.2.0.1 ... linking ... done.
Loading package process-1.1.0.2 ... linking ... done.
Loading package storable-complex-0.2.1 ... linking ... done.
Loading package hmatrix-0.13.1.0 ... can't load .so/.DLL for: gsl-0.dll (addDLL: could not load DLL)

{-# LANGUAGE NoMonomorphismRestriction #-}
import Data.List
import Data.Complex
import System.IO
import System.Directory
import Data.Ratio
 
--Операция с матрицами, с исп-ем соотв. элементов из обеих матриц и запись их в соотв. элементы 3-й.
mMOp :: (a -> b -> c) -> [[a]] -> [[b]] -> [[c]]
mMOp op = zipWith (zipWith op)
 
mPlus :: Num a => [[a]] -> [[a]] -> [[a]]
mPlus = mMOp (+)
 
mMinus :: Num a => [[a]] -> [[a]] -> [[a]]
mMinus = mMOp (-)
 
--vector scalar multiplication 
vSMul :: Num a => [a] -> [a] -> a
vSMul x y= sum$zipWith (*) x y
 
vNMul number vector = map (*number) vector
 
vPlus = zipWith (+)
 
vMinus = zipWith (-)
 
mMul :: Num t => [[t]] -> [[t]] -> [[t]]
mMul xm ym = [[vSMul xml yml| yml<-(transpose ym)]|xml<-xm]
 
 
vAbs :: Floating a=> [a]->a
vAbs a =  (sqrt . sum . map (^2)) a
 
--Убрать нный элемент из списка
dropNE n list = take (n-1) list ++ drop n list 
 
--drop a cross from matrix
dropXM a b matrix = map (dropNE b) (dropNE a matrix)
 
--Убрать и получить нный элемент списка
dropGetNE n l = (l!!(n-1),take (n-1) l ++ drop n l)
 
dropGetXM a b m = let (x,y) = (dropGetNE a m) in (x!!(b-1), map (dropNE b) y )
 
-- determenant of matrix
detM2 [[m]] = m
--detM m = sum [let (x,y) = dropGetXM a a m in (-x*detM y) | a<-[1..(length m)]]
detM2 m = sum [let (x,y) = dropGetXM a 1 m in ((-1)^a*x*detM2 y) | a<-[1..(length m)]]
 
--Наибыстрейший из общих реализован по определению:
detM3 [[m]] = m
detM3 m = sum [let (x,y) = dropGetN1XM a m in ((-1)^a*x*detM3 y) | a<-[1..(length m)]]
 
dropGetN1XM n m = let (x,y) = (dropGetNE n m) in (head x, map tail y )
 
--Наибыстрейший из здешних, использует операцию деления, приводит матрицы к треуг. виду:
detM [[m]] = m
detM m = 
    let
        (x,(y:ys)) =  let z = map (getdivBy1stWOut1) m in (\(x,y)->(x,reverse y)) (foldl' (\(w,x) (y,z) -> (w*y,(z:x))) (1,[]) z)
    in
    x*detM (map (\x -> zipWith (-) x y) ys) --map (\a -> zipWith (-) a x ) y
 
getdivBy1stWOut1 (l:ls) = (l,map (/l) ls)

detM1 [[m]] =  m
detM1 m = 
    let 
        k = findIndex (\x->head x /= 0) m
        f (mh:mt) = 
            let (prodm,listm) = getDivBy1stWOut1 mh
            in (\(x,y)->x*detM1 (reverse y)) 
                    (foldl'
                        (\(w1,w2) x->
                            (\(y1,y2)->(y1*w1,(y2:w2)))
                                (if 
                                    head x /= 0 
                                then 
                                    (\(x1,x2)->
                                        (x1,zipWith (-) x2 listm)
                                    ) (getDivBy1stWOut1 x)
                                else (1,tail x)
                                )
                        ) (prodm,[]) mt
                    )
    in
        (case k of 
        Nothing-> 0
        Just 0 -> f m
        Just x -> negate$f (replaceGet1Index (x-1) (head m) (tail m)))

(>>=)([(.),(.)]>>=(\(?)->((?)id id)"."))((.)(.)(++)[](\(.) -> ['(',(.),')']))