@@ -152,10 +152,10 @@ inline void TwistIFFTUInt(PolynomialInFD<P> &res, const Polynomial<P> &a)
152
152
if constexpr (std::is_same_v<typename P::T, uint32_t >)
153
153
fftplvl1.execute_reverse_uint (res.data (), a.data ());
154
154
// if constexpr (std::is_same_v<typename P::T, uint64_t>)
155
- // fftplvl1.execute_reverse_torus64(res.data(), a.data());
155
+ // fftplvl1.execute_reverse_torus64(res.data(), a.data());
156
156
}
157
157
// else if constexpr (std::is_same_v<typename P::T, uint64_t>)
158
- // fftplvl2.execute_reverse_torus64(res.data(), a.data());
158
+ // fftplvl2.execute_reverse_torus64(res.data(), a.data());
159
159
else
160
160
static_assert (false_v<typename P::T>, " Undefined TwistIFFT!" );
161
161
}
@@ -316,8 +316,7 @@ inline void PolyMul(Polynomial<P> &res, const Polynomial<P> &a,
316
316
}
317
317
318
318
template <class P >
319
- inline void PolyMulRescaleUnsigned (Polynomial<P> &res,
320
- const Polynomial<P> &a,
319
+ inline void PolyMulRescaleUnsigned (Polynomial<P> &res, const Polynomial<P> &a,
321
320
const Polynomial<P> &b)
322
321
{
323
322
// if constexpr (std::is_same_v<typename P::T, uint32_t>) {
@@ -350,16 +349,20 @@ inline void PolyMulNaive(Polynomial<P> &res, const Polynomial<P> &a,
350
349
}
351
350
352
351
template <class P >
353
- inline void PolyMulNaieveRescaleUnsigned (Polynomial<P> &res, const Polynomial<P> &a,
354
- const Polynomial<P> &b)
352
+ inline void PolyMulNaieveRescaleUnsigned (Polynomial<P> &res,
353
+ const Polynomial<P> &a,
354
+ const Polynomial<P> &b)
355
355
{
356
356
for (int i = 0 ; i < P::n; i++) {
357
357
__int128_t ri = 0 ;
358
358
for (int j = 0 ; j <= i; j++)
359
- ri += static_cast <__int128_t >(a[j]) * static_cast <__int128_t >(b[i - j]);
359
+ ri += static_cast <__int128_t >(a[j]) *
360
+ static_cast <__int128_t >(b[i - j]);
360
361
for (int j = i + 1 ; j < P::n; j++)
361
- ri -= static_cast <__int128_t >(a[j]) * static_cast <__int128_t >(b[P::n + i - j]);
362
- // res[i] = static_cast<typename P::T>((ri) >> (std::numeric_limits<typename P::T>::digits - 3));
362
+ ri -= static_cast <__int128_t >(a[j]) *
363
+ static_cast <__int128_t >(b[P::n + i - j]);
364
+ // res[i] = static_cast<typename P::T>((ri) >>
365
+ // (std::numeric_limits<typename P::T>::digits - 3));
363
366
res[i] = static_cast <typename P::T>((ri) >> 29 );
364
367
}
365
368
}
0 commit comments