[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.72 and 1.75

version 1.72, 2007/09/19 05:42:59 version 1.75, 2008/01/06 00:36:32
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.71 2007/09/17 12:47:45 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/dp.c,v 1.74 2007/10/21 07:47:59 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_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 94  void Pdp_set_weight(),Pdp_set_top_weight();
Line 95  void Pdp_set_weight(),Pdp_set_top_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();
 void Pnd_gr(),Pnd_gr_trace(),Pnd_f4(),Pnd_f4_trace();  void Pnd_gr(),Pnd_gr_trace(),Pnd_f4(),Pnd_f4_trace();
 void Pnd_gr_postproc();  void Pnd_gr_postproc(), Pnd_weyl_gr_postproc();
 void Pnd_weyl_gr(),Pnd_weyl_gr_trace();  void Pnd_weyl_gr(),Pnd_weyl_gr_trace();
 void Pnd_nf();  void Pnd_nf();
 void Pdp_initial_term();  void Pdp_initial_term();
 void Pdp_order();  void Pdp_order();
 void Pdp_inv_or_split();  void Pdp_inv_or_split();
   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();
Line 137  struct ftab dp_tab[] = {
Line 139  struct ftab dp_tab[] = {
         {"dp_nf_f",Pdp_nf_f,4},          {"dp_nf_f",Pdp_nf_f,4},
         {"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_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 156  struct ftab dp_tab[] = {
Line 159  struct ftab dp_tab[] = {
         {"nd_gr_trace",Pnd_gr_trace,5},          {"nd_gr_trace",Pnd_gr_trace,5},
         {"nd_f4_trace",Pnd_f4_trace,5},          {"nd_f4_trace",Pnd_f4_trace,5},
         {"nd_gr_postproc",Pnd_gr_postproc,5},          {"nd_gr_postproc",Pnd_gr_postproc,5},
           {"nd_weyl_gr_postproc",Pnd_weyl_gr_postproc,5},
         {"nd_weyl_gr",Pnd_weyl_gr,4},          {"nd_weyl_gr",Pnd_weyl_gr,4},
         {"nd_weyl_gr_trace",Pnd_weyl_gr_trace,5},          {"nd_weyl_gr_trace",Pnd_weyl_gr_trace,5},
         {"nd_nf",Pnd_nf,5},          {"nd_nf",Pnd_nf,5},
Line 255  struct ftab dp_supp_tab[] = {
Line 259  struct ftab dp_supp_tab[] = {
         {"dp_tdiv",Pdp_tdiv,2},          {"dp_tdiv",Pdp_tdiv,2},
         {"dp_minp",Pdp_minp,2},          {"dp_minp",Pdp_minp,2},
         {"dp_compute_last_w",Pdp_compute_last_w,5},          {"dp_compute_last_w",Pdp_compute_last_w,5},
           {"dp_compute_last_t",Pdp_compute_last_t,5},
         {"dp_compute_essential_df",Pdp_compute_essential_df,2},          {"dp_compute_essential_df",Pdp_compute_essential_df,2},
   
         {0,0,0}          {0,0,0}
 };  };
   
 NODE compute_last_w(NODE g,NODE gh,int n,int **v,int row1,int **m1,int row2,int **m2);  NODE compute_last_w(NODE g,NODE gh,int n,int **v,int row1,int **m1,int row2,int **m2);
   Q compute_last_t(NODE g,NODE gh,Q t,VECT w1,VECT w2,NODE *homo,VECT *wp);
   
   void Pdp_compute_last_t(NODE arg,LIST *rp)
   {
           NODE g,gh,homo,n;
           LIST hlist;
           VECT v1,v2,w;
           Q t;
   
           g = (NODE)BDY((LIST)ARG0(arg));
           gh = (NODE)BDY((LIST)ARG1(arg));
           t = (Q)ARG2(arg);
           v1 = (VECT)ARG3(arg);
           v2 = (VECT)ARG4(arg);
           t = compute_last_t(g,gh,t,v1,v2,&homo,&w);
           MKLIST(hlist,homo);
           n = mknode(3,t,w,hlist);
           MKLIST(*rp,n);
   }
   
 void Pdp_compute_last_w(NODE arg,LIST *rp)  void Pdp_compute_last_w(NODE arg,LIST *rp)
 {  {
         NODE g,gh,r;          NODE g,gh,r;
Line 975  LIST *rp;
Line 999  LIST *rp;
         MKLIST(*rp,n);          MKLIST(*rp,n);
 }  }
   
   DP *dp_true_nf_and_quotient_marked (NODE b,DP g,DP *ps,DP *hps,DP *rp,P *dnp);
   
   void Pdp_true_nf_and_quotient_marked(arg,rp)
   NODE arg;
   LIST *rp;
   {
           NODE b,n;
           DP *ps,*hps;
           DP g;
           DP nm;
           VECT quo;
           P dn;
           int full;
   
           do_weyl = 0; dp_fcoeffs = 0;
           asir_assert(ARG0(arg),O_LIST,"dp_true_nf_and_quotient_marked");
           asir_assert(ARG1(arg),O_DP,"dp_true_nf_and_quotient_marked");
           asir_assert(ARG2(arg),O_VECT,"dp_true_nf_and_quotient_marked");
           asir_assert(ARG3(arg),O_VECT,"dp_true_nf_and_quotient_marked");
           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));
                   NEWVECT(quo); quo->len = ((VECT)ARG2(arg))->len;
                   quo->body = (pointer *)dp_true_nf_and_quotient_marked(b,g,ps,hps,&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 2074  LIST *rp;
Line 2130  LIST *rp;
         nd_gr_postproc(f,v,m,ord,do_check,rp);          nd_gr_postproc(f,v,m,ord,do_check,rp);
 }  }
   
   void Pnd_weyl_gr_postproc(arg,rp)
   NODE arg;
   LIST *rp;
   {
           LIST f,v;
           int m,do_check;
           struct order_spec *ord;
   
           do_weyl = 1;
           asir_assert(ARG0(arg),O_LIST,"nd_gr");
           asir_assert(ARG1(arg),O_LIST,"nd_gr");
           asir_assert(ARG2(arg),O_N,"nd_gr");
           f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);
           f = remove_zero_from_list(f);
           if ( !BDY(f) ) {
                   *rp = f; return;
           }
           m = QTOS((Q)ARG2(arg));
           create_order_spec(0,ARG3(arg),&ord);
           do_check = ARG4(arg) ? 1 : 0;
           nd_gr_postproc(f,v,m,ord,do_check,rp);
   }
   
 void Pnd_gr_trace(arg,rp)  void Pnd_gr_trace(arg,rp)
 NODE arg;  NODE arg;
 LIST *rp;  LIST *rp;
Line 2398  VECT *rp;
Line 2477  VECT *rp;
                                 BDY(v)[i] = BDY(node);                                  BDY(v)[i] = BDY(node);
                 }                  }
                 for ( i = 0; i < v->len; i++ )                  for ( i = 0; i < v->len; i++ )
                         if ( !INT(BDY(v)[i]) || SGN((Q)BDY(v)[i]) < 0 )                          if ( !INT(BDY(v)[i]) || (BDY(v)[i] && SGN((Q)BDY(v)[i]) < 0) )
                                 error("dp_set_top_weight : each element must be a non-negative integer");                                  error("dp_set_top_weight : each element must be a non-negative integer");
                 current_top_weight_vector_obj = v;                  current_top_weight_vector_obj = v;
                 current_top_weight_vector = (N *)MALLOC(v->len*sizeof(N));                  current_top_weight_vector = (N *)MALLOC(v->len*sizeof(N));
                 for ( i = 0; i < v->len; i++ ) {                  for ( i = 0; i < v->len; i++ ) {
                         current_top_weight_vector[i] = NM((Q)BDY(v)[i]);                          current_top_weight_vector[i] = !BDY(v)[i]?0:NM((Q)BDY(v)[i]);
                 }                  }
                 *rp = current_top_weight_vector_obj;                  *rp = current_top_weight_vector_obj;
         }          }
Legend:
Removed from v.1.72  
changed lines
  Added in v.1.75
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>