[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.10 and 1.12

version 1.10, 2000/11/07 06:06:39 version 1.12, 2000/12/11 02:00:41
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.9 2000/08/22 05:04:05 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/engine/dist.c,v 1.11 2000/12/05 06:59:16 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
   
 #define NV(p) ((p)->nv)  
 #define C(p) ((p)->c)  
   
 #define ORD_REVGRADLEX 0  #define ORD_REVGRADLEX 0
 #define ORD_GRADLEX 1  #define ORD_GRADLEX 1
 #define ORD_LEX 2  #define ORD_LEX 2
Line 62 
Line 59 
 #define ORD_BGRADREV 7  #define ORD_BGRADREV 7
 #define ORD_BLEXREV 8  #define ORD_BLEXREV 8
 #define ORD_ELIM 9  #define ORD_ELIM 9
   #define ORD_WEYL_ELIM 10
   
 struct cdl {  
         P c;  
         DL d;  
 };  
   
 int (*cmpdl)()=cmpdl_revgradlex;  int (*cmpdl)()=cmpdl_revgradlex;
 int (*primitive_cmpdl[3])() = {cmpdl_revgradlex,cmpdl_gradlex,cmpdl_lex};  int (*primitive_cmpdl[3])() = {cmpdl_revgradlex,cmpdl_gradlex,cmpdl_lex};
   
Line 79  void comm_muld_tab(VL,int,struct cdl *,int,struct cdl 
Line 72  void comm_muld_tab(VL,int,struct cdl *,int,struct cdl 
   
 void mkwc(int,int,Q *);  void mkwc(int,int,Q *);
   
   int cmpdl_weyl_elim();
   
 int do_weyl;  int do_weyl;
   
 int dp_nelim,dp_fcoeffs;  int dp_nelim,dp_fcoeffs;
Line 149  struct order_spec *spec;
Line 144  struct order_spec *spec;
                                         cmpdl = cmpdl_blexrev; break;                                          cmpdl = cmpdl_blexrev; break;
                                 case ORD_ELIM:                                  case ORD_ELIM:
                                         cmpdl = cmpdl_elim; break;                                          cmpdl = cmpdl_elim; break;
                                   case ORD_WEYL_ELIM:
                                           cmpdl = cmpdl_weyl_elim; break;
                                 case ORD_LEX: default:                                  case ORD_LEX: default:
                                         cmpdl = cmpdl_lex; break;                                          cmpdl = cmpdl_lex; break;
                         }                          }
Line 835  DL *dr;
Line 832  DL *dr;
         }          }
 }  }
   
   /* d1 += d2 */
   
   void adddl_destructive(n,d1,d2)
   int n;
   DL d1,d2;
   {
           DL dt;
           int i;
   
           d1->td += d2->td;
           for ( i = 0; i < n; i++ )
                   d1->d[i] += d2->d[i];
   }
   
 int compd(vl,p1,p2)  int compd(vl,p1,p2)
 VL vl;  VL vl;
 DP p1,p2;  DP p1,p2;
Line 1061  DL d1,d2;
Line 1072  DL d1,d2;
                 return -1;                  return -1;
         else          else
                 return cmpdl_revgradlex(n,d1,d2);                  return cmpdl_revgradlex(n,d1,d2);
   }
   
   int cmpdl_weyl_elim(n,d1,d2)
   int n;
   DL d1,d2;
   {
           int e1,e2,i;
   
           for ( i = 1, e1 = 0, e2 = 0; i <= dp_nelim; i++ ) {
                   e1 += d1->d[n-i]; e2 += d2->d[n-i];
           }
           if ( e1 > e2 )
                   return 1;
           else if ( e1 < e2 )
                   return -1;
           else if ( d1->td > d2->td )
                   return 1;
           else if ( d1->td < d2->td )
                   return -1;
           else return -cmpdl_revlex(n,d1,d2);
 }  }
   
 int cmpdl_order_pair(n,d1,d2)  int cmpdl_order_pair(n,d1,d2)
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.12
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>