[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.78 and 1.80

version 1.78, 2009/09/09 08:13:24 version 1.80, 2010/01/19 06:17:22
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.77 2009/01/05 02:08:17 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/dp.c,v 1.79 2009/10/09 04:02:11 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "base.h"  #include "base.h"
Line 68  void Pdp_set_sugar();
Line 68  void Pdp_set_sugar();
 void Pdp_cri1(),Pdp_cri2(),Pdp_subd(),Pdp_mod(),Pdp_red_mod(),Pdp_tdiv();  void Pdp_cri1(),Pdp_cri2(),Pdp_subd(),Pdp_mod(),Pdp_red_mod(),Pdp_tdiv();
 void Pdp_prim(),Pdp_red_coef(),Pdp_mag(),Pdp_set_kara(),Pdp_rat();  void Pdp_prim(),Pdp_red_coef(),Pdp_mag(),Pdp_set_kara(),Pdp_rat();
 void Pdp_nf(),Pdp_true_nf(),Pdp_true_nf_marked(),Pdp_true_nf_marked_mod();  void Pdp_nf(),Pdp_true_nf(),Pdp_true_nf_marked(),Pdp_true_nf_marked_mod();
 void Pdp_true_nf_and_quotient_marked();  void Pdp_true_nf_and_quotient_marked(),Pdp_true_nf_and_quotient_marked_mod();
 void Pdp_nf_mod(),Pdp_true_nf_mod();  void Pdp_nf_mod(),Pdp_true_nf_mod();
 void Pdp_criB(),Pdp_nelim();  void Pdp_criB(),Pdp_nelim();
 void Pdp_minp(),Pdp_sp_mod();  void Pdp_minp(),Pdp_sp_mod();
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 140  struct ftab dp_tab[] = {
Line 141  struct ftab dp_tab[] = {
         {"dp_true_nf",Pdp_true_nf,4},          {"dp_true_nf",Pdp_true_nf,4},
         {"dp_true_nf_marked",Pdp_true_nf_marked,4},          {"dp_true_nf_marked",Pdp_true_nf_marked,4},
         {"dp_true_nf_and_quotient_marked",Pdp_true_nf_and_quotient_marked,4},          {"dp_true_nf_and_quotient_marked",Pdp_true_nf_and_quotient_marked,4},
           {"dp_true_nf_and_quotient_marked_mod",Pdp_true_nf_and_quotient_marked_mod,5},
         {"dp_true_nf_marked_mod",Pdp_true_nf_marked_mod,5},          {"dp_true_nf_marked_mod",Pdp_true_nf_marked_mod,5},
         {"dp_nf_mod",Pdp_nf_mod,5},          {"dp_nf_mod",Pdp_nf_mod,5},
         {"dp_true_nf_mod",Pdp_true_nf_mod,5},          {"dp_true_nf_mod",Pdp_true_nf_mod,5},
Line 245  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 1032  LIST *rp;
Line 1035  LIST *rp;
         MKLIST(*rp,n);          MKLIST(*rp,n);
 }  }
   
   DP *dp_true_nf_and_quotient_marked_mod (NODE b,DP g,DP *ps,DP *hps,int mod,DP *rp,P *dnp);
   
   void Pdp_true_nf_and_quotient_marked_mod(arg,rp)
   NODE arg;
   LIST *rp;
   {
           NODE b,n;
           DP *ps,*hps;
           DP g;
           DP nm;
           VECT quo;
           P dn;
           int full,mod;
   
           do_weyl = 0; dp_fcoeffs = 0;
           asir_assert(ARG0(arg),O_LIST,"dp_true_nf_and_quotient_marked_mod");
           asir_assert(ARG1(arg),O_DP,"dp_true_nf_and_quotient_marked_mod");
           asir_assert(ARG2(arg),O_VECT,"dp_true_nf_and_quotient_marked_mod");
           asir_assert(ARG3(arg),O_VECT,"dp_true_nf_and_quotient_marked_mod");
           asir_assert(ARG4(arg),O_N,"dp_true_nf_and_quotient_marked_mod");
           if ( !(g = (DP)ARG1(arg)) ) {
                   nm = 0; dn = (P)ONE;
           } else {
                   b = BDY((LIST)ARG0(arg));
                   ps = (DP *)BDY((VECT)ARG2(arg));
                   hps = (DP *)BDY((VECT)ARG3(arg));
                   mod = QTOS((Q)ARG4(arg));
                   NEWVECT(quo); quo->len = ((VECT)ARG2(arg))->len;
                   quo->body = (pointer *)dp_true_nf_and_quotient_marked_mod(b,g,ps,hps,mod,&nm,&dn);
           }
           n = mknode(3,nm,dn,quo);
           MKLIST(*rp,n);
   }
   
 void Pdp_true_nf_marked_mod(arg,rp)  void Pdp_true_nf_marked_mod(arg,rp)
 NODE arg;  NODE arg;
 LIST *rp;  LIST *rp;
Line 1224  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 2146  LIST *rp;
Line 2207  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 2217  LIST *rp;
Line 2279  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);          nd_gr(f,v,m,0,ord,rp);
           do_weyl = 0;
 }  }
   
 void Pnd_weyl_gr_trace(arg,rp)  void Pnd_weyl_gr_trace(arg,rp)
Line 2240  LIST *rp;
Line 2303  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)
Legend:
Removed from v.1.78  
changed lines
  Added in v.1.80
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>