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

Diff for /OpenXM_contrib2/asir2000/engine/dist.c between version 1.42 and 1.43

version 1.42, 2006/08/27 22:17:27 version 1.43, 2007/09/17 12:47:45
Line 45 
Line 45 
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,   * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.   * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *   *
  * $OpenXM: OpenXM_contrib2/asir2000/engine/dist.c,v 1.41 2006/08/08 02:39:17 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/engine/dist.c,v 1.42 2006/08/27 22:17:27 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
   
Line 65 
Line 65 
 #define ORD_HOMO_WW_DRL_ZIGZAG 13  #define ORD_HOMO_WW_DRL_ZIGZAG 13
   
 int cmpdl_drl_zigzag(), cmpdl_homo_ww_drl_zigzag();  int cmpdl_drl_zigzag(), cmpdl_homo_ww_drl_zigzag();
   int cmpdl_top_weight();
   
 int (*cmpdl)()=cmpdl_revgradlex;  int (*cmpdl)()=cmpdl_revgradlex;
   int (*cmpdl_tie_breaker)();
 int (*primitive_cmpdl[3])() = {cmpdl_revgradlex,cmpdl_gradlex,cmpdl_lex};  int (*primitive_cmpdl[3])() = {cmpdl_revgradlex,cmpdl_gradlex,cmpdl_lex};
   
 int do_weyl;  int do_weyl;
Line 118  int has_sfcoef_p(P f)
Line 120  int has_sfcoef_p(P f)
         }          }
 }  }
   
   extern N *current_top_weight_vector;
   static int current_top_weight_len;
   
 void initd(struct order_spec *spec)  void initd(struct order_spec *spec)
 {  {
           int len,i;
   
         switch ( spec->id ) {          switch ( spec->id ) {
                 case 3:                  case 3:
                         cmpdl = cmpdl_composite;                          cmpdl = cmpdl_composite;
Line 165  void initd(struct order_spec *spec)
Line 172  void initd(struct order_spec *spec)
                         }                          }
                         break;                          break;
         }          }
           if ( current_top_weight_vector ) {
                   cmpdl_tie_breaker = cmpdl;
                   cmpdl = cmpdl_top_weight;
                   for ( len = 0, i = 0; i < spec->nv; i++ )
                           if ( current_top_weight_vector[i] )
                                   len = MAX(PL(current_top_weight_vector[i]),len);
                   current_top_weight_len = len;
           }
         dp_current_spec = spec;          dp_current_spec = spec;
 }  }
   
Line 1638  int cmpdl_matrix(int n,DL d1,DL d2)
Line 1653  int cmpdl_matrix(int n,DL d1,DL d2)
                         return -1;                          return -1;
         }          }
         return 0;          return 0;
   }
   
   int cmpdl_top_weight(int n,DL d1,DL d2)
   {
           int *w;
           N sum,wm,wma,t;
           N *a;
           struct oN tn;
           int len,i,sgn,tsgn;
           int *t1,*t2;
   
           w = (int *)ALLOCA(n*sizeof(int));
           len = current_top_weight_len+3;
           t1 = d1->d; t2 = d2->d;
           for ( i = 0; i < n; i++ ) w[i] = t1[i]-t2[i];
           sum = (N)W_ALLOC(len); sgn = 0;
           wm = (N)W_ALLOC(len);
           wma = (N)W_ALLOC(len);
           a = current_top_weight_vector;
           for ( i = 0; i < n; i++ ) {
                   if ( !a[i] || !w[i] ) continue;
                   tn.p = 1;
                   if ( w[i] > 0 ) {
                           tn.b[0] = w[i]; tsgn = 1;
                   } else {
                           tn.b[0] = -w[i]; tsgn = -1;
                   }
                   _muln(a[i],&tn,wm);
                   if ( !sgn ) {
                           sgn = tsgn;
                           t = wm; wm = sum; sum = t;
                   } else if ( sgn == tsgn ) {
                           _addn(sum,wm,wma);
                           if ( !PL(wma) )
                                   sgn = 0;
                           t = wma; wma = sum; sum = t;
                   } else {
                           sgn *= _subn(sum,wm,wma);
                           t = wma; wma = sum; sum = t;
                   }
           }
           if ( sgn > 0 ) return 1;
           else if ( sgn < 0 ) return -1;
           else return (*cmpdl_tie_breaker)(n,d1,d2);
 }  }
   
 GeoBucket create_bucket()  GeoBucket create_bucket()
Legend:
Removed from v.1.42  
changed lines
  Added in v.1.43
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>