From 0c6450f69b1c9a0e99d0e34968342d05b8a99882 Mon Sep 17 00:00:00 2001 From: David Wells Date: Sat, 29 Aug 2015 13:33:54 -0400 Subject: [PATCH] Compute powers of -1 by checking exponent parity. --- source/base/polynomials_bernstein.cc | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/source/base/polynomials_bernstein.cc b/source/base/polynomials_bernstein.cc index c9c3e8849b..942d3e9314 100644 --- a/source/base/polynomials_bernstein.cc +++ b/source/base/polynomials_bernstein.cc @@ -18,7 +18,7 @@ namespace Binomial std::vector coeff(n + 1, number(0.0)); for (unsigned int i = k; i < n + 1; ++i) { - coeff[i] = pow(number(-1), number(i - k)) + coeff[i] = ((i - k) % 2 == 0 ? 1 : -1) * boost::math::binomial_coefficient(n, i) * boost::math::binomial_coefficient(i, k); } -- 2.39.5