Dear all,
I am trying to create a simple parfile to help a PhD student. The goal is to set up puncture initial data with NRPyEllipticET and then use QuasiLocalMeasures/AHFinderDirect to compute properties of the punctures and compare with the expected results. No evolution thorn is expected/required and the simulation ends after the initial data has been set and the first round of analysis has finished.
The parfile attached partially works. The puncture initial data is set up and NRPyEllipticET's output is what I'd expect. However, the diagnostics are clearly wrong since: 1) AHFinderDirect fails to find the puncture (in this example I am setting the mass of one of the punctures to zero); and 2) the QuasiLocalMeasures diagnostics is never triggered, despite being scheduled.
Could someone provide some guidance on how to set up these thorns correctly, specifically SphericalSurfaces, QuasiLocalMeasures, and AHFinderDirect? Am I missing additional thorns?
Also, if someone has a working TwoPuncture example parfile instead, please share it with me as I'd be happy to adapt it!
Thanks in advance! Leo