// Set b := 2*b-a*b^2 mod 2^(intDsize*n2)
cl_UDS_mul_square(dest_LSDptr,n,b2_LSDptr); // b^2
cl_UDS_mul(b2_LSDptr,n2,a_LSDptr,n2,prod_LSDptr); // a*b^2
// Set b := 2*b-a*b^2 mod 2^(intDsize*n2)
cl_UDS_mul_square(dest_LSDptr,n,b2_LSDptr); // b^2
cl_UDS_mul(b2_LSDptr,n2,a_LSDptr,n2,prod_LSDptr); // a*b^2