I'm trying to implement a data structure where if I had the use of infinity for numerical comparison purposes, it would si开发者_JAVA技巧mply things greatly. Note this isn't maxBound/minBound, because a value can be <= maxbound, but all values would be < infinity.
No hope?
Well how about that! It turns out if you just type 1/0
it returns Infinity
! On ghci:
Prelude> 1/0
Infinity
Prelude> :t 1/0
1/0 :: (Fractional t) => t
Prelude> let inf=1/0
Prelude> filter (>=inf) [1..]
and then of course it runs forever, never finding a number bigger than infinity. (But see ephemient's comments below on the actual behavior of [1..]
)
infinity = read "Infinity"
Maybe you want a Maybe type?
data Infinite a = Infinite | Only a
then write a Num instance for Num a => Infinite a, with the numeric rules you need.
Try something like this. However, to get Num
operations (like +
or -
) you will need to define Num
instance for Infinitable a
type. Just like I've done it for Ord
class.
data Infinitable a = Regular a | NegativeInfinity | PositiveInfinity deriving (Eq, Show)
instance Ord a => Ord (Infinitable a) where
compare NegativeInfinity NegativeInfinity = EQ
compare PositiveInfinity PositiveInfinity = EQ
compare NegativeInfinity _ = LT
compare PositiveInfinity _ = GT
compare _ PositiveInfinity = LT
compare _ NegativeInfinity = GT
compare (Regular x) (Regular y) = compare x y
main =
let five = Regular 5
pinf = PositiveInfinity::Infinitable Integer
ninf = NegativeInfinity::Infinitable Integer
results = [(pinf > five), (ninf < pinf), (five > ninf)]
in
do putStrLn (show results)
λ: let infinity = (read "Infinity")::Double
λ: infinity > 1e100
True
λ: -infinity < -1e100
True
Take a look at my RangedSets library, which does exactly this in a very general way. I defined a "Boundary" type so that a value of type "Boundary a" is always either above or below any given "a". Boundaries can be "AboveAll", "BelowAll", "Above x" and "Below x".
You can work with the ieee754
package which offers an IEEE
typeclass. This typeclass has an infinity
member. The typeclass is implemented for Float
, Double
, CFloat
and CDouble
.
You thus install the ieee-754
package, and then obtain positive infinity (or negative infinity) with:
ghci> import Numeric.IEEE
ghci> infinity :: Double
Infinity
ghci> -infinity :: Double
-Infinity
For Float
s, Double
s, etc. it is implemented as [src]:
infinity = 1/0
If your use case is that you have boundary conditions that sometimes need to be checked, but sometimes not, you can solve it like this:
type Bound a = Maybe a
withinBounds :: (Num a, Ord a) => Bound a -> Bound a -> a -> Bool
withinBounds lo hi v = maybe True (<=v) lo && maybe True (v<=) hi
There is a more principled approach based on an idea from non-standard analysis. Given a totally ordered ring R of characteristic zero, you can consider the Laurent ring R[inf,1/inf] with the natural lexicographic total ordering. For example, you have:
for all x>0 in R,
.. -inf < -x < -d < -d^2 < .. < 0 < .. < d^2 < d < x < inf < inf^2 < ..
where d = 1/inf.
This way the Laurent ring R[inf,1/inf] is again a totally ordered Z-algebra, i.e. an instance of Num
, with other niceties you possibly want, including +/-infinity, +/-infinitesimal, second-order infinitesimals, etc.. But note that it's not Archimedian and induction will no longer work, which is a sort of second-order arithmetic. For implementation take a look at this example. As in the comment in the code this construction should work for other algebras, such as the list monad. You can think of lists where two elements are "infinitely close" "second-order infinitely far away" etc. (which leads to a generalization of rose trees.)
精彩评论