<div dir="ltr">Leonardo<div><br></div><div>In this case the initial guess for the horizon radius might be off. As a rule of thumb, the surface where \alpha = 0.3 is the horizon. Your guess of r = 1 might be off.</div><div><br></div><div>-erik</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jun 28, 2023 at 4:58 PM Leonardo Werneck <<a href="mailto:wernecklr@gmail.com">wernecklr@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:arial,sans-serif;font-size:small">He Erik,</div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:small">Thanks for the response!</div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:small">The puncture is placed at x=5 because of the nature of the grid used by the elliptic solver. However, when the solution is interpolated to the carpet grid the position of the puncture is shifted to x=0 (this is controlled by the parameter NRPyEllipticET::position_shift[0]; please see the attached plot of the lapse function at t=0).</div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:small">Cheers,</div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:small">Leo</div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:small"><img src="cid:ii_ljg763e90" alt="image.png" width="562" height="391"><br></div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:small"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jun 28, 2023 at 1:49 PM Erik Schnetter <<a href="mailto:schnetter@gmail.com" target="_blank">schnetter@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">Leonardo<br>
<br>
I've briefly looked at the parameter file. The horizon finder needs to<br>
know an initial guess for both position and shape of the horizon. You<br>
are setting this location to 0, but are placing the puncture at x=5.<br>
The horizon finder then fails. QuasiLocalMeasures then sees that there<br>
is no horizon and does not run.<br>
<br>
-erik<br>
<br>
On Wed, Jun 28, 2023 at 3:50 PM Leonardo Werneck <<a href="mailto:wernecklr@gmail.com" target="_blank">wernecklr@gmail.com</a>> wrote:<br>
><br>
> Dear all,<br>
><br>
> 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.<br>
><br>
> 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.<br>
><br>
> Could someone provide some guidance on how to set up these thorns correctly, specifically SphericalSurfaces, QuasiLocalMeasures, and AHFinderDirect? Am I missing additional thorns?<br>
><br>
> Also, if someone has a working TwoPuncture example parfile instead, please share it with me as I'd be happy to adapt it!<br>
><br>
> Thanks in advance!<br>
> Leo<br>
><br>
> --<br>
> ------------------------------<br>
> Dr. Leonardo R. Werneck, Ph.D.<br>
> _______________________________________________<br>
> Users mailing list<br>
> <a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
> <a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
<br>
<br>
<br>
-- <br>
Erik Schnetter <<a href="mailto:schnetter@gmail.com" target="_blank">schnetter@gmail.com</a>><br>
<a href="http://www.perimeterinstitute.ca/personal/eschnetter/" rel="noreferrer" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a><br>
</blockquote></div>
</blockquote></div><br clear="all"><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Erik Schnetter <<a href="mailto:schnetter@gmail.com" target="_blank">schnetter@gmail.com</a>></div><div><a href="http://www.perimeterinstitute.ca/personal/eschnetter/" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a><br></div><div><br></div></div></div></div></div>