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

Diff for /OpenXM_contrib2/asir2000/builtin/dp.c between version 1.23 and 1.26

version 1.23, 2001/11/19 01:40:04 version 1.26, 2003/01/04 09:06:15
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/builtin/dp.c,v 1.22 2001/11/19 00:57:10 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/dp.c,v 1.25 2002/06/06 01:18:05 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "base.h"  #include "base.h"
 #include "parse.h"  #include "parse.h"
   
 extern int dp_fcoeffs;  
 extern int dp_nelim;  extern int dp_nelim;
 extern int dp_order_pair_length;  extern int dp_order_pair_length;
 extern struct order_pair *dp_order_pair;  extern struct order_pair *dp_order_pair;
Line 85  void Pdp_weyl_gr_main(),Pdp_weyl_gr_mod_main(),Pdp_wey
Line 84  void Pdp_weyl_gr_main(),Pdp_weyl_gr_mod_main(),Pdp_wey
 void Pdp_weyl_f4_main(),Pdp_weyl_f4_mod_main(),Pdp_weyl_f4_f_main();  void Pdp_weyl_f4_main(),Pdp_weyl_f4_mod_main(),Pdp_weyl_f4_f_main();
 void Pdp_weyl_mul(),Pdp_weyl_mul_mod();  void Pdp_weyl_mul(),Pdp_weyl_mul_mod();
 void Pdp_weyl_set_weight();  void Pdp_weyl_set_weight();
   void Pdp_set_weight();
 void Pdp_nf_f(),Pdp_weyl_nf_f();  void Pdp_nf_f(),Pdp_weyl_nf_f();
 void Pdp_lnf_f();  void Pdp_lnf_f();
   
   LIST remove_zero_from_list(LIST);
   
 struct ftab dp_tab[] = {  struct ftab dp_tab[] = {
         /* content reduction */          /* content reduction */
         {"dp_ptozp",Pdp_ptozp,1},          {"dp_ptozp",Pdp_ptozp,1},
Line 118  struct ftab dp_tab[] = {
Line 120  struct ftab dp_tab[] = {
         /* Buchberger algorithm */          /* Buchberger algorithm */
         {"dp_gr_main",Pdp_gr_main,5},          {"dp_gr_main",Pdp_gr_main,5},
         {"dp_gr_mod_main",Pdp_gr_mod_main,5},          {"dp_gr_mod_main",Pdp_gr_mod_main,5},
         {"dp_gr_f_main",Pdp_gr_f_main,4},          {"dp_gr_f_main",Pdp_gr_f_main,5},
         {"dp_gr_checklist",Pdp_gr_checklist,2},          {"dp_gr_checklist",Pdp_gr_checklist,2},
   
         /* F4 algorithm */          /* F4 algorithm */
Line 151  struct ftab dp_tab[] = {
Line 153  struct ftab dp_tab[] = {
         {"dp_weyl_f4_mod_main",Pdp_weyl_f4_mod_main,4},          {"dp_weyl_f4_mod_main",Pdp_weyl_f4_mod_main,4},
   
         /* misc */          /* misc */
           {"dp_set_weight",Pdp_set_weight,-1},
         {"dp_weyl_set_weight",Pdp_weyl_set_weight,-1},          {"dp_weyl_set_weight",Pdp_weyl_set_weight,-1},
         {0,0,0},          {0,0,0},
 };  };
Line 320  DP *rp;
Line 323  DP *rp;
         n = v->len;          n = v->len;
         NEWDL(dl,n); d = dl->d;          NEWDL(dl,n); d = dl->d;
         for ( i = 0, td = 0; i < n; i++ ) {          for ( i = 0, td = 0; i < n; i++ ) {
                 d[i] = QTOS((Q)(v->body[i])); td += d[i];                  d[i] = QTOS((Q)(v->body[i])); td += MUL_WEIGHT(d[i],i);
         }          }
         dl->td = td;          dl->td = td;
         NEWMP(m); m->dl = dl; m->c = (P)ONE; NEXT(m) = 0;          NEWMP(m); m->dl = dl; m->c = (P)ONE; NEXT(m) = 0;
Line 929  DP *rp;
Line 932  DP *rp;
         n = p1->nv; d1 = BDY(p1)->dl; d2 = BDY(p2)->dl;          n = p1->nv; d1 = BDY(p1)->dl; d2 = BDY(p2)->dl;
         NEWDL(d,n);          NEWDL(d,n);
         for ( i = 0, td = 0; i < n; i++ ) {          for ( i = 0, td = 0; i < n; i++ ) {
                 d->d[i] = MAX(d1->d[i],d2->d[i]); td += d->d[i];                  d->d[i] = MAX(d1->d[i],d2->d[i]); td += MUL_WEIGHT(d->d[i],i);
         }          }
         d->td = td;          d->td = td;
         NEWMP(m); m->dl = d; m->c = (P)ONE; NEXT(m) = 0;          NEWMP(m); m->dl = d; m->c = (P)ONE; NEXT(m) = 0;
Line 1244  LIST *rp;
Line 1247  LIST *rp;
         asir_assert(ARG2(arg),O_N,"dp_gr_main");          asir_assert(ARG2(arg),O_N,"dp_gr_main");
         asir_assert(ARG3(arg),O_N,"dp_gr_main");          asir_assert(ARG3(arg),O_N,"dp_gr_main");
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);
           f = remove_zero_from_list(f);
           if ( !BDY(f) ) {
                   *rp = f; return;
           }
         homo = (Num)ARG2(arg);          homo = (Num)ARG2(arg);
         m = (Q)ARG3(arg);          m = (Q)ARG3(arg);
         if ( !m )          if ( !m )
Line 1262  LIST *rp;
Line 1269  LIST *rp;
 {  {
         LIST f,v;          LIST f,v;
         Num homo;          Num homo;
           int m,field,t;
         struct order_spec ord;          struct order_spec ord;
           NODE n;
   
         do_weyl = 0;          do_weyl = 0;
         asir_assert(ARG0(arg),O_LIST,"dp_gr_f_main");          asir_assert(ARG0(arg),O_LIST,"dp_gr_f_main");
         asir_assert(ARG1(arg),O_LIST,"dp_gr_f_main");          asir_assert(ARG1(arg),O_LIST,"dp_gr_f_main");
         asir_assert(ARG2(arg),O_N,"dp_gr_f_main");          asir_assert(ARG2(arg),O_N,"dp_gr_f_main");
           asir_assert(ARG3(arg),O_N,"dp_gr_f_main");
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);
           f = remove_zero_from_list(f);
           if ( !BDY(f) ) {
                   *rp = f; return;
           }
         homo = (Num)ARG2(arg);          homo = (Num)ARG2(arg);
         create_order_spec(ARG3(arg),&ord);          m = QTOS((Q)ARG3(arg));
         dp_gr_main(f,v,homo,0,1,&ord,rp);          if ( m )
                   error("dp_gr_f_main : trace lifting is not implemented yet");
           create_order_spec(ARG4(arg),&ord);
           field = 0;
           for ( n = BDY(f); n; n = NEXT(n) ) {
                   t = get_field_type(BDY(n));
                   if ( !t )
                           continue;
                   if ( t < 0 )
                           error("dp_gr_f_main : incosistent coefficients");
                   if ( !field )
                           field = t;
                   else if ( t != field )
                           error("dp_gr_f_main : incosistent coefficients");
           }
           dp_gr_main(f,v,homo,m?1:0,field,&ord,rp);
 }  }
   
 void Pdp_f4_main(arg,rp)  void Pdp_f4_main(arg,rp)
Line 1285  LIST *rp;
Line 1314  LIST *rp;
         asir_assert(ARG0(arg),O_LIST,"dp_f4_main");          asir_assert(ARG0(arg),O_LIST,"dp_f4_main");
         asir_assert(ARG1(arg),O_LIST,"dp_f4_main");          asir_assert(ARG1(arg),O_LIST,"dp_f4_main");
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);
           f = remove_zero_from_list(f);
           if ( !BDY(f) ) {
                   *rp = f; return;
           }
         create_order_spec(ARG2(arg),&ord);          create_order_spec(ARG2(arg),&ord);
         dp_f4_main(f,v,&ord,rp);          dp_f4_main(f,v,&ord,rp);
 }  }
