Fixed problem with glibc using an inaccurate implementation of

std::pow<complex> for computing powers of primitive root in ckks.cpp.
microsoft · Aug 14, 2019 · bdad6ab · bdad6ab
1 parent e1c63a8
commit bdad6ab
Showing 1 changed file with 2 additions and 3 deletions.
diff --git a/native/src/seal/ckks.cpp b/native/src/seal/ckks.cpp
@@ -62,11 +62,10 @@ namespace seal
 
         roots_ = allocate<complex<double>>(coeff_count, pool_);
         inv_roots_ = allocate<complex<double>>(coeff_count, pool_);
-        complex<double> psi{ cos((2 * PI_) / static_cast<double>(m)), 
-            sin((2 * PI_) / static_cast<double>(m)) };
+        double psi_arg = 2 * PI_ / static_cast<double>(m);
         for (size_t i = 0; i < coeff_count; i++)
         {
-            roots_[i] = pow(psi, static_cast<double>(reverse_bits(i, logn)));
+            roots_[i] = polar<double>(1.0, psi_arg * static_cast<double>(reverse_bits(i, logn)));
             inv_roots_[i] = 1.0 / roots_[i];
         }
     }