[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.81 and 1.82

version 1.81, 2010/04/16 07:13:42 version 1.82, 2010/05/01 02:17:49
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.80 2010/01/19 06:17:22 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/dp.c,v 1.81 2010/04/16 07:13:42 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "base.h"  #include "base.h"
Line 106  void Pdp_compute_last_w();
Line 106  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();  void Pdp_symb_add();
   void Pdp_mono_raddec();
   
 LIST dp_initial_term();  LIST dp_initial_term();
 LIST dp_order();  LIST dp_order();
Line 265  struct ftab dp_supp_tab[] = {
Line 266  struct ftab dp_supp_tab[] = {
         {"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_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},
           {"dp_mono_raddec",Pdp_mono_raddec,2},
   
         {0,0,0}          {0,0,0}
 };  };
Line 2629  VECT *rp;
Line 2631  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;
           }
   }
   
   NODE mono_raddec(NODE ideal);
   
   void Pdp_mono_raddec(NODE arg,LIST *rp)
   {
           NODE ideal,rd,t,t1,r,r1,u;
           VL vl0,vl;
           int nv,i,bpi;
           int *s;
           DP dp;
           P *v;
           LIST l;
   
           ideal = BDY((LIST)ARG0(arg));
           if ( !ideal ) *rp = (LIST)ARG0(arg);
           else {
                   t = BDY((LIST)ARG1(arg));
                   nv = length(t);
                   v = (P)MALLOC(nv*sizeof(P));
                   for ( vl0 = 0, i = 0; t; t = NEXT(t), i++ ) {
                           NEXTVL(vl0,vl); VR(vl) = VR((P)BDY(t));
                           MKV(VR(vl),v[i]);
                   }
                   if ( vl0 ) NEXT(vl) = 0;
                   for ( t = 0, r = ideal; r; r = NEXT(r) ) {
                           ptod(CO,vl0,BDY(r),&dp); MKNODE(t1,dp,t); t = t1;
                   }
                   rd = mono_raddec(t);
                   r = 0;
                   bpi = (sizeof(int)/sizeof(char))*8;
                   for ( u = rd; u; u = NEXT(u) ) {
                           s = (int *)BDY(u);
                           for ( i = nv-1, t = 0; i >= 0; i-- )
                                   if ( s[i/bpi]&(1<<(i%bpi)) ) {
                                           MKNODE(t1,v[i],t); t = t1;
                                   }
                           MKLIST(l,t); MKNODE(r1,l,r); r = r1;
                   }
                   MKLIST(*rp,r);
         }          }
 }  }
   
Legend:
Removed from v.1.81  
changed lines
  Added in v.1.82
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>