For verifying frsqrt, we could use the bigfloat library (python bindings to mpfr). It supports emulating different ieee 754-2008 fp formats. We will have to remember that ieee 754-2008 defines sqrt(-0) to be -0 and rsqrt(-0) to be -infinity (not +infinity). https://pypi.org/project/bigfloat/