search

alvincjk / EMRI_Kludge_Suite

25 stars 19 forks source link

Circular orbits #6

Open matthiasfabry opened 5 years ago

matthiasfabry commented 5 years ago

When starting a run with initial eccentricity 0 (ie a circular orbit) the code quickly spits out: Zero eccentricity? Use the circular code, punky.

I can trace this message to come from src file: IEKG.cc line 29

Indeed the reason that is given there is sound, but there are plenty of circular tools built in right? Is there a switch to make the code use these circular tools?

alvincjk commented 5 years ago

As stated in the readme, this is a known issue. There is indeed some legacy code containing circular tools, but I have not gotten round to interfacing them with the code that is being maintained.

In the meantime, I would suggest using an epsilon value (say 10^-3) as a workaround. If something more precise is required, the kludges are not quite suitable for the job to begin with. Although there is currently no publicly available alternative, some other projects in the Black Hole Perturbation Toolkit might eventually include functionality for waveform generation (e.g. https://github.com/BlackHolePerturbationToolkit/GremlinEq).