Maria
I see that you only set up one horizon. You need to set at least N_horizons to 2. All parameters that you set for horizon one (those with array index [1]), will also need to be set for the second horizons. It's also weird that the horizon is at x=0; in a binary system, I'd expect it to be away from the origin.
You probably want to set AHFinderDirect::which_surface_to_store_info[2] = 1 to store information about the second horizon in a different spherical surface.
-erik
On Tue, Apr 20, 2021 at 12:56 PM Hamilton, Maria babiuc@marshall.edu wrote:
Hello,
Could I please get some help in activating the horizon finder correctly in a bns run? This is what I have in the par file:
ActiveThorns = "AHFinderDirect" AHFinderDirect::find_after_individual_time[1] = 1500 AHFinderDirect::find_every = 64 AHFinderDirect::output_h_every = 128
AHFinderDirect::geometry_interpolator_name = "Lagrange polynomial interpolation" AHFinderDirect::geometry_interpolator_pars = "order=4" AHFinderDirect::surface_interpolator_name = "Lagrange polynomial interpolation" AHFinderDirect::surface_interpolator_pars = "order=4" AHFinderDirect::which_surface_to_store_info[1] = 0
AHFinderDirect::initial_guess_method [1] = "coordinate sphere" AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.3 AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0.0
AHFinderDirect::N_horizons = 1 AHFinderDirect::max_allowable_horizon_radius[1] = 3 AHFinderDirect::reset_horizon_after_not_finding[1] = "no"
AHFinderDirect::origin_x[1] = 0 AHFinderDirect::move_origins = "yes" AHFinderDirect::predict_origin_movement = "yes" AHFinderDirect::reshape_while_moving = "yes"
AHFinderDirect::track_origin_from_grid_scalar[1] = "yes" AHFinderDirect::track_origin_source_x[1] = "Hydro_Analysis::Hydro_Analysis_rho_max_loc[0]" AHFinderDirect::track_origin_source_y[1] = "Hydro_Analysis::Hydro_Analysis_rho_max_loc[1]" AHFinderDirect::track_origin_source_z[1] = "Hydro_Analysis::Hydro_Analysis_rho_max_loc[2]" AHFinderDirect::track_origin_from_grid_scalar[2] = "yes" AHFinderDirect::track_origin_source_x[2] = "Hydro_Analysis::Hydro_Analysis_rho_max_loc[0]" AHFinderDirect::track_origin_source_y[2] = "Hydro_Analysis::Hydro_Analysis_rho_max_loc[1]" AHFinderDirect::track_origin_source_z[2] = "Hydro_Analysis::Hydro_Analysis_rho_max_loc[2]"
And this is what happens when is activated:
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. Rank 0 with PID 22304 received signal 6 Writing backtrace to BNS_BSSN_Sly_myAHF_dx12/backtrace.0.txt
Primary job terminated normally, but 1 process returned a non-zero exit code. Per user-direction, the job has been aborted.
thanks, Maria _______________________ Maria C. Babiuc Hamilton, Ph.D. Professor, Department of Physics College of Science, Marshall University, 1 John Marshall Drive, Huntington, WV, 25755 Room S 257, Phone: (304)696-2754
Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users