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

Error message for KerrGeoIBSO at extended precisions. #55

Open
Kevin-Cunningham opened this issue Aug 30, 2023 · 0 comments
Open

Error message for KerrGeoIBSO at extended precisions. #55

Kevin-Cunningham opened this issue Aug 30, 2023 · 0 comments

Comments

@Kevin-Cunningham
Copy link

The line

orbit = KerrGeoOrbit[0.532, 1532, 0, .7`32]

Produces the error:
"FindRoot: The precision of the argument function ... is less than WorkingPrecision
(30.999999999999996`)."

where the ... is the full expression for IBSOPoly in SpecialOrbits.m

This does not happen for all pairs of a and x, but when it does occur the following small change should fix it in most cases.

The FindRoot is in KerrGeoIBSO where the WorkingPrecision is prec - 1 (e.g. 31 in the example above). Changing this to prec -1.1 or prec - 2 should solve the issue.

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