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

Diff for /OpenXM_contrib2/asir2000/engine/nd.c between version 1.216 and 1.218

version 1.216, 2013/12/20 04:35:34 version 1.218, 2014/02/24 01:45:28
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.215 2013/12/20 02:02:24 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.217 2014/02/03 02:43:05 noro Exp $ */
   
 #include "nd.h"  #include "nd.h"
   
Line 2703  int ndv_newps(int m,NDV a,NDV aq,int f4)
Line 2703  int ndv_newps(int m,NDV a,NDV aq,int f4)
     nd_ps[nd_psn] = a;      nd_ps[nd_psn] = a;
     if ( aq ) {      if ( aq ) {
         nd_ps_trace[nd_psn] = aq;          nd_ps_trace[nd_psn] = aq;
                 nd_ps_gz[nd_psn] = ndvtondvgz(aq);                  if ( !nd_vc ) nd_ps_gz[nd_psn] = ndvtondvgz(aq);
         register_hcf(aq);          register_hcf(aq);
         nd_bound[nd_psn] = ndv_compute_bound(aq);          nd_bound[nd_psn] = ndv_compute_bound(aq);
         SG(r) = SG(aq); ndl_copy(HDL(aq),DL(r));          SG(r) = SG(aq); ndl_copy(HDL(aq),DL(r));
Line 2711  int ndv_newps(int m,NDV a,NDV aq,int f4)
Line 2711  int ndv_newps(int m,NDV a,NDV aq,int f4)
         if ( !m ) register_hcf(a);          if ( !m ) register_hcf(a);
         nd_bound[nd_psn] = ndv_compute_bound(a);          nd_bound[nd_psn] = ndv_compute_bound(a);
         SG(r) = SG(a); ndl_copy(HDL(a),DL(r));          SG(r) = SG(a); ndl_copy(HDL(a),DL(r));
                 if ( !m ) nd_ps_gz[nd_psn] = ndvtondvgz(a);                  if ( !m && !nd_vc ) nd_ps_gz[nd_psn] = ndvtondvgz(a);
     }      }
     if ( nd_demand ) {      if ( nd_demand ) {
         if ( aq ) {          if ( aq ) {
Line 2789  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
Line 2789  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
         hc = HCU(w[i].p);          hc = HCU(w[i].p);
         if ( trace ) {          if ( trace ) {
             a = nd_ps_trace[i] = ndv_dup(0,w[i].p);              a = nd_ps_trace[i] = ndv_dup(0,w[i].p);
                         nd_ps_gz[i] = ndvtondvgz(a);                          if ( !nd_vc ) nd_ps_gz[i] = ndvtondvgz(a);
             if ( !dont_removecont) ndv_removecont(0,a);              if ( !dont_removecont) ndv_removecont(0,a);
             register_hcf(a);              register_hcf(a);
             am = nd_ps[i] = ndv_dup(mod,a);              am = nd_ps[i] = ndv_dup(mod,a);
Line 2799  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
Line 2799  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
             ndv_removecont(mod,am);              ndv_removecont(mod,am);
         } else {          } else {
             a = nd_ps[i] = ndv_dup(mod,w[i].p);              a = nd_ps[i] = ndv_dup(mod,w[i].p);
                         if ( !mod ) nd_ps_gz[i] = ndvtondvgz(a);                          if ( !mod && !nd_vc ) nd_ps_gz[i] = ndvtondvgz(a);
             if ( mod || !dont_removecont ) ndv_removecont(mod,a);              if ( mod || !dont_removecont ) ndv_removecont(mod,a);
             if ( !mod ) register_hcf(a);              if ( !mod ) register_hcf(a);
         }          }
Line 6068  NODE nd_f4(int m,int **indp)
Line 6068  NODE nd_f4(int m,int **indp)
                 node = BDY((LIST)BDY(tn));                  node = BDY((LIST)BDY(tn));
                             if ( QTOS((Q)ARG0(node)) == sugar ) break;                              if ( QTOS((Q)ARG0(node)) == sugar ) break;
             }              }
             if ( !tn ) error("nd_f4 : inconsistent non-zero list");                          if ( tn ) {
                         for ( t = l, ll0 = 0; t; t = NEXT(t) ) {                                  for ( t = l, ll0 = 0; t; t = NEXT(t) ) {
                 for ( tn = BDY((LIST)ARG1(node)); tn; tn = NEXT(tn) ) {                          for ( tn = BDY((LIST)ARG1(node)); tn; tn = NEXT(tn) ) {
                                   i1s = QTOS((Q)ARG0(BDY((LIST)BDY(tn))));                                          i1s = QTOS((Q)ARG0(BDY((LIST)BDY(tn))));
                                   i2s = QTOS((Q)ARG1(BDY((LIST)BDY(tn))));                                          i2s = QTOS((Q)ARG1(BDY((LIST)BDY(tn))));
                                   if ( t->i1 == i1s && t->i2 == i2s ) break;                                          if ( t->i1 == i1s && t->i2 == i2s ) break;
                                 }                                          }
                             if ( tn ) {                                  if ( tn ) {
                                     if ( !ll0 ) ll0 = t;                                          if ( !ll0 ) ll0 = t;
                                         else NEXT(ll) = t;                                                  else NEXT(ll) = t;
                                         ll = t;                                                  ll = t;
                                 }                                          }
             }                  }
                         if ( ll0 ) NEXT(ll) = 0;                                  if ( ll0 ) NEXT(ll) = 0;
                     l = ll0;                          l = ll0;
                           } else l = 0;
         }          }
         bucket = create_pbucket();          bucket = create_pbucket();
         stat = nd_sp_f4(m,0,l,bucket);          stat = nd_sp_f4(m,0,l,bucket);
