[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.79 and 1.81

version 1.79, 2009/10/09 04:02:11 version 1.81, 2010/04/16 07:13:42
Line 44 
Line 44 
  * 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.78 2009/09/09 08:13:24 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/dp.c,v 1.80 2010/01/19 06:17:22 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "base.h"  #include "base.h"
Line 105  void Pdp_compute_last_t();
Line 105  void Pdp_compute_last_t();
 void Pdp_compute_last_w();  void Pdp_compute_last_w();
 void Pdp_compute_essential_df();  void Pdp_compute_essential_df();
 void Pdp_get_denomlist();  void Pdp_get_denomlist();
   void Pdp_symb_add();
   
 LIST dp_initial_term();  LIST dp_initial_term();
 LIST dp_order();  LIST dp_order();
Line 246  struct ftab dp_supp_tab[] = {
Line 247  struct ftab dp_supp_tab[] = {
         {"dp_rest",Pdp_rest,1},          {"dp_rest",Pdp_rest,1},
         {"dp_initial_term",Pdp_initial_term,1},          {"dp_initial_term",Pdp_initial_term,1},
         {"dp_order",Pdp_order,1},          {"dp_order",Pdp_order,1},
           {"dp_symb_add",Pdp_symb_add,2},
   
         /* degree and size */          /* degree and size */
         {"dp_td",Pdp_td,1},          {"dp_td",Pdp_td,1},
Line 1259  DP *rp;
Line 1261  DP *rp;
         dp_subd(p1,p2,rp);          dp_subd(p1,p2,rp);
 }  }
   
   void Pdp_symb_add(arg,rp)
   NODE arg;
   DP *rp;
   {
           DP p1,p2,r;
           NODE s0;
           MP mp0,mp;
           int nv;
   
           p1 = (DP)ARG0(arg); p2 = (DP)ARG1(arg);
           asir_assert(p1,O_DP,"dp_symb_add");
           asir_assert(p2,O_DP,"dp_symb_add");
           if ( p1->nv != p2->nv )
                   error("dp_sumb_add : invalid input");
           nv = p1->nv;
           s0 = symb_merge(dp_dllist(p1),dp_dllist(p2),nv);
           for ( mp0 = 0; s0; s0 = NEXT(s0) ) {
                   NEXTMP(mp0,mp); mp->dl = (DL)BDY(s0); mp->c = (P)ONE;
           }
           NEXT(mp) = 0;
           MKDP(nv,mp0,r); r->sugar = MAX(p1->sugar,p2->sugar);
           *rp = r;
   }
   
 void Pdp_mul_trunc(arg,rp)  void Pdp_mul_trunc(arg,rp)
 NODE arg;  NODE arg;
 DP *rp;  DP *rp;
Line 2104  NODE arg;
Line 2130  NODE arg;
 LIST *rp;  LIST *rp;
 {  {
         LIST f,v;          LIST f,v;
         int m,homo;          int m,find;
           Obj homo;
         struct order_spec *ord;          struct order_spec *ord;
   
         do_weyl = 0;          do_weyl = 0;
Line 2118  LIST *rp;
Line 2145  LIST *rp;
         }          }
         m = QTOS((Q)ARG2(arg));          m = QTOS((Q)ARG2(arg));
         create_order_spec(0,ARG3(arg),&ord);          create_order_spec(0,ARG3(arg),&ord);
         nd_gr(f,v,m,1,ord,rp);          find = get_opt("homo",&homo);
           nd_gr(f,v,m,find&&homo,1,ord,rp);
 }  }
   
 void Pnd_gr(arg,rp)  void Pnd_gr(arg,rp)
Line 2126  NODE arg;
Line 2154  NODE arg;
 LIST *rp;  LIST *rp;
 {  {
         LIST f,v;          LIST f,v;
         int m,homo;          int m,find;
           Obj homo;
         struct order_spec *ord;          struct order_spec *ord;
   
         do_weyl = 0;          do_weyl = 0;
Line 2140  LIST *rp;
Line 2169  LIST *rp;
         }          }
         m = QTOS((Q)ARG2(arg));          m = QTOS((Q)ARG2(arg));
         create_order_spec(0,ARG3(arg),&ord);          create_order_spec(0,ARG3(arg),&ord);
         nd_gr(f,v,m,0,ord,rp);          find = get_opt("homo",&homo);
           nd_gr(f,v,m,find&&homo,0,ord,rp);
 }  }
   
 void Pnd_gr_postproc(arg,rp)  void Pnd_gr_postproc(arg,rp)
