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

Diff for /OpenXM/src/asir-contrib/testing/noro/mpmat.rr between version 1.1 and 1.2

version 1.1, 2009/12/25 06:43:58 version 1.2, 2010/01/15 06:09:09
Line 1 
Line 1 
   def printsys(F,V,Name) {
           D = map(dp_ptod,F,V);
           output(Name);
           print(length(V),0); print(" ",0); print(length(F));
           map(printpoly,D);
           output();
   }
   
   def printpoly(F)
   {
           for ( I = 0, T = F; T; T = dp_rest(T), I++ );
           print(I);
           for ( T = F; T; T = dp_rest(T) ) {
                   print(dp_hc(T),0); print(" ",0); print(dp_etov(dp_ht(T)));
           }
   }
   
   def allmat(F,V,O,Name) {
           printsys(F,V,Name);
           G = nd_gr_trace(F,V,1,1,O);
           N = length(V);
           T = ttttt;
           for ( I = 0; I < N; I++ ) {
                   NameI = Name+rtostr(I+1);
                   UI = minipoly_mat(G,V,O,V[I],T);
                   utabout(UI,NameI);
           }
   }
   
 def utabout(U,Name)  def utabout(U,Name)
 {  {
         output(Name);          output(Name);
Line 40  def minipoly_mat(G,V,O,V1,V0)
Line 69  def minipoly_mat(G,V,O,V1,V0)
                         TL = cons(V1^J,TL);                          TL = cons(V1^J,TL);
                 NF = gennf(G,TL,V,O,V1,1)[0];                  NF = gennf(G,TL,V,O,V1,1)[0];
                 R = tolex_main(V,O,NF,[MP],M,MB);                  R = tolex_main(V,O,NF,[MP],M,MB);
                 return [Mat,subst(ptozp(R[0]),V1,V0)];                  return [Mat,sqpart(subst(ptozp(R[0]),V1,V0))];
         }          }
 }  }
   
Line 71  def utabtomat(G,V,V1,MB)
Line 100  def utabtomat(G,V,V1,MB)
                         }                          }
         }          }
         return [M,LCM];          return [M,LCM];
   }
   
   def sqpart(F)
   {
           G = gcd(F,diff(F,var(F)));
           return sdiv(F,G);
 }  }
 end$  end$
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>