source: git/Tst/Manual/determinecenter.tst @ 4efdd1c

spielwiese
Last change on this file since 4efdd1c was 894057, checked in by Oleksandr Motsak <motsak@…>, 13 years ago
ADD: Tests from online manual (res+stat on mamawutz): short ones Tst/Manual/s.lst
  • Property mode set to 100644
File size: 1.3 KB
Line 
1LIB "tst.lib"; tst_init();
2LIB "resbinomial.lib";
3ring r = 0,(x(1..4)),dp;
4list flag=identifyvar();
5ideal J=x(1)^2-x(2)^2*x(3)^5, x(1)*x(3)^3+x(4)^6;
6list Lmb=1,list(0,0,0,0),list(0,0,0,0),list(0,0,0,0),iniD(4),iniD(4),list(0,0,0,0),-1;
7list L=data(J,2,4);
8list LL=determinecenter(L[1],L[2],2,4,0,0,Lmb,flag,0,-1); // Compute the first center
9LL[1];  // index of variables in the center
10LL[2];  // exponents of ideals J_4,J_3,J_2,J_1
11LL[3];  // list of orders of J_4,J_3,J_2,J_1
12LL[4];  // list of critical values
13LL[5];  // components of the resolution function t
14LL[6];  // list of D_4,D_3,D_2,D_1
15LL[7];  // list of H_4,H_3,H_2,H_1 (exceptional divisors)
16LL[8];  // list of all exceptional divisors acumulated
17LL[9];  // auxiliary invariant
18LL[10]; // intvec pointing out the last step where the function t has dropped
19ring r= 0,(x(1..4)),dp;
20list flag=identifyvar();
21ideal J=x(1)^3-x(2)^2*x(3)^5, x(1)*x(3)^3+x(4)^5;
22list Lmb=2,list(0,0,0,0),list(0,0,0,0),list(0,0,0,0),iniD(4),iniD(4),list(0,0,0,0),-1;
23list L2=data(J,2,4);
24list L3=determinecenter(L2[1],L2[2],2,4,0,0,Lmb,flag,0,-1); // Example with rational exponents in E-Coeff
25L3[1]; // index of variables in the center
26L3[2]; // exponents of ideals J_4,J_3,J_2,J_1
27L3[3]; // list of orders of J_4,J_3,J_2,J_1
28L3[4]; // list of critical values
29L3[5]; // components of the resolution function
30tst_status(1);$
Note: See TracBrowser for help on using the repository browser.