Issue #22604: Fix assertion error in debug mode when dividing a complex number by (nan+0j).
diff --git a/Objects/complexobject.c b/Objects/complexobject.c
index 5ee0c15..2659a23 100644
--- a/Objects/complexobject.c
+++ b/Objects/complexobject.c
@@ -96,7 +96,7 @@
const double abs_breal = b.real < 0 ? -b.real : b.real;
const double abs_bimag = b.imag < 0 ? -b.imag : b.imag;
- if (abs_breal >= abs_bimag) {
+ if (abs_breal >= abs_bimag) {
/* divide tops and bottom by b.real */
if (abs_breal == 0.0) {
errno = EDOM;
@@ -109,7 +109,7 @@
r.imag = (a.imag - a.real * ratio) / denom;
}
}
- else {
+ else if (abs_bimag >= abs_breal) {
/* divide tops and bottom by b.imag */
const double ratio = b.real / b.imag;
const double denom = b.real * ratio + b.imag;
@@ -117,6 +117,10 @@
r.real = (a.real * ratio + a.imag) / denom;
r.imag = (a.imag * ratio - a.real) / denom;
}
+ else {
+ /* At least one of b.real or b.imag is a NaN */
+ r.real = r.imag = Py_NAN;
+ }
return r;
}