reactivate NewtonHensel step

parent 8a7a3edb
......@@ -175,7 +175,7 @@ Sum(const EffectiveDivisor& D1, const EffectiveDivisor& D2) {
ZZ_pX new_f = D1.get_f()*D2.get_f();
ZZ_pX new_g = (D1.get_g()*a2*(D2.get_f()/gcd) +
D2.get_g()*a1*(D1.get_f()/gcd)) % (new_f/gcd);
// new_g = NewtonHenselStep(*D1.curve().get_pdefpol(), new_g, new_f/gcd);
new_g = NewtonHenselStep(*D1.curve().get_pdefpol(), new_g, new_f/gcd);
new_g = new_g % new_f;
assert(new_g % D1.get_f() == D1.get_g());
assert(new_g % D2.get_f() == D2.get_g());
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment