rounded

Correctly-rounded arbitrary-precision floating-point arithmetic

https://github.com/ekmett/rounded

Version on this page:0.1.0.1
LTS Haskell 13.30:0.1.0.1
Stackage Nightly 2019-01-07:0.1.0.1
Latest on Hackage:1.1.1

See all snapshots rounded appears in

BSD-3-Clause licensed by Edward A. Kmett, Daniel G. Peebles, Claude Heiland-Allen
Maintained by Claude Heiland-Allen
This version can be pinned in stack with:rounded-0.1.0.1@sha256:1c5403ec50f4163877a39e155dc353e20f008063853b94000dbaa5aa9f989f86,3231
Used by 1 package in nightly-2018-12-17(full list with versions):

rounded

Build Status

This package provides properly rounded floating point numbers of arbitrary precision. It does so by wrapping the GNU MPFR library.

Phantom types carry the information about the precision and rounding mode, letting you treat properly rounded floating point numbers as instances of Num or Floating, like any other numeric type in Haskell.

Unlike other attempts to port MPFR to Haskell, this library does not require you to cripple Integer performance or link your code in an unnatural way.

Usage

{-# LANGUAGE DataKinds #-}
{-# LANGUAGE ScopedTypeVariables #-}
import Numeric.Rounded
import Data.Proxy

To use a 53 bit significand (the same size as used by a Double), and round down intermediate results:

>>> pi :: Rounded TowardZero Double
3.141592653589793

We can also round away from zero, or use other rounding modes.

>>> pi :: Rounded AwayFromZero Double
3.1415926535897936

We can specify the significand size directly using type literals in GHC:

>>> kCatalan :: Rounded TowardZero 128
0.915965594177219015054603514932384110773

You can also specify a dynamic significand size at runtime:

>>> reifyPrecision 512 (\(_ :: Proxy p) -> show (logBase 10 2 :: Rounded TowardNearest p))
"0.3010299956639811952137388947244930267681898814621085413104274611271081892744245094869272521181861720406844771914309953790947678811335235059996923337046956"

or a dynamic rounding mode:

ghci> reifyRounding TowardZero (\(_ :: Proxy r) -> show (logBase 10 2 :: Rounded r 512))
"0.30102999566398119521373889472449302676818988146210854131042746112710818927442450948692725211818617204068447719143099537909476788113352350599969233370469556"

Contact Information

Please, feel free to contact me with questions, concerns, or bug fixes.

I can be reached as ekmett via github or as edwardk on the #haskell IRC channel on irc.freenode.net.

-Edward Kmett

Changes

Revision history for rounded

0.1.0.1 – 2018-10-31

  • License changed from LGPL to BSD3.

0.1 – 2018-10-30

  • Bindings rewritten as a straight FFI binding without a patched MPFR.
  • Added Simple module.

0.0.1

  • Repository initialized based on older work on a package called precision
  • Daniel Peebles was able to get a patched MPFR to work