Line 1322  LIST *rp;
Line 1355  LIST *rp;
         asir_assert(ARG1(arg),O_LIST,"dp_f4_mod_main");          asir_assert(ARG1(arg),O_LIST,"dp_f4_mod_main");
         asir_assert(ARG2(arg),O_N,"dp_f4_mod_main");          asir_assert(ARG2(arg),O_N,"dp_f4_mod_main");
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg); m = QTOS((Q)ARG2(arg));          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg); m = QTOS((Q)ARG2(arg));
           f = remove_zero_from_list(f);
           if ( !BDY(f) ) {
                   *rp = f; return;
           }
         if ( !m )          if ( !m )
                 error("dp_f4_mod_main : invalid argument");                  error("dp_f4_mod_main : invalid argument");
         create_order_spec(ARG3(arg),&ord);          create_order_spec(ARG3(arg),&ord);
Line 1343  LIST *rp;
Line 1380  LIST *rp;
         asir_assert(ARG2(arg),O_N,"dp_gr_mod_main");          asir_assert(ARG2(arg),O_N,"dp_gr_mod_main");
         asir_assert(ARG3(arg),O_N,"dp_gr_mod_main");          asir_assert(ARG3(arg),O_N,"dp_gr_mod_main");
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);
           f = remove_zero_from_list(f);
           if ( !BDY(f) ) {
                   *rp = f; return;
           }
         homo = (Num)ARG2(arg); m = QTOS((Q)ARG3(arg));          homo = (Num)ARG2(arg); m = QTOS((Q)ARG3(arg));
         if ( !m )          if ( !m )
                 error("dp_gr_mod_main : invalid argument");                  error("dp_gr_mod_main : invalid argument");
