version 1.3, 2000/08/22 05:04:22 |
version 1.6, 2003/12/12 03:08:29 |
|
|
* 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.2 2000/08/21 08:31:42 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/lib/primdec,v 1.5 2003/10/20 00:58:47 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 297 def primadec_main(F,REMS,H,VL) |
|
Line 301 def primadec_main(F,REMS,H,VL) |
|
def isolated(IP,F,GF,RA,H,SC,VL) |
def isolated(IP,F,GF,RA,H,SC,VL) |
{ |
{ |
T0 = newvect(4,time()); |
T0 = newvect(4,time()); |
if ( RA[0] == "begin" ) |
if ( RA[0] == "begin" ) { |
PD = primedec_main([RA[1]],VL); |
PD = primedec_main([RA[1]],VL); |
else if ( RA[0] == "ext" || RA[0] == "sep" ) { |
PD = map(dp_gr_main,PD,VL,0,1,PRIMEORD); /* XXX */ |
|
} else if ( RA[0] == "ext" || RA[0] == "sep" ) { |
if ( RA[0] == "sep" ) |
if ( RA[0] == "sep" ) |
T = prime_irred(idealsav(RA[1]),VL); |
T = prime_irred(idealsav(RA[1]),VL); |
else |
else |