Fermat (computer algebra system)
From Wikipedia, the free encyclopedia
| Fermat | |
| Developer: | Robert H. Lewis |
|---|---|
| Latest release: | 3.8.1 / March 13, 2007 |
| OS: | Mac OS X, Mac OS, Linux, Unix, Windows |
| Use: | computer algebra system |
| License: | shareware |
| Website: | http://www.bway.net/~lewis |
Fermat is a shareware program developed by Prof. Robert H. Lewis of Fordham University. It is a computer algebra system, in which items being computed can be integers (of arbitrary size), rational numbers, real numbers, complex numbers, modular numbers, finite field elements, multivariable polynomials, rational functions, or polynomials modulo other polynomials. The main areas of application are multivariate rational function arithmetic and matrix algebra over rings of multivariate polynomials or rational functions. Fermat does not do simplification of transcendental functions or symbolic integration.
A session with Fermat usually starts by choosing rational or modular "mode" to establish the ground field (or ground ring) F as Z or Z/n. On top of this may be attached any number of symbolic variables t1,t2,...,tn, thereby creating the polynomial ring F[t1,t2,...,tn] and its quotient field. Further, some polynomials p,q,... involving some of the ti can be chosen to mod out with, creating the quotient ring F(t1,t2,...) / (p,q,...). Finally, it is possible to allow Laurent polynomials, those with negative as well as positive exponents. Once the computational ring is established in this way, all computations are of elements of this ring. The computational ring can be changed later in the session.
The polynomial gcd procedures, which call each other in a highly recursive manner, are about 7000 lines of code.
Fermat has extensive built-in primitives for array and matrix manipulations, such as submatrix, sparse matrix, determinant, normalize, column reduce, row echelon, Smith form, and matrix inverse. It is consistently faster than some well known computer algebra systems, especially in multivariate polynomial gcd. It is also space efficient.
The basic data item in Fermat is a multivariate rational function or quolynomial. The numerator and denominator are polynomials with no common factor. Polynomials are implemented recursively as general linked lists, unlike some systems that implement polynomials as lists of monomials. To implement (most) finite fields, the user finds an irreducible monic polynomial in a symbolic variable, say p(t1), and commands Fermat to mod out by it. This may be continued recursively, q(t2,t1), etc. Low level data structures are set up to facilitate arithmetic and gcd over this newly created ground field. Two special fields, GF(28) and GF(216), are more efficiently implemented at the bit level.
To help implement the Dixon resultant technique, special features have been added to the determinant function. These provide a dramatic increase in the speed of resultant calculations with systems of polynomial equations that exhibit symmetry.
Fermat provides a complete programming language. Programs and data can be saved to an ordinary text file that can be examined as such, read during a later session, or read by some other software system.
[edit] History
In an earlier version, called FFermat, the basic number type is floating point numbers of 18 digits. That version allows for numerical computing techniques, has extensive graphics capabilities, no sophisticated polynomial gcd algorithms, and is available only for Mac OS.
Fermat was originally written in Pascal for a DEC Vax, then for Mac OS during 1985 - 1996. It was ported to Microsoft Windows in 1998. In 2003 it was translated into C and ported to Linux (Intel machines) and Unix (Sparc/Sun). It is about 88,000 lines of C code.
