version 1.13, 2004/02/09 23:37:12 |
version 1.15, 2006/02/24 01:15:56 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/lib/primdec_mod,v 1.12 2003/10/20 00:58:47 takayama Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/lib/primdec_mod,v 1.14 2004/07/30 02:24:11 noro Exp $ */ |
|
|
extern Hom,GBTime$ |
extern Hom,GBTime$ |
extern DIVLIST,INTIDEAL,ORIGINAL,ORIGINALDIMENSION,STOP,Trials,REM$ |
extern DIVLIST,INTIDEAL,ORIGINAL,ORIGINALDIMENSION,STOP,Trials,REM$ |
Line 1974 def contraction(P,V,W) |
|
Line 1974 def contraction(P,V,W) |
|
/* This procedure is called by zeroprimedecomposition. */ |
/* This procedure is called by zeroprimedecomposition. */ |
/* So, P is supposed to be a GB w.r.t. DRL. */ |
/* So, P is supposed to be a GB w.r.t. DRL. */ |
|
|
|
Ord0 = dp_ord(); |
Ord=0; |
Ord=0; |
YSet=setminus(W,V); |
YSet=setminus(W,V); |
|
|
Ord1 = [[Ord,length(V)],[0,length(YSet)]]; |
Ord1 = [[Ord,length(V)],[0,length(YSet)]]; |
GP1 = dp_gr_f_main(P,W,Hom,Ord1); |
W1 = append(V,YSet); |
|
GP1 = dp_gr_f_main(P,W1,Hom,Ord1); |
|
|
Factor = extcont_factor(GP1,V,Ord); |
Factor = extcont_factor(GP1,V,Ord); |
for ( F = 1, T = Factor; T != []; T = cdr(T) ) |
for ( F = 1, T = Factor; T != []; T = cdr(T) ) |
Line 1990 def contraction(P,V,W) |
|
Line 1992 def contraction(P,V,W) |
|
for ( T = G; T != []; T = cdr(T) ) |
for ( T = G; T != []; T = cdr(T) ) |
if ( !member(Vt,vars(car(T))) ) |
if ( !member(Vt,vars(car(T))) ) |
R = cons(car(T),R); |
R = cons(car(T),R); |
|
dp_ord(Ord0); |
return [R,F]; |
return [R,F]; |
} |
} |
|
|