Evolution of parameters

[Priti101]

Dear Joe,

If I set self.gr = True in Delaunay.py, I think 2.5 PN effect is turned on.
However, when I plot the evolution of a1 I do not get the smooth decay (as expected).
I tried it for the initial conditions set in your code.

Regards,
Priti

[joe-antognini]

Hi, @Priti101. By default a1 is 1 AU and the masses are 1 M_Sun, so the gravitational radiation will be exceedingly weak even over a Hubble time. Have you tried a value of a1 like 0.001? That should decay in perhaps ~10^5 yr.

[joshcrawford11]

Hi Joe,

I'm experiencing a similar issue as above. If I design a triple system where the inner binary will merge within a reasonable time frame (such as the conditions you mentioned above) and allow the system to evolve up to and including the merger, the evolve function consistently returns an error that it is attempting to access an index outside the array limits. Would you have any suggestions on how to fix this issue?

Kind regards,
Josh

[joe-antognini]

Unfortunately the TripleDelaunay code has a bug in it where the maxoutput parameter isn't getting used, so it will sometimes trigger an index error. I can try to fix it this weekend.

If you can make the test particle approximation I recommend using the TripleVectorial class instead since it will run much faster.