Return to pgcd.k CVS log

Up to [local] / OpenXM / src / k097 / debug

Annotation of OpenXM/src/k097/debug/pgcd.k, Revision 1.1.1.1

1.1       maekawa     1:
                      2:
                      3: /* polynomial gcd.k */
                      4:
                      5: class Pgcd extends Object {
                      6:  pzero = Poly("0");
                      7:  def gcd0(f,g,xx) {
                      8:    local tmp,r;
                      9:    if (Degree(g,xx) > Degree(f,xx)) {
                     10:       tmp = f; f = g; g = tmp;
                     11:    }
                     12:    while (!IsConstant(g)) {
                     13:      r = Reduction(f,[g]);
                     14:      Println([f,g,r[1]]);
                     15:      r = r[0];
                     16:      f = g;
                     17:      g = r;
                     18:    }
                     19:    if (g == pzero) return(f);
                     20:    else return(g);
                     21:  }
                     22:  def gcd(f,g,xx) {
                     23:    local tmp,r;
                     24:    if (Degree(g,xx) > Degree(f,xx)) {
                     25:       tmp = f; f = g; g = tmp;
                     26:    }
                     27:    while (!IsConstant(g)) {
                     28:      r = pseudoRemainder(f,g,xx);
                     29:      f = g;
                     30:      g = r;
                     31:    }
                     32:    if (g == pzero) return(f);
                     33:    else return(g);
                     34:  }
                     35:  def pseudoRemainder(f,g,xx) {
                     36:    local c,r;
                     37:    c = Init(g);
                     38:    Println([c,xx]);
                     39:    c = Coefficients(c,xx);
                     40:    c = c[1,0];
                     41:    r = Reduction(c^(Degree(f,xx)-Degree(g,xx)+1)*f,[g]);
                     42:    Println([f,g,r[1]]);
                     43:    return(r[0]);
                     44:  }
                     45: }
                     46:
                     47:
                     48:
                     49:
                     50: