Copyright | (c) Armando Santos 2019-2020 |
---|---|
Maintainer | armandoifsantos@gmail.com |
Stability | experimental |
Safe Haskell | None |
Language | Haskell2010 |
- Type safe matrix representation
- Constraint type aliases
- Type aliases
- Primitives
- Auxiliary type families
- Matrix construction and conversion
- Functor instance equivalent function
- Applicative/Monoidal instance equivalent functions
- Selective equivalent instance function
- Monad equivalent instance function
- Misc
- Biproduct approach
- Matrix composition and lifting
- Relation
- Matrix printing
The LAoP discipline generalises relations and functions treating them as Boolean matrices and in turn consider these as arrows.
LAoP is a library for algebraic (inductive) construction and manipulation of matrices in Haskell. See my Msc Thesis for the motivation behind the library, the underlying theory, and implementation details.
This module offers a newtype wrapper around Matrix
that
uses arbitrary types instead of canonical data types for the matrices
dimensions.
NOTE: If the types in the dimensions are custom they must need to
implement a Generic
instance.
Synopsis
- newtype Matrix e (cols :: Type) (rows :: Type) = M (Matrix e (Normalize cols) (Normalize rows))
- type Countable a = KnownNat (Count a)
- type CountableDimensions a b = (Countable a, Countable b)
- type CountableN a = KnownNat (Count (Normalize a))
- type CountableDimensionsN a b = (CountableN a, CountableN b)
- type FromListsN e a b = FromLists e (Normalize a) (Normalize b)
- type Liftable e a b = (Bounded a, Bounded b, Enum a, Enum b, Eq b, Num e, Ord e)
- type Trivial a = Normalize (Normalize a) ~ Normalize (Normalize (Normalize a))
- type TrivialP a b = Normalize (a, b) ~ Normalize (Normalize a, Normalize b)
- type Zero = Void
- type One = ()
- empty :: Matrix e Zero Zero
- one :: e -> Matrix e One One
- join :: Matrix e a rows -> Matrix e b rows -> Matrix e (Either a b) rows
- fork :: Matrix e cols a -> Matrix e cols b -> Matrix e cols (Either a b)
- type family FromNat (n :: Nat) :: Type where ...
- type family Count (d :: Type) :: Nat where ...
- type family Normalize (d :: Type) :: Type where ...
- class FromLists e cols rows
- fromLists :: FromListsN e cols rows => [[e]] -> Matrix e cols rows
- toLists :: Matrix e cols rows -> [[e]]
- toList :: Matrix e cols rows -> [e]
- matrixBuilder' :: (FromListsN e cols rows, CountableDimensionsN cols rows) => ((Int, Int) -> e) -> Matrix e cols rows
- matrixBuilder :: (FromListsN e a b, Enum a, Enum b, Bounded a, Bounded b, Eq a, CountableDimensionsN a b) => ((a, b) -> e) -> Matrix e a b
- row :: FromListsN e cols () => [e] -> Matrix e cols One
- col :: FromListsN e () rows => [e] -> Matrix e One rows
- zeros :: (Num e, FromListsN e cols rows, CountableDimensionsN cols rows) => Matrix e cols rows
- ones :: (Num e, FromListsN e cols rows, CountableDimensionsN cols rows) => Matrix e cols rows
- bang :: forall e cols. (Num e, Enum e, FromListsN e cols (), CountableN cols) => Matrix e cols One
- point :: (Bounded a, Enum a, Eq a, Num e, Ord e, CountableN a, FromListsN e a One) => a -> Matrix e One a
- constant :: (Num e, FromListsN e cols rows, CountableDimensionsN cols rows) => e -> Matrix e cols rows
- fmapM :: (Liftable e a b, CountableDimensionsN a b, FromListsN e b a) => (a -> b) -> Matrix e c a -> Matrix e c b
- bimapM :: (Liftable e a b, Liftable e c d, CountableDimensionsN a c, CountableDimensionsN b d, FromListsN e d c, FromListsN e b a) => (a -> b) -> (c -> d) -> Matrix e a c -> Matrix e b d
- unitM :: Num e => Matrix e () ()
- multM :: (CountableDimensionsN a b, CountableN (a, b), Num e, FromListsN e (a, b) a, FromListsN e (a, b) b, TrivialP a b) => Matrix e c a -> Matrix e c b -> Matrix e c (a, b)
- selectM :: (Num e, FromListsN e b b, CountableN b) => Matrix e cols (Either a b) -> Matrix e a b -> Matrix e cols b
- returnM :: forall e a. (Num e, Enum e, Enum a, FromListsN e () a, Countable a) => a -> Matrix e One a
- bindM :: Num e => Matrix e a b -> Matrix e b c -> Matrix e a c
- columns :: CountableN cols => Matrix e cols rows -> Int
- rows :: CountableN rows => Matrix e cols rows -> Int
- tr :: Matrix e cols rows -> Matrix e rows cols
- (.|) :: Num e => e -> Matrix e cols rows -> Matrix e cols rows
- (./) :: Fractional e => Matrix e cols rows -> e -> Matrix e cols rows
- cond :: (Trivial a, Trivial2 a, Trivial3 a, CountableN a, FromListsN e () a, FromListsN e a (), FromListsN e a a, Liftable e a Bool) => (a -> Bool) -> Matrix e a b -> Matrix e a b -> Matrix e a b
- abideJF :: Matrix e cols rows -> Matrix e cols rows
- abideFJ :: Matrix e cols rows -> Matrix e cols rows
- zipWithM :: (e -> f -> g) -> Matrix e a b -> Matrix f a b -> Matrix g a b
- (===) :: Matrix e cols a -> Matrix e cols b -> Matrix e cols (Either a b)
- p1 :: (Num e, CountableDimensionsN n m, FromListsN e n m, FromListsN e m m) => Matrix e (Either m n) m
- p2 :: (Num e, CountableDimensionsN n m, FromListsN e m n, FromListsN e n n) => Matrix e (Either m n) n
- (|||) :: Matrix e a rows -> Matrix e b rows -> Matrix e (Either a b) rows
- i1 :: (Num e, CountableDimensionsN n m, FromListsN e n m, FromListsN e m m) => Matrix e m (Either m n)
- i2 :: (Num e, CountableDimensionsN n m, FromListsN e m n, FromListsN e n n) => Matrix e n (Either m n)
- (-|-) :: (Num e, CountableDimensionsN j k, FromListsN e k k, FromListsN e j k, FromListsN e k j, FromListsN e j j) => Matrix e n k -> Matrix e m j -> Matrix e (Either n m) (Either k j)
- (><) :: forall e m p n q. (Num e, CountableDimensionsN m n, CountableDimensionsN p q, CountableDimensionsN (m, n) (p, q), FromListsN e (m, n) m, FromListsN e (m, n) n, FromListsN e (p, q) p, FromListsN e (p, q) q, TrivialP m n, TrivialP p q) => Matrix e m p -> Matrix e n q -> Matrix e (m, n) (p, q)
- fstM :: forall e m k. (Num e, CountableDimensionsN m k, CountableN (m, k), FromListsN e (m, k) m, TrivialP m k) => Matrix e (m, k) m
- sndM :: forall e m k. (Num e, CountableDimensionsN k m, CountableN (m, k), FromListsN e (m, k) k, TrivialP m k) => Matrix e (m, k) k
- kr :: forall e cols a b. (Num e, CountableDimensionsN a b, CountableN (a, b), FromListsN e (a, b) a, FromListsN e (a, b) b, TrivialP a b) => Matrix e cols a -> Matrix e cols b -> Matrix e cols (a, b)
- iden :: (Num e, FromListsN e a a, CountableN a) => Matrix e a a
- comp :: Num e => Matrix e cr rows -> Matrix e cols cr -> Matrix e cols rows
- fromF' :: (Liftable e a b, CountableDimensionsN cols rows, FromListsN e rows cols) => (a -> b) -> Matrix e cols rows
- fromF :: (Liftable e a b, CountableDimensionsN a b, FromListsN e b a) => (a -> b) -> Matrix e a b
- toRel :: (Liftable (Natural 0 1) a b, CountableDimensionsN a b, FromListsN (Natural 0 1) b a) => (a -> b -> Bool) -> Matrix (Natural 0 1) a b
- pretty :: (CountableDimensionsN cols rows, Show e) => Matrix e cols rows -> String
- prettyPrint :: (CountableDimensionsN cols rows, Show e) => Matrix e cols rows -> IO ()
Documentation
LAoP (Linear Algebra of Programming) Inductive Matrix definition.
LAoP generalises relations and functions treating them as Boolean matrices and in turn consider these as arrows. This library offers many of the combinators mentioned in the work of Macedo (2012) and Oliveira (2012).
This definition is a wrapper around Type
but
dimensions are arbitrary data types. Type inference might not
be as desired.
Type safe matrix representation
newtype Matrix e (cols :: Type) (rows :: Type) Source #
Instances
Num e => Category (Matrix e) Source # | It is possible to implement a constrained version of the category type class. |
Eq e => Eq (Matrix e cols rows) Source # | |
Num e => Num (Matrix e cols rows) Source # | |
Defined in LAoP.Matrix.Type (+) :: Matrix e cols rows -> Matrix e cols rows -> Matrix e cols rows # (-) :: Matrix e cols rows -> Matrix e cols rows -> Matrix e cols rows # (*) :: Matrix e cols rows -> Matrix e cols rows -> Matrix e cols rows # negate :: Matrix e cols rows -> Matrix e cols rows # abs :: Matrix e cols rows -> Matrix e cols rows # signum :: Matrix e cols rows -> Matrix e cols rows # fromInteger :: Integer -> Matrix e cols rows # | |
Ord e => Ord (Matrix e cols rows) Source # | |
Defined in LAoP.Matrix.Type compare :: Matrix e cols rows -> Matrix e cols rows -> Ordering # (<) :: Matrix e cols rows -> Matrix e cols rows -> Bool # (<=) :: Matrix e cols rows -> Matrix e cols rows -> Bool # (>) :: Matrix e cols rows -> Matrix e cols rows -> Bool # (>=) :: Matrix e cols rows -> Matrix e cols rows -> Bool # max :: Matrix e cols rows -> Matrix e cols rows -> Matrix e cols rows # min :: Matrix e cols rows -> Matrix e cols rows -> Matrix e cols rows # | |
Show e => Show (Matrix e cols rows) Source # | |
NFData e => NFData (Matrix e cols rows) Source # | |
Defined in LAoP.Matrix.Type | |
type Object (Matrix e) a Source # | |
Defined in LAoP.Matrix.Type |
Constraint type aliases
type Countable a = KnownNat (Count a) Source #
Constraint type synonyms to keep the type signatures less convoluted
type CountableDimensions a b = (Countable a, Countable b) Source #
type CountableDimensionsN a b = (CountableN a, CountableN b) Source #
Type aliases
Primitives
join :: Matrix e a rows -> Matrix e b rows -> Matrix e (Either a b) rows Source #
Matrix Join
constructor
fork :: Matrix e cols a -> Matrix e cols b -> Matrix e cols (Either a b) Source #
Matrix Fork
constructor
Auxiliary type families
type family FromNat (n :: Nat) :: Type where ... Source #
Type family that computes of a given type dimension from a given natural
Thanks to Li-Yao Xia this type family is super fast.
type family Count (d :: Type) :: Nat where ... Source #
Type family that computes the cardinality of a given type dimension.
It can also count the cardinality of custom types that implement the
Generic
instance.
Count (Natural n m) = (m - n) + 1 | |
Count (List a) = (^) 2 (Count a) | |
Count (Either a b) = (+) (Count a) (Count b) | |
Count (a, b) = * (Count a) (Count b) | |
Count (a -> b) = (^) (Count b) (Count a) | |
Count (M1 _ _ f p) = Count (f p) | |
Count (K1 _ _ _) = 1 | |
Count (V1 _) = 0 | |
Count (U1 _) = 1 | |
Count ((:*:) a b p) = Count (a p) * Count (b p) | |
Count ((:+:) a b p) = Count (a p) + Count (b p) | |
Count d = Count (Rep d R) |
type family Normalize (d :: Type) :: Type where ... Source #
Type family that normalizes the representation of a given data structure
Matrix construction and conversion
class FromLists e cols rows Source #
Type class for defining the fromList
conversion function.
Given that it is not possible to branch on types at the term level type classes are needed very much like an inductive definition but on types.
Instances
FromLists e () () Source # | |
Defined in LAoP.Matrix.Internal | |
FromLists e Void Void Source # | |
(FromLists e () a, FromLists e () b, Countable a) => FromLists e () (Either a b) Source # | |
FromLists e () rows => FromLists e () (Either () rows) Source # | |
(FromLists e a (), FromLists e b (), Countable a) => FromLists e (Either a b) () Source # | |
FromLists e cols () => FromLists e (Either () cols) () Source # | |
(FromLists e (Either a b) c, FromLists e (Either a b) d, Countable c) => FromLists e (Either a b) (Either c d) Source # | |
fromLists :: FromListsN e cols rows => [[e]] -> Matrix e cols rows Source #
Build a matrix out of a list of list of elements. Throws a runtime error if the dimensions do not match.
matrixBuilder' :: (FromListsN e cols rows, CountableDimensionsN cols rows) => ((Int, Int) -> e) -> Matrix e cols rows Source #
Matrix builder function. Constructs a matrix provided with a construction function that operates with indices.
matrixBuilder :: (FromListsN e a b, Enum a, Enum b, Bounded a, Bounded b, Eq a, CountableDimensionsN a b) => ((a, b) -> e) -> Matrix e a b Source #
Matrix builder function. Constructs a matrix provided with a construction function that operates with arbitrary types.
zeros :: (Num e, FromListsN e cols rows, CountableDimensionsN cols rows) => Matrix e cols rows Source #
The zero matrix. A matrix wholly filled with zeros.
ones :: (Num e, FromListsN e cols rows, CountableDimensionsN cols rows) => Matrix e cols rows Source #
The ones matrix. A matrix wholly filled with ones.
Also known as T (Top) matrix.
bang :: forall e cols. (Num e, Enum e, FromListsN e cols (), CountableN cols) => Matrix e cols One Source #
The T (Top) row vector matrix.
point :: (Bounded a, Enum a, Eq a, Num e, Ord e, CountableN a, FromListsN e a One) => a -> Matrix e One a Source #
Point constant relation
constant :: (Num e, FromListsN e cols rows, CountableDimensionsN cols rows) => e -> Matrix e cols rows Source #
The constant matrix constructor. A matrix wholly filled with a given value.
Functor instance equivalent function
fmapM :: (Liftable e a b, CountableDimensionsN a b, FromListsN e b a) => (a -> b) -> Matrix e c a -> Matrix e c b Source #
Functor instance equivalent function
bimapM :: (Liftable e a b, Liftable e c d, CountableDimensionsN a c, CountableDimensionsN b d, FromListsN e d c, FromListsN e b a) => (a -> b) -> (c -> d) -> Matrix e a c -> Matrix e b d Source #
Bifunctor equivalent function
Applicative/Monoidal instance equivalent functions
multM :: (CountableDimensionsN a b, CountableN (a, b), Num e, FromListsN e (a, b) a, FromListsN e (a, b) b, TrivialP a b) => Matrix e c a -> Matrix e c b -> Matrix e c (a, b) Source #
Applicative instance equivalent unit
function,
Selective equivalent instance function
selectM :: (Num e, FromListsN e b b, CountableN b) => Matrix e cols (Either a b) -> Matrix e a b -> Matrix e cols b Source #
Selective functors select
operator equivalent inspired by the
ArrowMonad solution presented in the paper.
Monad equivalent instance function
returnM :: forall e a. (Num e, Enum e, Enum a, FromListsN e () a, Countable a) => a -> Matrix e One a Source #
Monad instance equivalent return
function,
bindM :: Num e => Matrix e a b -> Matrix e b c -> Matrix e a c Source #
Monad instance equivalent '(>>=)' function,
Misc
Get dimensions
Matrix Transposition
Scalar multiplication/division of matrices
(.|) :: Num e => e -> Matrix e cols rows -> Matrix e cols rows infixl 7 Source #
Scalar multiplication of matrices.
(./) :: Fractional e => Matrix e cols rows -> e -> Matrix e cols rows infixl 7 Source #
Scalar multiplication of matrices.
McCarthy's Conditional
cond :: (Trivial a, Trivial2 a, Trivial3 a, CountableN a, FromListsN e () a, FromListsN e a (), FromListsN e a a, Liftable e a Bool) => (a -> Bool) -> Matrix e a b -> Matrix e a b -> Matrix e a b Source #
McCarthy's Conditional expresses probabilistic choice.
Matrix "abiding"
abideJF :: Matrix e cols rows -> Matrix e cols rows Source #
Matrix "abiding" followin the 'Join'-'Fork' abide law.
Law:
Join
(Fork
a c) (Fork
b d) ==Fork
(Join
a b) (Join
c d)
abideFJ :: Matrix e cols rows -> Matrix e cols rows Source #
Matrix "abiding" followin the 'Fork'-'Join' abide law.
Law:
Fork
(Join
a b) (Join
c d) ==Join
(Fork
a c) (Fork
b d)
Zip Matrices
zipWithM :: (e -> f -> g) -> Matrix e a b -> Matrix f a b -> Matrix g a b Source #
Zip two matrices with a given binary function
Biproduct approach
Fork
(===) :: Matrix e cols a -> Matrix e cols b -> Matrix e cols (Either a b) infixl 2 Source #
Matrix Fork
constructor
Projections
p1 :: (Num e, CountableDimensionsN n m, FromListsN e n m, FromListsN e m m) => Matrix e (Either m n) m Source #
Biproduct first component projection
p2 :: (Num e, CountableDimensionsN n m, FromListsN e m n, FromListsN e n n) => Matrix e (Either m n) n Source #
Biproduct second component projection
Join
(|||) :: Matrix e a rows -> Matrix e b rows -> Matrix e (Either a b) rows infixl 3 Source #
Matrix Join
constructor
Injections
i1 :: (Num e, CountableDimensionsN n m, FromListsN e n m, FromListsN e m m) => Matrix e m (Either m n) Source #
Biproduct first component injection
i2 :: (Num e, CountableDimensionsN n m, FromListsN e m n, FromListsN e n n) => Matrix e n (Either m n) Source #
Biproduct second component injection
Bifunctors
(-|-) :: (Num e, CountableDimensionsN j k, FromListsN e k k, FromListsN e j k, FromListsN e k j, FromListsN e j j) => Matrix e n k -> Matrix e m j -> Matrix e (Either n m) (Either k j) infixl 5 Source #
Matrix coproduct functor also known as matrix direct sum.
(><) :: forall e m p n q. (Num e, CountableDimensionsN m n, CountableDimensionsN p q, CountableDimensionsN (m, n) (p, q), FromListsN e (m, n) m, FromListsN e (m, n) n, FromListsN e (p, q) p, FromListsN e (p, q) q, TrivialP m n, TrivialP p q) => Matrix e m p -> Matrix e n q -> Matrix e (m, n) (p, q) infixl 4 Source #
Matrix product functor also known as Kronecker product
Applicative matrix combinators
Note that given the restrictions imposed it is not possible to implement the standard type classes present in standard Haskell. *** Matrix pairing projections
fstM :: forall e m k. (Num e, CountableDimensionsN m k, CountableN (m, k), FromListsN e (m, k) m, TrivialP m k) => Matrix e (m, k) m Source #
Khatri Rao product first component projection matrix.
sndM :: forall e m k. (Num e, CountableDimensionsN k m, CountableN (m, k), FromListsN e (m, k) k, TrivialP m k) => Matrix e (m, k) k Source #
Khatri Rao product second component projection matrix.
Matrix pairing
kr :: forall e cols a b. (Num e, CountableDimensionsN a b, CountableN (a, b), FromListsN e (a, b) a, FromListsN e (a, b) b, TrivialP a b) => Matrix e cols a -> Matrix e cols b -> Matrix e cols (a, b) Source #
Khatri Rao Matrix product also known as matrix pairing.
NOTE: That this is not a true categorical product, see for instance:
| fstM . kr a b == a kr a b ==> | | sndM . kr a b == b
Emphasis on the implication symbol.
Matrix composition and lifting
Arrow matrix combinators
Note that given the restrictions imposed it is not possible to implement the standard type classes present in standard Haskell.
iden :: (Num e, FromListsN e a a, CountableN a) => Matrix e a a Source #
iden matrix
comp :: Num e => Matrix e cr rows -> Matrix e cols cr -> Matrix e cols rows Source #
Matrix composition. Equivalent to matrix-matrix multiplication.
This definition takes advantage of divide-and-conquer and fusion laws from LAoP.
fromF' :: (Liftable e a b, CountableDimensionsN cols rows, FromListsN e rows cols) => (a -> b) -> Matrix e cols rows Source #
Lifts functions to matrices with arbitrary dimensions.
NOTE: Be careful to not ask for a matrix bigger than the cardinality of
types a
or b
allows.
fromF :: (Liftable e a b, CountableDimensionsN a b, FromListsN e b a) => (a -> b) -> Matrix e a b Source #
Lifts functions to matrices with dimensions matching a
and b
cardinality's.
Relation
toRel :: (Liftable (Natural 0 1) a b, CountableDimensionsN a b, FromListsN (Natural 0 1) b a) => (a -> b -> Bool) -> Matrix (Natural 0 1) a b Source #
Lifts relation functions to Boolean Matrix
Matrix printing
pretty :: (CountableDimensionsN cols rows, Show e) => Matrix e cols rows -> String Source #
Matrix pretty printer
prettyPrint :: (CountableDimensionsN cols rows, Show e) => Matrix e cols rows -> IO () Source #
Matrix pretty printer