<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! </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="🙂">🙂</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-family:"Segoe UI", "Segoe UI Web (West European)", "Segoe UI", -apple-system, system-ui, Roboto, "Helvetica Neue", sans-serif;font-size:15px;background-color:rgb(255, 255, 255);display:inline !important">track_origin</span>"?
Also, <span style="background-color:rgb(255, 255, 255);display:inline !important">
although we expect only one merged black hole, </span>would it be necessary to set "<span style="font-family:"Segoe UI", "Segoe UI Web (West European)", "Segoe UI", -apple-system, system-ui, Roboto, "Helvetica Neue", sans-serif;font-size:14px;background-color:rgb(255, 255, 255);display:inline !important">n_horizons
= 2"<span style="font-family:Calibri, Arial, Helvetica, sans-serif;font-size:16px;background-color:rgb(255, 255, 255);display:inline !important"> ?</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 </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="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 <bruno.giacomazzo@unimib.it><br>
<b>Sent:</b> Wednesday, April 21, 2021 5:12 AM<br>
<b>To:</b> Erik Schnetter <schnetter@cct.lsu.edu><br>
<b>Cc:</b> Hamilton, Maria <babiuc@marshall.edu>; users@einsteintoolkit.org <users@einsteintoolkit.org><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 class="x_gmail_quote">
<div dir="ltr" class="x_gmail_attr">Il giorno mar 20 apr 2021 alle ore 19:30 Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu">schnetter@cct.lsu.edu</a>> 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 <<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" 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 <<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" 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&data=04%7C01%7Cbabiuc%40marshall.edu%7Cd983cd2c5ed34965850f08d904a5a966%7C239ab2783bba4c78b41d8508a541e025%7C0%7C0%7C637545931848216589%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=Z5GyYbUuBLu4AjnpBBKUMcQTXAKmaklJ3fDB%2Fk%2BQUf8%3D&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: </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" 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>