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

Spurious underflow in cyl_bessel_k #1229

Open
jzmaddock opened this issue Dec 23, 2024 · 0 comments
Open

Spurious underflow in cyl_bessel_k #1229

jzmaddock opened this issue Dec 23, 2024 · 0 comments

Comments

@jzmaddock
Copy link
Collaborator

cyl_bessel_k(100, 746) returns zero when evaluating at double precision (with no extended long double type), but the result should be 9.49147e-66.

Underflow occurs in k0 and k1 prior to iteration. The problem is also present if v is perturbed to be a non-integer in which case CF2_ik underflows.

Reduced test case from #1228.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant