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
realToFrac
function goes throughRational
, which loses information asRational
cannot represent all floating point values.Making floating types an instance of
Ord
makes no sense.Ord
is for totally ordered data types, which floats are not. As a result, a number of library functions (such asmax
andsort
) produce nonsensical results.The
Enum
instance for floating types similarly makes little sense. WhilefromEnum
andtoEnum
functions useInt
instead ofInteger
, limiting their usefulness,pred
andsucc
can be defined in a much more useful way.Functions that should care about negative zeros, such as
signum
andabs
, 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
Double
toCDouble
in 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.