jrenaud90 / TidalPy

Software suite to perform solid-body tidal dissipation calculations for rocky and icy worlds
Other
17 stars 3 forks source link

FLoating Point Issues in Takeuchi Starting Conditions #41

Open jrenaud90 opened 9 months ago

jrenaud90 commented 9 months ago

An earlier version of TidalPy used the full definition of Takeuchi & Saito (1972)'s starting conditions in the radial function. However, testing revealed that there were issues matching the limiting values (as z -> 0) since the full solution has terms that are like: (Very Big) / (Very Big) ~ near 1. Errors would compound leaving the final answer larger or smaller than 1.

Current Solution: Have reverted back to using the limiting version of these equations.