[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.24 and 1.26

version 1.24, 2002/01/28 00:54:41 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.23 2001/11/19 01:40:04 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 89  void Pdp_set_weight();
Line 88  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 119  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 1246  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 1264  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 1287  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 1324  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 1345  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 1369  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 1396  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 1413  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 1431  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 1453  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 1511  VECT *rp;
Line 1570  VECT *rp;
                 for ( i = 0; i < n; i++ )                  for ( i = 0; i < n; i++ )
                         current_weyl_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.24  
changed lines
  Added in v.1.26
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>