added function congruence_groups_between_gamma0_and_gamma1