Changeset 2d2c8be in git
- Timestamp:
- May 25, 2000, 1:22:19 PM (24 years ago)
- Branches:
- (u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')
- Children:
- e762aa7102176256bbb9d4ba9809f356f0a1f5ea
- Parents:
- 5674d5213a7afa320d0912b0dc1cd45291072f6f
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/primdec.lib
r5674d5 r2d2c8be 1 // $Id: primdec.lib,v 1.7 1 2000-05-25 10:07:19 pfister Exp $1 // $Id: primdec.lib,v 1.72 2000-05-25 11:22:19 pfister Exp $ 2 2 /////////////////////////////////////////////////////////////////////////////// 3 3 // primdec.lib // … … 11 11 /////////////////////////////////////////////////////////////////////////////// 12 12 13 version="$Id: primdec.lib,v 1.7 1 2000-05-25 10:07:19 pfister Exp $";13 version="$Id: primdec.lib,v 1.72 2000-05-25 11:22:19 pfister Exp $"; 14 14 info=" 15 15 LIBRARY: primdec.lib PROCEDURES FOR PRIMARY DECOMPOSITION … … 2000 2000 list primary,indep,ltras; 2001 2001 intvec @vh,isat; 2002 int @wr,@k,@n,@m,@n1,@n2,@n3,homo,seri,keepdi ;2002 int @wr,@k,@n,@m,@n1,@n2,@n3,homo,seri,keepdi,@w; 2003 2003 ideal peek=i; 2004 2004 ideal ser,tras; … … 2025 2025 2026 2026 homo=homog(i); 2027 if((find(ordstr(basering),"w")!=0)||(find(ordstr(basering),"W")!=0))2028 {2029 homo=0;2030 }2031 2032 2027 if(homo==1) 2033 2028 { … … 2055 2050 return(primary); 2056 2051 } 2057 intvec @hilb=hilb(tras,1); 2052 for(n=1;n<=nvars(basering);n++) 2053 { 2054 @w[n]=ord(var(n)); 2055 } 2056 intvec @hilb=hilb(tras,1,@w); 2058 2057 intvec keephilb=@hilb; 2059 2058 } … … 2082 2081 if((ordstr(@P)[1]!="(C,lp)")&&(ordstr(@P)[3]!="(C,lp)")) 2083 2082 { 2084 ideal @j=std(fetch(@P,i),@hilb );2083 ideal @j=std(fetch(@P,i),@hilb,@w); 2085 2084 } 2086 2085 else … … 2279 2278 else 2280 2279 { 2281 ideal jwork=std(imap(gnir,@j),@hilb );2280 ideal jwork=std(imap(gnir,@j),@hilb,@w); 2282 2281 } 2283 2282 … … 2318 2317 if(homo==1) 2319 2318 { 2320 ideal @j=std(phi(@j),@hilb );2319 ideal @j=std(phi(@j),@hilb,@w); 2321 2320 } 2322 2321 else … … 2500 2499 if(homo==1) 2501 2500 { 2502 @hilb=hilb(jwork,1 );2501 @hilb=hilb(jwork,1,@w); 2503 2502 } 2504 2503 … … 2584 2583 if(homo==1) 2585 2584 { 2586 ideal @j=std(phi(jkeep),keephilb );2585 ideal @j=std(phi(jkeep),keephilb,@w); 2587 2586 } 2588 2587 else
Note: See TracChangeset
for help on using the changeset viewer.