BUG: special: Fix handling of powm1 overflow errors (#17855)

scipy · Jan 26, 2023 · e63cc01 · e63cc01
1 parent 0243a96
commit e63cc01
Showing 1 changed file with 6 additions and 1 deletion.
diff --git a/scipy/special/boost_special_functions.h b/scipy/special/boost_special_functions.h
@@ -90,7 +90,12 @@ Real powm1_wrap(Real x, Real y)
         z = NAN;
     } catch (const std::overflow_error& e) {
         sf_error("powm1", SF_ERROR_OVERFLOW, NULL);
-        z = INFINITY;
+        if (x > 0) {
+            z = INFINITY;
+        }
+        else {
+            z = -INFINITY;
+        }
     } catch (const std::underflow_error& e) {
         sf_error("powm1", SF_ERROR_UNDERFLOW, NULL);
         z = 0;