Skip to content

Errors for special orbits with numeric orbital parameters #60

@znasipak

Description

@znasipak

Using x = 1.0 rather than x=1 typically leads to errors. For example,

orbit = KerrGeoOrbit[0.9`40, 15.0`40, 0.0`40, 1.0`40]

produces the error:

FindRoot: The precision of the argument function \
(0``79.76399993315745+0``38.63109226687787 KerrGeodesics`\
SpecialOrbits`Private`p^2+<<78>> (-<<78>>-<<21>> <<37>>) <<37>>^\
5+(-4+<<37>>)^2 <<37>>^6+0.\
6560999999999999999999999999999999999999999999999999999998`39.\
39794000867204 KerrGeodesics`SpecialOrbits`Private`p^3 (0``37.\
95078197732982+1.`38.19382002601611 KerrGeodesics`SpecialOrbits`\
Private`p)) is less than WorkingPrecision (39.`)

while

orbit = KerrGeoOrbit[0.9`40, 15.0`40, 0.0`40, 1]

evaluates without error.

Metadata

Metadata

Assignees

No one assigned

    Labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions