invbase.rlg 9.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412
  1. Tue Apr 15 00:32:48 2008 run on win32
  2. % ***** Example 1 *****
  3. g:=invbase{4*x^2 + x*y^2 - z +1/4,
  4. 2*x + y^2*z + 1/2,
  5. x^2*z - 1/2*x - y^2};
  6. 3 2 3 2
  7. g := {8*x*y*z - 2*x*y*z + 4*y - 4*y*z + 16*x*y + 17*y*z - 4*y,
  8. 4 2 2 2
  9. 8*y - 8*x*z - 256*y + 2*x*z + 64*z - 96*x + 20*z - 9,
  10. 3
  11. 2*y *z + 4*x*y + y,
  12. 3 2 2 2
  13. 8*x*z - 2*x*z + 4*y - 4*z + 16*x + 17*z - 4,
  14. 3 3 2
  15. - 4*y*z - 8*y + 6*x*y*z + y*z - 36*x*y - 8*y,
  16. 2 2 2
  17. 4*x*y + 32*y - 8*z + 12*x - 2*z + 1,
  18. 2
  19. 2*y *z + 4*x + 1,
  20. 3 2 2
  21. - 4*z - 8*y + 6*x*z + z - 36*x - 8,
  22. 2 2 2
  23. 8*x - 16*y + 4*z - 6*x - z}
  24. h:=invlex g;
  25. 6 5 4 3 2
  26. h := {3976*x + 37104*z - 600*z + 2111*z + 122062*z + 232833*z - 680336*z
  27. + 288814,
  28. 2 6 5 4 3 2
  29. 1988*y - 76752*z + 1272*z - 4197*z - 251555*z - 481837*z + 1407741*z
  30. - 595666,
  31. 7 6 5 4 3 2
  32. 16*z - 8*z + z + 52*z + 75*z - 342*z + 266*z - 60}
  33. % ***** Example 2 *****
  34. on trinvbase$
  35. invtorder revgradlex,{x,y,z}$
  36. g:=invbase{x^3 + y^2 + z - 3,
  37. y^3 + z^2 + x - 3,
  38. z^3 + x^2 + y - 3};
  39. ---------- ORDER = 3 ----------
  40. ---------- ORDER = 4 ----------
  41. ---------- ORDER = 5 ----------
  42. ---------- ORDER = 6 ----------
  43. ---------- ORDER = 7 ----------
  44. reductions = 77 zeros = 11 maxord = 7 order = 7 length = 13
  45. D i m e n s i o n = 0
  46. N u m b e r o f s o l u t i o n s = 27
  47. 2 2 3 2 2 2 2 2 2 2 2 2
  48. g := {x *y *z - 3*x *y - x*y *z - x *z + x*y*z + x *y + 3*x*y + 3*x
  49. 2
  50. - 3*x*y + y + z - 3,
  51. 2 3 2 2 2 2 2
  52. x *y*z + x *y - 3*x *y - x*y*z + x*z + x + 3*x*y - 3*x,
  53. 2 3 2 2 2 2 2 2
  54. x*y *z - 3*x*y - y *z - x*z + y*z - x + x*y + 3*y + 3*x - 3*y,
  55. 2 3 2 2 2 2
  56. x *y + x *z - 3*x - y - z + 3,
  57. 2 3 2 2 2
  58. x *z + x *y - x*y - 3*x - x*z + 3*x,
  59. 3 2 2
  60. x*y*z + x*y - 3*x*y - y*z + z + x + 3*y - 3,
  61. 2 3 2 2 2 2
  62. y *z + x *y - 3*y - z - x + 3,
  63. 3 2 2
  64. x*y + x*z + x - 3*x,
  65. 3 2
  66. x*z + x*y - y - 3*x - z + 3,
  67. 3 2 2
  68. y*z + x *y + y - 3*y,
  69. 3 2
  70. x + y + z - 3,
  71. 3 2
  72. y + z + x - 3,
  73. 3 2
  74. z + x + y - 3}
  75. h:=invlex g;
  76. h := { - 412373224241856640945111992285148*x
  77. 26
  78. - 1449641911307232269543863070491*z
  79. 25
  80. - 2168612583844782211565651535007*z
  81. 24
  82. - 2847785553349083352614138977565*z
  83. 23
  84. + 35576725674692081471990149502410*z
  85. 22
  86. + 54428253744724168431241789131696*z
  87. 21
  88. + 72399213723404842594731673129040*z
  89. 20
  90. - 367271934803243933721304377312611*z
  91. 19
  92. - 577412401939211224792461395441215*z
  93. 18
  94. - 752437808233499373488146484648759*z
  95. 17
  96. + 2023265153056028087298524971059780*z
  97. 16
  98. + 3362763223678472034221124579531852*z
  99. 15
  100. + 4206754352383617663824252489277347*z
  101. 14
  102. - 6294684651757967009725536832231313*z
  103. 13
  104. - 11645937803380007452970955449190202*z
  105. 12
  106. - 13912359441969785881761771576274650*z
  107. 11
  108. + 10813944944367254864931915957111635*z
  109. 10
  110. + 24146769890624467199683669920316403*z
  111. 9
  112. + 28253894162862384778437975597863994*z
  113. 8
  114. - 9413195341759783675090699662838024*z
  115. 7
  116. - 28732526014615244592092156992897700*z
  117. 6
  118. - 34274801170918929476253738727746640*z
  119. 5
  120. + 3129736563440111416048255862484824*z
  121. 4
  122. + 17956474721641990844572020234799903*z
  123. 3
  124. + 21526113174342847360723274047268152*z
  125. 2
  126. + 795762450545743140366490379212137*z
  127. - 6078501600786528783018721470971548*z
  128. - 3909915395631179340911139035268300,
  129. - 412373224241856640945111992285148*y
  130. 26
  131. + 3680069960199680647552580014011*z
  132. 25
  133. + 4946533576928304373640222248439*z
  134. 24
  135. + 6522058320833813074018729716109*z
  136. 23
  137. - 91123955793021263648983859056246*z
  138. 22
  139. - 122860148727246593163920662895892*z
  140. 21
  141. - 161652285275223157884596590612424*z
  142. 20
  143. + 962753147411097965886678769071203*z
  144. 19
  145. + 1303906344577106971108666976068347*z
  146. 18
  147. + 1646174502798616879170351863301227*z
  148. 17
  149. - 5539016636709239326199213127901604*z
  150. 16
  151. - 7732787650045519336370661934943044*z
  152. 15
  153. - 9110016144563661988538140239320223*z
  154. 14
  155. + 18612337918090152097453612706753413*z
  156. 13
  157. + 27965492180063505085033283788513066*z
  158. 12
  159. + 30440317356106139389125602029763822*z
  160. 11
  161. - 36863224004805998790098755360970471*z
  162. 10
  163. - 62542906673581589636380853858043447*z
  164. 9
  165. - 64689461678563738668073440578715518*z
  166. 8
  167. + 42623160090556250860454187465583768*z
  168. 7
  169. + 83548043234053149543179359124170180*z
  170. 6
  171. + 85865493477306743665317502795142584*z
  172. 5
  173. - 27434780477528021937653276615015928*z
  174. 4
  175. - 61602505785524913541319871156904287*z
  176. 3
  177. - 62515628463318116801915981996829328*z
  178. 2
  179. + 5925778048881538700551831705942583*z
  180. + 24088990130824351149845277501309728*z
  181. + 15820742036151533576971241715895080,
  182. 27 24 21 19 18 17 16
  183. - z + 27*z - 317*z + 18*z + 2067*z + 50*z - 279*z
  184. 15 14 13 12 11 10
  185. - 8156*z - 645*z + 1674*z + 20359*z + 3044*z - 4645*z
  186. 9 8 7 6 5 4 3
  187. - 33644*z - 6288*z + 6388*z + 36936*z + 5925*z - 4957*z - 23187*z
  188. 2
  189. - 4063*z + 4342*z + 5352}
  190. % ***** Example 3 (limited by the degree bound) *****
  191. invtorder revgradlex,{x,z,y,t}$
  192. k:=5$
  193. on errcont$
  194. invbase{x^(k+1)-y^(k-1)*z*t,
  195. x*z^(k-1)-y**k,
  196. x^k*y-z^k*t};
  197. ---------- ORDER = 6 ----------
  198. ---------- ORDER = 7 ----------
  199. ---------- ORDER = 8 ----------
  200. ---------- ORDER = 9 ----------
  201. ---------- ORDER = 10 ----------
  202. ---------- ORDER = 11 ----------
  203. ---------- ORDER = 12 ----------
  204. ---------- ORDER = 13 ----------
  205. ---------- ORDER = 14 ----------
  206. ---------- ORDER = 15 ----------
  207. ---------- ORDER = 16 ----------
  208. ---------- ORDER = 17 ----------
  209. ---------- ORDER = 18 ----------
  210. ---------- ORDER = 19 ----------
  211. ---------- ORDER = 20 ----------
  212. ---------- ORDER = 21 ----------
  213. ***** Maximum degree bound exceeded.
  214. invtempbasis;
  215. 17 2 16
  216. { - t*z + x *y ,
  217. 13 3 11
  218. - t*z + x *y ,
  219. 9 4 6
  220. - t*z + x *y ,
  221. 4 6
  222. - t*y *z + x ,
  223. 5 5
  224. - t*z + x *y,
  225. 4 5
  226. x*z - y }
  227. end$
  228. Time for test: 28 ms