Description
Using the algorithms in Logar-Sturmfels and Fabianska-Quadrat, this package computes a free basis of a projective module over a polynomial ring with coefficients in the rationals, integers, or a finite field. It also provides methods to solve related problems involving completing a unimodular matrix to a square invertible matrix over a polynomial ring with coefficients in the rationals, integers, or a finite field, or a Laurent polynomial ring with coefficients in the rationals or a finite field.
For mathematical background and applications, see
- A. Fabianska. Algorithmic analysis of presentations of groups and modules. http://darwin.bth.rwth-aachen.de/opus/volltexte/2009/2950/, Jan 2009.
- T. Y. Lam. Serre's problem on projective modules. Springer Monographs in Mathematics. Springer-Verlag, Berlin, 2006.
- A. Logar and B. Sturmfels. Algorithms for the Quillen-Suslin theorem. J. Algebra, 145(1):231-239, 1992.
- A. Fabianska and A. Quadrat. Applications of the Quillen-Suslin theorem to multidimensional systems theory. Grobner bases in control theory and signal processing. Radon Series Comp. Appl. Math (3):23-106, 2007.
A first example computes a free basis for the kernel of a unimodular row.
i1 : R = QQ[x,y]
o1 = R
o1 : PolynomialRing
|
i2 : f = matrix{{x^2*y+1,x+y-2,2*x*y}}
o2 = | x2y+1 x+y-2 2xy |
1 3
o2 : Matrix R <-- R
|
i3 : P = ker f
o3 = image {3} | 0 2x+2y-4 2y2-4y |
{1} | 2xy -2x2y-2xy2+4xy-2 -2xy3+4xy2-2y |
{2} | -x-y+2 xy+y2-2x-4y+4 y3-4y2+4y+1 |
3
o3 : R-module, submodule of R
|
i4 : B = computeFreeBasis P
o4 = {3} | -2x-2y+4 -2y2+4y |
{1} | 2x2y+2xy2-4xy+2 2xy3-4xy2+2y |
{2} | -xy-y2+2x+4y-4 -y3+4y2-4y-1 |
3 2
o4 : Matrix R <-- R
|
i5 : image B == P
o5 = true
|