initial interface to MKL sparse solver

1 files changed, 32 insertions, 0 deletions

diff --git a/packages/sparse/src/Numeric/LinearAlgebra/Sparse.hs b/packages/sparse/src/Numeric/LinearAlgebra/Sparse.hs
new file mode 100644
index 0000000..ccf28b7
--- /dev/null
+++ b/packages/sparse/src/Numeric/LinearAlgebra/Sparse.hs
@@ -0,0 +1,32 @@
+{-# LANGUAGE ForeignFunctionInterface #-}
+{-# LANGUAGE RecordWildCards #-}
+
+
+
+module Numeric.LinearAlgebra.Sparse (
+    dss
+) where
+
+import Foreign.C.Types(CInt(..))
+import Data.Packed.Development
+import System.IO.Unsafe(unsafePerformIO)
+import Foreign(Ptr)
+import Numeric.HMatrix
+import Text.Printf
+import Numeric.LinearAlgebra.Util((~!~))
+
+
+type IV t = CInt -> Ptr CInt   -> t
+type  V t = CInt -> Ptr Double -> t
+type SMxV = V (IV (IV (V (V (IO CInt)))))
+
+dss :: CSR -> Vector Double -> Vector Double
+dss CSR{..} b = unsafePerformIO $ do
+    size b /= csrNRows ~!~ printf "dss: incorrect sizes: (%d,%d) x %d" csrNRows csrNCols (size b)
+    r <- createVector csrNCols
+    app5 c_dss vec csrVals vec csrCols vec csrRows vec b vec r "dss"
+    return r
+
+foreign import ccall unsafe "dss"
+  c_dss :: SMxV
+