<div dir="ltr">Maria,<div> I do not track the centers of the stars in the AHFinderDirect, but in case of a strong kick in the post-merger remnant that could be useful since the BH may form far from the origin.</div><div><br></div><div>I always looked for two horizons with two different initial guesses in order to be sure to find at least 1 solution for the AH. It is not necessary though.</div><div><br></div><div>Cheers,</div><div>Bruno</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mer 21 apr 2021 alle ore 22:13 Hamilton, Maria <<a href="mailto:babiuc@marshall.edu">babiuc@marshall.edu</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Bruno,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Awesome! <span style="background-color:rgb(255,255,255);display:inline">
Thank you for sharing! </span>I'll give it a try! </div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span id="gmail-m_6443084474636642736🙂">🙂</span><br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Just a quick question: so, there is no need to "<span style="color:rgb(32,31,30);font-size:15px;background-color:rgb(255,255,255);display:inline">track_origin</span>"?
Also, <span style="background-color:rgb(255,255,255);display:inline">
although we expect only one merged black hole, </span>would it be necessary to set "<span style="font-size:14px;background-color:rgb(255,255,255);display:inline">n_horizons
= 2"<span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;background-color:rgb(255,255,255);display:inline"> ?</span></span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Maria</div>
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div id="gmail-m_6443084474636642736Signature">
<div>
<div id="gmail-m_6443084474636642736divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif">
<div name="divtagdefaultwrapper">
<div>
<div dir="ltr">_______________________</div>
<div dir="ltr"><i><font size="2">Maria C. Babiuc Hamilton, Ph.D.</font></i></div>
<div dir="ltr"><i><font size="2">Professor, Department of Physics </font></i></div>
<div dir="ltr"><i style="font-size:12pt"><font size="2"><i style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"><font size="2"></font></i>College of Science, </font></i><i style="font-size:12pt"><font size="2">Marshall University,</font></i></div>
<div dir="ltr"><i><font size="2">1 John Marshall Drive, </font></i><i style="font-size:12pt"><font size="2">Huntington, WV, 25755</font></i></div>
<div dir="ltr"><i><font size="2"><i style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"><font size="2"><i style="font-size:12pt"><font size="2">Room S 257, <span id="gmail-m_6443084474636642736ms-rterangepaste-start"></span><i style="font-size:16px"><font size="2">Phone:
(304)696-2754</font></i></font></i></font></i><br>
</font></i></div>
<div dir="ltr"><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div id="gmail-m_6443084474636642736appendonsend"></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_6443084474636642736divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Bruno Giacomazzo <<a href="mailto:bruno.giacomazzo@unimib.it" target="_blank">bruno.giacomazzo@unimib.it</a>><br>
<b>Sent:</b> Wednesday, April 21, 2021 5:12 AM<br>
<b>To:</b> Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>><br>
<b>Cc:</b> Hamilton, Maria <<a href="mailto:babiuc@marshall.edu" target="_blank">babiuc@marshall.edu</a>>; <a href="mailto:users@einsteintoolkit.org" target="_blank">users@einsteintoolkit.org</a> <<a href="mailto:users@einsteintoolkit.org" target="_blank">users@einsteintoolkit.org</a>><br>
<b>Subject:</b> Re: [Users] Assertion `ox' failed</font>
<div> </div>
</div>
<div>
<div dir="ltr">x=0 is OK since in a BNS system the BH will form at the origin where the NSs merge (unless there is a kick). These are the parameters I use in my BNS runs:
<div><br>
</div>
<div>#------------------ahfinderdirect----------------------------------------------<br>
<br>
ahfinderdirect::find_after_individual_time[1] = 2000 <br>
ahfinderdirect::find_after_individual_time[2] = 2000 <br>
ahfinderdirect::find_every = 64<br>
ahfinderdirect::geometry_interpolator_name = "Lagrange polynomial interpolation"<br>
ahfinderdirect::geometry_interpolator_pars = "order=4"<br>
ahfinderdirect::initial_guess__coord_sphere__radius[1] = 1<br>
ahfinderdirect::initial_guess__coord_sphere__radius[2] = 3<br>
ahfinderdirect::mask_buffer_thickness = 0<br>
ahfinderdirect::mask_radius_multiplier = 0.6<br>
ahfinderdirect::mask_radius_offset = 0<br>
ahfinderdirect::max_allowable_horizon_radius[1] = 15<br>
ahfinderdirect::max_allowable_horizon_radius[2] = 15<br>
ahfinderdirect::max_allowable_theta_growth_iterations = 10<br>
ahfinderdirect::max_allowable_theta_nonshrink_iterations = 10<br>
ahfinderdirect::max_newton_iterations__initial = 50<br>
ahfinderdirect::max_newton_iterations__subsequent = 50<br>
ahfinderdirect::move_origins = yes<br>
ahfinderdirect::n_horizons = 2<br>
<br>
ahfinderdirect::origin_x[1] = 0<br>
ahfinderdirect::origin_x[2] = 0<br>
ahfinderdirect::output_h_every = 64<br>
ahfinderdirect::predict_origin_movement = yes<br>
ahfinderdirect::reset_horizon_after_not_finding[1] = no<br>
ahfinderdirect::reset_horizon_after_not_finding[2] = no<br>
ahfinderdirect::reshape_while_moving = yes<br>
<br>
ahfinderdirect::surface_interpolator_name = "Lagrange polynomial interpolation"<br>
ahfinderdirect::surface_interpolator_pars = "order=4"<br>
ahfinderdirect::verbose_level = "physics details"<br>
ahfinderdirect::which_surface_to_store_info[1] = 0<br>
ahfinderdirect::which_surface_to_store_info[2] = 1<br>
</div>
<div><br>
</div>
</div>
<br>
<div>
<div dir="ltr">Il giorno mar 20 apr 2021 alle ore 19:30 Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>> ha scritto:<br>
</div>
<blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Maria<br>
<br>
I see that you only set up one horizon. You need to set at least<br>
N_horizons to 2. All parameters that you set for horizon one (those<br>
with array index [1]), will also need to be set for the second<br>
horizons. It's also weird that the horizon is at x=0; in a binary<br>
system, I'd expect it to be away from the origin.<br>
<br>
You probably want to set<br>
AHFinderDirect::which_surface_to_store_info[2] = 1<br>
to store information about the second horizon in a different spherical surface.<br>
<br>
<br>
-erik<br>
<br>
On Tue, Apr 20, 2021 at 12:56 PM Hamilton, Maria <<a href="mailto:babiuc@marshall.edu" target="_blank">babiuc@marshall.edu</a>> wrote:<br>
><br>
> Hello,<br>
><br>
> Could I please get some help in activating the horizon finder correctly in a bns run?<br>
> This is what I have in the par file:<br>
><br>
> ActiveThorns = "AHFinderDirect"<br>
> AHFinderDirect::find_after_individual_time[1] = 1500<br>
> AHFinderDirect::find_every = 64<br>
> AHFinderDirect::output_h_every = 128<br>
><br>
> AHFinderDirect::geometry_interpolator_name = "Lagrange polynomial interpolation"<br>
> AHFinderDirect::geometry_interpolator_pars = "order=4"<br>
> AHFinderDirect::surface_interpolator_name = "Lagrange polynomial interpolation"<br>
> AHFinderDirect::surface_interpolator_pars = "order=4"<br>
> AHFinderDirect::which_surface_to_store_info[1] = 0<br>
><br>
> AHFinderDirect::initial_guess_method [1] = "coordinate sphere"<br>
> AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.3<br>
> AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0.0<br>
><br>
> AHFinderDirect::N_horizons = 1<br>
> AHFinderDirect::max_allowable_horizon_radius[1] = 3<br>
> AHFinderDirect::reset_horizon_after_not_finding[1] = "no"<br>
><br>
> AHFinderDirect::origin_x[1] = 0<br>
> AHFinderDirect::move_origins = "yes"<br>
> AHFinderDirect::predict_origin_movement = "yes"<br>
> AHFinderDirect::reshape_while_moving = "yes"<br>
><br>
> AHFinderDirect::track_origin_from_grid_scalar[1] = "yes"<br>
> AHFinderDirect::track_origin_source_x[1] = "Hydro_Analysis::Hydro_Analysis_rho_max_loc[0]"<br>
> AHFinderDirect::track_origin_source_y[1] = "Hydro_Analysis::Hydro_Analysis_rho_max_loc[1]"<br>
> AHFinderDirect::track_origin_source_z[1] = "Hydro_Analysis::Hydro_Analysis_rho_max_loc[2]"<br>
> AHFinderDirect::track_origin_from_grid_scalar[2] = "yes"<br>
> AHFinderDirect::track_origin_source_x[2] = "Hydro_Analysis::Hydro_Analysis_rho_max_loc[0]"<br>
> AHFinderDirect::track_origin_source_y[2] = "Hydro_Analysis::Hydro_Analysis_rho_max_loc[1]"<br>
> AHFinderDirect::track_origin_source_z[2] = "Hydro_Analysis::Hydro_Analysis_rho_max_loc[2]"<br>
><br>
><br>
> And this is what happens when is activated:<br>
><br>
> cactus_bns: /home/babiuc/ET_2020_11/Cactus/configs/bns/build/AHFinderDirect/driver/Newton.cc:109: void AHFinderDirect::track_origin(const cGH*, AHFinderDirect::patch_system&, AHFinderDirect::AH_data*, int, bool): Assertion `ox' failed.<br>
> Rank 0 with PID 22304 received signal 6<br>
> Writing backtrace to BNS_BSSN_Sly_myAHF_dx12/backtrace.0.txt<br>
> -------------------------------------------------------<br>
> Primary job terminated normally, but 1 process returned<br>
> a non-zero exit code. Per user-direction, the job has been aborted.<br>
><br>
><br>
> thanks,<br>
> Maria<br>
> _______________________<br>
> Maria C. Babiuc Hamilton, Ph.D.<br>
> Professor, Department of Physics<br>
> College of Science, Marshall University,<br>
> 1 John Marshall Drive, Huntington, WV, 25755<br>
> Room S 257, Phone: (304)696-2754<br>
><br>
> _______________________________________________<br>
> Users mailing list<br>
> <a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
> <a href="https://nam02.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.einsteintoolkit.org%2Fmailman%2Flistinfo%2Fusers&data=04%7C01%7Cbabiuc%40marshall.edu%7Cd983cd2c5ed34965850f08d904a5a966%7C239ab2783bba4c78b41d8508a541e025%7C0%7C0%7C637545931848206633%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=7iVP0HVSRR8t5%2FiDS%2FRtXDYabxLYqzcBnSbnBor6GLg%3D&reserved=0" rel="noreferrer" target="_blank">
http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
<br>
<br>
<br>
-- <br>
Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>><br>
<a href="https://nam02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.perimeterinstitute.ca%2Fpersonal%2Feschnetter%2F&data=04%7C01%7Cbabiuc%40marshall.edu%7Cd983cd2c5ed34965850f08d904a5a966%7C239ab2783bba4c78b41d8508a541e025%7C0%7C0%7C637545931848206633%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=Z44Hbr5b%2BX%2BFZ8%2BWG1%2FunVbkyF27bBmqLMMrFR8K%2FkI%3D&reserved=0" rel="noreferrer" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a><br>
_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
<a href="https://nam02.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.einsteintoolkit.org%2Fmailman%2Flistinfo%2Fusers&data=04%7C01%7Cbabiuc%40marshall.edu%7Cd983cd2c5ed34965850f08d904a5a966%7C239ab2783bba4c78b41d8508a541e025%7C0%7C0%7C637545931848216589%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=Z5GyYbUuBLu4AjnpBBKUMcQTXAKmaklJ3fDB%2Fk%2BQUf8%3D&reserved=0" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
</blockquote>
</div>
<br clear="all">
<div><br>
</div>
-- <br>
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr">
<p><font color="#000000">Pr<font face="arial, sans-serif">of. Bruno Giacomazzo<br>
Department of Physics<br>
University of Milano-Bicocca<br>
</font></font><span style="color:rgb(0,0,0)"><font face="arial, sans-serif">Piazza della Scienza 3<br>
</font></span><span style="color:rgb(0,0,0)"><font face="arial, sans-serif">20126 Milano<br>
</font></span><span style="font-size:12.8px;color:rgb(0,0,0)">Italy</span></p>
<p><font color="#000000"><span style="font-size:12.8px">email: </span><span style="font-size:12.8px"><a href="mailto:bruno.giacomazzo@unimib.it" target="_blank">bruno.giacomazzo@unimib.it</a><br>
</span><span style="font-size:12.8px">phone: (+39) 02 6448 2321</span></font><br>
<font color="#000000"><span style="font-size:12.8px">web: </span></font><span style="color:rgb(0,0,0);font-size:12.8px"><a href="https://nam02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.brunogiacomazzo.org%2F&data=04%7C01%7Cbabiuc%40marshall.edu%7Cd983cd2c5ed34965850f08d904a5a966%7C239ab2783bba4c78b41d8508a541e025%7C0%7C0%7C637545931848216589%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=cG4VZXEfloZMWMsP9ByZUahKbwmS5eACsoCvxW8HI1w%3D&reserved=0" style="font-size:12.8px" target="_blank">http://www.brunogiacomazzo.org</a></span></p>
<p><font color="#000000">----------------------------------------------------------------------<br>
<span style="font-size:12.8px">There are only 10 types of people in the world:<br>
</span><span style="font-size:12.8px">Those who understand binary, and those who don't<br>
</span><span style="font-size:12.8px">----------------------------------------------------------------------<br>
</span></font></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><p><font color="#000000">Pr<font face="arial, sans-serif">of. Bruno Giacomazzo<br>Department of Physics<br>University of Milano-Bicocca<br></font></font><span style="color:rgb(0,0,0)"><font face="arial, sans-serif">Piazza della Scienza 3<br></font></span><span style="color:rgb(0,0,0)"><font face="arial, sans-serif">20126 Milano<br></font></span><span style="font-size:12.8px;color:rgb(0,0,0)">Italy</span></p><p><font color="#000000"><span style="font-size:12.8px">email: </span><span style="font-size:12.8px"><a href="mailto:bruno.giacomazzo@unimib.it" target="_blank">bruno.giacomazzo@unimib.it</a><br></span><span style="font-size:12.8px">phone: (+39) 02 6448 2321</span></font><br><font color="#000000"><span style="font-size:12.8px">web: </span></font><span style="color:rgb(0,0,0);font-size:12.8px"><a href="http://www.brunogiacomazzo.org/" style="font-size:12.8px" target="_blank">http://www.brunogiacomazzo.org</a></span></p><p><font color="#000000">----------------------------------------------------------------------<br><span style="font-size:12.8px">There are only 10 types of people in the world:<br></span><span style="font-size:12.8px">Those who understand binary, and those who don't<br></span><span style="font-size:12.8px">----------------------------------------------------------------------<br></span></font></p></div></div></div></div></div></div></div>