Return +-Py_HUGE_VAL for tgamma(+-0) instead of risking FP exceptions by computing...

diff --git a/Modules/mathmodule.c b/Modules/mathmodule.c
index cebb4ff91ba58e6dc25b5196e28643809fe3c389..7e73bfed62b5a4c8de9ee9beca8e11baf938499b 100644 (file)
--- a/Modules/mathmodule.c
+++ b/Modules/mathmodule.c
@@ -239,7 +239,8 @@ m_tgamma(double x)
     }
     if (x == 0.0) {
         errno = EDOM;
-        return 1.0/x; /* tgamma(+-0.0) = +-inf, divide-by-zero */
+        /* tgamma(+-0.0) = +-inf, divide-by-zero */
+        return copysign(Py_HUGE_VAL, x);
     }
 
     /* integer arguments */