[BACK]Return to de CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / lib

Diff for /OpenXM_contrib2/asir2000/lib/de between version 1.1 and 1.5

version 1.1, 2005/08/02 07:16:42 version 1.5, 2005/08/09 05:34:10
Line 3  module de;
Line 3  module de;
 localf split_lexgb;  localf split_lexgb;
 localf sort_lex_dec,sort_lex_inc;  localf sort_lex_dec,sort_lex_inc;
 localf inverse_or_split, linear_dim;  localf inverse_or_split, linear_dim;
 localf sp_sqrt,calcb,dp_monic_mod,monic_gb;  localf dp_monic_mod,monic_gb;
   localf membership_test;
 localf dp_chrem,intdptoratdp,intdpltoratdpl;  localf dp_chrem,intdptoratdp,intdpltoratdpl;
 localf comp_by_ht,dp_gr_mod,gr_chrem;  localf comp_by_ht,dp_gr_mod,gr_chrem;
   localf construct_sqfrbasis;
   
 /*  /*
  * G : a 0-dim lex gb, reduced   * G : a 0-dim lex gb, reduced
Line 52  def inverse_or_split(V,Id,F)
Line 54  def inverse_or_split(V,Id,F)
         Ret = inv_or_split_dalg(DF);          Ret = inv_or_split_dalg(DF);
         if ( type(Ret) == 1 ) {          if ( type(Ret) == 1 ) {
                 /* Ret = 1/F */                  /* Ret = 1/F */
                 return Ret;                  Dp = dalgtodp(Ret);
                   return dp_dtop(Dp[0],V)/Dp[1];
         } else {          } else {
                 /* Ret = GB(Id:F) */                  /* Ret = GB(Id:F) */
                 /* compute GB(Id+<f>) */                  /* compute GB(Id+<f>) */
                 Gquo = append(map(ptozp,map(dp_dtop,Ret,V)),Id);                  Gquo = append(map(ptozp,map(dp_dtop,Ret,V)),Id);
                 Gquo = nd_gr(Gquo,V,0,2);                  /* inter-reduction */
                   Gquo = nd_gr_postproc(Gquo,V,0,2,0);
                 DTotal = linear_dim(Id,V,2);                  DTotal = linear_dim(Id,V,2);
                 Dquo = linear_dim(Gquo,V,2);                  Dquo = linear_dim(Gquo,V,2);
                 Drem = DTotal-Dquo;                  Drem = DTotal-Dquo;
Line 67  def inverse_or_split(V,Id,F)
Line 71  def inverse_or_split(V,Id,F)
         }          }
 }  }
   
   /* add F(X,V) to Id(B) */
   /* returns a list of splitted ideals */
   /* B should be a triangular basis */
   
   def construct_sqfrbasis(F,X,B,V)
   {
           if ( type(F) == 1 )
                   return [];
           B = sort_lex_dec(B,V);
           V1 = cons(X,V);
           F = nd_nf(F,reverse(B),cons(X,V),2,0);
           D = deg(F,X);
           H = coef(F,D,X);
           if ( type(H) == 2 ) {
                   Ret = inverse_or_split(V,B,H);
                   if ( type(Ret) == 4 ) {
                           /* H != 0 on Id_nz, H = 0 on Id_z */
                           B0=construct_sqfrbasis(F,X,Ret[0],V);
                           B1=construct_sqfrbasis(F,X,Ret[1],V);
                           return append(B0,B1);
                   } else
                           F = nd_nf(F*Ret,reverse(B),cons(X,V),2,0);
           }
           B1 = cons(F,B);
           /* F is monic */
           M = minipoly(B1,V1,2,X,zzz);
           S = sqfr(M); S = cdr(S);
           if ( length(S) == 1 && car(S)[1] == 1 )
                   return [cons(F,B)];
           else {
                   R = [];
                   for ( T = S; T != []; T = cdr(T) ) {
                           G = nd_gr_trace(cons(subst(car(T)[0],zzz,X),B1),V1,1,1,2);
                           R1 = split_lexgb(G,V1);
                           R = append(R1,R);
                   }
                   return R;
           }
   }
   
 def sort_lex_dec(B,V)  def sort_lex_dec(B,V)
 {  {
         dp_ord(2);          dp_ord(2);
Line 92  def linear_dim(G,V,Ord)
Line 136  def linear_dim(G,V,Ord)
         return length(MB);          return length(MB);
 }  }
   
   def membership_test(B,G,V,O)
   {
           B = map(ptozp,B);
           G = map(ptozp,G);
           for ( T = B; T != []; T = cdr(T) )
                   if ( nd_nf(car(T),G,V,O,0) ) return 0;
           return 1;
   }
   
 def gr_chrem(B,V,O,Dim)  def gr_chrem(B,V,O,Dim)
 {  {
         B = map(ptozp,B);          B = map(ptozp,B);
Line 105  def gr_chrem(B,V,O,Dim)
Line 158  def gr_chrem(B,V,O,Dim)
                 Mod *= P;                  Mod *= P;
                 if ( G != [] )                  if ( G != [] )
                         HS = HSM;                          HS = HSM;
                 R1 = intdpltoratdpl(G,Mod);                  M = idiv(isqrt(2*Mod),2);
                   R1 = intdpltoratdpl(G,Mod,M);
                 if ( R1 ) {                  if ( R1 ) {
                         if ( Found && R == R1 )                          if ( Found && R == R1
                                   && (GB=nd_gr_postproc(map(dp_dtop,R,V),V,0,O,1))
                                   && membership_test(B,GB,V,O) )
                                 break;                                  break;
                         else {                          else {
                                 R = R1; Found = 1;                                  R = R1; Found = 1;
                         }                          }
                 }                  }
         }          }
         return map(dp_dtop,R,V);          return GB;
 }  }
   
 def comp_by_ht(A,B)  def comp_by_ht(A,B)