Line 2181  LIST *rp;
Line 2211  LIST *rp;
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);
         f = remove_zero_from_list(f);          f = remove_zero_from_list(f);
         if ( !BDY(f) ) {          if ( !BDY(f) ) {
                 *rp = f; return;                  *rp = f; do_weyl = 0; return;
         }          }
         m = QTOS((Q)ARG2(arg));          m = QTOS((Q)ARG2(arg));
         create_order_spec(0,ARG3(arg),&ord);          create_order_spec(0,ARG3(arg),&ord);
         do_check = ARG4(arg) ? 1 : 0;          do_check = ARG4(arg) ? 1 : 0;
         nd_gr_postproc(f,v,m,ord,do_check,rp);          nd_gr_postproc(f,v,m,ord,do_check,rp);
           do_weyl = 0;
 }  }
   
 void Pnd_gr_trace(arg,rp)  void Pnd_gr_trace(arg,rp)
Line 2242  NODE arg;
Line 2273  NODE arg;
 LIST *rp;  LIST *rp;
 {  {
         LIST f,v;          LIST f,v;
         int m,homo;          int m,find;
           Obj homo;
         struct order_spec *ord;          struct order_spec *ord;
   
         do_weyl = 1;          do_weyl = 1;
Line 2252  LIST *rp;
Line 2284  LIST *rp;
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);
         f = remove_zero_from_list(f);          f = remove_zero_from_list(f);
         if ( !BDY(f) ) {          if ( !BDY(f) ) {
                 *rp = f; return;                  *rp = f; do_weyl = 0; return;
         }          }
         m = QTOS((Q)ARG2(arg));          m = QTOS((Q)ARG2(arg));
         create_order_spec(0,ARG3(arg),&ord);          create_order_spec(0,ARG3(arg),&ord);
         nd_gr(f,v,m,0,ord,rp);          find = get_opt("homo",&homo);
           nd_gr(f,v,m,find&&homo,0,ord,rp);
           do_weyl = 0;
 }  }
   
 void Pnd_weyl_gr_trace(arg,rp)  void Pnd_weyl_gr_trace(arg,rp)
Line 2275  LIST *rp;
Line 2309  LIST *rp;
         f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);          f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);
         f = remove_zero_from_list(f);          f = remove_zero_from_list(f);
         if ( !BDY(f) ) {          if ( !BDY(f) ) {
                 *rp = f; return;                  *rp = f; do_weyl = 0; return;
         }          }
         homo = QTOS((Q)ARG2(arg));          homo = QTOS((Q)ARG2(arg));
         m = QTOS((Q)ARG3(arg));          m = QTOS((Q)ARG3(arg));
         create_order_spec(0,ARG4(arg),&ord);          create_order_spec(0,ARG4(arg),&ord);
         nd_gr_trace(f,v,m,homo,0,ord,rp);          nd_gr_trace(f,v,m,homo,0,ord,rp);
           do_weyl = 0;
 }  }
   
 void Pnd_nf(arg,rp)  void Pnd_nf(arg,rp)
Line 2761  int dpv_hp(DPV p)
Line 2796  int dpv_hp(DPV p)
                         break;                          break;
         }          }
 }  }
   
   int get_opt(char *key0,Obj *r) {
      NODE tt,p;
      char *key;
   
      if ( current_option ) {
        for ( tt = current_option; tt; tt = NEXT(tt) ) {
          p = BDY((LIST)BDY(tt));
          key = BDY((STRING)BDY(p));
          /*  value = (Obj)BDY(NEXT(p)); */
          if ( !strcmp(key,key0) )  {
                *r = (Obj)BDY(NEXT(p));
                return 1;
              }
        }
      }
      return 0;
   }
   
Legend:
Removed from v.1.79  
changed lines
  Added in v.1.81
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>