Полной точности суммирования в Haskell


Я пытаюсь ускорить следующий код для суммирования парных с полной точностью, основанной на языке Python рецепт от Raymond Hettinger. Любые мысли о том, как сделать это быстрее?

-- Full precision summation based on http://code.activestate.com/recipes/393090/
msum :: [Double] -> Double
msum [] = 0
msum (x:xs) = sum $ foldl' (\ps x' -> reverse (partials x' ps)) [x] xs

partials :: Double -> [Double] -> [Double]
partials x = foldl' (\acc p -> case hilo x p of (hi, 0) -> hi:acc
                                                (hi, lo) -> hi:lo:acc) []

hilo :: Double -> Double -> (Double, Double)
hilo x y | abs x < abs y = hilo y x
         | otherwise = (hi, lo) where hi = x + y
                                      lo = y - (hi - x)

Это уже второй проект, можно увидеть предыдущую версию кода в этом суть.



553
6
задан 12 апреля 2011 в 11:04 Источник Поделиться
Комментарии
1 ответ

Вы могли бы написать хило так не хило :: двойной -> двойной -> [двойной], который возвращает [привет,Ло] для Ло /= 0 и [привет] для Ло = 0. Затем вы можете написать частичными качестве

partials x = foldl' (\acc p -> (hilo x p) ++ acc) -> hi:acc []

Кстати: как мог Ло стать ненулевой учитывая

lo = y - (hi - x) = y - (x + y - x) = y - y = 0

Закругление?

2
ответ дан 14 апреля 2011 в 02:04 Источник Поделиться