poulpy/spqlios/lib/test/testlib/negacyclic_polynomial.cpp
spqlios/lib/test/testlib/negacyclic_polynomial.cppRaw
1#include "negacyclic_polynomial_impl.h"
2
3// explicit instantiation
4EXPLICIT_INSTANTIATE_POLYNOMIAL(__int128_t);
5EXPLICIT_INSTANTIATE_POLYNOMIAL(int64_t);
6EXPLICIT_INSTANTIATE_POLYNOMIAL(double);
7
8double infty_dist(const rnx_f64& a, const rnx_f64& b) {
9 const uint64_t nn = a.nn();
10 const double* aa = a.data();
11 const double* bb = b.data();
12 double res = 0.;
13 for (uint64_t i = 0; i < nn; ++i) {
14 double d = fabs(aa[i] - bb[i]);
15 if (d > res) res = d;
16 }
17 return res;
18}