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

Add features to KerrSpinSeq.__call__ (solve at new value of a) #5

Closed duetosymmetry closed 5 years ago

duetosymmetry commented 5 years ago

Interpolating was added by 3a9d14ff3a22498e57a6447785c26395e3217287 . Right now there is no checking that the value of a makes sense. Ideally we have a param that whether or not to extend the sequence up to the value of interest in the event that the requested a is greater than the highest a computed so far.

duetosymmetry commented 5 years ago

Add a parameter named something like interp_only which controls whether or not to simply interpolate. Right now, we interpolate for a starting guess, and then polish the root (this would correspond to interp_only=False. However we could just interpolate (interp_only=True). To get the decomposition coefficients when interpolating, we can just get the appropriate eigenvector of the angular matrix evaluated at the corresponding QNM frequency, and simply assume that it's close enough that we don't have to polish the root. Alternatively, maybe people don't even want the eigenvector, in which case we can skip this step.

duetosymmetry commented 5 years ago

Added options interp_only and resolve_if_found for more control in 322ad38d9ec343f66c3fe1a493ebf97653cf54cd .