[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.56 and 1.57

version 1.56, 2021/12/05 22:41:04 version 1.57, 2021/12/07 08:48:34
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2018/engine/nd.c,v 1.55 2021/11/29 09:19:33 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2018/engine/nd.c,v 1.56 2021/12/05 22:41:04 noro Exp $ */
   
 #include "nd.h"  #include "nd.h"
   
Line 2541  int comp_hn(P a, P b)
Line 2541  int comp_hn(P a, P b)
   DCP dc;    DCP dc;
   
   subp(CO,a,b,&s);    subp(CO,a,b,&s);
   if ( !s ) return 99999999; /* XXX */    if ( !s ) return -1;
   else if ( OID(s) == 1 ) return 0;    else if ( OID(s) == 1 ) return 0;
   else {    else {
     for ( dc = DC(s); NEXT(dc); dc = NEXT(dc) );      for ( dc = DC(s); NEXT(dc); dc = NEXT(dc) );
Line 2704  get_eg(&eg2); add_eg(&eg_update,&eg1,&eg2);
Line 2704  get_eg(&eg2); add_eg(&eg_update,&eg1,&eg2);
       g = update_base(g,nh);        g = update_base(g,nh);
       FREENDP(l);        FREENDP(l);
       if ( nd_hpdata ) {        if ( nd_hpdata ) {
         int dg;          int dg,sugar0;
   
         update_hpdata(&current_hpdata,nh);          update_hpdata(&current_hpdata,nh);
         dg = comp_hn(final_hpdata.hn,current_hpdata.hn);          dg = comp_hn(final_hpdata.hn,current_hpdata.hn);
           if ( dg < 0 ) {
              fprintf(asir_out,"We found a gb\n");
              d = 0;
           }
         if ( dg > sugar ) {          if ( dg > sugar ) {
           printexpr(CO,(Obj)current_hpdata.hn); printf("\n");             // printexpr(CO,(Obj)current_hpdata.hn);
              fprintf(asir_out,"\n");
         }          }
         while ( d && dg > sugar ) {          sugar0 = sugar;
           while ( d && dg > sugar0 ) {
           if ( DP_Print ) fprintf(asir_out,"sugar=%d done.\n",sugar);            if ( DP_Print ) fprintf(asir_out,"sugar=%d done.\n",sugar);
           d = nd_remove_same_sugar(d,sugar);            d = nd_remove_same_sugar(d,sugar);
           sugar++;            sugar0++;
         }          }
       }        }
     } else {      } else {
Line 2726  get_eg(&eg2); add_eg(&eg_update,&eg1,&eg2);
Line 2732  get_eg(&eg2); add_eg(&eg_update,&eg1,&eg2);
         MKNODE(t,list,nd_alltracelist);          MKNODE(t,list,nd_alltracelist);
         nd_alltracelist = t; nd_tracelist = 0;          nd_alltracelist = t; nd_tracelist = 0;
       }        }
       if ( DP_Print ) { printf("."); fflush(stdout); }        if ( DP_Print ) { fprintf(asir_out,"."); fflush(asir_out); }
         FREENDP(l);          FREENDP(l);
     }      }
   }    }
   conv_ilist(nd_demand,0,g,indp);    conv_ilist(nd_demand,0,g,indp);
   if ( !checkonly && DP_Print ) {    if ( !checkonly && DP_Print ) {
     printf("\nnd_gb done. Nnd_add=%d,Npairs=%d, Nnfnz=%d,Nnfz=%d,",Nnd_add,Npairs,Nnfnz,Nnfz);      fprintf(asir_out,"\nnd_gb done. Nnd_add=%d,Npairs=%d, Nnfnz=%d,Nnfz=%d,",Nnd_add,Npairs,Nnfnz,Nnfz);
     printf("Nremoved=%d\n",NcriB+NcriMF+Ncri2);      fprintf(asir_out,"Nremoved=%d\n",NcriB+NcriMF+Ncri2);
     fflush(asir_out);      fflush(asir_out);
   }    }
   if ( DP_Print ) {    if ( DP_Print ) {
Line 3505  again:
Line 3511  again:
         d = update_pairs(d,g,nh,0);          d = update_pairs(d,g,nh,0);
         g = update_base(g,nh);          g = update_base(g,nh);
         if ( nd_hpdata ) {          if ( nd_hpdata ) {
           int dg;            int dg,sugar0;
   
           update_hpdata(&current_hpdata,nh);            update_hpdata(&current_hpdata,nh);
           dg = comp_hn(final_hpdata.hn,current_hpdata.hn);            dg = comp_hn(final_hpdata.hn,current_hpdata.hn);
             if ( dg < 0 ) {
                fprintf(asir_out,"We found a gb\n");
                d = 0;
             }
           if ( dg > sugar ) {            if ( dg > sugar ) {
             printexpr(CO,(Obj)current_hpdata.hn); printf("\n");              // printexpr(CO,(Obj)current_hpdata.hn);
               fprintf(asir_out,"\n");
           }            }
           while ( d && dg > sugar ) {            sugar0 = sugar;
             while ( d && dg > sugar0 ) {
             if ( DP_Print ) fprintf(asir_out,"sugar=%d done.\n",sugar);              if ( DP_Print ) fprintf(asir_out,"sugar=%d done.\n",sugar);
             d = nd_remove_same_sugar(d,sugar);              d = nd_remove_same_sugar(d,sugar);
             sugar++;              sugar0++;
           }            }
         }          }
       } else {        } else {
         if ( DP_Print ) { printf("*"); fflush(stdout); }          if ( DP_Print ) { fprintf(asir_out,"*"); fflush(asir_out); }
       }        }
     } else {      } else {
       if ( DP_Print ) { printf("."); fflush(stdout); }        if ( DP_Print ) { fprintf(asir_out,"."); fflush(asir_out); }
     }      }
     FREENDP(l);      FREENDP(l);
   }    }
Line 3534  again:
Line 3546  again:
     }      }
   }    }
   conv_ilist(nd_demand,1,g,indp);    conv_ilist(nd_demand,1,g,indp);
   if ( DP_Print ) { printf("\nnd_gb_trace done.\n"); fflush(stdout); }    if ( DP_Print ) { fprintf(asir_out,"\nnd_gb_trace done.\n"); fflush(asir_out); }
   return g;    return g;
 }  }
   
Legend:
Removed from v.1.56  
changed lines
  Added in v.1.57
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>