Thanks to visit codestin.com
Credit goes to github.com

Skip to content

gh-101678: refactor the math module to use special functions from c11 #101679

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

Merged
merged 13 commits into from
Feb 9, 2023
Merged

gh-101678: refactor the math module to use special functions from c11 #101679

Changes from 1 commit
Commits
Show all changes
13 commits
Select commit Hold shift + click to select a range
d202b5c
gh-101678: refactor the math module to use special functions from c11
skirpichev Feb 8, 2023
7d14105
Use libc's lgamma/tgamma instead of custom implementations
skirpichev Feb 8, 2023
468096b
Revert "Use libc's lgamma/tgamma instead of custom implementations"
skirpichev Feb 8, 2023
5b75ef8
Comment on custom implementations of tgamma/lgamma
skirpichev Feb 8, 2023
538810d
UPD commit msg: gh-101678: refactor the math module to use erf/erfc/l…
skirpichev Feb 8, 2023
fa76552
Grammar fix
skirpichev Feb 9, 2023
42829cb
Restore braces in m_log2 per PEP7
skirpichev Feb 9, 2023
27334d8
Merge branch 'main' into math-c11
skirpichev Feb 9, 2023
c097545
Merge branch 'main' into math-c11
skirpichev Feb 9, 2023
47970e2
Merge math_1_to_whatever() and math_1()
skirpichev Feb 9, 2023
bebbc79
Provide context for log1p workaround
skirpichev Feb 9, 2023
ee5e689
Revert "Merge math_1_to_whatever() and math_1()"
skirpichev Feb 9, 2023
06e9c8c
Fix issues prefixes
skirpichev Feb 9, 2023
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
UPD commit msg: gh-101678: refactor the math module to use erf/erfc/l… 
…og2 from c11
  • Loading branch information
@skirpichev
skirpichev committed Feb 8, 2023
commit 538810dd92b95e8973b0c621a957741fa7d73d2f
21 changes: 1 addition & 20 deletions Modules/mathmodule.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -641,27 +641,8 @@ m_log2(double x)
}
}

if (x > 0.0) {
#ifdef HAVE_LOG2
if (x > 0.0)
return log2(x);
#else
double m;
int e;
m = frexp(x, &e);
/* We want log2(m * 2**e) == log(m) / log(2) + e. Care is needed when
* x is just greater than 1.0: in that case e is 1, log(m) is negative,
* and we get significant cancellation error from the addition of
* log(m) / log(2) to e. The slight rewrite of the expression below
* avoids this problem.
*/
if (x >= 1.0) {
return log(2.0 * m) / log(2.0) + (e - 1);
}
else {
return log(m) / log(2.0) + e;
}
#endif
}
else if (x == 0.0) {
errno = EDOM;
return -Py_HUGE_VAL; /* log2(0) = -inf, divide-by-zero */
Expand Down