[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.58 and 1.60

version 1.58, 2003/09/05 07:00:37 version 1.60, 2003/09/06 01:33:26
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.57 2003/09/05 05:02:53 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.59 2003/09/05 13:20:14 noro Exp $ */
   
 #include "ca.h"  #include "ca.h"
 #include "inline.h"  #include "inline.h"
Line 136  extern int *current_weyl_weight_vector;
Line 136  extern int *current_weyl_weight_vector;
 #define LEN(a) ((a)->len)  #define LEN(a) ((a)->len)
 #define LCM(a) ((a)->lcm)  #define LCM(a) ((a)->lcm)
 #define GET_EXP(d,a) (((d)[nd_epos[a].i]>>nd_epos[a].s)&nd_mask0)  #define GET_EXP(d,a) (((d)[nd_epos[a].i]>>nd_epos[a].s)&nd_mask0)
   #define GET_EXP_MASK(d,a,m) ((((d)[nd_epos[a].i]&(m)[nd_epos[a].i])>>nd_epos[a].s)&nd_mask0)
 #define PUT_EXP(r,a,e) ((r)[nd_epos[a].i] |= ((e)<<nd_epos[a].s))  #define PUT_EXP(r,a,e) ((r)[nd_epos[a].i] |= ((e)<<nd_epos[a].s))
 #define XOR_EXP(r,a,e) ((r)[nd_epos[a].i] ^= ((e)<<nd_epos[a].s))  #define XOR_EXP(r,a,e) ((r)[nd_epos[a].i] ^= ((e)<<nd_epos[a].s))
   
Line 557  int ndl_weight(unsigned int *d)
Line 558  int ndl_weight(unsigned int *d)
         unsigned int t,u;          unsigned int t,u;
         int i,j;          int i,j;
   
         for ( t = 0, i = nd_exporigin; i < nd_wpd; i++ ) {          if ( current_dl_weight_vector )
                 u = d[i];                  for ( i = 0, t = 0; i < nd_nvar; i++ ) {
                 for ( j = 0; j < nd_epw; j++, u>>=nd_bpe )                          u = GET_EXP(d,i);
                         t += (u&nd_mask0);                          t += MUL_WEIGHT(u,i);
         }                  }
           else
                   for ( t = 0, i = nd_exporigin; i < nd_wpd; i++ ) {
                           u = d[i];
                           for ( j = 0; j < nd_epw; j++, u>>=nd_bpe )
                                   t += (u&nd_mask0);
                   }
         return t;          return t;
 }  }
   
Line 572  int ndl_weight_mask(unsigned int *d,int index)
Line 579  int ndl_weight_mask(unsigned int *d,int index)
         int i,j;          int i,j;
   
         mask = nd_blockmask->mask[index];          mask = nd_blockmask->mask[index];
         for ( t = 0, i = nd_exporigin; i < nd_wpd; i++ ) {          if ( current_dl_weight_vector )
                 u = d[i]&mask[i];                  for ( i = 0, t = 0; i < nd_nvar; i++ ) {
                 for ( j = 0; j < nd_epw; j++, u>>=nd_bpe )                          u = GET_EXP_MASK(d,i,mask);
                         t += (u&nd_mask0);                          t += MUL_WEIGHT(u,i);
         }                  }
           else
                   for ( t = 0, i = nd_exporigin; i < nd_wpd; i++ ) {
                           u = d[i]&mask[i];
                           for ( j = 0; j < nd_epw; j++, u>>=nd_bpe )
                                   t += (u&nd_mask0);
                   }
         return t;          return t;
 }  }
   
Line 3249  EPOS nd_create_epos(struct order_spec *ord)
Line 3262  EPOS nd_create_epos(struct order_spec *ord)
         }          }
         return epos;          return epos;
 }  }
   
   /* external interface */
   
   void nd_nf_p(P f,LIST g,LIST v,int m,struct order_spec *ord,P *rp)
   {
           NODE t;
           ND nd,nf;
           VL vv;
           DP d,r,dm;
           int stat;
   
           pltovl(v,&vv);
           nd_nvar = length(vv);
           nd_init_ord(ord);
           initd(ord);
           for ( t = BDY(g); NEXT(t); t = NEXT(t) ) {
                   ptod(CO,vv,(P)BDY(t),&d); BDY(t) = (pointer)d;
           }
           ptod(CO,vv,(P)BDY(t),&d); BDY(t) = (pointer)d;
           NEWNODE(NEXT(t)); ptod(CO,vv,f,&d); BDY(NEXT(t)) = d; NEXT(NEXT(t)) = 0;
           nd_setup(m,0,BDY(g));
           nd_psn--;
           nd_scale=2;
           while ( 1 ) {
                   nd = (pointer)ndvtond(m,nd_ps[nd_psn]);
                   stat = nd_nf(m,nd,nd_ps,1,&nf);
                   if ( !stat ) {
                           nd_psn++;
                           nd_reconstruct(m,0,0);
                           nd_psn--;
                   } else
                           break;
           }
           r = (DP)ndtodp(m,nf);
           if ( m ) _dtop_mod(CO,vv,r,rp);
           else dtop(CO,vv,r,rp);
   }
   

Legend:
Removed from v.1.58  
changed lines
  Added in v.1.60

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>