1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
|
{-# LANGUAGE ForeignFunctionInterface #-}
-----------------------------------------------------------------------------
{- |
Module : Numeric.GSL.Special.Internal
Copyright : (c) Alberto Ruiz 2007
License : GPL-style
Maintainer : Alberto Ruiz (aruiz at um dot es)
Stability : provisional
Portability : uses ffi
Support for Special functions.
<http://www.gnu.org/software/gsl/manual/html_node/Special-Functions.html#Special-Functions>
-}
-----------------------------------------------------------------------------
#include <gsl/gsl_sf_result.h>
#let alignment t = "%lu", (unsigned long)offsetof(struct {char x__; t (y__); }, y__)
module Numeric.GSL.Special.Internal (
createSFR,
create2SFR,
createSFR_E10,
Precision(..),
Gsl_mode_t,
Size_t,
precCode
)
where
import Foreign.Storable
import Foreign.Ptr
import Foreign.Marshal
import System.IO.Unsafe(unsafePerformIO)
import Numeric.LinearAlgebra.Devel(check,(//))
import Foreign.C.Types
data Precision = PrecDouble | PrecSingle | PrecApprox
precCode :: Precision -> Int
precCode PrecDouble = 0
precCode PrecSingle = 1
precCode PrecApprox = 2
type Gsl_mode_t = Int
type Size_t = CSize
---------------------------------------------------
data Gsl_sf_result = SF Double Double
deriving (Show)
instance Storable Gsl_sf_result where
sizeOf _ = #size gsl_sf_result
alignment _ = #alignment gsl_sf_result
peek ptr = do
val <- (#peek gsl_sf_result, val) ptr
err <- (#peek gsl_sf_result, err) ptr
return (SF val err)
poke ptr (SF val err) = do
(#poke gsl_sf_result, val) ptr val
(#poke gsl_sf_result, err) ptr err
data Gsl_sf_result_e10 = SFE Double Double CInt
deriving (Show)
instance Storable Gsl_sf_result_e10 where
sizeOf _ = #size gsl_sf_result_e10
alignment _ = #alignment gsl_sf_result_e10
peek ptr = do
val <- (#peek gsl_sf_result_e10, val) ptr
err <- (#peek gsl_sf_result_e10, err) ptr
e10 <- (#peek gsl_sf_result_e10, e10) ptr
return (SFE val err e10)
poke ptr (SFE val err e10) = do
(#poke gsl_sf_result_e10, val) ptr val
(#poke gsl_sf_result_e10, err) ptr err
(#poke gsl_sf_result_e10, e10) ptr e10
----------------------------------------------------------------
-- | access to one sf_result
createSFR :: String -> (Ptr a -> IO CInt) -> (Double, Double)
createSFR s f = unsafePerformIO $ do
p <- malloc :: IO (Ptr Gsl_sf_result)
f (castPtr p) // check s
SF val err <- peek p
free p
return (val,err)
----------------------------------------------------------------
-- | access to two sf_result's
create2SFR :: String -> (Ptr a -> Ptr a -> IO CInt) -> ((Double, Double),(Double, Double))
create2SFR s f = unsafePerformIO $ do
p1 <- malloc :: IO (Ptr Gsl_sf_result)
p2 <- malloc :: IO (Ptr Gsl_sf_result)
f (castPtr p1) (castPtr p2) // check s
SF val1 err1 <- peek p1
SF val2 err2 <- peek p2
free p1
free p2
return ((val1,err1),(val2,err2))
---------------------------------------------------------------------
-- the sf_result_e10 contains two doubles and the exponent
-- | access to sf_result_e10
createSFR_E10 :: String -> (Ptr a -> IO CInt) -> (Double, Int, Double)
createSFR_E10 s f = unsafePerformIO $ do
p <- malloc :: IO (Ptr Gsl_sf_result_e10)
f (castPtr p) // check s
SFE val err expo <- peek p
free p
return (val, fromIntegral expo,err)
|
