version 1.4, 2001/01/15 05:18:20 |
version 1.7, 2006/02/24 01:15:56 |
|
|
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* |
* |
* $OpenXM: OpenXM_contrib2/asir2000/lib/primdec,v 1.3 2000/08/22 05:04:22 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/lib/primdec,v 1.6 2003/12/12 03:08:29 takayama Exp $ |
*/ |
*/ |
/* Primary decomposition & Radical decomposition program */ |
/* Primary decomposition & Radical decomposition program */ |
/* written by T.Shimoyama, Fujitsu Lab. Date: 1995.10.12 */ |
/* written by T.Shimoyama, Fujitsu Lab. Date: 1995.10.12 */ |
Line 64 COMMONCHECK if 1 : redundancy check by intersection (i |
|
Line 64 COMMONCHECK if 1 : redundancy check by intersection (i |
|
SELECTFLAG selection strategy of separators (0 -- 3) |
SELECTFLAG selection strategy of separators (0 -- 3) |
*/ |
*/ |
|
|
if (vtype(minipoly) != 3) load("gr")$$ |
if (!module_definedp("gr")) load("gr")$ else {}$ |
|
module primdec $ |
|
/* Empty for now. It will be used in a future. */ |
|
endmodule$ |
|
|
|
|
#define GR(R,F,V,O) T2=newvect(4,time());R=dp_gr_main(F,V,0,0,O);GRTIME+=newvect(4,time())-T2; |
#define GR(R,F,V,O) T2=newvect(4,time());R=dp_gr_main(F,V,0,0,O);GRTIME+=newvect(4,time())-T2; |
#define HGRM(R,F,V,O) T2=newvect(4,time());R=dp_gr_main(F,V,1,1,O);GRTIME+=newvect(4,time())-T2; |
#define HGRM(R,F,V,O) T2=newvect(4,time());R=dp_gr_main(F,V,1,1,O);GRTIME+=newvect(4,time())-T2; |
#define NF(R,IN,F,G,O) T2=newvect(4,time());R=dp_nf(IN,F,G,O);NFTIME+=newvect(4,time())-T2; |
#define NF(R,IN,F,G,O) T2=newvect(4,time());R=dp_nf(IN,F,G,O);NFTIME+=newvect(4,time())-T2; |
Line 1384 def zraddec(F,X) |
|
Line 1388 def zraddec(F,X) |
|
/* contraction from V to X */ |
/* contraction from V to X */ |
def radcont(Q,V,X) |
def radcont(Q,V,X) |
{ |
{ |
dp_ord(irem(PRIMEORD,3)); |
|
for (R=[],I=length(Q)-1;I>=0;I--) { |
for (R=[],I=length(Q)-1;I>=0;I--) { |
|
dp_ord(irem(PRIMEORD,3)); |
G=Q[I]; |
G=Q[I]; |
for (E=1,J=0;J<length(G);J++) |
for (E=1,J=0;J<length(G);J++) |
E = lcm(E,dp_hc(dp_ptod(G[J],V))); |
E = lcm(E,dp_hc(dp_ptod(G[J],V))); |