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
|
{-# LANGUAGE NoImplicitPrelude #-}
module Compiler.RealFloat where
import Compiler.Error
import Compiler.Num
import Compiler.RealFrac
import Compiler.Floating
import Data.Bool
import Data.Eq
import Data.Ord
import Data.Tuple
import Data.Function
import Foreign.Pair
class (RealFrac a, Floating a) => RealFloat a where
floatRadix :: a -> Integer
floatDigits :: a -> Int
floatRange :: a -> (Int, Int)
decodeFloat :: a -> (Integer, Int)
encodeFloat :: Integer -> Int -> a
exponent :: a -> Int
significand :: a -> a
scaleFloat :: Int -> a -> a
isNaN :: a -> Bool
isInfinite :: a -> Bool
isDenormalized :: a -> Bool
isNegativeZero :: a -> Bool
isIEEE :: a -> Bool
atan2 :: a -> a -> a
instance RealFloat Double where
floatRadix _ = 2
floatDigits _ = 53
floatRange _ = (-1021,1024)
decodeFloat x = decodeDouble x
encodeFloat sig exp = encodeDouble sig exp
exponent x = case decodeFloat x of (m,n) -> if m == 0 then 0 else n + floatDigits x
significand x = case decodeFloat x of (m,_) -> encodeFloat m (negate (floatDigits x))
scaleFloat 0 x = x
scaleFloat k x
| isFix = x
| otherwise = case decodeFloat x of
(m,n) -> encodeFloat m (n + clamp bd k)
where bd = snd (floatRange x) - fst (floatRange x) + 4*floatDigits x
isFix = x == 0 || (not $ isDoubleFinite x)
clamp :: Int -> Int -> Int
clamp bd k = max (-bd) (min bd k)
isNaN x = isDoubleNaN x
isInfinite x = isDoubleInfinite x
isDenormalized x = isDoubleDenormalized x
isNegativeZero x = isDoubleNegativeZero x
isIEEE _ = True
atan2 x y = atan2_double x y
foreign import bpcall "Real:" isDoubleNaN :: Double -> Bool
foreign import bpcall "Real:" isDoubleInfinite :: Double -> Bool
foreign import bpcall "Real:" isDoubleFinite :: Double -> Bool
foreign import bpcall "Real:" isDoubleDenormalized :: Double -> Bool
foreign import bpcall "Real:" isDoubleNegativeZero :: Double -> Bool
foreign import bpcall "Real:" atan2_double :: Double -> Double -> Double
foreign import bpcall "Real:" decodeDoubleRaw :: Double -> EPair Integer Int
foreign import bpcall "Real:" encodeDouble :: Integer -> Int -> Double
decodeDouble = pair_from_c . decodeDoubleRaw
|
