void PowerfulDCRT::dcrtToPowerful(Vec<ZZ>& out, const DoubleCRT& dcrt) const
{
const IndexSet& set = dcrt.getIndexSet();
if (empty(set)) { // sanity check
clear(out);
return;
}
zz_pBak bak; bak.save(); // backup NTL's current modulus
ZZ product = conv<ZZ>(1L);
for (long i = set.first(); i <= set.last(); i = set.next(i)) {
pConvVec[i].restoreModulus();
zz_pX oneRowPoly;
long newPrime = dcrt.getOneRow(oneRowPoly,i);
HyperCube<zz_p> oneRowPwrfl(indexes.shortSig);
pConvVec[i].polyToPowerful(oneRowPwrfl, oneRowPoly);
if (i == set.first()) // just copy
conv(out, oneRowPwrfl.getData());
else // CRT
intVecCRT(out, product, oneRowPwrfl.getData(), newPrime); // in NumbTh
product *= newPrime;
}
}
