next | previous | forward | backward | up | top | index | toc | Macaulay2 website
Macaulay2Doc :: Matrix // Matrix

Matrix // Matrix -- factor a map through another

Synopsis

Description

If f is a matrix, and g is a matrix or Gröbner basis, then quotient(f,g) is an alternate notation for f//g.

If either f or g is a ring element, then it is taken to be a scalar matrix acting on H. If both are ring elements, then the result is also a ring element. If g is a MonomialIdeal, then it is taken to be the matrix of generators of g. Finally, if g is a GroebnerBasis object, then the Gröbner basis as so far computed is used. In these latter two cases, no Gröbner bases will be computed.

The resulting matrix h is such that f - g*h is the reduction of f modulo a Gröbner basis for the image of g.

If the remainder f - g*h is zero, then the quotient f//g satisfies the equation f === g * (f//g) and the quotient g\\f satisfies the equation f === g * (g\\f).

One common use is the following. If an ideal contains 1, then we may write 1 in terms of the generators of the ideal. First we make an ideal.

i1 : A = ZZ/101[x,y,z]

o1 = A

o1 : PolynomialRing
i2 : F = x^4 - y*z*(1-x)^2 - z - y^3

      4    2       3
o2 = x  - x y*z - y  + 2x*y*z - y*z - z

o2 : A
i3 : I = ideal(F,diff(x,F),diff(y,F),diff(z,F))

             4    2       3                       3                     2      2                2
o3 = ideal (x  - x y*z - y  + 2x*y*z - y*z - z, 4x  - 2x*y*z + 2y*z, - x z - 3y  + 2x*z - z, - x y + 2x*y - y - 1)

o3 : Ideal of A
Transposing the (row) matrix of generators of the ideal puts the generators on separate lines and shows the degrees.
i4 : transpose gens I

o4 = {-4} | x4-x2yz-y3+2xyz-yz-z |
     {-3} | 4x3-2xyz+2yz         |
     {-3} | -x2z-3y2+2xz-z       |
     {-3} | -x2y+2xy-y-1         |

             4       1
o4 : Matrix A  <--- A
Next we test whether 1 is in the ideal.
i5 : 1 % I

o5 = 0

o5 : A
We see that 1 is in the ideal. Now we represent 1 in terms of the generators of I.
i6 : h = 1 // gens I

o6 = {4} | -46x3y2-4x2y3+11x4z+4x3yz-2x2y2z+22x3z2+50x2yz2+8x3y+36x2y2-xy3+14x3z+48xy2z+11x2z2-26xyz2-12x3-49x2y-7xy2+5y3-44x2z+1
     {3} | -39x4y2+x3y3-28x5z-x4yz-50x3y2z+45x4z2+38x3yz2-2x4y-9x3y2-25x2y3+47x4z-12x2y2z-28x3z2-44x2yz2+3x4-13x3y+8x2y2+24xy3+11
     {3} | 49x3y3+35x2y4+30x4yz-35x3y2z-33x2y3z-41x3yz2+17x2y2z2+31x3y2-12x2y3+34xy4+29x3yz-16xy3z+30x2yz2-25xy2z2+4x3y+50x2y2+36
     {3} | -26x3y2z-33x2y3z+15x4z2+33x3yz2+34x2y2z2+30x3z3-42x2yz3-35x3yz-29x2y2z+15xy3z+20x3z2+2x2yz2-9xy2z2+26x2z3-38xyz3+25x3y
     ----------------------------------------------------------------------------------------------------------------------------
     2xyz-16y2z-13xz2-46yz2+19x2-28xy-y2+22xz+9yz-15z2-49x-11y-35z+28                                                            
     x3z-32x2yz+41xy2z-22x2z2+4xyz2-30x3+7x2y-12xy2+45x2z-17xyz+21y2z+29xz2+2yz2-13x2+28xy-20y2-17xz+26yz-7x-19y-47z             
     xy3+32y4-19x2yz-4xy2z+39y3z+38xyz2+49y2z2-40x2y+43xy2+34y3-26x2z+27xyz-3y2z+49xz2+5yz2-2x2+50xy-30y2+15xz-38yz-40z2+35x-43y+
     -13x3z-25x2yz-50xy2z-37y3z-27x2z2-43xyz2+19y2z2-15xz3+yz3+x2y-14x2z+49xyz+7y2z+15xz2-39yz2-46z3-2xy+47xz+20yz+9z2+y+6z-1    
     ----------------------------------------------------------------------------------------------------------------------------
            |
            |
     26z-34 |
            |

             4       1
o6 : Matrix A  <--- A
i7 : gens I * h

o7 = | 1 |

             1       1
o7 : Matrix A  <--- A

See also