search

duetosymmetry / qnm

Python package for computing Kerr quasinormal mode frequencies, separation constants, and spherical-spheroidal mixing coefficients
https://qnm.readthedocs.io
MIT License
39 stars 19 forks source link

Automatically set cf_tol based on desired (omega) tolerance #22

Open duetosymmetry opened 5 years ago

duetosymmetry commented 5 years ago

The tolerance for computed continued fractions, cf_tol, should be related to the desired tolerance on omega, tol, by the gradient of the "objective function" (continued fraction equation). We can estimate it with a finite difference derivative.

If cf_tol ends up being smaller than DBL_EPSILON, we probably won't be able to polish omega to the desired tolerance.