[BACK]Return to nd.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2018 / engine

Diff for /OpenXM_contrib2/asir2018/engine/nd.c between version 1.63 and 1.64

version 1.63, 2022/01/13 08:15:02 version 1.64, 2022/03/09 07:50:51
Line 1 
Line 1 
 /* $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);
Legend:
Removed from v.1.63  
changed lines
  Added in v.1.64
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>