search

alvincjk / EMRI_Kludge_Suite

25 stars 19 forks source link

Equatorial orbits #7

Open matthiasfabry opened 5 years ago

matthiasfabry commented 5 years ago

When I launch a run with the example parameters, only changing the inclination to be equatorial: iota = 0.0, The code hangs in a for loop that I traced to be in the source file: GKInsp.cc line 473 This for loop has no specified range and indeed runs forever when iota = 0 in the parameter file. Does this for loop need a specified range or what is the cause for this infinite loop?

alvincjk commented 5 years ago

This is a known issue and appears to be a problem only for the numerical kludge, which is currently not being maintained. As with the zero-eccentricity case, I would suggest using an epsilon value as an interim workaround.