Line 1367  LIST *rp;
Line 1408  LIST *rp;
         asir_assert(ARG2(arg),O_N,"dp_weyl_gr_main");          asir_assert(ARG2(arg),O_N,"dp_weyl_gr_main");
         asir_assert(ARG3(arg),O_N,"dp_weyl_gr_main");          asir_assert(ARG3(arg),O_N,"dp_weyl_gr_main");
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);
           f = remove_zero_from_list(f);
           if ( !BDY(f) ) {
                   *rp = f; return;
           }
         homo = (Num)ARG2(arg);          homo = (Num)ARG2(arg);
         m = (Q)ARG3(arg);          m = (Q)ARG3(arg);
         if ( !m )          if ( !m )
Line 1394  LIST *rp;
Line 1439  LIST *rp;
         asir_assert(ARG2(arg),O_N,"dp_weyl_gr_main");          asir_assert(ARG2(arg),O_N,"dp_weyl_gr_main");
         asir_assert(ARG3(arg),O_N,"dp_weyl_gr_main");          asir_assert(ARG3(arg),O_N,"dp_weyl_gr_main");
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);
           f = remove_zero_from_list(f);
           if ( !BDY(f) ) {
                   *rp = f; return;
           }
         homo = (Num)ARG2(arg);          homo = (Num)ARG2(arg);
         create_order_spec(ARG3(arg),&ord);          create_order_spec(ARG3(arg),&ord);
         do_weyl = 1;          do_weyl = 1;
Line 1411  LIST *rp;
Line 1460  LIST *rp;
         asir_assert(ARG0(arg),O_LIST,"dp_weyl_f4_main");          asir_assert(ARG0(arg),O_LIST,"dp_weyl_f4_main");
         asir_assert(ARG1(arg),O_LIST,"dp_weyl_f4_main");          asir_assert(ARG1(arg),O_LIST,"dp_weyl_f4_main");
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);
           f = remove_zero_from_list(f);
           if ( !BDY(f) ) {
                   *rp = f; return;
           }
         create_order_spec(ARG2(arg),&ord);          create_order_spec(ARG2(arg),&ord);
         do_weyl = 1;          do_weyl = 1;
         dp_f4_main(f,v,&ord,rp);          dp_f4_main(f,v,&ord,rp);
Line 1429  LIST *rp;
Line 1482  LIST *rp;
         asir_assert(ARG1(arg),O_LIST,"dp_weyl_f4_main");          asir_assert(ARG1(arg),O_LIST,"dp_weyl_f4_main");
         asir_assert(ARG2(arg),O_N,"dp_f4_main");          asir_assert(ARG2(arg),O_N,"dp_f4_main");
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg); m = QTOS((Q)ARG2(arg));          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg); m = QTOS((Q)ARG2(arg));
           f = remove_zero_from_list(f);
           if ( !BDY(f) ) {
                   *rp = f; return;
           }
         if ( !m )          if ( !m )
                 error("dp_weyl_f4_mod_main : invalid argument");                  error("dp_weyl_f4_mod_main : invalid argument");
         create_order_spec(ARG3(arg),&ord);          create_order_spec(ARG3(arg),&ord);
