Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Remove remnants of support for non-IEEE 754 systems from cmathmodule.c #121268

Closed
skirpichev opened this issue Jul 2, 2024 · 3 comments
Closed
Labels
type-feature A feature request or enhancement

Comments

@skirpichev
Copy link
Member

skirpichev commented Jul 2, 2024

Feature or enhancement

Proposal:

Proposed patch:

diff --git a/Modules/cmathmodule.c b/Modules/cmathmodule.c
index bf86a211bc..591442334e 100644
--- a/Modules/cmathmodule.c
+++ b/Modules/cmathmodule.c
@@ -185,15 +185,8 @@ cmath_acos_impl(PyObject *module, Py_complex z)
     if (fabs(z.real) > CM_LARGE_DOUBLE || fabs(z.imag) > CM_LARGE_DOUBLE) {
         /* avoid unnecessary overflow for large arguments */
         r.real = atan2(fabs(z.imag), z.real);
-        /* split into cases to make sure that the branch cut has the
-           correct continuity on systems with unsigned zeros */
-        if (z.real < 0.) {
-            r.imag = -copysign(log(hypot(z.real/2., z.imag/2.)) +
-                               M_LN2*2., z.imag);
-        } else {
-            r.imag = copysign(log(hypot(z.real/2., z.imag/2.)) +
-                              M_LN2*2., -z.imag);
-        }
+        r.imag = -copysign(log(hypot(z.real/2., z.imag/2.)) +
+                           M_LN2*2., z.imag);
     } else {
         s1.real = 1.-z.real;
         s1.imag = -z.imag;
@@ -386,11 +379,7 @@ cmath_atanh_impl(PyObject *module, Py_complex z)
         */
         h = hypot(z.real/2., z.imag/2.);  /* safe from overflow */
         r.real = z.real/4./h/h;
-        /* the two negations in the next line cancel each other out
-           except when working with unsigned zeros: they're there to
-           ensure that the branch cut has the correct continuity on
-           systems that don't support signed zeros */
-        r.imag = -copysign(Py_MATH_PI/2., -z.imag);
+        r.imag = copysign(Py_MATH_PI/2., z.imag);
         errno = 0;
     } else if (z.real == 1. && ay < CM_SQRT_DBL_MIN) {
         /* C99 standard says:  atanh(1+/-0.) should be inf +/- 0i */

Removing this legacy was proposed in #102067, but before merging #31790. Maybe now that change does make sense?

Has this already been discussed elsewhere?

No response given

Links to previous discussion of this feature:

#102067 (comment)

Linked PRs

@skirpichev skirpichev added the type-feature A feature request or enhancement label Jul 2, 2024
skirpichev added a commit to skirpichev/cpython that referenced this issue Aug 6, 2024
as now building Python now requires support of IEEE 754 floating point
numbers.
@serhiy-storchaka
Copy link
Member

For reference, the original code was added in 6f34109 (ported to Python 3 in 53876d9).

vstinner pushed a commit that referenced this issue Aug 9, 2024
)

As now building Python now requires support of IEEE 754 floating point
numbers.
@vstinner
Copy link
Member

vstinner commented Aug 9, 2024

Fixed by change b6e745a.

@vstinner vstinner closed this as completed Aug 9, 2024
@serhiy-storchaka
Copy link
Member

I would like to hear opinions of @tiran and @mdickinson about this change.

blhsing pushed a commit to blhsing/cpython that referenced this issue Aug 22, 2024
…python#122716)

As now building Python now requires support of IEEE 754 floating point
numbers.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
type-feature A feature request or enhancement
Projects
None yet
Development

No branches or pull requests

3 participants