numeric-prelude-0.4.2: An experimental alternative hierarchy of numeric type classes

Safe HaskellSafe-Inferred
LanguageHaskell98

Algebra.IntegralDomain

Contents

Synopsis

Class

class C a => C a where Source

IntegralDomain corresponds to a commutative ring, where a mod b picks a canonical element of the equivalence class of a in the ideal generated by b. div and mod satisfy the laws

                        a * b === b * a
(a `div` b) * b + (a `mod` b) === a
              (a+k*b) `mod` b === a `mod` b
                    0 `mod` b === 0

Typical examples of IntegralDomain include integers and polynomials over a field. Note that for a field, there is a canonical instance defined by the above rules; e.g.,

instance IntegralDomain.C Rational where
    divMod a b =
       if isZero b
         then (undefined,a)
         else (a\/b,0)

It shall be noted, that div, mod, divMod have a parameter order which is unfortunate for partial application. But it is adapted to mathematical conventions, where the operators are used in infix notation.

Minimal definition: divMod or (div and mod)

Minimal complete definition

divMod | div, mod

Methods

div, mod :: a -> a -> a infixl 7 `div`, `mod` Source

divMod :: a -> a -> (a, a) Source

Instances

C Int 
C Int8 
C Int16 
C Int32 
C Int64 
C Integer 
C Word 
C Word8 
C Word16 
C Word32 
C Word64 
C T 
(Ord a, C a) => C (T a) 
(C a, C a) => C (T a) 
Integral a => C (T a) 
(C a, C a) => C (T a)

The C instance is intensionally built from the C structure of the polynomial coefficients. If we would use Integral.C a superclass, then the Euclidean algorithm could not determine the greatest common divisor of e.g. [1,1] and [2].

C a => C (T a) 
C a => C (T a) 
(Ord a, C a, C a) => C (T a)

divMod is implemented in terms of divModStrict. If it is needed we could also provide a function that accesses the divisor first in a lazy way and then uses a strict divisor for subsequent rounds of the subtraction loop. This way we can handle the cases "dividend smaller than divisor" and "dividend greater than divisor" in a lazy and efficient way. However changing the way of operation within one number is also not nice.

Derived functions

divModZero :: (C a, C a) => a -> a -> (a, a) Source

Allows division by zero. If the divisor is zero, then the dividend is returned as remainder.

divides :: (C a, C a) => a -> a -> Bool Source

sameResidueClass :: (C a, C a) => a -> a -> a -> Bool Source

divChecked :: (C a, C a) => a -> a -> a Source

Returns the result of the division, if divisible. Otherwise undefined.

safeDiv :: (C a, C a) => a -> a -> a Source

Deprecated: use divChecked instead

Returns the result of the division, if divisible. Otherwise undefined.

even :: (C a, C a) => a -> Bool Source

odd :: (C a, C a) => a -> Bool Source

divUp :: C a => a -> a -> a Source

divUp n m is similar to div but it rounds up the quotient, such that divUp n m * m = roundUp n m.

roundDown :: C a => a -> a -> a Source

roundDown n m rounds n down to the next multiple of m. That is, roundDown n m is the greatest multiple of m that is at most n. The parameter order is consistent with div and friends, but maybe not useful for partial application.

roundUp :: C a => a -> a -> a Source

roundUp n m rounds n up to the next multiple of m. That is, roundUp n m is the greatest multiple of m that is at most n.

Algorithms

decomposeVarPositional :: (C a, C a) => [a] -> a -> [a] Source

decomposeVarPositional [b0,b1,b2,...] x decomposes x into a positional representation with mixed bases x0 + b0*(x1 + b1*(x2 + b2*x3)) E.g. decomposeVarPositional (repeat 10) 123 == [3,2,1]

decomposeVarPositionalInf :: C a => [a] -> a -> [a] Source

Properties

propInverse :: (Eq a, C a, C a) => a -> a -> Property Source

propMultipleDiv :: (Eq a, C a, C a) => a -> a -> Property Source

propMultipleMod :: (Eq a, C a, C a) => a -> a -> Property Source

propProjectAddition :: (Eq a, C a, C a) => a -> a -> a -> Property Source

propProjectMultiplication :: (Eq a, C a, C a) => a -> a -> a -> Property Source

propUniqueRepresentative :: (Eq a, C a, C a) => a -> a -> a -> Property Source

propSameResidueClass :: (Eq a, C a, C a) => a -> a -> a -> Property Source