[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.61

version 1.56, 2021/12/05 22:41:04 version 1.61, 2021/12/18 04:43:59
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.60 2021/12/17 23:53:13 noro Exp $ */
   
 #include "nd.h"  #include "nd.h"
   
Line 2505  void setup_hpdata(HPDATA final,HPDATA current)
Line 2505  void setup_hpdata(HPDATA final,HPDATA current)
   
   final->n = n = nd_nvar;    final->n = n = nd_nvar;
   final->hn = (P)BDY(nd_hpdata);    final->hn = (P)BDY(nd_hpdata);
   #if 0
   if ( NEXT(nd_hpdata) != 0 && (weight=(LIST)BDY(NEXT(nd_hpdata))) != 0 ) {    if ( NEXT(nd_hpdata) != 0 && (weight=(LIST)BDY(NEXT(nd_hpdata))) != 0 ) {
     wlen = length(BDY(weight));      wlen = length(BDY(weight));
     if ( n != wlen )      if ( n != wlen )
Line 2514  void setup_hpdata(HPDATA final,HPDATA current)
Line 2515  void setup_hpdata(HPDATA final,HPDATA current)
       w[i] = ZTOS((Z)BDY(nd));        w[i] = ZTOS((Z)BDY(nd));
   } else    } else
     w = 0;      w = 0;
   #else
     w = current_dl_weight_vector;
   #endif
   MKVECT(x,n);    MKVECT(x,n);
   for ( i = 0; i < n; i++ ) {    for ( i = 0; i < n; i++ ) {
     NEWDL(dl,n); dl->d[i] = 1; dl->td = 1; BDY(x)[i] = dl;      NEWDL(dl,n); dl->d[i] = 1; dl->td = 1; BDY(x)[i] = dl;
Line 2541  int comp_hn(P a, P b)
Line 2545  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 2584  ND_pairs nd_remove_same_sugar( ND_pairs d, int sugar)
Line 2588  ND_pairs nd_remove_same_sugar( ND_pairs d, int sugar)
 {  {
   struct oND_pairs root;    struct oND_pairs root;
   ND_pairs prev,cur;    ND_pairs prev,cur;
     int i;
   
   root.next = d;    root.next = d;
   prev = &root; cur = d;    prev = &root; cur = d;
     i = 0;
   while ( cur ) {    while ( cur ) {
     if ( SG(cur) == sugar )      if ( SG(cur) == sugar ) {
       prev->next = cur->next;        prev->next = cur->next;
     else        i++;
       } else
       prev = cur;        prev = cur;
     cur = cur->next;      cur = cur->next;
   }    }
     if ( DP_Print ) fprintf(asir_out,"[%d]",i);
   return root.next;    return root.next;
 }  }
   
Line 2704  get_eg(&eg2); add_eg(&eg_update,&eg1,&eg2);
Line 2712  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 > sugar ) {          if ( dg < 0 ) {
           printexpr(CO,(Obj)current_hpdata.hn); printf("\n");             int d_len;
              for ( d_len = 0; d; d = d->next, d_len++);
              fprintf(asir_out,"[%d] We found a gb\n",d_len);
         }          }
         while ( d && dg > sugar ) {          sugar0 = sugar;
           if ( DP_Print ) fprintf(asir_out,"sugar=%d done.\n",sugar);          while ( d && dg > sugar0 ) {
           d = nd_remove_same_sugar(d,sugar);            d = nd_remove_same_sugar(d,sugar0);
           sugar++;            sugar0++;
         }          }
       }        }
     } else {      } else {
Line 2726  get_eg(&eg2); add_eg(&eg_update,&eg1,&eg2);
Line 2736  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. Nbase=%d,Nnd_add=%d,Npairs=%d, Nnfnz=%d,Nnfz=%d,",nd_psn,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 3187  get_eg(&eg2); add_eg(&eg_remove,&eg1,&eg2);
Line 3197  get_eg(&eg2); add_eg(&eg_remove,&eg1,&eg2);
  g = conv_ilist_s(nd_demand,0,indp);   g = conv_ilist_s(nd_demand,0,indp);
  if ( DP_Print ) {   if ( DP_Print ) {
    printf("\ndlen=%d,nd_sba done. nd_add=%d,Nsyz=%d,Nsamesig=%d,Nnominimal=%d\n",dlen,Nnd_add,Nsyz,Nsamesig,Nnominimal);     printf("\ndlen=%d,nd_sba done. nd_add=%d,Nsyz=%d,Nsamesig=%d,Nnominimal=%d\n",dlen,Nnd_add,Nsyz,Nsamesig,Nnominimal);
    printf("Nnfnz=%d,Nnfz=%d,Nnfsingular=%d\n",Nnfnz,Nnfz,Nnfs);     printf("Nbase=%d,Nnfnz=%d,Nnfz=%d,Nnfsingular=%d\n",Nnfnz,Nnfz,Nnfs,nd_psn);
    fflush(stdout);     fflush(stdout);
    if ( nd_sba_redundant_check )     if ( nd_sba_redundant_check )
    printf("Nredundant=%d\n",Nredundant);     printf("Nredundant=%d\n",Nredundant);
Line 3505  again:
Line 3515  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 > sugar ) {            if ( dg < 0 ) {
             printexpr(CO,(Obj)current_hpdata.hn); printf("\n");               int d_len;
                for ( d_len = 0; d; d = d->next, d_len++);
                fprintf(asir_out,"[%d] We found a gb\n",d_len);
           }            }
           while ( d && dg > sugar ) {            sugar0 = sugar;
             if ( DP_Print ) fprintf(asir_out,"sugar=%d done.\n",sugar);            while ( d && dg > sugar0 ) {
             d = nd_remove_same_sugar(d,sugar);              d = nd_remove_same_sugar(d,sugar0);
             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. Nbase=%d\n",nd_psn); fflush(asir_out); }
   return g;    return g;
 }  }
   
Legend:
Removed from v.1.56  
changed lines
  Added in v.1.61
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>