{-# LANGUAGE TypeOperators #-}
{-# LANGUAGE TypeFamilies #-}
module Internal.Devel where
import Control.Monad ( when )
import Foreign.C.Types ( CInt )
import Foreign.Ptr(Ptr)
import Control.Exception as E ( SomeException, catch )
import Internal.Vector(Vector,avec)
import Foreign.Storable(Storable)
(//) :: x -> (x -> y) -> y
infixl 0 //
(//) = flip ($)
errorCode :: CInt -> String
errorCode 2000 = "bad size"
errorCode 2001 = "bad function code"
errorCode 2002 = "memory problem"
errorCode 2003 = "bad file"
errorCode 2004 = "singular"
errorCode 2005 = "didn't converge"
errorCode 2006 = "the input matrix is not positive definite"
errorCode 2007 = "not yet supported in this OS"
errorCode n = "code "++show n
foreign import ccall unsafe "asm_finit" finit :: IO ()
check :: String -> IO CInt -> IO ()
check msg f = do
err <- f
when (err/=0) $ error (msg++": "++errorCode err)
return ()
infixl 0 #|
(#|) :: IO CInt -> String -> IO ()
(#|) = flip check
mbCatch :: IO x -> IO (Maybe x)
mbCatch act = E.catch (Just `fmap` act) f
where f :: SomeException -> IO (Maybe x)
f _ = return Nothing
type CM b r = CInt -> CInt -> Ptr b -> r
type CV b r = CInt -> Ptr b -> r
type OM b r = CInt -> CInt -> CInt -> CInt -> Ptr b -> r
type CIdxs r = CV CInt r
type Ok = IO CInt
infixr 5 :>, ::>, ..>
type (:>) t r = CV t r
type (::>) t r = OM t r
type (..>) t r = CM t r
class TransArray c
where
type Trans c b
type TransRaw c b
apply :: c -> (b -> IO r) -> (Trans c b) -> IO r
applyRaw :: c -> (b -> IO r) -> (TransRaw c b) -> IO r
infixl 1 `apply`, `applyRaw`
instance Storable t => TransArray (Vector t)
where
type Trans (Vector t) b = CInt -> Ptr t -> b
type TransRaw (Vector t) b = CInt -> Ptr t -> b
apply = avec
{-# INLINE apply #-}
applyRaw = avec
{-# INLINE applyRaw #-}