Line 1451  LIST *rp;
Line 1508  LIST *rp;
         asir_assert(ARG2(arg),O_N,"dp_weyl_gr_mod_main");          asir_assert(ARG2(arg),O_N,"dp_weyl_gr_mod_main");
         asir_assert(ARG3(arg),O_N,"dp_weyl_gr_mod_main");          asir_assert(ARG3(arg),O_N,"dp_weyl_gr_mod_main");
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);
           f = remove_zero_from_list(f);
           if ( !BDY(f) ) {
                   *rp = f; return;
           }
         homo = (Num)ARG2(arg); m = QTOS((Q)ARG3(arg));          homo = (Num)ARG2(arg); m = QTOS((Q)ARG3(arg));
         if ( !m )          if ( !m )
                 error("dp_weyl_gr_mod_main : invalid argument");                  error("dp_weyl_gr_mod_main : invalid argument");
Line 1460  LIST *rp;
Line 1521  LIST *rp;
         do_weyl = 0;          do_weyl = 0;
 }  }
   
 static VECT current_weight_vector_obj;  static VECT current_dl_weight_vector_obj;
 int *current_weight_vector;  int *current_dl_weight_vector;
   
   void Pdp_set_weight(arg,rp)
   NODE arg;
   VECT *rp;
   {
           VECT v;
           int i,n;
   
           if ( !arg )
                   *rp = current_dl_weight_vector_obj;
           else if ( !ARG0(arg) ) {
                   current_dl_weight_vector_obj = 0;
                   current_dl_weight_vector = 0;
                   *rp = 0;
           } else {
                   asir_assert(ARG0(arg),O_VECT,"dp_set_weight");
                   v = (VECT)ARG0(arg);
                   current_dl_weight_vector_obj = v;
                   n = v->len;
                   current_dl_weight_vector = (int *)CALLOC(n,sizeof(int));
                   for ( i = 0; i < n; i++ )
                           current_dl_weight_vector[i] = QTOS((Q)v->body[i]);
                   *rp = v;
           }
   }
   
   static VECT current_weyl_weight_vector_obj;
   int *current_weyl_weight_vector;
   
 void Pdp_weyl_set_weight(arg,rp)  void Pdp_weyl_set_weight(arg,rp)
 NODE arg;  NODE arg;
 VECT *rp;  VECT *rp;
Line 1471  VECT *rp;
Line 1560  VECT *rp;
         int i,n;          int i,n;
   
         if ( !arg )          if ( !arg )
                 *rp = current_weight_vector_obj;                  *rp = current_weyl_weight_vector_obj;
         else {          else {
                 asir_assert(ARG0(arg),O_VECT,"dp_weyl_set_weight");                  asir_assert(ARG0(arg),O_VECT,"dp_weyl_set_weight");
                 v = (VECT)ARG0(arg);                  v = (VECT)ARG0(arg);
                 current_weight_vector_obj = v;                  current_weyl_weight_vector_obj = v;
                 n = v->len;                  n = v->len;
                 current_weight_vector = (int *)CALLOC(n,sizeof(int));                  current_weyl_weight_vector = (int *)CALLOC(n,sizeof(int));
                 for ( i = 0; i < n; i++ )                  for ( i = 0; i < n; i++ )
                         current_weight_vector[i] = QTOS((Q)v->body[i]);                          current_weyl_weight_vector[i] = QTOS((Q)v->body[i]);
                 *rp = v;                  *rp = v;
           }
   }
   
   LIST remove_zero_from_list(LIST l)
   {
           NODE n,r0,r;
           LIST rl;
   
           asir_assert(l,O_LIST,"remove_zero_from_list");
           n = BDY(l);
           for ( r0 = 0; n; n = NEXT(n) )
                   if ( BDY(n) ) {
                           NEXTNODE(r0,r);
                           BDY(r) = BDY(n);
                   }
           if ( r0 )
                   NEXT(r) = 0;
           MKLIST(rl,r0);
           return rl;
   }
   
   int get_field_type(P p)
   {
           int type,t;
           DCP dc;
   
           if ( !p )
                   return 0;
           else if ( NUM(p) )
                   return NID((Num)p);
           else {
                   type = 0;
                   for ( dc = DC(p); dc; dc = NEXT(dc) ) {
                           t = get_field_type(COEF(dc));
                           if ( !t )
                                   continue;
                           if ( t < 0 )
                                   return t;
                           if ( !type )
                                   type = t;
                           else if ( t != type )
                                   return -1;
                   }
                   return type;
         }          }
 }  }
Legend:
Removed from v.1.23  
changed lines
  Added in v.1.26
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>