Skip to content

gh-117999: fixed small nonnegative integer powers of complex numbers #118000

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.

Sign up for GitHub

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

Draft
skirpichev wants to merge 13 commits into
base: main
Choose a base branch
from
Draft

gh-117999: fixed small nonnegative integer powers of complex numbers #118000

Changes from 1 commit
Commits
Show all changes
13 commits
Select commit Hold shift + click to select a range
70d21c7
gh-117999: fixed small nonnegative integer powers of complex numbers
skirpichev Apr 17, 2024
8e0c482
Polishing of c_powu(): exclude unreachible cases and magic numbers
skirpichev May 29, 2024
6fb1d12
Use _Py_c_pow() to test implementation
skirpichev May 30, 2024
87c420a
Merge branch 'main' into fix-c_powi-117999
skirpichev May 31, 2024
389441d
Apply suggestions from code review
skirpichev Aug 19, 2024
9f24ad1
Apply suggestions from code review
skirpichev Aug 19, 2024
807bcc6
Address review: subTest and asserts
skirpichev Aug 19, 2024
b3ec9c6
Merge branch 'main' into fix-c_powi-117999
skirpichev Sep 17, 2025
8cc2a2b
rename news
skirpichev Sep 17, 2025
9e9b454
take tests from #124243
skirpichev Sep 17, 2025
7f25e09
Use new _Py_rc_quot()
skirpichev Sep 17, 2025
99c7be1
run _Py_ADJUST_ERANGE2 only for finite input
skirpichev Sep 17, 2025
a86abc1
inline c_1
skirpichev Sep 17, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
run _Py_ADJUST_ERANGE2 only for finite input
  • Loading branch information
@skirpichev
skirpichev committed Sep 17, 2025
commit 99c7be1e10f3fa2d7bd032354e763427206aef39
14 changes: 11 additions & 3 deletions Objects/complexobject.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -333,7 +333,11 @@ _Py_c_pow(Py_complex a, Py_complex b)
r.real = len*cos(phase);
r.imag = len*sin(phase);

_Py_ADJUST_ERANGE2(r.real, r.imag);
if (isfinite(a.real) && isfinite(a.imag)
&& isfinite(b.real) && isfinite(b.imag))
{
_Py_ADJUST_ERANGE2(r.real, r.imag);
}
}
return r;
}
Expand Down Expand Up @@ -753,9 +757,13 @@ complex_pow(PyObject *v, PyObject *w, PyObject *z)
errno = 0;
// Check whether the exponent has a small integer value, and if so use
// a faster and more accurate algorithm.
if (b.imag == 0.0 && b.real == floor(b.real) && fabs(b.real) <= INT_EXP_CUTOFF) {
if (b.imag == 0.0 && b.real == floor(b.real)
&& fabs(b.real) <= INT_EXP_CUTOFF)
{
p = c_powi(a, (long)b.real);
_Py_ADJUST_ERANGE2(p.real, p.imag);
if (isfinite(a.real) && isfinite(a.imag)) {
_Py_ADJUST_ERANGE2(p.real, p.imag);
}
}
else {
p = _Py_c_pow(a, b);
Expand Down