Changeset a39a07 in git
- Timestamp:
- May 13, 1998, 12:54:07 PM (25 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')
- Children:
- edf0e9681a3d776622b3fefcdd261407d1a25780
- Parents:
- 546de146b5808486dd0067fe60ba45a307acd96b
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/primdec.lib
r546de1 ra39a07 1 // $Id: primdec.lib,v 1.1 5 1998-05-13 07:38:35obachman Exp $1 // $Id: primdec.lib,v 1.16 1998-05-13 10:54:07 obachman Exp $ 2 2 /////////////////////////////////////////////////////// 3 3 // primdec.lib … … 11 11 ////////////////////////////////////////////////////// 12 12 13 version="$Id: primdec.lib,v 1.1 5 1998-05-13 07:38:35obachman Exp $";13 version="$Id: primdec.lib,v 1.16 1998-05-13 10:54:07 obachman Exp $"; 14 14 info=" 15 15 LIBRARY: primdec.lib: PROCEDURE FOR PRIMARY DECOMPOSITION (I) … … 844 844 int nva = nvars(basering); 845 845 int @k,@s,@n,@k1,zz; 846 list primary,lres ,lres1,act,@lh,@wh;846 list primary,lres0,lres1,act,@lh,@wh; 847 847 map phi,psi,phi1,psi1; 848 848 ideal jmap,jmap1,jmap2,helpprim,@qh,@qht,ser1; … … 1082 1082 1083 1083 1084 kill lres ;1085 list lres ;1084 kill lres0; 1085 list lres0; 1086 1086 if(size(@lh)==2) 1087 1087 { 1088 1088 helpprim=@lh[2]; 1089 lres [1]=primary[2*@k-1];1089 lres0[1]=primary[2*@k-1]; 1090 1090 ser1=psi(helpprim); 1091 lres [2]=psi1(ser1);1092 if(size(reduce(lres [2],lres[1],1))==0)1091 lres0[2]=psi1(ser1); 1092 if(size(reduce(lres0[2],lres0[1],1))==0) 1093 1093 { 1094 1094 primary[2*@k]=primary[2*@k-1]; … … 1105 1105 @f=act[1][@n]^act[2][@n]; 1106 1106 ser1=psi(@f); 1107 lres [2*@n-1]=interred(primary[2*@k-1]+psi1(ser1));1107 lres0[2*@n-1]=interred(primary[2*@k-1]+psi1(ser1)); 1108 1108 helpprim=@lh[2*@n]; 1109 1109 ser1=psi(helpprim); 1110 lres [2*@n]=psi1(ser1);1110 lres0[2*@n]=psi1(ser1); 1111 1111 } 1112 1112 } … … 1114 1114 { 1115 1115 lres1=psi(@lh); 1116 lres =psi1(lres1);1116 lres0=psi1(lres1); 1117 1117 } 1118 1118 } … … 1130 1130 execute(@ri); 1131 1131 list @lvec; 1132 list @lr=imap(@P,lres );1132 list @lr=imap(@P,lres0); 1133 1133 ideal @lr1; 1134 1134 … … 1196 1196 kill @lvec; 1197 1197 setring @P; 1198 lres =imap(@Phelp1,@lr);1198 lres0=imap(@Phelp1,@lr); 1199 1199 kill @Phelp1; 1200 for(@n=1;@n<=size(lres );@n++)1200 for(@n=1;@n<=size(lres0);@n++) 1201 1201 { 1202 lres [@n]=clearSB(lres[@n]);1203 attrib(lres [@n],"isSB",1);1202 lres0[@n]=clearSB(lres0[@n]); 1203 attrib(lres0[@n],"isSB",1); 1204 1204 } 1205 1205 1206 primary[2*@k-1]=lres [1];1207 primary[2*@k]=lres [2];1206 primary[2*@k-1]=lres0[1]; 1207 primary[2*@k]=lres0[2]; 1208 1208 @s=size(primary)/2; 1209 for(@n=1;@n<=size(lres )/2-1;@n++)1209 for(@n=1;@n<=size(lres0)/2-1;@n++) 1210 1210 { 1211 primary[2*@s+2*@n-1]=lres [2*@n+1];1212 primary[2*@s+2*@n]=lres [2*@n+2];1211 primary[2*@s+2*@n-1]=lres0[2*@n+1]; 1212 primary[2*@s+2*@n]=lres0[2*@n+2]; 1213 1213 } 1214 1214 @k--;
Note: See TracChangeset
for help on using the changeset viewer.