diff --git a/AHFinderDirect/interface.ccl b/AHFinderDirect/interface.ccl index 7891212f..8826d96e 100644 --- a/AHFinderDirect/interface.ccl +++ b/AHFinderDirect/interface.ccl @@ -148,3 +148,20 @@ FUNCTION sf_IdFromName (CCTK_INT IN id, CCTK_POINTER_TO_CONST IN name) USES FUNCTION sf_IdFromName ################################################################################ + +# CarpetX Interpolator +CCTK_INT FUNCTION DriverInterpolate( + CCTK_POINTER_TO_CONST IN cctkGH, + CCTK_INT IN N_dims, + CCTK_INT IN local_interp_handle, + CCTK_INT IN param_table_handle, + CCTK_INT IN coord_system_handle, + CCTK_INT IN N_interp_points, + CCTK_INT IN interp_coords_type, + CCTK_POINTER_TO_CONST ARRAY IN interp_coords, + CCTK_INT IN N_input_arrays, + CCTK_INT ARRAY IN input_array_indices, + CCTK_INT IN N_output_arrays, + CCTK_INT ARRAY IN output_array_types, + CCTK_POINTER ARRAY IN output_arrays) +REQUIRES FUNCTION DriverInterpolate diff --git a/AHFinderDirect/src/gr/expansion.cc b/AHFinderDirect/src/gr/expansion.cc index c3ca3c07..8f4b4121 100644 --- a/AHFinderDirect/src/gr/expansion.cc +++ b/AHFinderDirect/src/gr/expansion.cc @@ -826,7 +826,7 @@ sleep(1); CCTK_Barrier(cgi.GH); #endif -status = CCTK_InterpGridArrays(cgi.GH, N_GRID_DIMS, +status = DriverInterpolate(cgi.GH, N_GRID_DIMS, gi.operator_handle, gi.param_table_handle, cgi.coord_system_handle, N_interp_points, diff --git a/PunctureTracker/interface.ccl b/PunctureTracker/interface.ccl index abfeced8..dd58112e 100644 --- a/PunctureTracker/interface.ccl +++ b/PunctureTracker/interface.ccl @@ -17,3 +17,19 @@ CCTK_REAL pt_vel[10] TYPE=scalar { pt_vel_t pt_vel_x pt_vel_y pt_vel_z } "Velocity of punctures" + +CCTK_INT FUNCTION DriverInterpolate( + CCTK_POINTER_TO_CONST IN cctkGH, + CCTK_INT IN N_dims, + CCTK_INT IN local_interp_handle, + CCTK_INT IN param_table_handle, + CCTK_INT IN coord_system_handle, + CCTK_INT IN N_interp_points, + CCTK_INT IN interp_coords_type, + CCTK_POINTER_TO_CONST ARRAY IN interp_coords, + CCTK_INT IN N_input_arrays, + CCTK_INT ARRAY IN input_array_indices, + CCTK_INT IN N_output_arrays, + CCTK_INT ARRAY IN output_array_types, + CCTK_POINTER ARRAY IN output_arrays) +REQUIRES FUNCTION DriverInterpolate diff --git a/PunctureTracker/src/puncture_tracker.cc b/PunctureTracker/src/puncture_tracker.cc index 89110129..2ad9e9a5 100644 --- a/PunctureTracker/src/puncture_tracker.cc +++ b/PunctureTracker/src/puncture_tracker.cc @@ -106,7 +106,7 @@ extern "C" void PunctureTracker_Track(CCTK_ARGUMENTS) { // Interpolation operator const int operator_handle = - CCTK_InterpHandle("Lagrange polynomial interpolation"); + CCTK_InterpHandle("CarpetX"); if (operator_handle < 0) { CCTK_WARN(CCTK_WARN_ALERT, "Can't get interpolation handle"); return; @@ -169,7 +169,7 @@ extern "C" void PunctureTracker_Track(CCTK_ARGUMENTS) { // Interpolate int ierr; - ierr = CCTK_InterpGridArrays( + ierr = DriverInterpolate( cctkGH, dim, operator_handle, param_table_handle, coordsys_handle, num_points, CCTK_VARIABLE_REAL, interp_coords, num_vars, input_array_indices, num_vars, output_array_type_codes, output_arrays);