Opened 8 years ago
Closed 8 years ago
#656 closed bug (fixed)
primdecSY() failing example using given ordering
Reported by: | Owned by: | somebody | |
---|---|---|---|
Priority: | minor | Milestone: | 4-1-0 and higher |
Component: | dontKnow | Version: | 4-0-1 |
Keywords: | primdecSY bug given ordering | Cc: |
Description
LIB("primdec.lib"); ring rng = (0),(xy,xg,xt,xb),dp; ideal I = -11658/12181*xg*xt*xb-8219/674*xg*xb^2,307/3906*xy^2+20884/14685*xg^3+21229/26747*xg,-3493/11608*xy+13440/12053*xg*xt-1161/3359*xb-1377/857; def L1 = primdecSY(I,0); L1; def pc1 = primdecSY(L1[1][1],0); pc1; ASSUME(0, size(pc1)==1);
output
> L1; [1]: [1]: _[1]=xg*xb _[2]=10055169859*xy*xb+11549693016*xb^2+53690981544*xb _[3]=449104590197760*xg*xt-121194962310527*xy-139208449921848*xb-647137400549832 _[4]=1357996201*xb^3+12625815918*xb^2+29346773481*xb _[5]=103896688728*xg^3+5742065065*xy^2+57985102890*xg _[6]=2098625879096419232773276*xy*xg^2+429797962984283189542400*xy^2*xt+11205905511528825226115616*xg^2+1171250393220930037520505*xy+1345336048645961456490120*xb+6254054792474858515369080 _[7]=2052314018137716618924676595712000*xy^2*xt^2+2704279380778581227514567627811201*xy^2*xg+28879753670842996554294113888414832*xy*xg+5592798960856239544961320746374400*xy*xt+6424069608359399396642533059225600*xt*xb+77103735843339848393418416082209856*xg+29863529905251970129424864884070400*xt [2]: _[1]=xb _[2]=133701872640*xg*xt-36080667553*xy-192657755448 _[3]=1535082492565853405604721165825046536141096374288038267682816*xg^3+314384443120329561947410092241853885818408775425768449638400*xy*xg*xt-1693390411685174109220974568908148419775496470839017211878400*xg*xt+3964456979190256096145244292549837861047821651456631049925*xy+856734871590571619633443179410872467138359797800021466030080*xg+2440091446518952175400104854424672347144136799890452464665880 _[4]=6657713477159283485799806016*xy*xg^2+1363497762568380888521318400*xy^2*xt+1289320312116364717296307200*xy*xg*xt-347934824189465841842579065*xy^2+35549789455565591473748307456*xg^2+6884505587957174818395955200*xg*xt+9920230493078625888567896640 _[5]=610989585631948978542624768000*xy^2*xt^2+805084662333605605095137727839*xy^2*xg+8597723629307233270108399490448*xy*xg+1665019041636272564739899001600*xy*xt+22954372087924932537486875880384*xg+8890601341247981580656405145600*xt [2]: [1]: _[1]=7857492*xt+100115639*xb _[2]=3746866085456823230720*xg*xb+79357370566272284857*xy+91152440132776773768*xb+423739745643425045112 _[3]=103896688728*xg^3+5742065065*xy^2+57985102890*xg _[4]=20886268138853856146840277901684997569863143674774863590400*xy*xg^2-54501462218930437302891222070322252537074908154485944320000*xy^2*xb-508113784347294309608339828466045022416077454072809849760*xy*xg+111525141085758123710004679965217905427847720423890010726400*xg^2+11669060015336161926085582048397509637398665483504074257869*xy-2713144402351818576426257823065438195618917080244794428160*xg+13403459399721804569264388316552623544383235579366981001256*xb+62308573072832287290877129470579180667139228348369871871704 _[5]=271091799117902125196825419532000853142481631257841412654120960000*xy^2*xb^2+2200331786955787997491814271413008662504358903187930401698996720*xy^2*xg-107057795452692278426615122248670556896011888304132064073057711*xy^2+23497956776608591613710245993669723625507907730833541150461303040*xy*xg-58165201293212597306087438898136237517802037804925695590875166996*xy*xb-66669182351079448138049006889405612983221456502367585775715379168*xb^2-1143300053682728791407222688677341944610102726687661440346667152*xy+62735308369026012653634178404475477892010309480055824566930488320*xg-310581202797055851747341904423820676162076317166454322973266519136*xb-3052405028572953954462014780673096142026997446395465543083540416 [2]: _[1]=7857492*xt+100115639*xb _[2]=3746866085456823230720*xg*xb+79357370566272284857*xy+91152440132776773768*xb+423739745643425045112 _[3]=6266262831745919295412345101681645443828466604657514464419200474697044099032764210327265765315000751218143911630318940441818491557679252194742763520000000*xg^3-16351436489650916382705947286442895808309998954377523365215307354974763417650335488003733605603111661726566670037918262265983692996932866897084416000000000*xy*xg*xb-397791995153754243583308662289819853953812606737385683962294584416376446639455737093753721288677280271205225679579402466423069278700690032977510400000000*xy*xb+86950131941917509670971307576280861690104906074178941997831808312413960622478947347965805809829017426051525900910023587081936667261501316912951940363759360*xg*xb-7638302804406911717626165505214370163380384598773347326357358447483445743871264221308094021201926282946317502758192232294431829459162573623327282309609*xy+3497223053814684694328568762711662799403219028660966329993378570274206859237510907874447018383949663884338760085523554540771112800008193617533337600000000*xg+2115292224383397640250237201450304414500502461359728912824843356851512502310252015935210316917108233138749920289549203322429189346841002760913611864852984*xb+9833344974814671687433933536220189382748543490083682948281215994650776241031911280503783999911885547495756155732066551028792217248323935360488180373832456 _[4]=74189026970079562018671715440124392248037239763696704666553139200*xy*xg^2-193591809871827421960860159572444016479811853360325489549967360000*xy^2*xb-183642336586383923066068094546173740005043203830181350233678990080*xy*xg*xb-3889483270487625438206666361727350441260481622877047528503847123*xy^2-1804844551463238057855684271023158289782983088584416081220308480*xy*xg+396142654343384971807874703461092330539175639636237664820687667200*xg^2-4467586166611744598085906292265301958795155032352137380900801752*xy*xb-980583863089226723229895261078737917164033571177576363364092897280*xg*xb-87815246755029481447771197879950646534506252898308023700339424*xy-9637219148084972470757736247903943611871783177049361558956871680*xg+23754431861918883529098289287448006046911419544477076852864417856*xb+110427070305648748161143492430149262377963918774112603690194311104 _[5]=271091799117902125196825419532000853142481631257841412654120960000*xy^2*xb^2+2200331786955787997491814271413008662504358903187930401698996720*xy^2*xg-107057795452692278426615122248670556896011888304132064073057711*xy^2+23497956776608591613710245993669723625507907730833541150461303040*xy*xg-58165201293212597306087438898136237517802037804925695590875166996*xy*xb-66669182351079448138049006889405612983221456502367585775715379168*xb^2-1143300053682728791407222688677341944610102726687661440346667152*xy+62735308369026012653634178404475477892010309480055824566930488320*xg-310581202797055851747341904423820676162076317166454322973266519136*xb-3052405028572953954462014780673096142026997446395465543083540416 > def pc1 = primdecSY(L1[1][1],0); > pc1; [1]: [1]: _[1]=xg _[2]=10055169859*xy+11549693016*xb+53690981544 _[3]=1357996201*xb^2+12625815918*xb+29346773481 [2]: _[1]=36851*xb+171309 _[2]=xg _[3]=xy [2]: [1]: _[1]=xb _[2]=133701872640*xg*xt-36080667553*xy-192657755448 _[3]=103896688728*xg^3+5742065065*xy^2+57985102890*xg _[4]=624776980975415073764*xy*xg^2+127954142001870553600*xy^2*xt+3336083808136000365024*xg^2+348690203400098254695*xy+1861879962034789674120 _[5]=610989585631948978542624768000*xy^2*xt^2+805084662333605605095137727839*xy^2*xg+8597723629307233270108399490448*xy*xg+1665019041636272564739899001600*xy*xt+22954372087924932537486875880384*xg+8890601341247981580656405145600*xt [2]: _[1]=xb _[2]=133701872640*xg*xt-36080667553*xy-192657755448 _[3]=59789615854769376822075615326106048*xg^3+12244895747238632274508365455155200*xy*xg*xt+19297908218203249487956064497536000*xg^2*xt-5207716220729950935108011368282825*xy*xg-21794442752409714892396673910374400*xg*xt-11762857586443827578091842031557760*xy+5561465383766446940633621772038040*xg+31404709141396113092158123693678080 _[4]=624776980975415073764*xy*xg^2+127954142001870553600*xy^2*xt+3336083808136000365024*xg^2+348690203400098254695*xy+1861879962034789674120 _[5]=610989585631948978542624768000*xy^2*xt^2+805084662333605605095137727839*xy^2*xg+8597723629307233270108399490448*xy*xg+1665019041636272564739899001600*xy*xt+22954372087924932537486875880384*xg+8890601341247981580656405145600*xt > . ASSUME(0, size(pc1)==1); ? ASSUME failed:ASSUME(0, size(pc1)==1);
Note: See
TracTickets for help on using
tickets.
fixed by mlee with 8eec294f01ff794d6970f9b1dc8693a5ca131146 (see tr #659)