added default argument ``limit``  to ``sqrt_mod``