[BACK]Return to pgcd.k CVS log [TXT][DIR] Up to [local] / OpenXM / src / k097 / debug
File: [local] / OpenXM / src / k097 / debug / pgcd.k (download)

Revision 1.1.1.1 (vendor branch), Fri Oct 8 02:12:16 1999 UTC (24 years, 8 months ago) by maekawa
Branch: OpenXM, MAIN
CVS Tags: maekawa-ipv6, R_1_3_1-2, RELEASE_20000124, RELEASE_1_3_1_13b, RELEASE_1_2_3_12, RELEASE_1_2_3, RELEASE_1_2_2_KNOPPIX_b, RELEASE_1_2_2_KNOPPIX, RELEASE_1_2_2, RELEASE_1_2_1, RELEASE_1_1_3, RELEASE_1_1_2, KNOPPIX_2006, HEAD, DEB_REL_1_2_3-9, ALPHA
Changes since 1.1: +0 -0 lines 

o import OpenXM sources


/* polynomial gcd.k */

class Pgcd extends Object {
 pzero = Poly("0");
 def gcd0(f,g,xx) {
   local tmp,r;
   if (Degree(g,xx) > Degree(f,xx)) {
      tmp = f; f = g; g = tmp;
   }
   while (!IsConstant(g)) {
     r = Reduction(f,[g]);
     Println([f,g,r[1]]);
     r = r[0];
     f = g;
     g = r;
   }
   if (g == pzero) return(f);
   else return(g);
 }
 def gcd(f,g,xx) {
   local tmp,r;
   if (Degree(g,xx) > Degree(f,xx)) {
      tmp = f; f = g; g = tmp;
   }
   while (!IsConstant(g)) {
     r = pseudoRemainder(f,g,xx);
     f = g;
     g = r;
   }
   if (g == pzero) return(f);
   else return(g);
 }
 def pseudoRemainder(f,g,xx) {
   local c,r;
   c = Init(g);
   Println([c,xx]);
   c = Coefficients(c,xx);
   c = c[1,0];
   r = Reduction(c^(Degree(f,xx)-Degree(g,xx)+1)*f,[g]);
   Println([f,g,r[1]]);
   return(r[0]);
 }
}