allow p == 1 in dgs_bern_mp_init() as well
Add -n option to bench_gauss to pick number of trials
Make sigma and c static in dgs_disc_gauss_mp_init