looking at mpn/generic/mul_fft.c its not entirely clear how to call a "mpn_mul_mod_2^k+1"(mp_ptr dest,mp_srcptr src1,mp_srcptr src2,mp_size_t k) // dest= src1*src2 % 2^(k*GMP_NUMB_BITS)+1 any hints ? thanks jason