netlib-carray (empty) → 0.0

raw patch · 4 files changed

+195/−0 lines, 4 filesdep +basedep +carraydep +storable-complexsetup-changed

Dependencies added: base, carray, storable-complex, transformers

Files

LICENSE +27/−0
Setup.lhs +3/−0
netlib-carray.cabal +42/−0
src/Numeric/Netlib/CArray/Utility.hs +123/−0

+ LICENSE view

@@ -0,0 +1,27 @@+Copyright (c) Henning Thielemann 2017++All rights reserved.++Redistribution and use in source and binary forms, with or without+modification, are permitted provided that the following conditions+are met:+1. Redistributions of source code must retain the above copyright+   notice, this list of conditions and the following disclaimer.+2. Redistributions in binary form must reproduce the above copyright+   notice, this list of conditions and the following disclaimer in the+   documentation and/or other materials provided with the distribution.+3. Neither the name of the author nor the names of his contributors+   may be used to endorse or promote products derived from this software+   without specific prior written permission.++THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND+ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE+IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE+ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHORS OR CONTRIBUTORS BE LIABLE+FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL+DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS+OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)+HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT+LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY+OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF+SUCH DAMAGE.

+ Setup.lhs view

@@ -0,0 +1,3 @@+#! /usr/bin/env runhaskell+> import Distribution.Simple+> main = defaultMain

+ netlib-carray.cabal view

@@ -0,0 +1,42 @@+Name:             netlib-carray+Version:          0.0+License:          BSD3+License-File:     LICENSE+Author:           Henning Thielemann <haskell@henning-thielemann.de>+Maintainer:       Henning Thielemann <haskell@henning-thielemann.de>+Homepage:         http://hub.darcs.net/thielema/netlib-carray/+Category:         Math+Synopsis:         Helper modules for CArray wrappers to BLAS and LAPACK+Description:+  Netlib is a collection of packages for efficient numeric linear algebra.+  Most prominent parts of Netlib are BLAS and LAPACK.+  These packages contain functions for matrix computations,+  solution of simultaneous linear equations and eigenvalue problems.+  .+  This package provides definitions shared by+  the packages @blas-carray@ and @lapack-carray@.+Tested-With:      GHC==7.4.2, GHC==7.8.4+Cabal-Version:    >=1.14+Build-Type:       Simple++Source-Repository this+  Tag:         0.0+  Type:        darcs+  Location:    http://hub.darcs.net/thielema/netlib-carray/++Source-Repository head+  Type:        darcs+  Location:    http://hub.darcs.net/thielema/netlib-carray/++Library+  Build-Depends:+    carray >=0.1.5 && <0.2,+    storable-complex >=0.2.2 && <0.3,+    transformers >=0.4 && <0.6,+    base >=4.5 && <5++  GHC-Options:      -Wall -fwarn-missing-import-lists+  Hs-Source-Dirs:   src+  Default-Language: Haskell98+  Exposed-Modules:+    Numeric.Netlib.CArray.Utility

+ src/Numeric/Netlib/CArray/Utility.hs view

