Changeset 119e27b in git
- Timestamp:
- Mar 25, 2011, 3:48:41 PM (12 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')
- Children:
- a3a116bf31680a590c2cc3c1e27c44b1b4fc7dea
- Parents:
- 787f10764ca8eabeccca61ebad635bf6a2b12c5a
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/resbinomial.lib
r787f10 r119e27b 169 169 ///////////////////////////////////////////////////////////////////////////////// 170 170 171 staticproc data(ideal K,int k,int n)171 proc data(ideal K,int k,int n) 172 172 "USAGE: data(K,k,n); K any ideal, k integer (!=0), n integer (!=0) 173 173 COMPUTE: Construcs a list with the coefficients and exponents of one ideal … … 1005 1005 list flag=identifyvar(); 1006 1006 ideal J=x(1)^2-x(2)^2*x(3)^5, x(1)*x(3)^3+x(4)^6; 1007 list Lmb=1, Resbinomial::list0(4),Resbinomial::list0(4),Resbinomial::list0(4),iniD(4),iniD(4),Resbinomial::list0(4),-1;1007 list 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; 1008 1008 list L=data(J,2,4); 1009 1009 list LL=determinecenter(L[1],L[2],2,4,0,0,Lmb,flag,0,-1); // Compute the first center … … 1022 1022 list flag=identifyvar(); 1023 1023 ideal J=x(1)^3-x(2)^2*x(3)^5, x(1)*x(3)^3+x(4)^5; 1024 list Lmb=2, Resbinomial::list0(4),Resbinomial::list0(4),Resbinomial::list0(4),iniD(4),iniD(4),Resbinomial::list0(4),-1;1024 list 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; 1025 1025 list L2=data(J,2,4); 1026 1026 list L3=determinecenter(L2[1],L2[2],2,4,0,0,Lmb,flag,0,-1); // Example with rational exponents in E-Coeff … … 1133 1133 list flag=identifyvar(); 1134 1134 ideal J=x(1)^3-x(3)^2*y(4)^2,x(1)*x(7)*y(2)-x(6)^3*x(5)*y(4)^3,x(5)^3-x(5)^3*y(2)^2; 1135 list Lmb=2,list 0(7),list0(7),list0(7),iniD(7),iniD(7),list0(7),-1;1135 list Lmb=2,list(0,0,0,0,0,0,0),list(0,0,0,0,0,0,0),list(0,0,0,0,0,0,0),iniD(7),iniD(7),list(0,0,0,0,0,0,0),-1; 1136 1136 list L=data(J,3,7); 1137 1137 list L2=determinecenter(L[1],L[2],2,7,0,0,Lmb,flag,0,-1); // Computing the center 1138 1138 module auxpath=[0,-1]; 1139 list infochart=0,0,0,L[2],L[1],flag,0,list 0(7),auxpath,list(),list();1139 list infochart=0,0,0,L[2],L[1],flag,0,list(0,0,0,0,0,0,0),auxpath,list(),list(); 1140 1140 list L3=Blowupcenter(L2[1],1,1,infochart,2,7,0); 1141 1141 L3[1]; // current chart (parent,Y,center,expJ,Coef,flag,Hhist,blwhist,path,hipercoef,hiperexp) with sons: [12],...,[16] … … 1369 1369 list flag=identifyvar(); 1370 1370 ideal J=x(1)^4*x(2)^2, x(3)^3; 1371 list Lmb=1,list 0(3),list0(3),list0(3),iniD(3),iniD(3),list0(3),-1;1371 list Lmb=1,list(0,0,0),list(0,0,0),list(3),iniD(3),iniD(3),list(0,0,0),-1; 1372 1372 list L=data(J,2,3); 1373 1373 list LL=determinecenter(L[1],L[2],3,3,0,0,Lmb,flag,0,-1); // Calculate the center 1374 1374 module auxpath=[0,-1]; 1375 list infochart=0,0,0,L[2],L[1],flag,0,list 0(3),auxpath,list(),list();1375 list infochart=0,0,0,L[2],L[1],flag,0,list(0,0,0),auxpath,list(),list(); 1376 1376 list L3=Blowupcenter(LL[1],1,1,infochart,3,3,0); // blowing-up and looking to the x(3) chart 1377 1377 calculateI(L3[2][1][5],L3[2][1][4],3,3,3,L3[2][1][3],3,iniD(3)); // (I_3)
Note: See TracChangeset
for help on using the changeset viewer.