version 1.63, 2022/01/13 08:15:02 |
version 1.64, 2022/03/09 07:50:51 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2018/engine/nd.c,v 1.62 2021/12/19 08:23:28 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2018/engine/nd.c,v 1.63 2022/01/13 08:15:02 noro Exp $ */ |
|
|
#include "nd.h" |
#include "nd.h" |
|
|
Line 1303 INLINE int ndl_find_reducer_s(UINT *dg,SIG sig) |
|
Line 1303 INLINE int ndl_find_reducer_s(UINT *dg,SIG sig) |
|
quo->pos = nd_psh[i]->sig->pos; |
quo->pos = nd_psh[i]->sig->pos; |
_adddl(nd_nvar,DL(quo),nd_sba_hm[quo->pos],DL2(quo)); |
_adddl(nd_nvar,DL(quo),nd_sba_hm[quo->pos],DL2(quo)); |
ret = comp_sig(sig,quo); |
ret = comp_sig(sig,quo); |
|
// if ( ret >= 0 ) { singular = 0; break; } |
if ( ret > 0 ) { singular = 0; break; } |
if ( ret > 0 ) { singular = 0; break; } |
if ( ret == 0 ) { /* fprintf(asir_out,"s"); fflush(asir_out); */ singular = 1; } |
if ( ret == 0 ) { /* fprintf(asir_out,"s"); fflush(asir_out); */ singular = 1; } |
} |
} |
Line 2423 int do_diagonalize(int sugar,int m) |
|
Line 2424 int do_diagonalize(int sugar,int m) |
|
LIST l; |
LIST l; |
Z iq; |
Z iq; |
|
|
|
if ( diag_period == 0 ) return 1; |
for ( i = nd_psn-1; i >= 0 && SG(nd_psh[i]) == sugar; i-- ) { |
for ( i = nd_psn-1; i >= 0 && SG(nd_psh[i]) == sugar; i-- ) { |
if ( nd_gentrace ) { |
if ( nd_gentrace ) { |
/* Trace = [1,index,1,1] */ |
/* Trace = [1,index,1,1] */ |
Line 3344 int do_diagonalize_trace(int sugar,int m) |
|
Line 3346 int do_diagonalize_trace(int sugar,int m) |
|
Z iq; |
Z iq; |
P cont,cont1; |
P cont,cont1; |
|
|
|
if ( diag_period == 0 ) return 1; |
for ( i = nd_psn-1; i >= 0 && SG(nd_psh[i]) == sugar; i-- ) { |
for ( i = nd_psn-1; i >= 0 && SG(nd_psh[i]) == sugar; i-- ) { |
if ( nd_gentrace ) { |
if ( nd_gentrace ) { |
/* Trace = [1,index,1,1] */ |
/* Trace = [1,index,1,1] */ |
Line 5495 void nd_gr_trace(LIST f,LIST v,int trace,int homo,int |
|
Line 5498 void nd_gr_trace(LIST f,LIST v,int trace,int homo,int |
|
Obj obj; |
Obj obj; |
NumberField nf; |
NumberField nf; |
struct order_spec *ord1; |
struct order_spec *ord1; |
struct oEGT eg_check,eg0,eg1; |
struct oEGT eg_check,eg_gb,eg_intred,eg0,eg1; |
NODE tr,tl1,tl2,tl3,tl4; |
NODE tr,tl1,tl2,tl3,tl4; |
LIST l1,l2,l3,l4,l5; |
LIST l1,l2,l3,l4,l5; |
int *perm; |
int *perm; |
Line 5633 void nd_gr_trace(LIST f,LIST v,int trace,int homo,int |
|
Line 5636 void nd_gr_trace(LIST f,LIST v,int trace,int homo,int |
|
} |
} |
if ( MaxDeg > 0 ) nocheck = 1; |
if ( MaxDeg > 0 ) nocheck = 1; |
while ( 1 ) { |
while ( 1 ) { |
tl1 = tl2 = tl3 = tl4 = 0; |
get_eg(&eg0); |
|
tl1 = tl2 = tl3 = tl4 = 0; |
if ( Demand ) |
if ( Demand ) |
nd_demand = 1; |
nd_demand = 1; |
ret = ndv_setup(m,1,fd0,nd_gbblock?1:0,0,0); |
ret = ndv_setup(m,1,fd0,nd_gbblock?1:0,0,0); |
Line 5659 void nd_gr_trace(LIST f,LIST v,int trace,int homo,int |
|
Line 5663 void nd_gr_trace(LIST f,LIST v,int trace,int homo,int |
|
nd_init_ord(ord); |
nd_init_ord(ord); |
nd_setup_parameters(nvar,0); |
nd_setup_parameters(nvar,0); |
} |
} |
|
get_eg(&eg1); init_eg(&eg_gb); add_eg(&eg_gb,&eg0,&eg1); |
nd_demand = 0; |
nd_demand = 0; |
cand = ndv_reducebase(cand,perm); |
cand = ndv_reducebase(cand,perm); |
if ( nd_gentrace ) { tl1 = nd_alltracelist; nd_alltracelist = 0; } |
if ( nd_gentrace ) { tl1 = nd_alltracelist; nd_alltracelist = 0; } |
|
get_eg(&eg0); |
cand = ndv_reduceall(0,cand); |
cand = ndv_reduceall(0,cand); |
|
get_eg(&eg1); init_eg(&eg_intred); add_eg(&eg_intred,&eg0,&eg1); |
cbpe = nd_bpe; |
cbpe = nd_bpe; |
if ( nd_gentrace ) { tl2 = nd_alltracelist; nd_alltracelist = 0; } |
if ( nd_gentrace ) { tl2 = nd_alltracelist; nd_alltracelist = 0; } |
get_eg(&eg0); |
get_eg(&eg0); |
Line 5700 void nd_gr_trace(LIST f,LIST v,int trace,int homo,int |
|
Line 5707 void nd_gr_trace(LIST f,LIST v,int trace,int homo,int |
|
} |
} |
get_eg(&eg1); init_eg(&eg_check); add_eg(&eg_check,&eg0,&eg1); |
get_eg(&eg1); init_eg(&eg_check); add_eg(&eg_check,&eg0,&eg1); |
if ( DP_Print ) |
if ( DP_Print ) |
fprintf(asir_out,"check=%.3fsec\n",eg_check.exectime); |
fprintf(asir_out,"gb=%.3fsec,check=%.3fsec,intred=%.3fsec\n", |
|
eg_gb.exectime,eg_check.exectime,eg_intred.exectime); |
/* dp->p */ |
/* dp->p */ |
nd_bpe = cbpe; |
nd_bpe = cbpe; |
nd_setup_parameters(nd_nvar,0); |
nd_setup_parameters(nd_nvar,0); |