Merge pull request #1130 from boostorg/mp_fix

Correct ellint_1 logic for types with no infinity.
boostorg · May 4, 2024 · a53b013 · a53b013
2 parents 435d103 + 3dafc07
commit a53b013
Showing 1 changed file with 9 additions and 4 deletions.
diff --git a/include/boost/math/special_functions/ellint_1.hpp b/include/boost/math/special_functions/ellint_1.hpp
@@ -123,12 +123,17 @@ T ellint_f_imp(T phi, T k, const Policy& pol, T one_minus_k2)
           //
           T c = 1 / sinp;
           T c_minus_one = cosp / sinp;
-          T cross = fabs(c / (k * k));
           T arg2;
-          if ((cross > 0.9f) && (cross < 1.1f))
-             arg2 = c_minus_one + one_minus_k2;
+          if (k != 0)
+          {
+             T cross = fabs(c / (k * k));
+             if ((cross > 0.9f) && (cross < 1.1f))
+                arg2 = c_minus_one + one_minus_k2;
+             else
+                arg2 = c - k * k;
+          }
           else
-             arg2 = c - k * k;
+             arg2 = c;
           result = static_cast<T>(s * ellint_rf_imp(c_minus_one, arg2, c, pol));
        }
        else