| | LIB "modwalk.lib";
ring R1 = 0, (x,y,z,t), dp;
ideal I = 3x3+x2+1, 11y5+y3+2, 5z4+z2+4;
I = std(I);
ring R2 = 0, (x,y,z,t), lp;
ideal I = fetch(R1, I);
ideal J = modWalk(I);
J;
==> J[1]=x3+1/3x2+1/3
==> J[2]=z4+1/5z2+4/5
==> J[3]=y5+1/11y3+2/11
ring S1 = 0, (a,b,c,d), Dp;
ideal I = 5b2, ac2+9d3+3a2+5b, 2a2c+7abd+bcd+4a2, 2ad2+6b2d+7c3+8ad+4c;
I = std(I);
ring S2 = 0, (c,d,b,a), lp;
ideal I = fetch(S1,I);
// I is assumed to be a Dp-Groebner basis.
// We compute a lp-Groebner basis.
ideal J = modWalk(I,"Dp");
J;
==> J[1]=a25+16a24+96a23+256a22+256a21+256/9a20+1024/3a19+2048a18+65536/9a17+\
32768/3a16+16384/81a15+131072/81a14+1048576/81a13+1048576/27a12+1048576/9\
a11
==> J[2]=ba11+1522867351997104938459/91668001658017308797687087104a24+4293036\
9782248629690765/91668001658017308797687087104a23+80925218629630777478637\
/22917000414504327199421771776a22+7108535670237178684767/2864625051813040\
899927721472a21-3255817194541612658349/89519532869157528122741296a20+5380\
8965391546362724459/358078131476630112490965184a19+1534729815590907963215\
01/358078131476630112490965184a18-260815719913165309506063/44759766434578\
764061370648a17-1485276141860757031491027/89519532869157528122741296a16-4\
92332725360316960775/22379883217289382030685324a15+7423992361030571232440\
/16784912412967036523013993a14-17640364913371983121693/167849124129670365\
23013993a13-37723213977586186442564/5594970804322345507671331a12+92047580\
41857159721472414/5594970804322345507671331a11
==> J[3]=b2a6-63/2ba10+41087306587333357057895823883/924013456712814472680685\
83800832a24+93915562116924232413944264677/1320019223875449246686694054297\
6a23+1314662746341964624103002499857/30800448557093815756022861266944a22+\
2628268795042931967685617407557/23100336417820361817017145950208a21+32861\
7969148352577032114618159/2887542052227545227127143243776a20+175631724829\
284757906915538269/12993939235023953522072144596992a19+303812212296039043\
29090032467/206253003730538944794795945984a18+581332950565269518541458030\
9/6445406366579342024837373312a17+1392577308804410648719627124495/4060606\
01094498547564754518656a16+728489619081836101063651608791/135353533698166\
182521584839552a15+80182830319998353431517995037/913636352462621732020697\
666976a14+22597043807001043905240513127/32629869730807919000739202392a13+\
2699689128255025271541196980091/456818176231310866010348833488a12+6291773\
289016174025274735397/352483160672307766983293853a11+4/7a6
==> J[4]=b3a5+4/7ba5-5398327059462849163101023479/739210765370251578144548670\
406656a24+3500067651845908053488406611/92401345671281447268068583800832a2\
3+81936880336538263803253409291/46200672835640723634034291900416a22+85587\
506267677081700930548967/6600096119377246233433470271488a21+4349147060430\
84846670081811347/11550168208910180908508572975104a20+1021713731964491721\
016405426709/25987878470047907044144289193984a19+189087309367688338214840\
3875/812121202188997095129509037312a18+3873208423822126582196454031/10312\
6501865269472397397972992a17+104243440049859097996976327663/4060606010944\
98547564754518656a16+52530075262606982469983545909/5075757513681231844559\
4314832a15+6154556265260978917193662164647/365454540985048692808279066790\
4a14+138400848395446486358821423/5639730570756924271732701648a13+68577767\
53456717192310999393/38068181352609238834195736124a12+3439955547980759942\
57860510037/228409088115655433005174416744a11+9/2a10+18a9
==> J[5]=b4-63/2b3a4+4/7b2+3859043113737/128ba10-61261515/8ba9-525086793/32ba\
8+3969/16ba7-19845/4ba6+9/2ba5-317530772199391516703685862633925319890715\
5/9853890705461273608546303647569412096a24-889000025950798002413919214249\
17407047623059/17244308734557228814956031383246471168a23-4443358152931308\
1841158459544460024461570519/1437025727879769067913002615270539264a22-443\
94201884268758058296712245149507940575423/5388846479549134004673759807264\
52224a21-88794706503699382133468723144783169844707611/1077769295909826800\
934751961452904448a20-1763218978464215084468461907410740543335425/1732129\
22556936450150227993804931072a19-4436403904696107325808495191699219229446\
85/4210036312147760941151374849425408a18-21791665871565292166347154127869\
734967876143/33680290497182087529210998795403264a17-382244791239220889094\
021792022832016196677083/151561307237319393881449494579314688a16-25381669\
887869257097552146559952295821945667/6315054468221641411727062274138112a1\
5-3081988985408645474889243822837753592680251/487161344691383766047516232\
57636864a14-42991086686519981772739427292780808208480959/8525323532099215\
9058315340700864512a13-44791324614077928834488691100822888689474979/10656\
654415124019882289417587608064a12-250969820105414481651193596059252574282\
0133/197345452131926294116470696066816a11+771901337907/64a10-1701/4a9+437\
53689/4a8-19845/8a7-19845/2a6
==> J[6]=da-1323/800b3a4+63/20b2a4+2701701074331/1280ba10-61269453/320ba9+306\
291699/160ba8-27783/16ba6-189/200ba4-685883099436497069901143509465258298\
0621641/123173633818265920106828795594617651200a24-5486790468650205133145\
8584531975294116209709/61586816909132960053414397797308825600a23-60954530\
98260341081175657932145681344303369/1140496609428388149137303662913126400\
a22-54839253809480536502430419967744127237862429/384917605682081000333839\
9862331801600a21-6855741428155530311932580353074238311128479/481147007102\
601250417299982791475200a20-2878029541644786739384025698911163129301521/1\
732129225569364501502279938049310720a19-539637460347386518746735121876349\
5466944867/288688204261560750250379989674885120a18-1360851075145362976080\
2370233473935440353147/120286751775650312604324995697868800a17-2261317453\
3404843047583509526614944559001629/54129038299042640671946248064040960a16\
-19899168362053216404911620136404443756991/309379505595808417192193918976\
00a15-3443170203569948104100043669593479708510001/30447584043211485377969\
7645360230400a14-27203224432640761938740031635681832222207391/30447584043\
2114853779697645360230400a13-11007706754707314979975768772869745826795967\
3/152237920216057426889848822680115200a12-1380727908059280736173048219943\
2331439121223/6343246675669059453743700945004800a11-1929592445193/640a10+\
437570343/1600a9+164090367/100a8+3969/8a6+9/5a4
==> J[7]=db2+4/7d+81/160b3a4-189/200b3a3+9/5b2a3+250047/32ba10+4750893/320ba9\
-567/40ba8-567/10ba7+81/280ba4-27/50ba3-627430952592078879073202920720376\
153131/5173292620367168644486809414973941350400a24-2906859466655459274817\
4299538370243493/18476045072739888016024319339192647680a23-69048240138740\
571289193103124930690109/11975214398998075565941688460587827200a22+488085\
08639174871490992992380443024047/11547528170462430010015199586995404800a2\
1+5093069077149293873574182379116672811143/808326971932370100701063971089\
67833600a20+3286117544936595300844164893345508176887/36374713736956654531\
547878699035525120a19-57124126734867585149642894515011901057/189451634046\
6492423518118682241433600a18-296407013643924816106633846582502767991/2526\
021787288656564690824909655244800a17-193831832299525900077669479610197531\
4379/11367098042798954541108712093448601600a16+49045157132492520042165362\
0398209599/433131307834131784069071486566400a15+9565278519988924726047202\
0815175898220317/25575970596297647717494602210259353600a14-51800242230268\
556716933919846651112061/456713760648172280669546468040345600a13-15498981\
08498974500694541229046869490537/1598498162268602982343412638141209600a12\
-700479417253540849801700529315888259/26641636037810049705723543969020160\
a11+20611017/1600a10+21504771/320a9+243/50a8-1701/100a7+36/35a3
==> J[8]=d2
==> J[9]=ca-80796665/1042771968d2b2a+80796665/521385984d2ba-456425375/6069043\
2d2a-3048690575/30345216d2+91285075/182071296db4a+35/9216db4-1088785/5689\
728db3a2+91259683/91035648db3a+24027/2528768db3-8162999/79656192db2a3-653\
9063/79656192db2a2+25495277/45517824db2a+5095/1264384db2-5243/158048dba4-\
109809/632192dba3-173005/632192dba2+265684003/159312384dba+24027/4425344d\
b-7/256da5+5/96da3+5/24da2+35/32da+2645/2489256d+8617/1896576b5a-25921/15\
8048b5+7/256b4a2-2163/79024b4a+2163/19756b4+529/316096b3a3+7/128b3a2-1281\
1/59268b3a+15435/19756b3+311787/632192b2a4+529/158048b2a3+1/64b2a2-309/19\
756b2a+309/4939b2+19467/9878ba4+529/553168ba3+1/32ba2-1/8ba+1/2b+3/32a6+3\
/8a5+63/32a4+529/276584a3
==> J[10]=cb-1068739/27586560d2b2a+1068739/13793280d2ba-76214943/20230144d2a-\
508115055/10115072d2+25404981/101150720db4a-7203/1024db4+453789/790240db3\
a2+32184817/10115072db3a-25399101/12643840db3-324051/6321920db2a3+194523/\
790240db2a2+11921259/25287680db2a-9921275/1264384db2-27783/632192dba4-966\
217/3160960dba3-43211/395120dba2+1278655/2528768dba-31286843/3160960db-5/\
48da2-5/12da-172823/79024d+64827/3160960b5a-1351/148170b5+64827/1580480b4\
a-64827/395120b4+1323/1580480b3a3+9261/790240b3a-386/74085b3-583443/39512\
0b2a4-10887849/1580480b2a3+9261/395120b2a-9261/98780b2-583443/197560ba4-6\
222951/395120ba3-3/16a5-3/4a4-1555407/395120a3
==> J[11]=cd+3281/1008d2b2a-53819/8064d2ba+35/64d2b+5/96d2a2-155/768d2a+33775\
/8d2-7/192db5+7/256db4a+7/320db4-1/80db3a2+3/640db3a-23/96db3+441/128db2a\
3+1/64db2a+1/80db2+63/8dba3+1/32dba-1/8db+141/160da5+141/40da4+63/32da3-6\
3/80b4a2-63/40b3a2-9/20b2a2-9/10ba2
==> J[12]=c2+1/3cb2+5/3d+3a3
intvec w = 3,2,1,2;
ring S3 = 0, (c,d,b,a), (a(w),lp);
ideal I = fetch(S1,I);
// I is assumed to be a Dp-Groebner basis.
// We compute a (a(w),lp)-Groebner basis.
ideal J = modWalk(I,"Dp",w);
J;
==> J[1]=d2
==> J[2]=c2+3a3+1/3cb2+5/3d
==> J[3]=ca2+4ca+7/2b3+2b
==> J[4]=cda-6/7ba3+4/7c2-2/21cb3+1/7dba-10/21db
==> J[5]=db4a-192/49b2a3+128/49c2b-64/7cda-64/147cb4+4/63d2ba+2db3a-4db4+60/4\
9db2a-8db3+8/7dba-656/147db2-32/7db
==> J[6]=db6+333576da3-189/4b8+15876b4a2-5186640/2401b2a3-5085554404492727/51\
883209c2b+86846986315/726364926cd2+1815960/49cdb2-5085742263853031/148237\
74cda+8395442084/7203cb4-2304/49cb2a+4320/2401ca2+3806167/6174d2b3+369235\
87416091/242121642d2ba-43215/7db5+199982018/49db3a-4664160/7dba2-126b7+72\
/7b5a+31752b3a2+37367280/16807ba3-10171335148319374/51883209c2-1012272/49\
cdb+39190949752/16807cb3-4608/49cba-3806167/21609d2b2-129552/49db4+864216\
/2401db2a-117b6+144/7b4a+9072b2a2-20160cd+4663152/7cb2+17280/2401ca-13992\
486044/50421d2b-172796/49db3-4843728303706583/103766418dba-1296/7b5+288/4\
9b3a+18144ba2+9326304/7cb+185320d2-19529120/7203db2-1188/7b4+576/49b2a+20\
772144/16807db-20016/343b3-3312/49b2+8640/2401b
==> J[7]=cb3a-5/8da3-245/1024b8-7/32b6a+2/343b2a3-25430238177075283/512309629\
44c2b+338505036925/239077827072cd2+5/112cdb2-25430430033433795/1463741798\
4cda-108013/49392cb4+25/7cb2a-115/2744ca2+12978355/6096384d2b3+5544175817\
19895/717233481216d2ba-3/256db5-15431/2016db3a+5/4dba2-245/256b7-25/32b5a\
+1297475/115248ba3-25430430232582339/25615481472c2+5/84cdb-3240415/103723\
2cb3+16/7cba-12978355/21337344d2b2-47/2688db4+53/24696db2a-91/256b6-5/8b4\
a-5/84cd-5/4cb2-1291/686ca-52493006345/37340352d2b+1/192db3-2543087867776\
3075/102461925888dba+65/32b5-1/14b3a-5/2cb-25/72d2-667/98784db2+85/64b4-2\
/7b2a+6494431/1037232db-17/98b3+3/14ba+93/112b2-1291/1372b
==> J[8]=cb5-4/7db2a2+96/7b2a3-64/7c2b-4/441cdb2+32cda+116/21cb4-8/7dba2-8/44\
1cdb+32/7cb3+14db4-16/7db2a+16/7cb2+28db3+16/7cb+328/21db2+16db
==> J[9]=cdb3+18da3+2cdb2+63cb4+441/2db3a-36dba2+126cb3+36cb2+126dba+72cb+10d\
2
==> J[10]=ba4-2/3c2a+1/9cb3a+14/3cda-1/6dba2+49/12db3+5/9dba+7/3db
==> J[11]=b3a3+2744/5c2d+382/15c2b2-1127/15cdb3+7357/45cdba+129659/45cb5+2352\
/5cb3a+906269/90db4a-8232/5db2a2-1029/10b6a+236/5c2b+161/5cdb2+5474/45cda\
+28812/5cb4+1680cb2a+55174/405d2ba-98/15db5-2549/90db3a-735/2b5a-152/15c2\
-224/15cdb+8232/5cb3+5376/5cba+2156/45db4+259588/45db2a+2058/5b6-294b4a-2\
8cd+16464/5cb2-4032/5ca+1069/9db3-898/45dba+1470b5-168/5b3a+1232/45db2+11\
76b4-672/5b2a+3136/45db+672/5b3+504/5ba+2688/5b2-2016/5b
==> J[12]=b5a2+275/3c2d+32/9c2b2-235/18cdb3+28204/1323cdba+931/2cb5+64cb3a+17\
/1134d2b4+34/200120949d2b2a-7/18db6+175759/108db4a-117310/441db2a2-14b6a+\
4b4a2+16480960/3176523c2b-340/200120949cd2+194/63cdb2+4999816/453789cda+1\
799/2cb4+1600/7cb2a-17/567d2b3+8819321851/400241898d2ba-11/6db5-86053/756\
db3a+7930/441dba2-50b5a+32/7b3a2-12215296/3176523c2-124/21cdb+203cb3+1024\
/7cba+34/3969d2b2+1079/189db4+1232881/1323db2a+56b6-40b4a+16/7b2a2-400/63\
cd+514cb2-768/7ca-68/3969d2b+406/27db3-210473813/3176523dba+200b5-32/7b3a\
+16/7ba2-36cb+4484/1323db2+160b4-128/7b2a+12160/1323db+128/7b3+96/7ba+512\
/7b2-384/7b
==> J[13]=b7a+2931231/490c2d+20584/15435c2b2+1583507/8820cdb3+865156/5145cdba\
+758/5cb5+3904/245cb3a+45280/9261d2b2a+4/63db6+8178503/15435db4a-3032/35d\
b2a2-52/35b6a-1008b4a2+8917192/108045c2b+20176/15435cdb2+3119812/5145cda-\
8711357/140cb4+26496/343cb2a+147043/19845d2ba+863872/2205db5-3003452331/1\
3720db3a+1751209/49dba2-4b7-772/49b5a-2016b3a2+17258032/108045c2+2892256/\
2205cdb-8747741/70cb3+17152/245cba+53656/315db4+520568/1715db2a+208/35b6-\
736/49b4a-576b2a2+3946912/3087cd-8711357/245cb2-71168/1715ca+100/441d2b+7\
09600/3087db3-9004850681/72030dba+3088/49b5-432/1715b3a-1152ba2-3501522/4\
9cb+1502048/15435db2+2944/49b4-1984/245b2a+7552/2205db+1728/1715b3+8896/1\
715ba+7936/245b2-35584/1715b
==> J[14]=b9+5198192/5145c2d+101720954583616588/49029632505c2b2-135402014770/\
22880495169cd2b-865832/5145cdb3+50860863659839238/7004233215cdba+384/35cb\
5+3072/1715cb3a+2944/16807cba2-5191342/583443d2b4-222073876339030/6864148\
5507d2b2a-4/441db6+1380928/36015db4a-1536/245db2a2+4b8-96/245b6a-16607680\
/352947b2a3+203439862414515992/49029632505c2b+7936/36015cdb2-3065728/1200\
5cda-167877614888/15882615cb4+16384/2401cb2a+10382684/4084101d2b3+87904/1\
38915d2ba-664/15435db5-1334472908/36015db3a+2074976/343dba2+36/7b7-512/34\
3b5a-110595584/756315c2+256/108045cdb-5185904/245cb3+546304/84035cba+1679\
77620304/28588707d2b2+304/1715db4+50861936373471494/49029632505db2a+1504/\
245b6-512/343b4a-1024/21609cd-10353376/1715cb2-16384/12005ca-640/3087d2b+\
1056/2401db3-5347108784/252105dba+3728/343b5-6656/12005b3a-4149952/343cb-\
413536192/15882615db2+18944/2401b4-8704/12005b2a+4096/15435db+42304/12005\
b3+2048/12005ba+273152/84035b2-8192/12005b
==> J[15]=dba3-2/3c2d+1/9cdb3-1/6d2ba+5/9d2b
==> J[16]=db3a2-147/2cdb3-10cdba+63b4a2-32c2b-132cda-49/2db5+4db3a+126b3a2-64\
c2-82cdb-21/2db4+36b2a2-80cd-14db3-16dba+72ba2-6db2
==> J[17]=a5+1/9cb2a2-4/3c2a-7/6cb3+5/9da2-2/3cb
==> J[18]=da4-4/21c3+11/63cdba+8/63c2b-1/63db3a+10/63cdb+32/7cba+5/9d2a+2/63d\
b2a-b4a+4b4-4/7b2a+16/7b2
|
