From 69b0f9c46f09dc9fc5c25e706efb8a7da51d2480 Mon Sep 17 00:00:00 2001 From: Alberto Ruiz Date: Wed, 14 Mar 2012 12:27:09 +0100 Subject: new module Numeric.LinearAlgebra.Real --- examples/Real.hs | 113 ------------------------------------------------------- 1 file changed, 113 deletions(-) delete mode 100644 examples/Real.hs (limited to 'examples') diff --git a/examples/Real.hs b/examples/Real.hs deleted file mode 100644 index 9083b87..0000000 --- a/examples/Real.hs +++ /dev/null @@ -1,113 +0,0 @@ - --- Alternative interface and utilities for creation of real arrays, useful to work in interactive mode. - -module Real( - module Numeric.LinearAlgebra, - (<>), (*>), (<*), (<\>), (\>), - vector, - eye, - zeros, ones, - diagl, - row, - col, - (#),(&), (//), blocks, - rand, randn -) where - -import Numeric.LinearAlgebra hiding ((<>), (<\>)) -import System.Random(randomIO) - -infixl 7 <> --- | Matrix product ('multiply') -(<>) :: Field t => Matrix t -> Matrix t -> Matrix t -(<>) = multiply - -infixl 7 *> --- | matrix x vector -(*>) :: Field t => Matrix t -> Vector t -> Vector t -m *> v = flatten $ m <> (asColumn v) - -infixl 7 <* --- | vector x matrix -(<*) :: Field t => Vector t -> Matrix t -> Vector t -v <* m = flatten $ (asRow v) <> m - - --- | Least squares solution of a linear system for several right-hand sides, similar to the \\ operator of Matlab\/Octave. (\<\\\>) = 'linearSolveSVD'. -(<\>) :: (Field a) => Matrix a -> Matrix a -> Matrix a -infixl 7 <\> -(<\>) = linearSolveSVD - --- | Least squares solution of a linear system for a single right-hand side. See '(\<\\\>)'. -(\>) :: (Field a) => Matrix a -> Vector a -> Vector a -infixl 7 \> -m \> v = flatten (m <\> reshape 1 v) - --- | Pseudorandom matrix with uniform elements between 0 and 1. -randm :: RandDist - -> Int -- ^ rows - -> Int -- ^ columns - -> IO (Matrix Double) -randm d r c = do - seed <- randomIO - return (reshape c $ randomVector seed d (r*c)) - --- | Pseudorandom matrix with uniform elements between 0 and 1. -rand :: Int -> Int -> IO (Matrix Double) -rand = randm Uniform - --- | Pseudorandom matrix with normal elements -randn :: Int -> Int -> IO (Matrix Double) -randn = randm Gaussian - --- | Real identity matrix. -eye :: Int -> Matrix Double -eye = ident - --- | Create a real vector from a list. -vector :: [Double] -> Vector Double -vector = fromList - --- | Create a real diagonal matrix from a list. -diagl :: [Double] -> Matrix Double -diagl = diag . vector - --- | Create a matrix or zeros. -zeros :: Int -- ^ rows - -> Int -- ^ columns - -> Matrix Double -zeros r c = konst 0 (r,c) - --- | Create a matrix or ones. -ones :: Int -- ^ rows - -> Int -- ^ columns - -> Matrix Double -ones r c = konst 1 (r,c) - --- | Concatenation of real vectors. -infixl 9 # -(#) :: Vector Double -> Vector Double -> Vector Double -a # b = join [a,b] - --- | Horizontal concatenation of real matrices. -infixl 8 & -(&) :: Matrix Double -> Matrix Double -> Matrix Double -a & b = fromBlocks [[a,b]] - --- | Vertical concatenation of real matrices. -infixl 7 // -(//) :: Matrix Double -> Matrix Double -> Matrix Double -a // b = fromBlocks [[a],[b]] - --- | Real block matrix from a rectangular list of lists. -blocks :: [[Matrix Double]] -> Matrix Double -blocks = fromBlocks - --- | A real matrix with a single row, create from a list of elements. -row :: [Double] -> Matrix Double -row = asRow . vector - --- | A real matrix with a single column, created from a list of elements. -col :: [Double] -> Matrix Double -col = asColumn . vector - -- cgit v1.2.3