Changeset 06f16f in git for Singular/LIB/primdec.lib
- Timestamp:
- May 12, 2000, 2:20:26 PM (24 years ago)
- Branches:
- (u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')
- Children:
- d18870bf915431bbbed064042d13fe3ecb292569
- Parents:
- 38c1655b00c7eaba87af5a330f66679d70e33495
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/primdec.lib
r38c165 r06f16f 1 // $Id: primdec.lib,v 1.6 1 2000-05-11 11:22:40 pfister Exp $1 // $Id: primdec.lib,v 1.62 2000-05-12 12:20:26 krueger Exp $ 2 2 /////////////////////////////////////////////////////////////////////////////// 3 3 // primdec.lib // … … 11 11 /////////////////////////////////////////////////////////////////////////////// 12 12 13 version="$Id: primdec.lib,v 1.6 1 2000-05-11 11:22:40 pfister Exp $";13 version="$Id: primdec.lib,v 1.62 2000-05-12 12:20:26 krueger Exp $"; 14 14 info=" 15 15 LIBRARY: primdec.lib PROCEDURES FOR PRIMARY DECOMPOSITION … … 85 85 if(ordstr(basering)[1,2]!="dp") 86 86 { 87 execute "ring @Phelp=("+charstr(@P)+"),("+varstr(@P)+"),(C,dp);";87 execute("ring @Phelp=("+charstr(@P)+"),("+varstr(@P)+"),(C,dp);"); 88 88 ideal inew=std(imap(@P,id)); 89 89 ideal @h=imap(@P,h); … … 1182 1182 def bsr= basering; 1183 1183 string @id=string(h); 1184 execute "ring @r=0,("+@pa+","+varstr(bsr)+"),(C,dp);";1185 execute "ideal @i="+@id+";";1184 execute("ring @r=0,("+@pa+","+varstr(bsr)+"),(C,dp);"); 1185 execute("ideal @i="+@id+";"); 1186 1186 poly @p=lcmP(@i); 1187 1187 string @ps=string(@p); 1188 1188 setring bsr; 1189 execute "poly @p="+@ps+";";1189 execute("poly @p="+@ps+";"); 1190 1190 return(@p); 1191 1191 } … … 1507 1507 list l= prepareQuotientring(3); 1508 1508 l; 1509 execute l[1];1510 execute l[2];1509 execute(l[1]); 1510 execute(l[2]); 1511 1511 basering; 1512 1512 phi; … … 1562 1562 i=qr[1]; 1563 1563 1564 execute 1565 +ordstr(basering)+");" ;1564 execute("ring gnir = ("+charstr(basering)+"),("+varstr(basering)+"),(" 1565 +ordstr(basering)+");"); 1566 1566 1567 1567 … … 1597 1597 setring @P; 1598 1598 // kill gnir; 1599 execute 1600 ("+varstr(basering)+"),(C,lp);" ;1599 execute("ring gnir1 = ("+charstr(basering)+"), 1600 ("+varstr(basering)+"),(C,lp);"); 1601 1601 ideal i=fetch(@P,i); 1602 1602 list @pr=facstd(i); … … 1684 1684 def P=basering; 1685 1685 1686 execute "ring ir = ("+charstr(basering)+"),("+varstr(basering)+"),(C,lp);";1686 execute("ring ir = ("+charstr(basering)+"),("+varstr(basering)+"),(C,lp);"); 1687 1687 list l=fetch(P,li); 1688 1688 list @erg; … … 2046 2046 //---------------------------------------------------------------- 2047 2047 2048 execute "ring gnir = ("+charstr(basering)+"),("+varstr(basering)+"),(C,lp);";2048 execute("ring gnir = ("+charstr(basering)+"),("+varstr(basering)+"),(C,lp);"); 2049 2049 option(redSB); 2050 2050 … … 2198 2198 if(ordstr(@P)[1]=="w") 2199 2199 { 2200 execute "ring @Phelp=("+charstr(gnir)+"),("+varstr(gnir)+"),("+ordstr(@P)+");";2200 execute("ring @Phelp=("+charstr(gnir)+"),("+varstr(gnir)+"),("+ordstr(@P)+");"); 2201 2201 } 2202 2202 else 2203 2203 { 2204 execute "ring @Phelp=("+charstr(gnir)+"),("+varstr(gnir)+"),(C,dp);";2204 execute("ring @Phelp=("+charstr(gnir)+"),("+varstr(gnir)+"),(C,dp);"); 2205 2205 } 2206 2206 … … 2239 2239 //change the ring 2240 2240 { 2241 execute 2242 +ordstr(basering)+");" ;2241 execute("ring gnir1 = ("+charstr(basering)+"),("+varstr(basering)+"),(" 2242 +ordstr(basering)+");"); 2243 2243 ideal @j=fetch(gnir,@j); 2244 2244 attrib(@j,"isSB",1); … … 2249 2249 { 2250 2250 @va=string(maxideal(1)); 2251 execute 2252 +indep[@m][2]+");" ;2253 execute "map phi=gnir,"+@va+";";2251 execute("ring gnir1 = ("+charstr(basering)+"),("+indep[@m][1]+"),(" 2252 +indep[@m][2]+");"); 2253 execute("map phi=gnir,"+@va+";"); 2254 2254 if(homo==1) 2255 2255 { … … 2296 2296 //--------------------------------------------------------------------- 2297 2297 2298 execute quotring;2298 execute(quotring); 2299 2299 2300 2300 // @j considered in the quotientring … … 2507 2507 //change the ring 2508 2508 { 2509 execute 2510 varstr(basering)+"),("+ordstr(basering)+");" ;2509 execute("ring gnir1 = ("+charstr(basering)+"),("+ 2510 varstr(basering)+"),("+ordstr(basering)+");"); 2511 2511 ideal @j=fetch(gnir,jkeep); 2512 2512 attrib(@j,"isSB",1); … … 2515 2515 { 2516 2516 @va=string(maxideal(1)); 2517 execute 2518 restindep[@m][1]+"),(" +restindep[@m][2]+");" ;2519 execute "map phi=gnir,"+@va+";";2517 execute("ring gnir1 = ("+charstr(basering)+"),("+ 2518 restindep[@m][1]+"),(" +restindep[@m][2]+");"); 2519 execute("map phi=gnir,"+@va+";"); 2520 2520 if(homo==1) 2521 2521 { … … 2560 2560 //------------------------------------------------------------------ 2561 2561 2562 execute quotring;2562 execute(quotring); 2563 2563 2564 2564 // @j considered in the quotientring … … 2845 2845 //change the ring 2846 2846 { 2847 execute 2848 +ordstr(basering)+");" ;2847 execute("ring gnir1 = ("+charstr(basering)+"),("+varstr(basering)+"),(" 2848 +ordstr(basering)+");"); 2849 2849 ideal @j=fetch(@P,@j); 2850 2850 attrib(@j,"isSB",1); … … 2853 2853 { 2854 2854 @va=string(maxideal(1)); 2855 execute 2856 +indep[@m][2]+");" ;2857 execute "map phi=@P,"+@va+";";2855 execute("ring gnir1 = ("+charstr(basering)+"),("+indep[@m][1]+"),(" 2856 +indep[@m][2]+");"); 2857 execute("map phi=@P,"+@va+";"); 2858 2858 if(homo==1) 2859 2859 { … … 2897 2897 //------------------------------------------------------------------------ 2898 2898 2899 execute quotring;2899 execute(quotring); 2900 2900 2901 2901 // @j considered in the quotientring … … 3366 3366 def oldring=basering; 3367 3367 string n=system("neworder",PS); 3368 execute "ring r=("+charstr(oldring)+"),("+n+"),dp;";3368 execute("ring r=("+charstr(oldring)+"),("+n+"),dp;"); 3369 3369 ideal PS=imap(oldring,PS); 3370 3370 matrix m=char_series(PS); // We compute an irreducible … … 3963 3963 3964 3964 U[size(U)]=")"; // we compute the extractor of I (w.r.t. U) 3965 execute "ring RAU="+charstr(basering)+",("+A+U+",(dp("+string(a)+"),dp);";3965 execute("ring RAU="+charstr(basering)+",("+A+U+",(dp("+string(a)+"),dp);"); 3966 3966 ideal I=imap(R,SI); 3967 3967 //I=std(I,hv); // the standard basis in (R[U])[A] 3968 3968 I=std(I); // the standard basis in (R[U])[A] 3969 3969 A[size(A)]=")"; 3970 execute "ring Rloc=("+charstr(basering)+","+U+",("+A+",dp;";3970 execute("ring Rloc=("+charstr(basering)+","+U+",("+A+",dp;"); 3971 3971 ideal I=imap(RAU,I); 3972 3972 //"std in lokalisierung:"+newline,I; … … 4393 4393 if(dim(pr[1])>1) 4394 4394 { 4395 execute 4396 ("+varstr(basering)+"),(C,lp);" ;4395 execute("ring gnir = ("+charstr(basering)+"), 4396 ("+varstr(basering)+"),(C,lp);"); 4397 4397 ideal i=fetch(@P,i); 4398 4398 list @pr=facstd(i);
Note: See TracChangeset
for help on using the changeset viewer.