GenericPolynomial

A univariate, sparse, generic polynomial arithmetic library. That is, a polynomial in only one indeterminate, X, that only tracks terms with non-zero coefficients. This generic implementation has been tested and supports performing arithmetic on numeric types such as BigInteger, Complex, Decimal, Double, BigComplex, BigDecimal, BigRational, Int32, Int64 and more.

All arithmetic is done symbolically. That means the result of an arithmetic operation on two polynomials is another polynomial, not the result of evaluating those two polynomials and performing arithmetic on the results.

Generic Arithmetic Type Polynomial

• All classes and methods support a generic type T, which can be any data type that has arithmetic operator overloads. Explicitly, the following types have tests exercising all the standard arithmetic operations (Addition, Subtraction, Multiplication, Division, Exponentiation, Modulus, Square Root, Equality and Comparison (where applicable--Complex numbers are not orderable)) and so are well supported and come with verifiable proof that they work:

  • Every .NET CLR value type (Byte, Int16, Int32, Int64, UInt16, UInt32, UInt64, Float, Double, Decimal)
  • System.Numerics.BigInteger
  • System.Numerics.Complex
  • ExtendedNumerics.BigComplex (My arbitrary-precision complex number type library)
  • ExtendedNumerics.BigDecimal (My arbitrary-precision floating-point number type library.)
  • ExtendedNumerics.BigRational(My arbitrary precision rational number type library.)
  • ExtendedNumerics.Fraction (My arbitrary-precision rational number type library.)
  • Any type that has the Addition, Multiplication and Exponentiation operator overloads (at a minimum).

• Supports symbolic univariate polynomial (generic) arithmetic including:

  • Addition
  • Subtraction
  • Multiplication
  • Division
  • Modulus
  • Exponentiation
  • GCD of polynomials
  • Derivative
  • Integral
  • Reciprocal
  • Irreducibility checking
  • Polynomial evaluation by assigning to the invariant (X in this case) a value.

Polynomial Rings over a Finite Field

• Polynomial.Field supports addition, multiplication, division/modulus and inverse of a polynomial ring over a finite field. These operations do not support Complex, BigComplex, BigDecimal, or BigRational types.

  • What this effectively means in less-technical terms is that the polynomial arithmetic is performed in the usual way, but the result is then taken modulus two things:
    • Modulo an integer: All coefficients are reduced modulo an integer.
    • Modulo a polynomial: The whole polynomial is reduced modulo another, smaller, polynomial. This notion works much the same as regular modulus; The modulus polynomial, let's call it g, is declared to be equivalent to zero, and so every multiple of g is reduced to zero. You can think of it this way (although this is not how it's actually carried out): From a large polynomial, g is repeatedly subtracted from that polynomial until it can't subtract g anymore without crossing zero. The result is a polynomial that lies between 0 and g. Just like regular modulus, the result is always less than your modulus, or zero if the polynomial was a multiple of the modulus.
    • Effectively forms a quotient ring

• You can instantiate a polynomial in various ways:

  • From a string
    • This is the most massively-useful way and is the quickest way to start working with a particular polynomial you had in mind.
  • From its roots (Not all types supported)
    • Build a polynomial that has, as its roots, all of the numbers in the supplied array. If you want multiplicity of roots, include that number in the array multiple times.
  • From the base-m expansion of a number
    • Given a large number and a radix (base), call it m, a polynomial will be generated that is that number represented in the number base m.

• Other methods of interest that are related to, but not necessarily performed on a polynomial (Not all types supported):

  • Euler's Criterion
  • Legendre Symbol and Legendre Symbol Search
  • Tonelli-Shanks
  • Chinese Remainder Theorem
  • Polynomial evaluation by assigning to the invariant (X in this case) a value.

Other polynomial projects & numeric types

I've written a number of other polynomial implementations and numeric types catering to various specific scenarios. Depending on what you're trying to do, another implementation of this same library might be more appropriate. All of my polynomial projects should have feature parity, where appropriate¹.

• GenericArithmetic - A core math library. Its a class of static methods that allows you to perform arithmetic on an arbitrary numeric type, represented by the generic type T, who's concrete type is decided by the caller. This is implemented using System.Linq.Expressions and reflection to resolve the type's static overloadable operator methods at runtime, so it works on all the .NET numeric types automagically, as well as any custom numeric type, provided it overloads the numeric operators and standard method names for other common functions (Min, Max, Abs, Sqrt, Parse, Sign, Log, Round, etc.). Every generic arithmetic class listed below takes a dependency on this class.

• Polynomial - The original. A univariate polynomial that uses System.Numerics.BigInteger as the indeterminate type.

• GenericPolynomial - A univariate polynomial library that allows the indeterminate to be of an arbitrary type, as long as said type implements operator overloading. This is implemented dynamically, at run time, calling the operator overload methods using Linq.Expressions and reflection.

• CSharp11Preview.GenericMath.Polynomial - A univariate polynomial library that allows the indeterminate to be of an arbitrary type, but this version is implemented using C# 11's new Generic Math via static virtual members in interfaces.

• MultivariatePolynomial - A multivariate polynomial (meaning more than one indeterminate, e.g. 2XY^2) which uses BigInteger as the type for the indeterminates.
• GenericMultivariatePolynomial - A multivariate polynomial that allows the indeterminates to be of [the same] arbitrary type. GenericMultivariatePolynomial is to MultivariatePolynomial what GenericPolynomial is to Polynomial, and indeed is implemented using the same strategy as GenericPolynomial (i.e. dynamic calling of the operator overload methods at runtime using Linq.Expressions and reflection).

• ComplexPolynomial - A univariate polynomial library that has System.Numerics.Complex type indeterminates.
• ComplexMultivariatePolynomial - A multivariate polynomial library that has System.Numerics.Complex indeterminates.

• BigDecimal - An arbitrary precision, base-10 floating point number class.
• BigRational - Encodes a numeric value as an Integer + Fraction
• BigComplex - Essentially the same thing as System.Numerics.Complex but uses a System.Numerics.BigInteger type for the real and imaginary parts instead of a double.

• IntervalArithmetic. Instead of representing a value as a single number, interval arithmetic represents each value as a mathematical interval, or range of possibilities, [a,b], and allows the standard arithmetic operations to be performed upon them too, adjusting or scaling the underlying interval range as appropriate. See Wikipedia's article on Interval Arithmetic for further information.
• GNFS - A C# reference implementation of the General Number Field Sieve algorithm for the purpose of better understanding the General Number Field Sieve algorithm.

Footnotes

1. For example, the ComplexPolynomial implementation may be missing certain operations (namely: Irreducibility), because such a notion does not make sense or is ill defined in the context of complex numbers). ↩