<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body 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 !important">
Thank you for sharing!&nbsp;</span>I'll give it a try!&nbsp;</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span id="🙂">🙂</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 &quot;<span style="color:rgb(32, 31, 30);font-family:&quot;Segoe UI&quot;, &quot;Segoe UI Web (West European)&quot;, &quot;Segoe UI&quot;, -apple-system, system-ui, Roboto, &quot;Helvetica Neue&quot;, sans-serif;font-size:15px;background-color:rgb(255, 255, 255);display:inline !important">track_origin</span>&quot;?
 Also, <span style="background-color:rgb(255, 255, 255);display:inline !important">
although we expect only one merged black hole,&nbsp;</span>would it be necessary to set &quot;<span style="font-family:&quot;Segoe UI&quot;, &quot;Segoe UI Web (West European)&quot;, &quot;Segoe UI&quot;, -apple-system, system-ui, Roboto, &quot;Helvetica Neue&quot;, sans-serif;font-size:14px;background-color:rgb(255, 255, 255);display:inline !important">n_horizons
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 2&quot;<span style="font-family:Calibri, Arial, Helvetica, sans-serif;font-size:16px;background-color:rgb(255, 255, 255);display:inline !important">&nbsp;?</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="Signature">
<div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<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&nbsp;</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,&nbsp;</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,&nbsp;</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&nbsp;257,&nbsp;<span id="ms-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="appendonsend"></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Bruno Giacomazzo &lt;bruno.giacomazzo@unimib.it&gt;<br>
<b>Sent:</b> Wednesday, April 21, 2021 5:12 AM<br>
<b>To:</b> Erik Schnetter &lt;schnetter@cct.lsu.edu&gt;<br>
<b>Cc:</b> Hamilton, Maria &lt;babiuc@marshall.edu&gt;; users@einsteintoolkit.org &lt;users@einsteintoolkit.org&gt;<br>
<b>Subject:</b> Re: [Users] Assertion `ox' failed</font>
<div>&nbsp;</div>
</div>
<div>
<div dir="ltr">x=0 is OK since in a BNS system the BH will form at the origin&nbsp;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&nbsp;<br>
ahfinderdirect::find_after_individual_time[2] = 2000&nbsp;<br>
ahfinderdirect::find_every &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 64<br>
ahfinderdirect::geometry_interpolator_name = &quot;Lagrange polynomial interpolation&quot;<br>
ahfinderdirect::geometry_interpolator_pars = &quot;order=4&quot;<br>
ahfinderdirect::initial_guess__coord_sphere__radius[1] = 1<br>
ahfinderdirect::initial_guess__coord_sphere__radius[2] = 3<br>
ahfinderdirect::mask_buffer_thickness &nbsp; &nbsp;= 0<br>
ahfinderdirect::mask_radius_multiplier &nbsp; = 0.6<br>
ahfinderdirect::mask_radius_offset &nbsp; &nbsp; &nbsp; = 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 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = yes<br>
ahfinderdirect::n_horizons &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 2<br>
<br>
ahfinderdirect::origin_x[1] &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 0<br>
ahfinderdirect::origin_x[2] &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 0<br>
ahfinderdirect::output_h_every &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 64<br>
ahfinderdirect::predict_origin_movement &nbsp;= yes<br>
ahfinderdirect::reset_horizon_after_not_finding[1] = no<br>
ahfinderdirect::reset_horizon_after_not_finding[2] = no<br>
ahfinderdirect::reshape_while_moving &nbsp; &nbsp; = yes<br>
<br>
ahfinderdirect::surface_interpolator_name = &quot;Lagrange polynomial interpolation&quot;<br>
ahfinderdirect::surface_interpolator_pars = &quot;order=4&quot;<br>
ahfinderdirect::verbose_level &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &quot;physics details&quot;<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 class="x_gmail_quote">
<div dir="ltr" class="x_gmail_attr">Il giorno mar 20 apr 2021 alle ore 19:30 Erik Schnetter &lt;<a href="mailto:schnetter@cct.lsu.edu">schnetter@cct.lsu.edu</a>&gt; ha scritto:<br>
</div>
<blockquote class="x_gmail_quote" 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 &lt;<a href="mailto:babiuc@marshall.edu" target="_blank">babiuc@marshall.edu</a>&gt; wrote:<br>
&gt;<br>
&gt; Hello,<br>
&gt;<br>
&gt; Could I please get some help in activating the horizon finder correctly in a bns run?<br>
&gt; This is what I have in the par file:<br>
&gt;<br>
&gt; ActiveThorns = &quot;AHFinderDirect&quot;<br>
&gt; AHFinderDirect::find_after_individual_time[1] = 1500<br>
&gt; AHFinderDirect::find_every = 64<br>
&gt; AHFinderDirect::output_h_every = 128<br>
&gt;<br>
&gt; AHFinderDirect::geometry_interpolator_name = &quot;Lagrange polynomial interpolation&quot;<br>
&gt; AHFinderDirect::geometry_interpolator_pars = &quot;order=4&quot;<br>
&gt; AHFinderDirect::surface_interpolator_name = &quot;Lagrange polynomial interpolation&quot;<br>
&gt; AHFinderDirect::surface_interpolator_pars = &quot;order=4&quot;<br>
&gt; AHFinderDirect::which_surface_to_store_info[1] = 0<br>
&gt;<br>
&gt; AHFinderDirect::initial_guess_method [1] = &quot;coordinate sphere&quot;<br>
&gt; AHFinderDirect::initial_guess__coord_sphere__radius[1]&nbsp; &nbsp;= 0.3<br>
&gt; AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0.0<br>
&gt;<br>
&gt; AHFinderDirect::N_horizons&nbsp; = 1<br>
&gt; AHFinderDirect::max_allowable_horizon_radius[1] = 3<br>
&gt; AHFinderDirect::reset_horizon_after_not_finding[1] = &quot;no&quot;<br>
&gt;<br>
&gt; AHFinderDirect::origin_x[1] = 0<br>
&gt; AHFinderDirect::move_origins = &quot;yes&quot;<br>
&gt; AHFinderDirect::predict_origin_movement = &quot;yes&quot;<br>
&gt; AHFinderDirect::reshape_while_moving = &quot;yes&quot;<br>
&gt;<br>
&gt; AHFinderDirect::track_origin_from_grid_scalar[1] = &quot;yes&quot;<br>
&gt; AHFinderDirect::track_origin_source_x[1] = &quot;Hydro_Analysis::Hydro_Analysis_rho_max_loc[0]&quot;<br>
&gt; AHFinderDirect::track_origin_source_y[1] = &quot;Hydro_Analysis::Hydro_Analysis_rho_max_loc[1]&quot;<br>
&gt; AHFinderDirect::track_origin_source_z[1] = &quot;Hydro_Analysis::Hydro_Analysis_rho_max_loc[2]&quot;<br>
&gt; AHFinderDirect::track_origin_from_grid_scalar[2] = &quot;yes&quot;<br>
&gt; AHFinderDirect::track_origin_source_x[2] = &quot;Hydro_Analysis::Hydro_Analysis_rho_max_loc[0]&quot;<br>
&gt; AHFinderDirect::track_origin_source_y[2] = &quot;Hydro_Analysis::Hydro_Analysis_rho_max_loc[1]&quot;<br>
&gt; AHFinderDirect::track_origin_source_z[2] = &quot;Hydro_Analysis::Hydro_Analysis_rho_max_loc[2]&quot;<br>
&gt;<br>
&gt;<br>
&gt; And this is what happens when is activated:<br>
&gt;<br>
&gt; 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&amp;, AHFinderDirect::AH_data*, int, bool): Assertion `ox' failed.<br>
&gt; Rank 0 with PID 22304 received signal 6<br>
&gt; Writing backtrace to BNS_BSSN_Sly_myAHF_dx12/backtrace.0.txt<br>
&gt; -------------------------------------------------------<br>
&gt; Primary job&nbsp; terminated normally, but 1 process returned<br>
&gt; a non-zero exit code. Per user-direction, the job has been aborted.<br>
&gt;<br>
&gt;<br>
&gt; thanks,<br>
&gt; Maria<br>
&gt; _______________________<br>
&gt; Maria C. Babiuc Hamilton, Ph.D.<br>
&gt; Professor, Department of Physics<br>
&gt; College of Science, Marshall University,<br>
&gt; 1 John Marshall Drive, Huntington, WV, 25755<br>
&gt; Room S 257, Phone: (304)696-2754<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; Users mailing list<br>
&gt; <a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
&gt; <a href="https://nam02.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.einsteintoolkit.org%2Fmailman%2Flistinfo%2Fusers&amp;data=04%7C01%7Cbabiuc%40marshall.edu%7Cd983cd2c5ed34965850f08d904a5a966%7C239ab2783bba4c78b41d8508a541e025%7C0%7C0%7C637545931848206633%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&amp;sdata=7iVP0HVSRR8t5%2FiDS%2FRtXDYabxLYqzcBnSbnBor6GLg%3D&amp;reserved=0" originalsrc="http://lists.einsteintoolkit.org/mailman/listinfo/users" shash="YxD9v1NLBJI2YmfUyqWIi3WVLgaqJ2l1Zv+rqUApdT3DHvg8ydK0+mJUlKeZ80oq+RtgblDeM7+za7EkoiByEatPcusA764ET7oQ6C1ZWjwPdjqRwpxxs1X9R781hohwNAqQ8wj2xDGQda/xiOCY+Krs0+vEoaOjLSPgfVaSg5w=" rel="noreferrer" target="_blank">
http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
<br>
<br>
<br>
-- <br>
Erik Schnetter &lt;<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>&gt;<br>
<a href="https://nam02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.perimeterinstitute.ca%2Fpersonal%2Feschnetter%2F&amp;data=04%7C01%7Cbabiuc%40marshall.edu%7Cd983cd2c5ed34965850f08d904a5a966%7C239ab2783bba4c78b41d8508a541e025%7C0%7C0%7C637545931848206633%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&amp;sdata=Z44Hbr5b%2BX%2BFZ8%2BWG1%2FunVbkyF27bBmqLMMrFR8K%2FkI%3D&amp;reserved=0" originalsrc="http://www.perimeterinstitute.ca/personal/eschnetter/" shash="bTtTkIhn+E6GMvUs+MqqeVLkM9jjwhfg0baQ5uOCgh+t2eSlevNw8cUyMMhbZEuXqPqM6v1CHUpBQIEHw3mug3qL2wYmh6ZAiiNZszb2gIBSf4tg+scb96fcbBMLk9vZyiwT3yBSZ38Q9SVX8YcOmiZd6ocjApTjikvv2yEdR6o=" 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&amp;data=04%7C01%7Cbabiuc%40marshall.edu%7Cd983cd2c5ed34965850f08d904a5a966%7C239ab2783bba4c78b41d8508a541e025%7C0%7C0%7C637545931848216589%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&amp;sdata=Z5GyYbUuBLu4AjnpBBKUMcQTXAKmaklJ3fDB%2Fk%2BQUf8%3D&amp;reserved=0" originalsrc="http://lists.einsteintoolkit.org/mailman/listinfo/users" shash="tmyToORbdPvkzlTKEc5RTUwfyZxlIFhyuzewStzhHs1r7oWBFZqmqnOMlM/c/7KvuR30HCKKHSSSZ5cwfxVLFfW724b79UxPzIIpLRmEjDoiIHnvMHgFflazoCL2DvMCFWlpCl+PSnpmvMx+Jo/R48TO1ItXLUueoiNHRUJy9hU=" 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" class="x_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:&nbsp;</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:&nbsp;</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&amp;data=04%7C01%7Cbabiuc%40marshall.edu%7Cd983cd2c5ed34965850f08d904a5a966%7C239ab2783bba4c78b41d8508a541e025%7C0%7C0%7C637545931848216589%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&amp;sdata=cG4VZXEfloZMWMsP9ByZUahKbwmS5eACsoCvxW8HI1w%3D&amp;reserved=0" originalsrc="http://www.brunogiacomazzo.org/" shash="XdnuihBCOwCnrRVtOhrHCgDqiq9CCRgec5WEGJuYu7iEE0vO0t19mKoEIxjeOLgxlgn8VjgG6/BPvzfYjq9ewfakWuL1FT6rizmLEtEotv1/3AeyZMYbgvqrgUawOWNCZW+sLdbH9zekXuBVJ2JemEFINGNeDvtnEm2I25yndi4=" target="_blank" style="font-size:12.8px">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>
</body>
</html>