Alternative floating point support for GHC.
A replacement for the standard Haskell floating point types and supporting functions. There are a number of shortcomings which I feel severely hinder Haskell's utility for numerical computation. These shortcomings include
There is no way to sanely convert between Haskell's floating types -- not even between Double and CDouble. The implementation of the
realToFracfunction goes throughRational, which loses information asRationalcannot represent all floating point values.Making floating types an instance of
Ordmakes no sense.Ordis for totally ordered data types, which floats are not. As a result, a number of library functions (such asmaxandsort) produce nonsensical results.The
Enuminstance for floating types similarly makes little sense. WhilefromEnumandtoEnumfunctions useIntinstead ofInteger, limiting their usefulness,predandsucccan be defined in a much more useful way.Functions that should care about negative zeros, such as
signumandabs, do not.Some functions, such as
floor, have nonsensical behaviour for non-finite input.The selection of floating point library functions pales in comparison to C. This problem is made worse since, as noted above, it is impossible to convert losslessly from
DoubletoCDoublein order to use the FFI.There is no mechanism for handling rounding modes or exceptions.
This package is intended to address all of the above issues, and more. Also provided, for convenience, is an alternative to the standard Prelude which includes features from this library and the non-overlapping parts of the standard Prelude.