@@ -0,0 +1,123 @@+module Numeric.Netlib.CArray.Utility where++import qualified Data.Array.CArray as CArray+import Data.Array.IOCArray (IOCArray, withIOCArray)+import Data.Array.CArray (CArray, withCArray, Ix)++import qualified Foreign.Marshal.Utils as Marshal+import qualified Foreign.Marshal.Array as Array+import qualified Foreign.Marshal.Alloc as Alloc+import qualified Foreign.C.String as CStr+import qualified Foreign.C.Types as C+import Foreign.Storable.Complex ()+import Foreign.Storable (Storable, peek)+import Foreign.Ptr (Ptr)++import Control.Monad.Trans.Cont (ContT(ContT))+import Control.Monad.IO.Class (liftIO)+import Control.Monad (when)+import Control.Applicative ((<$>))++import Data.Complex (Complex)+++type FortranIO r = ContT r IO++run :: FortranIO r (IO a) -> FortranIO r a+run act = act >>= liftIO++runChecked :: String -> FortranIO r (Ptr C.CInt -> IO a) -> FortranIO r a+runChecked name act = do+   info <- alloca+   a <- run $ fmap ($info) act+   liftIO $ check name (peek info)+   return a++check :: String -> IO C.CInt -> IO ()+check msg f = do+   err <- f+   when (err/=0) $ error $ msg ++ ": " ++ show err++assert :: String -> Bool -> IO ()+assert msg success = when (not success) $ error $ "assertion failed: " ++ msg++ignore :: String -> Int -> IO ()+ignore _msg _dim = return ()+++newArray :: (Ix i, Storable e) => (i, i) -> IO (CArray i e)+newArray bnds = CArray.createCArray bnds (\_ -> return ())++newArray1 :: (Storable e) => Int -> IO (CArray Int e)+newArray1 m = newArray (0, m-1)++newArray2 :: (Storable e) => Int -> Int -> IO (CArray (Int,Int) e)+newArray2 m n = newArray ((0,0), (m-1,n-1))++newArray3 :: (Storable e) => Int -> Int -> Int -> IO (CArray (Int,Int,Int) e)+newArray3 m n k = newArray ((0,0,0), (m-1,n-1,k-1))+++sizes1 :: (Ix i) => (i,i) -> Int+sizes1 = CArray.rangeSize++sizes2 :: (Ix i, Ix j) => ((i,j),(i,j)) -> (Int,Int)+sizes2 ((i0,j0), (i1,j1)) =+   (CArray.rangeSize (i0,i1), CArray.rangeSize (j0,j1))++sizes3 :: (Ix i, Ix j, Ix k) => ((i,j,k),(i,j,k)) -> (Int,Int,Int)+sizes3 ((i0,j0,k0), (i1,j1,k1)) =+   (CArray.rangeSize (i0,i1),+    CArray.rangeSize (j0,j1),+    CArray.rangeSize (k0,k1))+++cint :: Int -> FortranIO r (Ptr C.CInt)+cint = ContT . Marshal.with . fromIntegral++range :: (Int,Int) -> FortranIO r (Ptr C.CInt)+range = cint . CArray.rangeSize++alloca :: (Storable a) => FortranIO r (Ptr a)+alloca = ContT Alloc.alloca++allocaArray :: (Storable a) => Int -> FortranIO r (Ptr a)+allocaArray = ContT . Array.allocaArray++bool :: Bool -> FortranIO r (Ptr Bool)+bool = ContT . Marshal.with++char :: Char -> FortranIO r (Ptr C.CChar)+char = ContT . Marshal.with . CStr.castCharToCChar++string :: String -> FortranIO r (Ptr C.CChar)+string = ContT . CStr.withCString++float :: Float -> FortranIO r (Ptr Float)+float = ContT . Marshal.with++double :: Double -> FortranIO r (Ptr Double)+double = ContT . Marshal.with++complexFloat :: Complex Float -> FortranIO r (Ptr (Complex Float))+complexFloat = ContT . Marshal.with++complexDouble :: Complex Double -> FortranIO r (Ptr (Complex Double))+complexDouble = ContT . Marshal.with+++array :: (Storable a) => CArray i a -> FortranIO r (Ptr a)+array = ContT . withCArray++arrayBounds :: (Storable a, Ix i) => CArray i a -> FortranIO r (Ptr a, (i,i))+arrayBounds v = flip (,) (CArray.bounds v) <$> array v++ioarray :: (Storable a) => IOCArray i a -> FortranIO r (Ptr a)+ioarray = ContT . withIOCArray+++unzipBounds :: ((i,j),(i,j)) -> ((i,i), (j,j))+unzipBounds ((i0,j0), (i1,j1)) = ((i0,i1), (j0,j1))++(^!) :: (Num a) => a -> Int -> a+x^!n = x^n