File: Alternative.hs

package info (click to toggle)
haskell-integer-conversion 0.1.0.1-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 84 kB
  • sloc: haskell: 294; makefile: 3
file content (111 lines) | stat: -rw-r--r-- 3,406 bytes parent folder | download
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
{-# LANGUAGE BangPatterns #-}
-- | A sub-quadratic algorithm for conversion of digits into 'Integer'.
-- Pairs of adjacent radix @b@ digits are combined into a single radix @b^2@ digit.
-- This process is repeated until we are left with a single digit.
-- This algorithm performs well only on large inputs,
-- so we use the simple algorithm for smaller inputs.
--
-- This implementation is taken from aeson-2.1.
module Alternative (
    byteStringToInteger,
    textToInteger,
    stringToInteger,
) where

import Data.Char (ord)
import Data.Word (Word8)

import qualified Data.ByteString as BS
import qualified Data.List       as L
import qualified Data.Text       as T

byteStringToInteger :: BS.ByteString -> Integer
byteStringToInteger bs
    -- here (and similarly in 'textToInteger') it could make sense
    -- to do first loop directly on 'ByteString' (or 'Text'),
    -- but as this is already a slow path, we opt rather for a simpler implementation.
    | l > 40    = valInteger' 10 l [ fromWord8 w | w <- BS.unpack bs ]
    | otherwise = byteStringToIntegerSimple bs
  where
    !l = BS.length bs

byteStringToIntegerSimple :: BS.ByteString -> Integer
byteStringToIntegerSimple = BS.foldl' step 0 where
  step a b = a * 10 + fromWord8 b

textToInteger :: T.Text -> Integer
textToInteger bs
    | l > 40    = valInteger' 10 l [ fromChar w | w <- T.unpack bs ]
    | otherwise = textToIntegerSimple bs
  where
    !l = T.length bs

textToIntegerSimple :: T.Text -> Integer
textToIntegerSimple = T.foldl' step 0 where
  step a b = a * 10 + fromChar b

stringToInteger :: String -> Integer
stringToInteger s
    | l > 40    = valInteger' 10 l (map fromChar s)
    | otherwise = stringToIntegerSimple s
  where
    !l = length s

stringToIntegerSimple :: String -> Integer
stringToIntegerSimple = L.foldl' step 0 where
  step a b = a * 10 + fromChar b

fromChar :: Char -> Integer
fromChar c = toInteger (ord c - 48 :: Int)
{-# INLINE fromChar #-}

fromWord8 :: Word8 -> Integer
fromWord8 w = toInteger (fromIntegral w - 48 :: Int)
{-# INLINE fromWord8 #-}

-- | A sub-quadratic algorithm.
--
-- Call 'valInteger'' directly if you know length of @digits@ in advance.
-- valInteger :: Integer -> [Integer] -> Integer
-- valInteger base ds = valInteger' base (length ds) ds

-- | A sub-quadratic algorithm implementation.
valInteger'
    :: Integer    -- ^ base
    -> Int        -- ^ length of digits
    -> [Integer]  -- ^ digits
    -> Integer
valInteger' = go
  where
    go :: Integer -> Int -> [Integer] -> Integer
    go _ _ []  = 0
    go _ _ [d] = d
    go b l ds
        | l > 40 = b' `seq` go b' l' (combine b ds')
        | otherwise = valIntegerSimple b ds
      where
        -- ensure that we have an even number of digits
        -- before we call combine:
        ds' = if even l then ds else 0 : ds
        b' = b * b
        l' = (l + 1) `quot` 2

    combine b (d1 : d2 : ds) = d `seq` (d : combine b ds)
      where
        d = d1 * b + d2
    combine _ []  = []
    combine _ [_] = errorWithoutStackTrace "this should not happen"

-- | The following algorithm is only linear for types whose Num operations
-- are in constant time.
--
-- We export this (mostly) for testing purposes.
--
valIntegerSimple :: Integer -> [Integer] -> Integer
valIntegerSimple base = go 0
  where
    go r [] = r
    go r (d : ds) = r' `seq` go r' ds
      where
        r' = r * base + fromIntegral d