Line 128  def comp_by_ht(A,B)
Line 184  def comp_by_ht(A,B)
                 return 0;                  return 0;
 }  }
   
 def intdpltoratdpl(G,Mod)  def intdpltoratdpl(G,Mod,M)
 {  {
         R = [];          for ( R = []; G != []; G = cdr(G) ) {
         M = calcb(Mod);  
         for ( ; G != []; G = cdr(G) ) {  
                 T = intdptoratdp(car(G),Mod,M);                  T = intdptoratdp(car(G),Mod,M);
                 if ( !T )                  if ( !T )
                         return 0;                          return 0;
Line 167  def dp_chrem(G,HS,Mod,GM,HSM,P)
Line 221  def dp_chrem(G,HS,Mod,GM,HSM,P)
                 return [];                  return [];
         R = [];          R = [];
         M1 = inv(Mod,P);          M1 = inv(Mod,P);
           ModM1 = Mod*M1;
         for ( ; G != []; G = cdr(G), GM = cdr(GM) ) {          for ( ; G != []; G = cdr(G), GM = cdr(GM) ) {
                 E = car(G); EM = car(GM);                  E = car(G); EM = car(GM);
                 E1 = E+(EM-E)*Mod*M1;                  E1 = E+(EM-E)*ModM1;
                 R = cons(E1,R);                  R = cons(E1,R);
         }          }
         return reverse(R);          return reverse(R);
Line 190  def dp_monic_mod(F,P)
Line 245  def dp_monic_mod(F,P)
         return dp_rat(FP/dp_hc(FP));          return dp_rat(FP/dp_hc(FP));
 }  }
   
 def calcb(M) {  
         N = 2*M;  
         T = sp_sqrt(N);  
         if ( T^2 <= N && N < (T+1)^2 )  
                 return idiv(T,2);  
         else  
                 error("afo");  
 }  
   
 def sp_sqrt(A) {  
         for ( J = 0, T = A; T >= 2^27; J++ ) {  
                 T = idiv(T,2^27)+1;  
         }  
         for ( I = 0; T >= 2; I++ ) {  
                 S = idiv(T,2);  
                 if ( T = S+S )  
                         T = S;  
                 else  
                         T = S+1;  
         }  
         X = (2^27)^idiv(J,2)*2^idiv(I,2);  
         while ( 1 ) {  
                 if ( (Y=X^2) < A )  
                         X += X;  
                 else if ( Y == A )  
                         return X;  
                 else  
                         break;  
         }  
         while ( 1 )  
                 if ( (Y = X^2) <= A )  
                         return X;  
                 else  
                         X = idiv(A + Y,2*X);  
 }  
   
 endmodule;  endmodule;
 end$  
 end$  end$
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.5
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>