diff options
author | Alberto Ruiz <aruiz@um.es> | 2010-02-07 09:29:50 +0000 |
---|---|---|
committer | Alberto Ruiz <aruiz@um.es> | 2010-02-07 09:29:50 +0000 |
commit | aef0333b5180ea79e539bd53194f1dfed20b7db5 (patch) | |
tree | 6ece2434ecacab194331120bd47d09ab04ae0f4f /examples/HMatrixReal.hs | |
parent | 634683fcfab73a0bd830fef03fb9f4603ba837b6 (diff) |
added odeSolve
Diffstat (limited to 'examples/HMatrixReal.hs')
-rw-r--r-- | examples/HMatrixReal.hs | 104 |
1 files changed, 0 insertions, 104 deletions
diff --git a/examples/HMatrixReal.hs b/examples/HMatrixReal.hs deleted file mode 100644 index 0edff10..0000000 --- a/examples/HMatrixReal.hs +++ /dev/null | |||
@@ -1,104 +0,0 @@ | |||
1 | |||
2 | -- Alternative interface and utilities for creation of real arrays, useful to work in interactive mode. | ||
3 | |||
4 | module HMatrixReal( | ||
5 | module Numeric.LinearAlgebra, | ||
6 | (<>), (*>), (<*), (<\>), (\>), | ||
7 | vector, | ||
8 | eye, | ||
9 | zeros, ones, | ||
10 | diagl, | ||
11 | row, | ||
12 | col, | ||
13 | (#),(&), (//), blocks, | ||
14 | rand | ||
15 | ) where | ||
16 | |||
17 | import Numeric.LinearAlgebra hiding ((<>), (<|>), (<->), (<\>), (.*), (*/)) | ||
18 | import System.Random(randomIO) | ||
19 | |||
20 | infixl 7 <> | ||
21 | -- | Matrix product ('multiply') | ||
22 | (<>) :: Field t => Matrix t -> Matrix t -> Matrix t | ||
23 | (<>) = multiply | ||
24 | |||
25 | infixl 7 *> | ||
26 | -- | matrix x vector | ||
27 | (*>) :: Field t => Matrix t -> Vector t -> Vector t | ||
28 | m *> v = flatten $ m <> (asColumn v) | ||
29 | |||
30 | infixl 7 <* | ||
31 | -- | vector x matrix | ||
32 | (<*) :: Field t => Vector t -> Matrix t -> Vector t | ||
33 | v <* m = flatten $ (asRow v) <> m | ||
34 | |||
35 | |||
36 | -- | Least squares solution of a linear system for several right-hand sides, similar to the \\ operator of Matlab\/Octave. (\<\\\>) = 'linearSolveSVD'. | ||
37 | (<\>) :: (Field a) => Matrix a -> Matrix a -> Matrix a | ||
38 | infixl 7 <\> | ||
39 | (<\>) = linearSolveSVD | ||
40 | |||
41 | -- | Least squares solution of a linear system for a single right-hand side. See '(\<\\\>)'. | ||
42 | (\>) :: (Field a) => Matrix a -> Vector a -> Vector a | ||
43 | infixl 7 \> | ||
44 | m \> v = flatten (m <\> reshape 1 v) | ||
45 | |||
46 | -- | Pseudorandom matrix with uniform elements between 0 and 1. | ||
47 | rand :: Int -- ^ rows | ||
48 | -> Int -- ^ columns | ||
49 | -> IO (Matrix Double) | ||
50 | rand r c = do | ||
51 | seed <- randomIO | ||
52 | return (reshape c $ randomVector seed Uniform (r*c)) | ||
53 | |||
54 | -- | Real identity matrix. | ||
55 | eye :: Int -> Matrix Double | ||
56 | eye = ident | ||
57 | |||
58 | -- | Create a real vector from a list. | ||
59 | vector :: [Double] -> Vector Double | ||
60 | vector = fromList | ||
61 | |||
62 | -- | Create a real diagonal matrix from a list. | ||
63 | diagl :: [Double] -> Matrix Double | ||
64 | diagl = diag . vector | ||
65 | |||
66 | -- | Create a matrix or zeros. | ||
67 | zeros :: Int -- ^ rows | ||
68 | -> Int -- ^ columns | ||
69 | -> Matrix Double | ||
70 | zeros r c = reshape c (constant 0 (r*c)) | ||
71 | |||
72 | -- | Create a matrix or ones. | ||
73 | ones :: Int -- ^ rows | ||
74 | -> Int -- ^ columns | ||
75 | -> Matrix Double | ||
76 | ones r c = reshape c (constant 1 (r*c)) | ||
77 | |||
78 | -- | Concatenation of real vectors. | ||
79 | infixl 9 # | ||
80 | (#) :: Vector Double -> Vector Double -> Vector Double | ||
81 | a # b = join [a,b] | ||
82 | |||
83 | -- | Horizontal concatenation of real matrices. | ||
84 | infixl 8 & | ||
85 | (&) :: Matrix Double -> Matrix Double -> Matrix Double | ||
86 | a & b = fromBlocks [[a,b]] | ||
87 | |||
88 | -- | Vertical concatenation of real matrices. | ||
89 | infixl 7 // | ||
90 | (//) :: Matrix Double -> Matrix Double -> Matrix Double | ||
91 | a // b = fromBlocks [[a],[b]] | ||
92 | |||
93 | -- | Real block matrix from a rectangular list of lists. | ||
94 | blocks :: [[Matrix Double]] -> Matrix Double | ||
95 | blocks = fromBlocks | ||
96 | |||
97 | -- | A real matrix with a single row, create from a list of elements. | ||
98 | row :: [Double] -> Matrix Double | ||
99 | row = asRow . vector | ||
100 | |||
101 | -- | A real matrix with a single column, created from a list of elements. | ||
102 | col :: [Double] -> Matrix Double | ||
103 | col = asColumn . vector | ||
104 | |||