adapting lapack wrappers

1 files changed, 14 insertions, 0 deletions

diff --git a/lib/Data/Packed/Internal/Matrix.hs b/lib/Data/Packed/Internal/Matrix.hs
index ec6657a..db53cd1 100644
--- a/lib/Data/Packed/Internal/Matrix.hs
+++ b/lib/Data/Packed/Internal/Matrix.hs
@@ -160,6 +160,10 @@ toRows m = toRows' 0 where
     toRows' k | k == r*c  = []
               | otherwise = subVector k c v : toRows' (k+c)
 
+-----------------------------------------------------------------
+
+liftMatrix f m = m { dat = f dat, tdat = f tdat } -- check sizes
+
 ------------------------------------------------------------------
 
 dotL a b = sum (zipWith (*) a b)
@@ -208,6 +212,16 @@ multiplyD order a b
 
 ----------------------------------------------------------------------
 
+outer u v = dat (multiply RowMajor r c)
+    where r = matrixFromVector RowMajor 1 u
+          c = matrixFromVector RowMajor (dim v) v
+
+dot u v = dat (multiply RowMajor r c) `at` 0
+    where r = matrixFromVector RowMajor (dim u) u
+          c = matrixFromVector RowMajor 1 v
+
+----------------------------------------------------------------------
+
 -- | extraction of a submatrix of a real matrix
 subMatrixR :: (Int,Int) -- ^ (r0,c0) starting position 
            -> (Int,Int) -- ^ (rt,ct) dimensions of submatrix