From 2b5883454c1b50afa76a80799c9cb4782b5f045f Mon Sep 17 00:00:00 2001 From: Leonardo de Moura Date: Sat, 12 Jan 2013 19:49:50 -0800 Subject: [PATCH] Add support for prem_gcd in square_free Signed-off-by: Leonardo de Moura --- src/math/realclosure/realclosure.cpp | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/src/math/realclosure/realclosure.cpp b/src/math/realclosure/realclosure.cpp index 37b438136..e4e174cc1 100644 --- a/src/math/realclosure/realclosure.cpp +++ b/src/math/realclosure/realclosure.cpp @@ -3475,12 +3475,17 @@ namespace realclosure { value_ref_buffer p_prime(*this); value_ref_buffer g(*this); derivative(sz, p, p_prime); - gcd(sz, p, p_prime.size(), p_prime.c_ptr(), g); + if (m_use_prem) + prem_gcd(sz, p, p_prime.size(), p_prime.c_ptr(), g); + else + gcd(sz, p, p_prime.size(), p_prime.c_ptr(), g); if (g.size() <= 1) { r.append(sz, p); } else { div(sz, p, g.size(), g.c_ptr(), r); + if (m_use_prem) + normalize_int_coeffs(r); } } }