|
|
D.5.5.15 calculateI
Procedure from library resbinomial.lib (see resbinomial_lib).
- Usage:
- calculateI(Coef,expJ,c,n,Y,a,b,D);
Coef, expJ, D lists, c, b numbers, n,Y integers, a intvec
- Return:
- ideal I, non monomial part of J
Example:
| | LIB "resbinomial.lib";
ring r = 0,(x(1..3)),dp;
list flag=identifyvar();
ideal J=x(1)^4*x(2)^2, x(3)^3;
list Lmb=1,list(0,0,0),list(0,0,0),list(3),iniD(3),iniD(3),list(0,0,0),-1;
list L=data(J,2,3);
list LL=determinecenter(L[1],L[2],3,3,0,0,Lmb,flag,0,-1); // Calculate the center
module auxpath=[0,-1];
list infochart=0,0,0,L[2],L[1],flag,0,list(0,0,0),auxpath,list(),list();
list L3=Blowupcenter(LL[1],1,1,infochart,3,3,0); // blowing-up and looking to the x(3) chart
calculateI(L3[2][1][5],L3[2][1][4],3,3,3,L3[2][1][3],3,iniD(3)); // (I_3)
==> [1]:
==> [1]:
==> [1]:
==> 4
==> [2]:
==> 2
==> [3]:
==> 3
==> [2]:
==> [1]:
==> [1]:
==> 0
==> [2]:
==> 0
==> [3]:
==> 0
// looking to the x(1) chart
calculateI(L3[2][2][5],L3[2][2][4],3,3,1,L3[2][2][3],3,iniD(3)); // (I_3)
==> [1]:
==> [1]:
==> [1]:
==> 3
==> [2]:
==> 2
==> [3]:
==> 0
==> [2]:
==> [1]:
==> [1]:
==> 0
==> [2]:
==> 0
==> [3]:
==> 3
|
|
