[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.70 and 1.73

version 1.70, 2007/09/15 10:17:08 version 1.73, 2007/10/14 02:32:21
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.69 2007/09/07 00:45:50 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/dp.c,v 1.72 2007/09/19 05:42: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 90  void Pdp_weyl_gr_main(),Pdp_weyl_gr_mod_main(),Pdp_wey
Line 91  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_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();
Line 102  void Pdp_order();
Line 103  void Pdp_order();
 void Pdp_inv_or_split();  void Pdp_inv_or_split();
 void Pdp_compute_last_w();  void Pdp_compute_last_w();
 void Pdp_compute_essential_df();  void Pdp_compute_essential_df();
   void Pdp_get_denomlist();
   
 LIST dp_initial_term();  LIST dp_initial_term();
 LIST dp_order();  LIST dp_order();
Line 136  struct ftab dp_tab[] = {
Line 138  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 191  struct ftab dp_tab[] = {
Line 194  struct ftab dp_tab[] = {
         /* misc */          /* misc */
         {"dp_inv_or_split",Pdp_inv_or_split,3},          {"dp_inv_or_split",Pdp_inv_or_split,3},
         {"dp_set_weight",Pdp_set_weight,-1},          {"dp_set_weight",Pdp_set_weight,-1},
           {"dp_set_top_weight",Pdp_set_top_weight,-1},
         {"dp_weyl_set_weight",Pdp_weyl_set_weight,-1},          {"dp_weyl_set_weight",Pdp_weyl_set_weight,-1},
   
           {"dp_get_denomlist",Pdp_get_denomlist,0},
         {0,0,0},          {0,0,0},
 };  };
   
Line 971  LIST *rp;
Line 977  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 2362  VECT *rp;
Line 2400  VECT *rp;
                         current_dl_weight_vector[i] = QTOS((Q)v->body[i]);                          current_dl_weight_vector[i] = QTOS((Q)v->body[i]);
                 *rp = v;                  *rp = v;
         }          }
   }
   
   VECT current_top_weight_vector_obj;
   N *current_top_weight_vector;
   
   void Pdp_set_top_weight(arg,rp)
   NODE arg;
   VECT *rp;
   {
           VECT v;
           int i,n;
           NODE node;
   
           if ( !arg )
                   *rp = current_top_weight_vector_obj;
           else if ( !ARG0(arg) ) {
                   current_top_weight_vector = 0;
                   current_top_weight_vector_obj = 0;
                   *rp = 0;
           } else {
                   if ( OID(ARG0(arg)) != O_VECT && OID(ARG0(arg)) != O_LIST )
                           error("dp_set_top_weight : invalid argument");
                   if ( OID(ARG0(arg)) == O_VECT )
                           v = (VECT)ARG0(arg);
                   else {
                           node = (NODE)BDY((LIST)ARG0(arg));
                           n = length(node);
                           MKVECT(v,n);
                           for ( i = 0; i < n; i++, node = NEXT(node) )
                                   BDY(v)[i] = BDY(node);
                   }
                   for ( i = 0; i < v->len; i++ )
                           if ( !INT(BDY(v)[i]) || SGN((Q)BDY(v)[i]) < 0 )
                                   error("dp_set_top_weight : each element must be a non-negative integer");
                   current_top_weight_vector_obj = v;
                   current_top_weight_vector = (N *)MALLOC(v->len*sizeof(N));
                   for ( i = 0; i < v->len; i++ ) {
                           current_top_weight_vector[i] = NM((Q)BDY(v)[i]);
                   }
                   *rp = current_top_weight_vector_obj;
           }
   }
   
   LIST get_denomlist();
   
   void Pdp_get_denomlist(LIST *rp)
   {
           *rp = get_denomlist();
 }  }
   
 static VECT current_weyl_weight_vector_obj;  static VECT current_weyl_weight_vector_obj;

Legend:
Removed from v.1.70  
changed lines
  Added in v.1.73

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>