<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="gmail_quote"><div dir="ltr" class="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="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="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@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>><br>
<a href="http://www.perimeterinstitute.ca/personal/eschnetter/" 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="http://lists.einsteintoolkit.org/mailman/listinfo/users" 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="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>