From a0ca9827e820145f17a053081c4f3ece06de87d1 Mon Sep 17 00:00:00 2001 From: dos-reis Date: Thu, 15 May 2008 09:40:44 +0000 Subject: * algebra/nregset.spad.pamphlet (recip$NormalizationPackage): Fix thinko. --- src/ChangeLog | 5 +++++ src/algebra/nregset.spad.pamphlet | 12 ++++++++---- 2 files changed, 13 insertions(+), 4 deletions(-) (limited to 'src') diff --git a/src/ChangeLog b/src/ChangeLog index f46e4df2..b9ff33a9 100644 --- a/src/ChangeLog +++ b/src/ChangeLog @@ -1,3 +1,8 @@ +2008-05-14 Gabriel Dos Reis + + * algebra/nregset.spad.pamphlet (recip$NormalizationPackage): Fix + thinko. + 2008-05-13 Gabriel Dos Reis * interp/c-util.boot (warnLiteral): Just warn, really. diff --git a/src/algebra/nregset.spad.pamphlet b/src/algebra/nregset.spad.pamphlet index a1ad85fe..9405dfbd 100644 --- a/src/algebra/nregset.spad.pamphlet +++ b/src/algebra/nregset.spad.pamphlet @@ -160,13 +160,17 @@ NormalizationPackage(R,E,V,P,TS): Exports == Implementation where -- ASSUME mvar(p) is algebraic w.r.t. ts v := mvar(p) ts_v := select(ts,v)::P + d : P + n : P if mdeg(p) < mdeg(ts_v) then - hesrg: Record (gcd : P, coef2 : P) := halfExtendedSubResultantGcd2(ts_v,p)$P - d: P := hesrg.gcd; n: P := hesrg.coef2 + hesrg2 := halfExtendedSubResultantGcd2(ts_v,p)$P + d := hesrg2.gcd + n := hesrg2.coef2 else - hesrg: Record (gcd : P, coef1 : P) := halfExtendedSubResultantGcd1(p,ts_v)$P - d: P := hesrg.gcd; n: P := hesrg.coef1 + hesrg1 := halfExtendedSubResultantGcd1(p,ts_v)$P + d := hesrg1.gcd + n := hesrg1.coef1 g := gcd(n,d) (n, d) := ((n exquo g)::P, (d exquo g)::P) remn, remd: Record(rnum:R,polnum:P,den:R) -- cgit v1.2.3