We should also have a TFM implementation of the `sqrtmod_prime()` function as done for LTM in https://github.com/libtom/libtommath/issues/33