In ``sqrt_mod`` fixed bug, added bound on root, added``all_roots`` argument;
in ``sqrt_mod_iter`` returned roots always in ZZ;
in ``_sqrt_mod_prime_power`` sorted result;
added tests