|
|
D.15.1.25 germWithSemiNormalizedNNB
Procedure from library arnold.lib (see arnold_lib).
- Usage:
- germWithSemiNormalizedNNB(N); N NormalForm as given by @ref(determineNormalizedNNBGerm)
- Return:
- a polynomial, returned as a Poly, that is stable equivalent to N.phi.sourcegerm.value with a normalized nondegenerate Newton boundary up to scalar multiplication
of each of its variables, if possible, or an ERROR message otherwise
Example:
| | LIB "arnold.lib";
ring R = 0,(x,y),ds;
poly g = (x^2+y^2)^2+5*x^(10)+y^(11);
poly phix = x+y^2+x^2+x*y+x^2*y+x*y^3;
poly phiy = y+y^2+2*x^2+x*y+y*x^2+y^2*x+x*y^4;
map phi = R,phix,phiy;
g = phi(g);
Poly F = makePoly(g);
NormalForm N = determineNormalForm(F);
N=determineGermWithSemiNormalizedNNB(N);
==> // coefficients: QQ[a]/(...) considered as a field
==> // number of vars : 2
==> // block 1 : ordering ds
==> // : names x y
==> // block 2 : ordering C
germWithSemiNormalizedNNB(N);
==> x^2*y^2+(-24576/27341*a^3+74880/27341*a^2-180138/27341*a+6394015/874912)*\
x^3*y^2+(24576/27341*a^3-74880/27341*a^2+180138/27341*a-2019455/874912)*x\
^2*y^3+(-26624/27341*a^3+81120/27341*a^2-390299/54682*a+14948265/3499648)\
*x^4*y^2+27/2*x^3*y^3+(26624/27341*a^3-81120/27341*a^2+390299/54682*a-215\
10105/3499648)*x^2*y^4+(-1024/1439*a^3+3120/1439*a^2-18511/5756*a+653965/\
368384)*x^6*y+(-5120/27341*a^3+15600/27341*a^2-150115/109364*a-862983/699\
9296)*x^5*y^2+(-86528/27341*a^3+263640/27341*a^2-5073887/218728*a+3380317\
41/13998592)*x^4*y^3+(64512/27341*a^3-196560/27341*a^2+1891449/109364*a-5\
2688771/6999296)*x^3*y^4+(-2560/27341*a^3+7800/27341*a^2-150115/218728*a-\
2175351/13998592)*x^2*y^5+(-1024/1439*a^3+3120/1439*a^2-18511/5756*a+6539\
65/368384)*x*y^6+(-2560/1439*a^3+574076/136705*a^2-608021/109364*a+203080\
399/139985920)*x^7*y+(-78848/27341*a^3+343848/27341*a^2-583339/27341*a+22\
6849033/13998592)*x^6*y^2+(-62336/27341*a^3+189930/27341*a^2-14621201/874\
912*a+590152643/55994368)*x^5*y^3+(-7040/27341*a^3+21450/27341*a^2-165126\
5/874912*a+408058939/55994368)*x^4*y^4+(-10240/27341*a^3+31200/27341*a^2-\
150115/54682*a-1464485/3499648)*x^3*y^5+(-68352/27341*a^3+5508/1439*a^2-1\
367961/437456*a-152450143/27997184)*x^2*y^6+(-2560/1439*a^3+907924/136705\
*a^2-287631/27341*a+1039453101/139985920)*x*y^7+(-139872/136705*a^3+19924\
9/54682*a^2-20013101/17498240*a+2178742623/223977472)*x^8*y+(-1284512/136\
705*a^3+1833153/54682*a^2-914550071/17498240*a+7924644043/223977472)*x^7*\
y^2+(-6048/27341*a^3+17693/2878*a^2-57708587/3499648*a+3729008129/2239774\
72)*x^6*y^3+(-17376/27341*a^3+105885/54682*a^2-16302489/3499648*a-6813108\
53/223977472)*x^5*y^4+(-154112/27341*a^3+469560/27341*a^2-9036923/218728*\
a+46510223/3499648)*x^4*y^5+(90016/27341*a^3-847847/54682*a^2+110241579/3\
499648*a-7432121633/223977472)*x^3*y^6+(-1214272/136705*a^3+606120/27341*\
a^2-250349693/8749120*a+29675405/2947072)*x^2*y^7+(-498432/136705*a^3+289\
342/27341*a^2-9968963/546820*a+1147813827/55994368)*x*y^8+(-721168/136705\
*a^3+391615/27341*a^2-954977673/34996480*a+7909138069/223977472)*x^9*y+(3\
31072/27341*a^3-6508557/273410*a^2+2291864/27341*a-43434291217/1119887360\
)*x^8*y^2+(-125024/136705*a^3+371901/27341*a^2-497791607/17498240*a+12570\
15587/55994368)*x^7*y^3+(-142904/27341*a^3+3339385/218728*a^2-563526729/1\
3998592*a+34094288135/895909888)*x^6*y^4+(-4584/1439*a^3+111735/11512*a^2\
-17203179/736768*a-378059247/47153152)*x^5*y^5+(-1840/27341*a^3-121475/10\
9364*a^2-24116685/6999296*a+4622048811/447954944)*x^4*y^6+(-453344/136705\
*a^3-143212/27341*a^2+272645353/17498240*a-2511559963/111988736)*x^3*y^7+\
(-632152/27341*a^3+59015733/1093640*a^2-1689748861/13998592*a+41316213788\
3/4479549440)*x^2*y^8+(-60528/136705*a^3+63146/27341*a^2+124799097/349964\
80*a+3451660149/223977472)*x*y^9+(7702/27341*a^3-12826181/4374560*a^2+361\
975373/55994368*a-178947898291/17918197760)*x^10*y+(-421364/136705*a^3+15\
215871/1093640*a^2+1731555959/139985920*a+154463363943/2239774720)*x^9*y^\
2+(3247828/136705*a^3-132412347/2187280*a^2+17482422267/139985920*a-87917\
3436357/8959098880)*x^8*y^3+(115974/27341*a^3-15483443/874912*a^2+1741375\
401/55994368*a+12521896607/3583639552)*x^7*y^4+(3676/27341*a^3-5057415/43\
7456*a^2+188868501/27997184*a+81896510659/1791819776)*x^6*y^5+(54926/2734\
1*a^3-265195/46048*a^2-958583543/55994368*a+87746088817/3583639552)*x^5*y\
^6+(9004/1439*a^3-13741961/437456*a^2+1618851427/27997184*a-46738331987/1\
791819776)*x^4*y^7+(-876158/136705*a^3-3343627/4374560*a^2-13328418249/27\
9971840*a-6095946383/943063040)*x^3*y^8+(-3625296/136705*a^3+151465779/21\
87280*a^2-9908389557/69992960*a+1666153135399/8959098880)*x^2*y^9+(12584/\
27341*a^3-39423/546820*a^2-58186327/13998592*a+45332313/34996480)*x*y^10+\
x^22+y^22
==>
|
|