3 files changed, 80 insertions, 10 deletions

diff --git a/examples/data.txt b/examples/data.txt
index 7031462..2b9bfea 100644
--- a/examples/data.txt
+++ b/examples/data.txt
@@ -1,10 +1,10 @@
- 1    1.1
- 2    3.9
- 3    9.2
- 4   15.8
- 5   25.3
- 6   35.7
- 7   49.4
- 8   63.6
- 9   81.5
-10   99.5
\ No newline at end of file
+ 0.9    1.1
+ 2.1    3.9
+ 3.1    9.2
+ 4.0   51.8
+ 4.9   25.3
+ 6.1   35.7
+ 7.0   49.4
+ 7.9    3.6
+ 9.1   81.5
+10.2   99.5
\ No newline at end of file
diff --git a/examples/speed.hs b/examples/speed.hs
new file mode 100644
index 0000000..a937f31
--- /dev/null
+++ b/examples/speed.hs
@@ -0,0 +1,50 @@
+import Numeric.LinearAlgebra
+import System
+import Data.List(foldl1')
+import System.CPUTime
+import Text.Printf
+import Debug.Trace
+
+debug x = trace (show x) x
+
+timing act = do
+    t0 <- getCPUTime
+    act
+    t1 <- getCPUTime
+    printf "%.2f CPU seconds\n" $ (fromIntegral ((t1 - t0) `div` (10^10)) / 100 :: Double)
+
+op a b = trans $ (trans a) <> (trans b)
+
+op2 a b = trans $ (trans a) + (trans b)
+
+rot' :: Double -> Matrix Double
+rot' a = matrix [[ c,0,s],
+                 [ 0,1,0],
+                 [-s,0,c]]
+    where c = cos a
+          s = sin a
+          matrix = fromLists
+
+rot :: Double -> Matrix Double
+rot a = (3><3) [ c,0,s
+               , 0,1,0
+               ,-s,0,c ]
+    where c = cos a
+          s = sin a
+
+
+fun n r = foldl1' (<>) (map r [0,delta..1]) where delta = 1 /(fromIntegral n-1)
+
+
+main = do
+    args <- getArgs
+    let [p,n,d] = map read args
+    let ms = replicate n ((ident d :: Matrix Double))
+    let mz = replicate n (diag (constant (0::Double) d))
+    timing $ case p of
+        0 -> print $ foldl1' (<>) ms
+        1 -> print $ foldl1' (<>) (map trans ms)
+        2 -> print $ foldl1' op ms
+        3 -> print $ foldl1' op2 mz
+        4 -> print $ fun n rot'
+        5 -> print $ fun n rot
diff --git a/examples/speed.m b/examples/speed.m
new file mode 100644
index 0000000..10dbf78
--- /dev/null
+++ b/examples/speed.m
@@ -0,0 +1,20 @@
+#! /usr/bin/octave -qf
+1; % measuring Octave computing times
+
+function r = rot(a)
+    c = cos(a);
+    s = sin(a);
+    r = [ c , 0, s;
+          0,  1, 0;
+         -s,  0, c];
+end
+
+t0=time();
+x = linspace(0,1,1E5);
+ac = eye(3);
+for a=x
+    ac = rot(a)*ac;
+end
+
+disp(ac);
+disp(time()-t0)