[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.35 and 1.38

version 1.35, 2003/08/20 07:04:31 version 1.38, 2003/08/20 08:56:29
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.34 2003/08/20 06:06:07 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.37 2003/08/20 08:38:20 noro Exp $ */
   
 #include "ca.h"  #include "ca.h"
 #include "inline.h"  #include "inline.h"
Line 372  void ndl_dehomogenize(unsigned int *d)
Line 372  void ndl_dehomogenize(unsigned int *d)
                         h = (d[1]>>((nd_epw-1)*nd_bpe))&nd_mask0;                          h = (d[1]>>((nd_epw-1)*nd_bpe))&nd_mask0;
                         for ( i = 1; i <= nd_wpd; i++ )                          for ( i = 1; i <= nd_wpd; i++ )
                                 d[i] = ((d[i]<<nd_bpe)&mask)                                  d[i] = ((d[i]<<nd_bpe)&mask)
                                         |(i+1<nd_wpd?((d[i+1]>>((nd_epw-1)*nd_bpe))&nd_mask0):0);                                          |(i+1<=nd_wpd?((d[i+1]>>((nd_epw-1)*nd_bpe))&nd_mask0):0);
                         TD(d) -= h;                          TD(d) -= h;
                 }                  }
         } else          } else
Line 505  INLINE void ndl_copy(unsigned int *d1,unsigned int *d2
Line 505  INLINE void ndl_copy(unsigned int *d1,unsigned int *d2
         }          }
 }  }
   
 /* XXX : TD is not added */  
   
 INLINE void ndl_add(unsigned int *d1,unsigned int *d2,unsigned int *d)  INLINE void ndl_add(unsigned int *d1,unsigned int *d2,unsigned int *d)
 {  {
         int i;          int i;
   
           TD(d) = TD(d1)+TD(d2);
         switch ( nd_wpd ) {          switch ( nd_wpd ) {
                 case 1:                  case 1:
                         d[1] = d1[1]+d2[1];                          d[1] = d1[1]+d2[1];
Line 1515  again:
Line 1514  again:
                                 printf("+"); fflush(stdout);                                  printf("+"); fflush(stdout);
                                 nh = nd_newps_trace(m,nf,nfq);                                  nh = nd_newps_trace(m,nf,nfq);
                                 /* failure; m|HC(nfq) */                                  /* failure; m|HC(nfq) */
                                 if ( nf < 0 ) return 0;                                  if ( nh < 0 ) return 0;
                                 d = update_pairs(d,g,nh);                                  d = update_pairs(d,g,nh);
                                 g = update_base(g,nh);                                  g = update_base(g,nh);
                         } else {                          } else {
Line 2472  void nd_setup_parameters() {
Line 2471  void nd_setup_parameters() {
   
 ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d)  ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d)
 {  {
         int i,obpe,oadv;          int i,obpe,oadv,h;
         NM prev_nm_free_list;          NM prev_nm_free_list;
         RHist mr0,mr;          RHist mr0,mr;
         RHist r;          RHist r;
           RHist *old_red;
         ND_pairs s0,s,t,prev_ndp_free_list;          ND_pairs s0,s,t,prev_ndp_free_list;
   
         obpe = nd_bpe;          obpe = nd_bpe;
Line 2504  ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d)
Line 2504  ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d)
                 SG(s) = SG(t);                  SG(s) = SG(t);
                 ndl_dup(obpe,LCM(t),LCM(s));                  ndl_dup(obpe,LCM(t),LCM(s));
         }          }
   
           old_red = (RHist *)ALLOCA(REDTAB_LEN*sizeof(RHist));
         for ( i = 0; i < REDTAB_LEN; i++ ) {          for ( i = 0; i < REDTAB_LEN; i++ ) {
                 for ( mr0 = 0, r = nd_red[i]; r; r = NEXT(r) ) {                  old_red[i] = nd_red[i];
                         NEXTRHist(mr0,mr);                  nd_red[i] = 0;
           }
           for ( i = 0; i < REDTAB_LEN; i++ )
                   for ( r = old_red[i]; r; r = NEXT(r) ) {
                           NEWRHist(mr);
                         mr->index = r->index;                          mr->index = r->index;
                         SG(mr) = SG(r);                          SG(mr) = SG(r);
                         ndl_dup(obpe,DL(r),DL(mr));                          ndl_dup(obpe,DL(r),DL(mr));
                           h = ndl_hash_value(DL(mr));
                           NEXT(mr) = nd_red[h];
                           nd_red[h] = mr;
                 }                  }
                 if ( mr0 ) NEXT(mr) = 0;          for ( i = 0; i < REDTAB_LEN; i++ ) old_red[i] = 0;
                 nd_red[i] = mr0;          old_red = 0;
         }  
         for ( i = 0; i < nd_psn; i++ ) {          for ( i = 0; i < nd_psn; i++ ) {
                 NEWRHist(r); SG(r) = SG(nd_psh[i]);                  NEWRHist(r); SG(r) = SG(nd_psh[i]);
                 ndl_dup(obpe,DL(nd_psh[i]),DL(r));                  ndl_dup(obpe,DL(nd_psh[i]),DL(r));
Line 2528  ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d)
Line 2536  ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d)
   
 void nd_reconstruct_direct(int mod,NDV *ps,int len)  void nd_reconstruct_direct(int mod,NDV *ps,int len)
 {  {
         int i,obpe,oadv;          int i,obpe,oadv,h;
         NM prev_nm_free_list;          NM prev_nm_free_list;
         RHist mr0,mr;          RHist mr0,mr;
         RHist r;          RHist r;
           RHist *old_red;
         ND_pairs s0,s,t,prev_ndp_free_list;          ND_pairs s0,s,t,prev_ndp_free_list;
   
         obpe = nd_bpe;          obpe = nd_bpe;
Line 2549  void nd_reconstruct_direct(int mod,NDV *ps,int len)
Line 2558  void nd_reconstruct_direct(int mod,NDV *ps,int len)
         _nm_free_list = 0; _ndp_free_list = 0;          _nm_free_list = 0; _ndp_free_list = 0;
         if ( mod != 0 )          if ( mod != 0 )
                 for ( i = len-1; i >= 0; i-- ) ndv_realloc(ps[i],obpe,oadv);                  for ( i = len-1; i >= 0; i-- ) ndv_realloc(ps[i],obpe,oadv);
           old_red = (RHist *)ALLOCA(REDTAB_LEN*sizeof(RHist));
         for ( i = 0; i < REDTAB_LEN; i++ ) {          for ( i = 0; i < REDTAB_LEN; i++ ) {
                 for ( mr0 = 0, r = nd_red[i]; r; r = NEXT(r) ) {                  old_red[i] = nd_red[i];
                         NEXTRHist(mr0,mr);                  nd_red[i] = 0;
           }
           for ( i = 0; i < REDTAB_LEN; i++ )
                   for ( r = old_red[i]; r; r = NEXT(r) ) {
                           NEWRHist(mr);
                         mr->index = r->index;                          mr->index = r->index;
                         SG(mr) = SG(r);                          SG(mr) = SG(r);
                         ndl_dup(obpe,DL(r),DL(mr));                          ndl_dup(obpe,DL(r),DL(mr));
                           h = ndl_hash_value(DL(mr));
                           NEXT(mr) = nd_red[h];
                           nd_red[h] = mr;
                 }                  }
                 if ( mr0 ) NEXT(mr) = 0;          for ( i = 0; i < REDTAB_LEN; i++ ) old_red[i] = 0;
                 nd_red[i] = mr0;          old_red = 0;
         }  
         prev_nm_free_list = 0;          prev_nm_free_list = 0;
         prev_ndp_free_list = 0;          prev_ndp_free_list = 0;
         GC_gcollect();          GC_gcollect();
Line 2706  ND ndv_mul_nm(int mod,NDV p,NM m0)
Line 2722  ND ndv_mul_nm(int mod,NDV p,NM m0)
                                 c1 = CM(m);                                  c1 = CM(m);
                                 DMAR(c1,c,0,mod,c2);                                  DMAR(c1,c,0,mod,c2);
                                 CM(mr) = c2;                                  CM(mr) = c2;
                                 TD(DL(mr)) = TD(DL(m))+td;  
                                 ndl_add(DL(m),d,DL(mr));                                  ndl_add(DL(m),d,DL(mr));
                         }                          }
                 } else {                  } else {
Line 2714  ND ndv_mul_nm(int mod,NDV p,NM m0)
Line 2729  ND ndv_mul_nm(int mod,NDV p,NM m0)
                         for ( i = 0; i < len; i++, NMV_ADV(m) ) {                          for ( i = 0; i < len; i++, NMV_ADV(m) ) {
                                 NEXTNM(mr0,mr);                                  NEXTNM(mr0,mr);
                                 mulq(CQ(m),q,&CQ(mr));                                  mulq(CQ(m),q,&CQ(mr));
                                 TD(DL(mr)) = TD(DL(m))+td;  
                                 ndl_add(DL(m),d,DL(mr));                                  ndl_add(DL(m),d,DL(mr));
                         }                          }
                 }                  }
Legend:
Removed from v.1.35  
changed lines
  Added in v.1.38
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>