add X9.31 prime generation routines from 0.9.8 branch