Line 7947  MAT nd_btog(LIST f,LIST v,int mod,struct order_spec *o
Line 7948  MAT nd_btog(LIST f,LIST v,int mod,struct order_spec *o
         pi0 = QTOS((Q)ARG0(pi)); pi1 = QTOS((Q)ARG1(pi));          pi0 = QTOS((Q)ARG0(pi)); pi1 = QTOS((Q)ARG1(pi));
         p[pi0] = c = (ND *)MALLOC(nb*sizeof(ND));          p[pi0] = c = (ND *)MALLOC(nb*sizeof(ND));
         ptomp(mod,(P)ARG2(pi),&inv);          ptomp(mod,(P)ARG2(pi),&inv);
           ((MQ)inv)->cont = invm(((MQ)inv)->cont,mod);
         u = ptond(CO,vv,(P)ONE);          u = ptond(CO,vv,(P)ONE);
         HCM(u) = ((MQ)inv)->cont;          HCM(u) = ((MQ)inv)->cont;
         c[pi1] = u;          c[pi1] = u;
Line 8018  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
Line 8020  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
     pi = BDY((LIST)BDY(t));      pi = BDY((LIST)BDY(t));
         pi0 = QTOS((Q)ARG0(pi)); pi1 = QTOS((Q)ARG1(pi));          pi0 = QTOS((Q)ARG0(pi)); pi1 = QTOS((Q)ARG1(pi));
         if ( pi1 == pos ) {          if ( pi1 == pos ) {
           ptomp(mod,(P)ARG2(pi),&inv);            ptomp(mod,(P)ARG2(pi),&inv);
             ((MQ)inv)->cont = invm(((MQ)inv)->cont,mod);
           u = ptond(CO,vv,(P)ONE);            u = ptond(CO,vv,(P)ONE);
           HCM(u) = ((MQ)inv)->cont;            HCM(u) = ((MQ)inv)->cont;
           p[pi0] = u;            p[pi0] = u;
Legend:
Removed from v.1.216  
changed lines
  Added in v.1.218
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>