Use minors based hnf in pernet stein if full column rank and dimension is less than a cutoff
use the right dimension to bound the cutoff between dixon and normal solving
several improvements to rref_multi_mod, changed to use the new _nmod_mat_rref too