numeric-prelude-0.4.3.3: An experimental alternative hierarchy of numeric type classes
Safe HaskellNone
LanguageHaskell98

MathObj.Polynomial.Core

Description

This module implements polynomial functions on plain lists. We use such functions in order to implement methods of other datatypes.

The module organization differs from that of ResidueClass: Here the Polynomial module exports the type that fits to the NumericPrelude type classes, whereas in ResidueClass the sub-modules export various flavors of them.

Synopsis

Documentation

horner :: C a => a -> [a] -> a Source #

Horner's scheme for evaluating a polynomial in a ring.

hornerCoeffVector :: C a v => a -> [v] -> v Source #

Horner's scheme for evaluating a polynomial in a module.

hornerArgVector :: (C a v, C v) => v -> [a] -> v Source #

normalize :: C a => [a] -> [a] Source #

It's also helpful to put a polynomial in canonical form. normalize strips leading coefficients that are zero.

shift :: C a => [a] -> [a] Source #

Multiply by the variable, used internally.

unShift :: [a] -> [a] Source #

equal :: (Eq a, C a) => [a] -> [a] -> Bool Source #

add :: C a => [a] -> [a] -> [a] Source #

sub :: C a => [a] -> [a] -> [a] Source #

negate :: C a => [a] -> [a] Source #

scale :: C a => a -> [a] -> [a] Source #

collinear :: (Eq a, C a) => [a] -> [a] -> Bool Source #

tensorProduct :: C a => [a] -> [a] -> [[a]] Source #

\(QC.NonEmpty xs) (QC.NonEmpty ys) -> PolyCore.tensorProduct xs ys == List.transpose (PolyCore.tensorProduct ys (intPoly xs))

tensorProductAlt :: C a => [a] -> [a] -> [[a]] Source #

mul :: C a => [a] -> [a] -> [a] Source #

mul is fast if the second argument is a short polynomial, ** relies on that fact.

mulShear :: C a => [a] -> [a] -> [a] Source #

\xs ys  ->  PolyCore.equal (intPoly $ PolyCore.mul xs ys) (PolyCore.mulShear xs ys)

mulShearTranspose :: C a => [a] -> [a] -> [a] Source #

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

\x y -> case (PolyCore.normalize x, PolyCore.normalize y) of (nx, ny) -> not (null (ratioPoly ny)) ==> mapSnd PolyCore.normalize (PolyCore.divMod nx ny) == mapPair (PolyCore.normalize, PolyCore.normalize) (PolyCore.divMod x y)
\x y -> not (isZero (ratioPoly y)) ==> let z = fst $ PolyCore.divMod (Poly.coeffs x) y in  PolyCore.normalize z == z
\x y -> case PolyCore.normalize $ ratioPoly y of ny -> not (null ny) ==> List.length (snd $ PolyCore.divMod x y) < List.length ny

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

The modulus will always have one element less than the divisor. This means that the modulus will be denormalized in some cases, e.g. mod [2,1,1] [1,1,1] == [1,0] instead of [1].

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

progression :: C a => [a] Source #

differentiate :: C a => [a] -> [a] Source #

integrate :: C a => a -> [a] -> [a] Source #

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

Integrates if it is possible to represent the integrated polynomial in the given ring. Otherwise undefined coefficients occur.

mulLinearFactor :: C a => a -> [a] -> [a] Source #

alternate :: C a => [a] -> [a] Source #

dilate :: C a => a -> [a] -> [a] Source #

shrink :: C a => a -> [a] -> [a] Source #