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

FP32 log1pf_u10 and FP64 log1p_u10 are infinitely inaccurate for large inputs #473

Open
ebavier opened this issue Oct 25, 2023 · 0 comments
Labels

Comments

@ebavier
Copy link

ebavier commented Oct 25, 2023

On x86 targets and possibly others, the log1pf_u10 and log1p_u10 functions are infinitely inaccurate for large input values. Currently

  • log1pf_u10 (x > 1e38) = inf, and
  • log1p_u10 (x > 1e307) = inf

No input value except x = inf should result in an inf result for log1p.

Relevant source:

if (d > 1e+38) r = SLEEF_INFINITYf;

if (d > 1e+307) r = SLEEF_INFINITY;

@blapie blapie added the algo label Nov 2, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants