[BACK]Return to ndbf.rr CVS log [TXT][DIR] Up to [local] / OpenXM / src / asir-contrib / testing / noro

Diff for /OpenXM/src/asir-contrib/testing/noro/ndbf.rr between version 1.9 and 1.10

version 1.9, 2009/11/12 01:39:54 version 1.10, 2010/04/28 05:58:43
Line 12 
Line 12 
   
 if (!module_definedp("gr")) load("gr")$ else{ }$  if (!module_definedp("gr")) load("gr")$ else{ }$
 if (!module_definedp("primdec")) load("primdec")$ else{ }$  if (!module_definedp("primdec")) load("primdec")$ else{ }$
   if (!module_definedp("newsyz")) load("noro_module_syz.rr")$ else{ }$
   /* Empty for now. It will be used in a future. */    /* Empty for now. It will be used in a future. */
   
 /* toplevel */  /* toplevel */
Line 21  module ndbf$
Line 22  module ndbf$
 /* bfunction */  /* bfunction */
   
 localf bfunction, in_ww, in_ww_main, ann, ann_n$  localf bfunction, in_ww, in_ww_main, ann, ann_n$
 localf ann0, psi, ww_weight, compare_first, generic_bfct$  localf ann0, ann_fa, psi, ww_weight, compare_first, generic_bfct$
 localf generic_bfct_1, initial_part, bfct, indicial1, bfct_via_gbfct$  localf generic_bfct_1, initial_part, bfct, indicial1, bfct_via_gbfct$
 localf bfct_via_gbfct_weight, bfct_via_gbfct_weight_1, bfct_via_gbfct_weight_2$  localf bfct_via_gbfct_weight, bfct_via_gbfct_weight_1, bfct_via_gbfct_weight_2$
 localf weyl_minipolym, weyl_minipoly, weyl_nf, weyl_nf_quo_check$  localf weyl_minipolym, weyl_minipoly, weyl_nf, weyl_nf_quo_check$
Line 469  def ann0(F)
Line 470  def ann0(F)
         return [Min,map(subst,Ann,s,Min,TMP_S,s,TMP_DS,ds)];          return [Min,map(subst,Ann,s,Min,TMP_S,s,TMP_DS,ds)];
 }  }
   
   /*
    * For a polynomial F and a scalar A,
    * compute generators of Ann(F^A).
    */
   
   def ann_fa(F,A)
   {
           if ( type(Syz=getopt(syz)) == -1 ) Syz = 0;
   
           F = subst(F,s,TMP_S);
           Ann = ann(F);
           Bf = bfunction(F);
   
           FList = cdr(fctr(Bf));
           for ( T = FList, Min = 0; T != []; T = cdr(T) ) {
                   LF = car(car(T));
                   Root = -coef(LF,0)/coef(LF,1);
                   if ( dn(Root) == 1 && Root < Min )
                           Min = Root;
           }
           D = A-Min;
           if ( dn(D) != 1 || D <= 0 )
                   return map(ptozp,map(subst,Ann,s,A,TMP_S,s,TMP_DS,ds));
   
           V = vars(F);
           for ( I = length(V)-1, DV = []; I >= 0; I-- )
                   DV = cons(strtov("d"+rtostr(V[I])),DV);
           VDV = append(V,DV);
           R = map(subst,Ann,s,Min,TMP_S,s,TMP_DS,ds);
           F = ptozp(F);
   
           if ( Syz ) {
                   /* syzygy method */
                   S = newsyz.module_syz(cons(F^D,R),VDV,1,0|weyl=1);
                   B = car(S);
                   for ( R = []; B != []; B = cdr(B) )
                           if ( H = car(car(B)) )
                                   R = cons(H,R);
           } else {
                   /* colon method */
                   for ( I = 0; I < D; I++ )
                           R = weyl_ideal_quotient(R,F,VDV);
           }
           return R;
   }
   
 def psi0(F,T,DT)  def psi0(F,T,DT)
 {  {
         D = dp_ptod(F,[T,DT]);          D = dp_ptod(F,[T,DT]);
Line 1323  def weyl_ideal_quotient(B,F,VDV)
Line 1370  def weyl_ideal_quotient(B,F,VDV)
         O1 = [[0,1],[0,N+1]];          O1 = [[0,1],[0,N+1]];
         GJ = nd_weyl_gr(J,VDV1,0,O1);          GJ = nd_weyl_gr(J,VDV1,0,O1);
         R = elimination(GJ,VDV);          R = elimination(GJ,VDV);
         return R;          return map(weyl_divide_by_right,R,F,VDV,0);
 }  }
   
 def bf_strat(F)  def bf_strat(F)
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>