Back to Forum | View unanswered posts | View active topics
Topic review - Triangularization of underdetermined systems |
Author |
Message |
|
|
Post subject: |
Triangularization of underdetermined systems |
|
|
Is there a Singular library analogous to "triang.lib" which triangularizes underdetermined systems (ideals with dim > 0)? For example, > LIB "triang.lib"; ... > ring r=0,(x,y,z,w),lp; > ideal i = x^2*y+z^3+w*x*y, x*w+y^2*x,x+x^2*y+z*w; > ideal j = groebner(i); > j; j[1]=z11w-4z9w2-2z8w2+6z7w3+2z6w3-4z5w4+z5w3+z5+2z4w4+z3w5-2z2w5+zw8+zw5 j[2]=yzw18-4yzw15+6yzw12-4yzw9-yzw8+yzw6-2yzw5-yzw2+z10w7-2z10w4+z10w+z9w4+ z9w+z8w11-6z8w8+9z8w5-4z8w2-z7w8+z7w5-6z7w2+z6w15-6z6w12+15z6w9-16z6w6-2z6w5+ 6z6w3-2z6w2-z5w12+3z5w9-4z5w6+8z5w3-2z4w16+9z4w13-16z4w10+13z4w7+z4w6-4z4w4+ 2z4w3+z4-z3w16+3z3w13-2z3w10+z3w6-2z3w4+2z3w3+z3+z2w17-4z2w14+6z2w11-z2w10- 4z2w8+z2w5+z2w4-zw17+4zw14-6zw11+4zw8-zw5 j[3]=4yz2w2+yzw15-5yzw12+7yzw9-3yzw6-yzw5-yzw2+z10w4-3z10w+z9w+z8w8-7z8w5+ 12z8w2-z7w5+2z7w2+z6w12-7z6w9+18z6w6-18z6w3-2z6w2-z5w9+4z5w6-4z5w3-2z4w13+ 11z4w10-19z4w7+12z4w4+z4w3-3z4-z3w13+4z3w10-2z3w7-2z3w4+z3w3+z3+z2w14-5z2w11+ 7z2w8-z2w7-3z2w5+z2w4-zw14+5zw11-7zw8+3zw5 j[4]=yz3-yzw4-z6w+2z4w2+z3w2-z2w3+zw3 j[5]=y2zw+zw2 j[6]=xw3+x-yz3w+yzw2-z3+zw j[7]=xz3-xzw-y2z3-yzw3-z6+2z4w-z2w2 j[8]=xy+xw2-yz3+yzw j[9]=x2-xz3+xzw+xw2 > triangL(j); ? Error: ideal is not zero-dimensional. > dim(j); 2 To paraphrase my question, how can I reproduce the output given by Maple's gsolve? > ideal := [x^2*y+z^3+w*x*y, x*w+y^2*x,x+x^2*y+z*w]: > vars := [x,y,z,w]: > Groebner[gsolve](ideal,vars); 2 {[[z, y - y + 1, x - y + 1, w + y - 1], plex(w, x, y, z), {y - 1, y, y + 1}], [[x, z], plex(z, x, y, w), {}], [[ 2 3 3 5 3 2 8 5 5 2 5 2 3 -z - 2 z y - z y - y z + y - y z + y , -z + y + x - z y - z 4 2 6 3 3 5 6 7 2 9 6 2 4 7 - z y + z y - z y - z y + z x - 2 z y - z - z y + z y 5 4 2 6 3 3 4 5 2 7 4 2 - z y , -y z + y + x y - z y - y z + z x - 2 z y - z - z y 2 7 3 4 + z y - z y , 2 2 7 7 2 4 4 2 2 3 2 2 6 5 x + z x + z y + y - z y - 2 z y - 3 y z - z y - z - z , 2 w + y ], plex(w, x, y, z), {x, z}], [[z, y + 1, x - 1, w + 1], plex(w, x, y, z), {}]} email: jtshorwo@math.uwaterloo.caPosted in old Singular Forum on: 2003-05-27 14:32:34+02
Is there a Singular library analogous to "triang.lib" which triangularizes underdetermined systems (ideals with dim > 0)?
For example,
> LIB "triang.lib"; ... > ring r=0,(x,y,z,w),lp; > ideal i = x^2*y+z^3+w*x*y, x*w+y^2*x,x+x^2*y+z*w; > ideal j = groebner(i); > j; j[1]=z11w-4z9w2-2z8w2+6z7w3+2z6w3-4z5w4+z5w3+z5+2z4w4+z3w5-2z2w5+zw8+zw5 j[2]=yzw18-4yzw15+6yzw12-4yzw9-yzw8+yzw6-2yzw5-yzw2+z10w7-2z10w4+z10w+z9w4+ z9w+z8w11-6z8w8+9z8w5-4z8w2-z7w8+z7w5-6z7w2+z6w15-6z6w12+15z6w9-16z6w6-2z6w5+ 6z6w3-2z6w2-z5w12+3z5w9-4z5w6+8z5w3-2z4w16+9z4w13-16z4w10+13z4w7+z4w6-4z4w4+ 2z4w3+z4-z3w16+3z3w13-2z3w10+z3w6-2z3w4+2z3w3+z3+z2w17-4z2w14+6z2w11-z2w10- 4z2w8+z2w5+z2w4-zw17+4zw14-6zw11+4zw8-zw5 j[3]=4yz2w2+yzw15-5yzw12+7yzw9-3yzw6-yzw5-yzw2+z10w4-3z10w+z9w+z8w8-7z8w5+ 12z8w2-z7w5+2z7w2+z6w12-7z6w9+18z6w6-18z6w3-2z6w2-z5w9+4z5w6-4z5w3-2z4w13+ 11z4w10-19z4w7+12z4w4+z4w3-3z4-z3w13+4z3w10-2z3w7-2z3w4+z3w3+z3+z2w14-5z2w11+ 7z2w8-z2w7-3z2w5+z2w4-zw14+5zw11-7zw8+3zw5 j[4]=yz3-yzw4-z6w+2z4w2+z3w2-z2w3+zw3 j[5]=y2zw+zw2 j[6]=xw3+x-yz3w+yzw2-z3+zw j[7]=xz3-xzw-y2z3-yzw3-z6+2z4w-z2w2 j[8]=xy+xw2-yz3+yzw j[9]=x2-xz3+xzw+xw2 > triangL(j); ? Error: ideal is not zero-dimensional. > dim(j); 2
To paraphrase my question, how can I reproduce the output given by Maple's gsolve?
> ideal := [x^2*y+z^3+w*x*y, x*w+y^2*x,x+x^2*y+z*w]: > vars := [x,y,z,w]: > Groebner[gsolve](ideal,vars); 2 {[[z, y - y + 1, x - y + 1, w + y - 1], plex(w, x, y, z), {y - 1, y, y + 1}],
[[x, z], plex(z, x, y, w), {}], [[
2 3 3 5 3 2 8 5 5 2 5 2 3 -z - 2 z y - z y - y z + y - y z + y , -z + y + x - z y - z
4 2 6 3 3 5 6 7 2 9 6 2 4 7 - z y + z y - z y - z y + z x - 2 z y - z - z y + z y
5 4 2 6 3 3 4 5 2 7 4 2 - z y , -y z + y + x y - z y - y z + z x - 2 z y - z - z y
2 7 3 4 + z y - z y ,
2 2 7 7 2 4 4 2 2 3 2 2 6 5 x + z x + z y + y - z y - 2 z y - 3 y z - z y - z - z ,
2 w + y ], plex(w, x, y, z), {x, z}],
[[z, y + 1, x - 1, w + 1], plex(w, x, y, z), {}]}
email: jtshorwo@math.uwaterloo.ca Posted in old Singular Forum on: 2003-05-27 14:32:34+02
|
|
|
|
Posted: Thu Aug 11, 2005 5:32 pm |
|
|
|
|
|
It is currently Fri May 13, 2022 10:56 am
|
|