summaryrefslogtreecommitdiff
path: root/lib/Numeric
diff options
context:
space:
mode:
Diffstat (limited to 'lib/Numeric')
-rw-r--r--lib/Numeric/GSL.hs2
-rw-r--r--lib/Numeric/GSL/Differentiation.hs4
-rw-r--r--lib/Numeric/GSL/Fitting.hs2
-rw-r--r--lib/Numeric/GSL/Fourier.hs2
-rw-r--r--lib/Numeric/GSL/Integration.hs12
-rw-r--r--lib/Numeric/GSL/Internal.hs10
-rw-r--r--lib/Numeric/GSL/Minimization.hs4
-rw-r--r--lib/Numeric/GSL/ODE.hs2
-rw-r--r--lib/Numeric/GSL/Polynomials.hs2
-rw-r--r--lib/Numeric/GSL/Root.hs4
-rw-r--r--lib/Numeric/GSL/Vector.hs58
-rw-r--r--lib/Numeric/LinearAlgebra/LAPACK.hs64
12 files changed, 83 insertions, 83 deletions
diff --git a/lib/Numeric/GSL.hs b/lib/Numeric/GSL.hs
index e53022e..5f39a3e 100644
--- a/lib/Numeric/GSL.hs
+++ b/lib/Numeric/GSL.hs
@@ -40,4 +40,4 @@ import Data.Complex
40 40
41-- | This action removes the GSL default error handler (which aborts the program), so that 41-- | This action removes the GSL default error handler (which aborts the program), so that
42-- GSL errors can be handled by Haskell (using Control.Exception) and ghci doesn't abort. 42-- GSL errors can be handled by Haskell (using Control.Exception) and ghci doesn't abort.
43foreign import ccall "GSL/gsl-aux.h no_abort_on_error" setErrorHandlerOff :: IO () 43foreign import ccall unsafe "GSL/gsl-aux.h no_abort_on_error" setErrorHandlerOff :: IO ()
diff --git a/lib/Numeric/GSL/Differentiation.hs b/lib/Numeric/GSL/Differentiation.hs
index ab5eb50..d2a332c 100644
--- a/lib/Numeric/GSL/Differentiation.hs
+++ b/lib/Numeric/GSL/Differentiation.hs
@@ -48,7 +48,7 @@ derivGen c h f x = unsafePerformIO $ do
48 freeHaskellFunPtr fp 48 freeHaskellFunPtr fp
49 return result 49 return result
50 50
51foreign import ccall "gsl-aux.h deriv" 51foreign import ccall safe "gsl-aux.h deriv"
52 c_deriv :: CInt -> FunPtr (Double -> Ptr () -> Double) -> Double -> Double 52 c_deriv :: CInt -> FunPtr (Double -> Ptr () -> Double) -> Double -> Double
53 -> Ptr Double -> Ptr Double -> IO CInt 53 -> Ptr Double -> Ptr Double -> IO CInt
54 54
@@ -84,4 +84,4 @@ derivBackward = derivGen 2
84 84
85{- | conversion of Haskell functions into function pointers that can be used in the C side 85{- | conversion of Haskell functions into function pointers that can be used in the C side
86-} 86-}
87foreign import ccall "wrapper" mkfun:: (Double -> Ptr() -> Double) -> IO( FunPtr (Double -> Ptr() -> Double)) 87foreign import ccall safe "wrapper" mkfun:: (Double -> Ptr() -> Double) -> IO( FunPtr (Double -> Ptr() -> Double))
diff --git a/lib/Numeric/GSL/Fitting.hs b/lib/Numeric/GSL/Fitting.hs
index 79db23f..6343b76 100644
--- a/lib/Numeric/GSL/Fitting.hs
+++ b/lib/Numeric/GSL/Fitting.hs
@@ -90,7 +90,7 @@ nlFitGen m f jac xiv epsabs epsrel maxit = unsafePerformIO $ do
90 freeHaskellFunPtr jp 90 freeHaskellFunPtr jp
91 return (subVector 2 p sol, path) 91 return (subVector 2 p sol, path)
92 92
93foreign import ccall "nlfit" 93foreign import ccall safe "nlfit"
94 c_nlfit:: CInt -> FunPtr TVV -> FunPtr TVM -> Double -> Double -> CInt -> CInt -> TVM 94 c_nlfit:: CInt -> FunPtr TVV -> FunPtr TVM -> Double -> Double -> CInt -> CInt -> TVM
95 95
96------------------------------------------------------- 96-------------------------------------------------------
diff --git a/lib/Numeric/GSL/Fourier.hs b/lib/Numeric/GSL/Fourier.hs
index 71bc8a1..4ef19b3 100644
--- a/lib/Numeric/GSL/Fourier.hs
+++ b/lib/Numeric/GSL/Fourier.hs
@@ -30,7 +30,7 @@ genfft code v = unsafePerformIO $ do
30 app2 (c_fft code) vec v vec r "fft" 30 app2 (c_fft code) vec v vec r "fft"
31 return r 31 return r
32 32
33foreign import ccall "gsl-aux.h fft" c_fft :: CInt -> TCVCV 33foreign import ccall unsafe "gsl-aux.h fft" c_fft :: CInt -> TCVCV
34 34
35 35
36{- | Fast 1D Fourier transform of a 'Vector' @(@'Complex' 'Double'@)@ using /gsl_fft_complex_forward/. It uses the same scaling conventions as GNU Octave. 36{- | Fast 1D Fourier transform of a 'Vector' @(@'Complex' 'Double'@)@ using /gsl_fft_complex_forward/. It uses the same scaling conventions as GNU Octave.
diff --git a/lib/Numeric/GSL/Integration.hs b/lib/Numeric/GSL/Integration.hs
index 2330fc6..defbf80 100644
--- a/lib/Numeric/GSL/Integration.hs
+++ b/lib/Numeric/GSL/Integration.hs
@@ -32,7 +32,7 @@ import System.IO.Unsafe(unsafePerformIO)
32 32
33{- | conversion of Haskell functions into function pointers that can be used in the C side 33{- | conversion of Haskell functions into function pointers that can be used in the C side
34-} 34-}
35foreign import ccall "wrapper" mkfun:: (Double -> Ptr() -> Double) -> IO( FunPtr (Double -> Ptr() -> Double)) 35foreign import ccall safe "wrapper" mkfun:: (Double -> Ptr() -> Double) -> IO( FunPtr (Double -> Ptr() -> Double))
36 36
37-------------------------------------------------------------------- 37--------------------------------------------------------------------
38{- | Numerical integration using /gsl_integration_qags/ (adaptive integration with singularities). For example: 38{- | Numerical integration using /gsl_integration_qags/ (adaptive integration with singularities). For example:
@@ -63,7 +63,7 @@ integrateQAGS prec n f a b = unsafePerformIO $ do
63 freeHaskellFunPtr fp 63 freeHaskellFunPtr fp
64 return result 64 return result
65 65
66foreign import ccall "gsl-aux.h integrate_qags" 66foreign import ccall safe "gsl-aux.h integrate_qags"
67 c_integrate_qags :: FunPtr (Double-> Ptr() -> Double) -> Double -> Double -> Double -> CInt 67 c_integrate_qags :: FunPtr (Double-> Ptr() -> Double) -> Double -> Double -> Double -> CInt
68 -> Ptr Double -> Ptr Double -> IO CInt 68 -> Ptr Double -> Ptr Double -> IO CInt
69 69
@@ -94,7 +94,7 @@ integrateQNG prec f a b = unsafePerformIO $ do
94 freeHaskellFunPtr fp 94 freeHaskellFunPtr fp
95 return result 95 return result
96 96
97foreign import ccall "gsl-aux.h integrate_qng" 97foreign import ccall safe "gsl-aux.h integrate_qng"
98 c_integrate_qng :: FunPtr (Double-> Ptr() -> Double) -> Double -> Double -> Double 98 c_integrate_qng :: FunPtr (Double-> Ptr() -> Double) -> Double -> Double -> Double
99 -> Ptr Double -> Ptr Double -> IO CInt 99 -> Ptr Double -> Ptr Double -> IO CInt
100 100
@@ -126,7 +126,7 @@ integrateQAGI prec n f = unsafePerformIO $ do
126 freeHaskellFunPtr fp 126 freeHaskellFunPtr fp
127 return result 127 return result
128 128
129foreign import ccall "gsl-aux.h integrate_qagi" 129foreign import ccall safe "gsl-aux.h integrate_qagi"
130 c_integrate_qagi :: FunPtr (Double-> Ptr() -> Double) -> Double -> CInt 130 c_integrate_qagi :: FunPtr (Double-> Ptr() -> Double) -> Double -> CInt
131 -> Ptr Double -> Ptr Double -> IO CInt 131 -> Ptr Double -> Ptr Double -> IO CInt
132 132
@@ -159,7 +159,7 @@ integrateQAGIU prec n f a = unsafePerformIO $ do
159 freeHaskellFunPtr fp 159 freeHaskellFunPtr fp
160 return result 160 return result
161 161
162foreign import ccall "gsl-aux.h integrate_qagiu" 162foreign import ccall safe "gsl-aux.h integrate_qagiu"
163 c_integrate_qagiu :: FunPtr (Double-> Ptr() -> Double) -> Double -> Double -> CInt 163 c_integrate_qagiu :: FunPtr (Double-> Ptr() -> Double) -> Double -> Double -> CInt
164 -> Ptr Double -> Ptr Double -> IO CInt 164 -> Ptr Double -> Ptr Double -> IO CInt
165 165
@@ -192,7 +192,7 @@ integrateQAGIL prec n f b = unsafePerformIO $ do
192 freeHaskellFunPtr fp 192 freeHaskellFunPtr fp
193 return result 193 return result
194 194
195foreign import ccall "gsl-aux.h integrate_qagil" 195foreign import ccall safe "gsl-aux.h integrate_qagil"
196 c_integrate_qagil :: FunPtr (Double-> Ptr() -> Double) -> Double -> Double -> CInt 196 c_integrate_qagil :: FunPtr (Double-> Ptr() -> Double) -> Double -> Double -> CInt
197 -> Ptr Double -> Ptr Double -> IO CInt 197 -> Ptr Double -> Ptr Double -> IO CInt
198 198
diff --git a/lib/Numeric/GSL/Internal.hs b/lib/Numeric/GSL/Internal.hs
index 303e7e1..84417ce 100644
--- a/lib/Numeric/GSL/Internal.hs
+++ b/lib/Numeric/GSL/Internal.hs
@@ -26,14 +26,14 @@ iv f n p = f (createV (fromIntegral n) copy "iv") where
26 return 0 26 return 0
27 27
28-- | conversion of Haskell functions into function pointers that can be used in the C side 28-- | conversion of Haskell functions into function pointers that can be used in the C side
29foreign import ccall "wrapper" 29foreign import ccall safe "wrapper"
30 mkVecfun :: (CInt -> Ptr Double -> Double) 30 mkVecfun :: (CInt -> Ptr Double -> Double)
31 -> IO( FunPtr (CInt -> Ptr Double -> Double)) 31 -> IO( FunPtr (CInt -> Ptr Double -> Double))
32 32
33foreign import ccall "wrapper" 33foreign import ccall safe "wrapper"
34 mkVecVecfun :: TVV -> IO (FunPtr TVV) 34 mkVecVecfun :: TVV -> IO (FunPtr TVV)
35 35
36foreign import ccall "wrapper" 36foreign import ccall safe "wrapper"
37 mkDoubleVecVecfun :: (Double -> TVV) -> IO (FunPtr (Double -> TVV)) 37 mkDoubleVecVecfun :: (Double -> TVV) -> IO (FunPtr (Double -> TVV))
38 38
39aux_vTov :: (Vector Double -> Vector Double) -> TVV 39aux_vTov :: (Vector Double -> Vector Double) -> TVV
@@ -46,10 +46,10 @@ aux_vTov f n p nr r = g where
46 g = do unsafeWith v $ \p' -> copyArray r p' (fromIntegral nr) 46 g = do unsafeWith v $ \p' -> copyArray r p' (fromIntegral nr)
47 return 0 47 return 0
48 48
49foreign import ccall "wrapper" 49foreign import ccall safe "wrapper"
50 mkVecMatfun :: TVM -> IO (FunPtr TVM) 50 mkVecMatfun :: TVM -> IO (FunPtr TVM)
51 51
52foreign import ccall "wrapper" 52foreign import ccall safe "wrapper"
53 mkDoubleVecMatfun :: (Double -> TVM) -> IO (FunPtr (Double -> TVM)) 53 mkDoubleVecMatfun :: (Double -> TVM) -> IO (FunPtr (Double -> TVM))
54 54
55aux_vTom :: (Vector Double -> Matrix Double) -> TVM 55aux_vTom :: (Vector Double -> Matrix Double) -> TVM
diff --git a/lib/Numeric/GSL/Minimization.hs b/lib/Numeric/GSL/Minimization.hs
index 88acfac..af85135 100644
--- a/lib/Numeric/GSL/Minimization.hs
+++ b/lib/Numeric/GSL/Minimization.hs
@@ -122,7 +122,7 @@ minimizeV method eps maxit szv f xiv = unsafePerformIO $ do
122 return (sol, path) 122 return (sol, path)
123 123
124 124
125foreign import ccall "gsl-aux.h minimize" 125foreign import ccall safe "gsl-aux.h minimize"
126 c_minimize:: CInt -> FunPtr (CInt -> Ptr Double -> Double) -> Double -> CInt -> TVVM 126 c_minimize:: CInt -> FunPtr (CInt -> Ptr Double -> Double) -> Double -> CInt -> TVVM
127 127
128---------------------------------------------------------------------------------- 128----------------------------------------------------------------------------------
@@ -179,7 +179,7 @@ minimizeVD method eps maxit istep tol f df xiv = unsafePerformIO $ do
179 freeHaskellFunPtr dfp 179 freeHaskellFunPtr dfp
180 return (sol,path) 180 return (sol,path)
181 181
182foreign import ccall "gsl-aux.h minimizeD" 182foreign import ccall safe "gsl-aux.h minimizeD"
183 c_minimizeD :: CInt 183 c_minimizeD :: CInt
184 -> FunPtr (CInt -> Ptr Double -> Double) 184 -> FunPtr (CInt -> Ptr Double -> Double)
185 -> FunPtr TVV 185 -> FunPtr TVV
diff --git a/lib/Numeric/GSL/ODE.hs b/lib/Numeric/GSL/ODE.hs
index c243f4b..3994210 100644
--- a/lib/Numeric/GSL/ODE.hs
+++ b/lib/Numeric/GSL/ODE.hs
@@ -116,7 +116,7 @@ odeSolveV' method mbjac h epsAbs epsRel f xiv ts = unsafePerformIO $ do
116 freeHaskellFunPtr fp 116 freeHaskellFunPtr fp
117 return sol 117 return sol
118 118
119foreign import ccall "ode" 119foreign import ccall safe "ode"
120 ode_c :: CInt -> Double -> Double -> Double -> FunPtr (Double -> TVV) -> FunPtr (Double -> TVM) -> TVVM 120 ode_c :: CInt -> Double -> Double -> Double -> FunPtr (Double -> TVV) -> FunPtr (Double -> TVM) -> TVVM
121 121
122------------------------------------------------------- 122-------------------------------------------------------
diff --git a/lib/Numeric/GSL/Polynomials.hs b/lib/Numeric/GSL/Polynomials.hs
index b21b1b6..694c003 100644
--- a/lib/Numeric/GSL/Polynomials.hs
+++ b/lib/Numeric/GSL/Polynomials.hs
@@ -55,4 +55,4 @@ polySolve' v | dim v > 1 = unsafePerformIO $ do
55 return r 55 return r
56 | otherwise = error "polySolve on a polynomial of degree zero" 56 | otherwise = error "polySolve on a polynomial of degree zero"
57 57
58foreign import ccall "gsl-aux.h polySolve" c_polySolve:: TVCV 58foreign import ccall unsafe "gsl-aux.h polySolve" c_polySolve:: TVCV
diff --git a/lib/Numeric/GSL/Root.hs b/lib/Numeric/GSL/Root.hs
index 41f7fb3..cd2982a 100644
--- a/lib/Numeric/GSL/Root.hs
+++ b/lib/Numeric/GSL/Root.hs
@@ -91,7 +91,7 @@ rootGen m f xi epsabs maxit = unsafePerformIO $ do
91 return (take n $ drop 1 sol, path) 91 return (take n $ drop 1 sol, path)
92 92
93 93
94foreign import ccall "root" 94foreign import ccall safe "root"
95 c_root:: CInt -> FunPtr TVV -> Double -> CInt -> TVM 95 c_root:: CInt -> FunPtr TVV -> Double -> CInt -> TVM
96 96
97------------------------------------------------------------------------- 97-------------------------------------------------------------------------
@@ -130,7 +130,7 @@ rootJGen m f jac xi epsabs maxit = unsafePerformIO $ do
130 return (take n $ drop 1 sol, path) 130 return (take n $ drop 1 sol, path)
131 131
132 132
133foreign import ccall "rootj" 133foreign import ccall safe "rootj"
134 c_rootj:: CInt -> FunPtr TVV -> FunPtr TVM -> Double -> CInt -> TVM 134 c_rootj:: CInt -> FunPtr TVV -> FunPtr TVM -> Double -> CInt -> TVM
135 135
136------------------------------------------------------- 136-------------------------------------------------------
diff --git a/lib/Numeric/GSL/Vector.hs b/lib/Numeric/GSL/Vector.hs
index b8c6d04..db34041 100644
--- a/lib/Numeric/GSL/Vector.hs
+++ b/lib/Numeric/GSL/Vector.hs
@@ -109,10 +109,10 @@ sumC x = unsafePerformIO $ do
109 app2 c_sumC vec x vec r "sumC" 109 app2 c_sumC vec x vec r "sumC"
110 return $ r @> 0 110 return $ r @> 0
111 111
112foreign import ccall safe "gsl-aux.h sumF" c_sumF :: TFF 112foreign import ccall unsafe "gsl-aux.h sumF" c_sumF :: TFF
113foreign import ccall safe "gsl-aux.h sumR" c_sumR :: TVV 113foreign import ccall unsafe "gsl-aux.h sumR" c_sumR :: TVV
114foreign import ccall safe "gsl-aux.h sumQ" c_sumQ :: TQVQV 114foreign import ccall unsafe "gsl-aux.h sumQ" c_sumQ :: TQVQV
115foreign import ccall safe "gsl-aux.h sumC" c_sumC :: TCVCV 115foreign import ccall unsafe "gsl-aux.h sumC" c_sumC :: TCVCV
116 116
117-- | product of elements 117-- | product of elements
118prodF :: Vector Float -> Float 118prodF :: Vector Float -> Float
@@ -142,10 +142,10 @@ prodC x = unsafePerformIO $ do
142 app2 c_prodC vec x vec r "prodC" 142 app2 c_prodC vec x vec r "prodC"
143 return $ r @> 0 143 return $ r @> 0
144 144
145foreign import ccall safe "gsl-aux.h prodF" c_prodF :: TFF 145foreign import ccall unsafe "gsl-aux.h prodF" c_prodF :: TFF
146foreign import ccall safe "gsl-aux.h prodR" c_prodR :: TVV 146foreign import ccall unsafe "gsl-aux.h prodR" c_prodR :: TVV
147foreign import ccall safe "gsl-aux.h prodQ" c_prodQ :: TQVQV 147foreign import ccall unsafe "gsl-aux.h prodQ" c_prodQ :: TQVQV
148foreign import ccall safe "gsl-aux.h prodC" c_prodC :: TCVCV 148foreign import ccall unsafe "gsl-aux.h prodC" c_prodC :: TCVCV
149 149
150-- | dot product 150-- | dot product
151dotF :: Vector Float -> Vector Float -> Float 151dotF :: Vector Float -> Vector Float -> Float
@@ -175,10 +175,10 @@ dotC x y = unsafePerformIO $ do
175 app3 c_dotC vec x vec y vec r "dotC" 175 app3 c_dotC vec x vec y vec r "dotC"
176 return $ r @> 0 176 return $ r @> 0
177 177
178foreign import ccall safe "gsl-aux.h dotF" c_dotF :: TFFF 178foreign import ccall unsafe "gsl-aux.h dotF" c_dotF :: TFFF
179foreign import ccall safe "gsl-aux.h dotR" c_dotR :: TVVV 179foreign import ccall unsafe "gsl-aux.h dotR" c_dotR :: TVVV
180foreign import ccall safe "gsl-aux.h dotQ" c_dotQ :: TQVQVQV 180foreign import ccall unsafe "gsl-aux.h dotQ" c_dotQ :: TQVQVQV
181foreign import ccall safe "gsl-aux.h dotC" c_dotC :: TCVCVCV 181foreign import ccall unsafe "gsl-aux.h dotC" c_dotC :: TCVCVCV
182 182
183------------------------------------------------------------------ 183------------------------------------------------------------------
184 184
@@ -210,25 +210,25 @@ vectorZipAux fun code u v = unsafePerformIO $ do
210toScalarR :: FunCodeS -> Vector Double -> Double 210toScalarR :: FunCodeS -> Vector Double -> Double
211toScalarR oper = toScalarAux c_toScalarR (fromei oper) 211toScalarR oper = toScalarAux c_toScalarR (fromei oper)
212 212
213foreign import ccall safe "gsl-aux.h toScalarR" c_toScalarR :: CInt -> TVV 213foreign import ccall unsafe "gsl-aux.h toScalarR" c_toScalarR :: CInt -> TVV
214 214
215-- | obtains different functions of a vector: norm1, norm2, max, min, posmax, posmin, etc. 215-- | obtains different functions of a vector: norm1, norm2, max, min, posmax, posmin, etc.
216toScalarF :: FunCodeS -> Vector Float -> Float 216toScalarF :: FunCodeS -> Vector Float -> Float
217toScalarF oper = toScalarAux c_toScalarF (fromei oper) 217toScalarF oper = toScalarAux c_toScalarF (fromei oper)
218 218
219foreign import ccall safe "gsl-aux.h toScalarF" c_toScalarF :: CInt -> TFF 219foreign import ccall unsafe "gsl-aux.h toScalarF" c_toScalarF :: CInt -> TFF
220 220
221-- | obtains different functions of a vector: only norm1, norm2 221-- | obtains different functions of a vector: only norm1, norm2
222toScalarC :: FunCodeS -> Vector (Complex Double) -> Double 222toScalarC :: FunCodeS -> Vector (Complex Double) -> Double
223toScalarC oper = toScalarAux c_toScalarC (fromei oper) 223toScalarC oper = toScalarAux c_toScalarC (fromei oper)
224 224
225foreign import ccall safe "gsl-aux.h toScalarC" c_toScalarC :: CInt -> TCVV 225foreign import ccall unsafe "gsl-aux.h toScalarC" c_toScalarC :: CInt -> TCVV
226 226
227-- | obtains different functions of a vector: only norm1, norm2 227-- | obtains different functions of a vector: only norm1, norm2
228toScalarQ :: FunCodeS -> Vector (Complex Float) -> Float 228toScalarQ :: FunCodeS -> Vector (Complex Float) -> Float
229toScalarQ oper = toScalarAux c_toScalarQ (fromei oper) 229toScalarQ oper = toScalarAux c_toScalarQ (fromei oper)
230 230
231foreign import ccall safe "gsl-aux.h toScalarQ" c_toScalarQ :: CInt -> TQVF 231foreign import ccall unsafe "gsl-aux.h toScalarQ" c_toScalarQ :: CInt -> TQVF
232 232
233------------------------------------------------------------------ 233------------------------------------------------------------------
234 234
@@ -236,25 +236,25 @@ foreign import ccall safe "gsl-aux.h toScalarQ" c_toScalarQ :: CInt -> TQVF
236vectorMapR :: FunCodeV -> Vector Double -> Vector Double 236vectorMapR :: FunCodeV -> Vector Double -> Vector Double
237vectorMapR = vectorMapAux c_vectorMapR 237vectorMapR = vectorMapAux c_vectorMapR
238 238
239foreign import ccall safe "gsl-aux.h mapR" c_vectorMapR :: CInt -> TVV 239foreign import ccall unsafe "gsl-aux.h mapR" c_vectorMapR :: CInt -> TVV
240 240
241-- | map of complex vectors with given function 241-- | map of complex vectors with given function
242vectorMapC :: FunCodeV -> Vector (Complex Double) -> Vector (Complex Double) 242vectorMapC :: FunCodeV -> Vector (Complex Double) -> Vector (Complex Double)
243vectorMapC oper = vectorMapAux c_vectorMapC (fromei oper) 243vectorMapC oper = vectorMapAux c_vectorMapC (fromei oper)
244 244
245foreign import ccall safe "gsl-aux.h mapC" c_vectorMapC :: CInt -> TCVCV 245foreign import ccall unsafe "gsl-aux.h mapC" c_vectorMapC :: CInt -> TCVCV
246 246
247-- | map of real vectors with given function 247-- | map of real vectors with given function
248vectorMapF :: FunCodeV -> Vector Float -> Vector Float 248vectorMapF :: FunCodeV -> Vector Float -> Vector Float
249vectorMapF = vectorMapAux c_vectorMapF 249vectorMapF = vectorMapAux c_vectorMapF
250 250
251foreign import ccall safe "gsl-aux.h mapF" c_vectorMapF :: CInt -> TFF 251foreign import ccall unsafe "gsl-aux.h mapF" c_vectorMapF :: CInt -> TFF
252 252
253-- | map of real vectors with given function 253-- | map of real vectors with given function
254vectorMapQ :: FunCodeV -> Vector (Complex Float) -> Vector (Complex Float) 254vectorMapQ :: FunCodeV -> Vector (Complex Float) -> Vector (Complex Float)
255vectorMapQ = vectorMapAux c_vectorMapQ 255vectorMapQ = vectorMapAux c_vectorMapQ
256 256
257foreign import ccall safe "gsl-aux.h mapQ" c_vectorMapQ :: CInt -> TQVQV 257foreign import ccall unsafe "gsl-aux.h mapQ" c_vectorMapQ :: CInt -> TQVQV
258 258
259------------------------------------------------------------------- 259-------------------------------------------------------------------
260 260
@@ -262,25 +262,25 @@ foreign import ccall safe "gsl-aux.h mapQ" c_vectorMapQ :: CInt -> TQVQV
262vectorMapValR :: FunCodeSV -> Double -> Vector Double -> Vector Double 262vectorMapValR :: FunCodeSV -> Double -> Vector Double -> Vector Double
263vectorMapValR oper = vectorMapValAux c_vectorMapValR (fromei oper) 263vectorMapValR oper = vectorMapValAux c_vectorMapValR (fromei oper)
264 264
265foreign import ccall safe "gsl-aux.h mapValR" c_vectorMapValR :: CInt -> Ptr Double -> TVV 265foreign import ccall unsafe "gsl-aux.h mapValR" c_vectorMapValR :: CInt -> Ptr Double -> TVV
266 266
267-- | map of complex vectors with given function 267-- | map of complex vectors with given function
268vectorMapValC :: FunCodeSV -> Complex Double -> Vector (Complex Double) -> Vector (Complex Double) 268vectorMapValC :: FunCodeSV -> Complex Double -> Vector (Complex Double) -> Vector (Complex Double)
269vectorMapValC = vectorMapValAux c_vectorMapValC 269vectorMapValC = vectorMapValAux c_vectorMapValC
270 270
271foreign import ccall safe "gsl-aux.h mapValC" c_vectorMapValC :: CInt -> Ptr (Complex Double) -> TCVCV 271foreign import ccall unsafe "gsl-aux.h mapValC" c_vectorMapValC :: CInt -> Ptr (Complex Double) -> TCVCV
272 272
273-- | map of real vectors with given function 273-- | map of real vectors with given function
274vectorMapValF :: FunCodeSV -> Float -> Vector Float -> Vector Float 274vectorMapValF :: FunCodeSV -> Float -> Vector Float -> Vector Float
275vectorMapValF oper = vectorMapValAux c_vectorMapValF (fromei oper) 275vectorMapValF oper = vectorMapValAux c_vectorMapValF (fromei oper)
276 276
277foreign import ccall safe "gsl-aux.h mapValF" c_vectorMapValF :: CInt -> Ptr Float -> TFF 277foreign import ccall unsafe "gsl-aux.h mapValF" c_vectorMapValF :: CInt -> Ptr Float -> TFF
278 278
279-- | map of complex vectors with given function 279-- | map of complex vectors with given function
280vectorMapValQ :: FunCodeSV -> Complex Float -> Vector (Complex Float) -> Vector (Complex Float) 280vectorMapValQ :: FunCodeSV -> Complex Float -> Vector (Complex Float) -> Vector (Complex Float)
281vectorMapValQ oper = vectorMapValAux c_vectorMapValQ (fromei oper) 281vectorMapValQ oper = vectorMapValAux c_vectorMapValQ (fromei oper)
282 282
283foreign import ccall safe "gsl-aux.h mapValQ" c_vectorMapValQ :: CInt -> Ptr (Complex Float) -> TQVQV 283foreign import ccall unsafe "gsl-aux.h mapValQ" c_vectorMapValQ :: CInt -> Ptr (Complex Float) -> TQVQV
284 284
285------------------------------------------------------------------- 285-------------------------------------------------------------------
286 286
@@ -288,25 +288,25 @@ foreign import ccall safe "gsl-aux.h mapValQ" c_vectorMapValQ :: CInt -> Ptr (Co
288vectorZipR :: FunCodeVV -> Vector Double -> Vector Double -> Vector Double 288vectorZipR :: FunCodeVV -> Vector Double -> Vector Double -> Vector Double
289vectorZipR = vectorZipAux c_vectorZipR 289vectorZipR = vectorZipAux c_vectorZipR
290 290
291foreign import ccall safe "gsl-aux.h zipR" c_vectorZipR :: CInt -> TVVV 291foreign import ccall unsafe "gsl-aux.h zipR" c_vectorZipR :: CInt -> TVVV
292 292
293-- | elementwise operation on complex vectors 293-- | elementwise operation on complex vectors
294vectorZipC :: FunCodeVV -> Vector (Complex Double) -> Vector (Complex Double) -> Vector (Complex Double) 294vectorZipC :: FunCodeVV -> Vector (Complex Double) -> Vector (Complex Double) -> Vector (Complex Double)
295vectorZipC = vectorZipAux c_vectorZipC 295vectorZipC = vectorZipAux c_vectorZipC
296 296
297foreign import ccall safe "gsl-aux.h zipC" c_vectorZipC :: CInt -> TCVCVCV 297foreign import ccall unsafe "gsl-aux.h zipC" c_vectorZipC :: CInt -> TCVCVCV
298 298
299-- | elementwise operation on real vectors 299-- | elementwise operation on real vectors
300vectorZipF :: FunCodeVV -> Vector Float -> Vector Float -> Vector Float 300vectorZipF :: FunCodeVV -> Vector Float -> Vector Float -> Vector Float
301vectorZipF = vectorZipAux c_vectorZipF 301vectorZipF = vectorZipAux c_vectorZipF
302 302
303foreign import ccall safe "gsl-aux.h zipF" c_vectorZipF :: CInt -> TFFF 303foreign import ccall unsafe "gsl-aux.h zipF" c_vectorZipF :: CInt -> TFFF
304 304
305-- | elementwise operation on complex vectors 305-- | elementwise operation on complex vectors
306vectorZipQ :: FunCodeVV -> Vector (Complex Float) -> Vector (Complex Float) -> Vector (Complex Float) 306vectorZipQ :: FunCodeVV -> Vector (Complex Float) -> Vector (Complex Float) -> Vector (Complex Float)
307vectorZipQ = vectorZipAux c_vectorZipQ 307vectorZipQ = vectorZipAux c_vectorZipQ
308 308
309foreign import ccall safe "gsl-aux.h zipQ" c_vectorZipQ :: CInt -> TQVQVQV 309foreign import ccall unsafe "gsl-aux.h zipQ" c_vectorZipQ :: CInt -> TQVQVQV
310 310
311----------------------------------------------------------------------- 311-----------------------------------------------------------------------
312 312
@@ -324,4 +324,4 @@ randomVector seed dist n = unsafePerformIO $ do
324 app1 (c_random_vector (fi seed) ((fi.fromEnum) dist)) vec r "randomVector" 324 app1 (c_random_vector (fi seed) ((fi.fromEnum) dist)) vec r "randomVector"
325 return r 325 return r
326 326
327foreign import ccall safe "random_vector" c_random_vector :: CInt -> CInt -> TV 327foreign import ccall unsafe "random_vector" c_random_vector :: CInt -> CInt -> TV
diff --git a/lib/Numeric/LinearAlgebra/LAPACK.hs b/lib/Numeric/LinearAlgebra/LAPACK.hs
index 7ad58e1..9bc6669 100644
--- a/lib/Numeric/LinearAlgebra/LAPACK.hs
+++ b/lib/Numeric/LinearAlgebra/LAPACK.hs
@@ -53,10 +53,10 @@ import System.IO.Unsafe(unsafePerformIO)
53 53
54----------------------------------------------------------------------------------- 54-----------------------------------------------------------------------------------
55 55
56foreign import ccall "multiplyR" dgemmc :: CInt -> CInt -> TMMM 56foreign import ccall unsafe "multiplyR" dgemmc :: CInt -> CInt -> TMMM
57foreign import ccall "multiplyC" zgemmc :: CInt -> CInt -> TCMCMCM 57foreign import ccall unsafe "multiplyC" zgemmc :: CInt -> CInt -> TCMCMCM
58foreign import ccall "multiplyF" sgemmc :: CInt -> CInt -> TFMFMFM 58foreign import ccall unsafe "multiplyF" sgemmc :: CInt -> CInt -> TFMFMFM
59foreign import ccall "multiplyQ" cgemmc :: CInt -> CInt -> TQMQMQM 59foreign import ccall unsafe "multiplyQ" cgemmc :: CInt -> CInt -> TQMQMQM
60 60
61isT Matrix{order = ColumnMajor} = 0 61isT Matrix{order = ColumnMajor} = 0
62isT Matrix{order = RowMajor} = 1 62isT Matrix{order = RowMajor} = 1
@@ -88,10 +88,10 @@ multiplyQ :: Matrix (Complex Float) -> Matrix (Complex Float) -> Matrix (Complex
88multiplyQ a b = multiplyAux cgemmc "cgemmc" a b 88multiplyQ a b = multiplyAux cgemmc "cgemmc" a b
89 89
90----------------------------------------------------------------------------- 90-----------------------------------------------------------------------------
91foreign import ccall "svd_l_R" dgesvd :: TMMVM 91foreign import ccall unsafe "svd_l_R" dgesvd :: TMMVM
92foreign import ccall "svd_l_C" zgesvd :: TCMCMVCM 92foreign import ccall unsafe "svd_l_C" zgesvd :: TCMCMVCM
93foreign import ccall "svd_l_Rdd" dgesdd :: TMMVM 93foreign import ccall unsafe "svd_l_Rdd" dgesdd :: TMMVM
94foreign import ccall "svd_l_Cdd" zgesdd :: TCMCMVCM 94foreign import ccall unsafe "svd_l_Cdd" zgesdd :: TCMCMVCM
95 95
96-- | Full SVD of a real matrix using LAPACK's /dgesvd/. 96-- | Full SVD of a real matrix using LAPACK's /dgesvd/.
97svdR :: Matrix Double -> (Matrix Double, Vector Double, Matrix Double) 97svdR :: Matrix Double -> (Matrix Double, Vector Double, Matrix Double)
@@ -211,10 +211,10 @@ leftSVAux f st x = unsafePerformIO $ do
211 211
212----------------------------------------------------------------------------- 212-----------------------------------------------------------------------------
213 213
214foreign import ccall "LAPACK/lapack-aux.h eig_l_R" dgeev :: TMMCVM 214foreign import ccall unsafe "LAPACK/lapack-aux.h eig_l_R" dgeev :: TMMCVM
215foreign import ccall "LAPACK/lapack-aux.h eig_l_C" zgeev :: TCMCMCVCM 215foreign import ccall unsafe "LAPACK/lapack-aux.h eig_l_C" zgeev :: TCMCMCVCM
216foreign import ccall "LAPACK/lapack-aux.h eig_l_S" dsyev :: CInt -> TMVM 216foreign import ccall unsafe "LAPACK/lapack-aux.h eig_l_S" dsyev :: CInt -> TMVM
217foreign import ccall "LAPACK/lapack-aux.h eig_l_H" zheev :: CInt -> TCMVCM 217foreign import ccall unsafe "LAPACK/lapack-aux.h eig_l_H" zheev :: CInt -> TCMVCM
218 218
219eigAux f st m = unsafePerformIO $ do 219eigAux f st m = unsafePerformIO $ do
220 l <- createVector r 220 l <- createVector r
@@ -325,10 +325,10 @@ eigOnlyH = vrev . fst. eigSHAux (zheev 1) "eigH'" . fmat
325vrev = flatten . flipud . reshape 1 325vrev = flatten . flipud . reshape 1
326 326
327----------------------------------------------------------------------------- 327-----------------------------------------------------------------------------
328foreign import ccall "linearSolveR_l" dgesv :: TMMM 328foreign import ccall unsafe "linearSolveR_l" dgesv :: TMMM
329foreign import ccall "linearSolveC_l" zgesv :: TCMCMCM 329foreign import ccall unsafe "linearSolveC_l" zgesv :: TCMCMCM
330foreign import ccall "cholSolveR_l" dpotrs :: TMMM 330foreign import ccall unsafe "cholSolveR_l" dpotrs :: TMMM
331foreign import ccall "cholSolveC_l" zpotrs :: TCMCMCM 331foreign import ccall unsafe "cholSolveC_l" zpotrs :: TCMCMCM
332 332
333linearSolveSQAux f st a b 333linearSolveSQAux f st a b
334 | n1==n2 && n1==r = unsafePerformIO $ do 334 | n1==n2 && n1==r = unsafePerformIO $ do
@@ -359,10 +359,10 @@ cholSolveC :: Matrix (Complex Double) -> Matrix (Complex Double) -> Matrix (Comp
359cholSolveC a b = linearSolveSQAux zpotrs "cholSolveC" (fmat a) (fmat b) 359cholSolveC a b = linearSolveSQAux zpotrs "cholSolveC" (fmat a) (fmat b)
360 360
361----------------------------------------------------------------------------------- 361-----------------------------------------------------------------------------------
362foreign import ccall "LAPACK/lapack-aux.h linearSolveLSR_l" dgels :: TMMM 362foreign import ccall unsafe "LAPACK/lapack-aux.h linearSolveLSR_l" dgels :: TMMM
363foreign import ccall "LAPACK/lapack-aux.h linearSolveLSC_l" zgels :: TCMCMCM 363foreign import ccall unsafe "LAPACK/lapack-aux.h linearSolveLSC_l" zgels :: TCMCMCM
364foreign import ccall "LAPACK/lapack-aux.h linearSolveSVDR_l" dgelss :: Double -> TMMM 364foreign import ccall unsafe "LAPACK/lapack-aux.h linearSolveSVDR_l" dgelss :: Double -> TMMM
365foreign import ccall "LAPACK/lapack-aux.h linearSolveSVDC_l" zgelss :: Double -> TCMCMCM 365foreign import ccall unsafe "LAPACK/lapack-aux.h linearSolveSVDC_l" zgelss :: Double -> TCMCMCM
366 366
367linearSolveAux f st a b = unsafePerformIO $ do 367linearSolveAux f st a b = unsafePerformIO $ do
368 r <- createMatrix ColumnMajor (max m n) nrhs 368 r <- createMatrix ColumnMajor (max m n) nrhs
@@ -401,8 +401,8 @@ linearSolveSVDC (Just rcond) a b = subMatrix (0,0) (cols a, cols b) $
401linearSolveSVDC Nothing a b = linearSolveSVDC (Just (-1)) (fmat a) (fmat b) 401linearSolveSVDC Nothing a b = linearSolveSVDC (Just (-1)) (fmat a) (fmat b)
402 402
403----------------------------------------------------------------------------------- 403-----------------------------------------------------------------------------------
404foreign import ccall "LAPACK/lapack-aux.h chol_l_H" zpotrf :: TCMCM 404foreign import ccall unsafe "LAPACK/lapack-aux.h chol_l_H" zpotrf :: TCMCM
405foreign import ccall "LAPACK/lapack-aux.h chol_l_S" dpotrf :: TMM 405foreign import ccall unsafe "LAPACK/lapack-aux.h chol_l_S" dpotrf :: TMM
406 406
407cholAux f st a = do 407cholAux f st a = do
408 r <- createMatrix ColumnMajor n n 408 r <- createMatrix ColumnMajor n n
@@ -427,8 +427,8 @@ mbCholS :: Matrix Double -> Maybe (Matrix Double)
427mbCholS = unsafePerformIO . mbCatch . cholAux dpotrf "cholS" . fmat 427mbCholS = unsafePerformIO . mbCatch . cholAux dpotrf "cholS" . fmat
428 428
429----------------------------------------------------------------------------------- 429-----------------------------------------------------------------------------------
430foreign import ccall "LAPACK/lapack-aux.h qr_l_R" dgeqr2 :: TMVM 430foreign import ccall unsafe "LAPACK/lapack-aux.h qr_l_R" dgeqr2 :: TMVM
431foreign import ccall "LAPACK/lapack-aux.h qr_l_C" zgeqr2 :: TCMCVCM 431foreign import ccall unsafe "LAPACK/lapack-aux.h qr_l_C" zgeqr2 :: TCMCVCM
432 432
433-- | QR factorization of a real matrix, using LAPACK's /dgeqr2/. 433-- | QR factorization of a real matrix, using LAPACK's /dgeqr2/.
434qrR :: Matrix Double -> (Matrix Double, Vector Double) 434qrR :: Matrix Double -> (Matrix Double, Vector Double)
@@ -448,8 +448,8 @@ qrAux f st a = unsafePerformIO $ do
448 mn = min m n 448 mn = min m n
449 449
450----------------------------------------------------------------------------------- 450-----------------------------------------------------------------------------------
451foreign import ccall "LAPACK/lapack-aux.h hess_l_R" dgehrd :: TMVM 451foreign import ccall unsafe "LAPACK/lapack-aux.h hess_l_R" dgehrd :: TMVM
452foreign import ccall "LAPACK/lapack-aux.h hess_l_C" zgehrd :: TCMCVCM 452foreign import ccall unsafe "LAPACK/lapack-aux.h hess_l_C" zgehrd :: TCMCVCM
453 453
454-- | Hessenberg factorization of a square real matrix, using LAPACK's /dgehrd/. 454-- | Hessenberg factorization of a square real matrix, using LAPACK's /dgehrd/.
455hessR :: Matrix Double -> (Matrix Double, Vector Double) 455hessR :: Matrix Double -> (Matrix Double, Vector Double)
@@ -469,8 +469,8 @@ hessAux f st a = unsafePerformIO $ do
469 mn = min m n 469 mn = min m n
470 470
471----------------------------------------------------------------------------------- 471-----------------------------------------------------------------------------------
472foreign import ccall "LAPACK/lapack-aux.h schur_l_R" dgees :: TMMM 472foreign import ccall unsafe "LAPACK/lapack-aux.h schur_l_R" dgees :: TMMM
473foreign import ccall "LAPACK/lapack-aux.h schur_l_C" zgees :: TCMCMCM 473foreign import ccall unsafe "LAPACK/lapack-aux.h schur_l_C" zgees :: TCMCMCM
474 474
475-- | Schur factorization of a square real matrix, using LAPACK's /dgees/. 475-- | Schur factorization of a square real matrix, using LAPACK's /dgees/.
476schurR :: Matrix Double -> (Matrix Double, Matrix Double) 476schurR :: Matrix Double -> (Matrix Double, Matrix Double)
@@ -488,8 +488,8 @@ schurAux f st a = unsafePerformIO $ do
488 where n = rows a 488 where n = rows a
489 489
490----------------------------------------------------------------------------------- 490-----------------------------------------------------------------------------------
491foreign import ccall "LAPACK/lapack-aux.h lu_l_R" dgetrf :: TMVM 491foreign import ccall unsafe "LAPACK/lapack-aux.h lu_l_R" dgetrf :: TMVM
492foreign import ccall "LAPACK/lapack-aux.h lu_l_C" zgetrf :: TCMVCM 492foreign import ccall unsafe "LAPACK/lapack-aux.h lu_l_C" zgetrf :: TCMVCM
493 493
494-- | LU factorization of a general real matrix, using LAPACK's /dgetrf/. 494-- | LU factorization of a general real matrix, using LAPACK's /dgetrf/.
495luR :: Matrix Double -> (Matrix Double, [Int]) 495luR :: Matrix Double -> (Matrix Double, [Int])
@@ -511,8 +511,8 @@ luAux f st a = unsafePerformIO $ do
511type TW a = CInt -> PD -> a 511type TW a = CInt -> PD -> a
512type TQ a = CInt -> CInt -> PC -> a 512type TQ a = CInt -> CInt -> PC -> a
513 513
514foreign import ccall "LAPACK/lapack-aux.h luS_l_R" dgetrs :: TMVMM 514foreign import ccall unsafe "LAPACK/lapack-aux.h luS_l_R" dgetrs :: TMVMM
515foreign import ccall "LAPACK/lapack-aux.h luS_l_C" zgetrs :: TQ (TW (TQ (TQ (IO CInt)))) 515foreign import ccall unsafe "LAPACK/lapack-aux.h luS_l_C" zgetrs :: TQ (TW (TQ (TQ (IO CInt))))
516 516
517-- | Solve a real linear system from a precomputed LU decomposition ('luR'), using LAPACK's /dgetrs/. 517-- | Solve a real linear system from a precomputed LU decomposition ('luR'), using LAPACK's /dgetrs/.
518lusR :: Matrix Double -> [Int] -> Matrix Double -> Matrix Double 518lusR :: Matrix Double -> [Int] -> Matrix Double -> Matrix Double