diff options
Diffstat (limited to 'src/gr')
-rw-r--r-- | src/gr/driver.cc | 4 | ||||
-rw-r--r-- | src/gr/horizon_function.cc | 3 |
2 files changed, 7 insertions, 0 deletions
diff --git a/src/gr/driver.cc b/src/gr/driver.cc index 58eeaf6..f63bcd4 100644 --- a/src/gr/driver.cc +++ b/src/gr/driver.cc @@ -152,7 +152,11 @@ patch_system ps(origin_x, origin_y, origin_z, if (STRING_EQUAL(method, "horizon")) then { ps.read_ghosted_gridfn(ghosted_gfns::gfn__h, "h.dat", false); + const int timer_handle = CCTK_TimerCreateI(); + CCTK_TimerStartI(timer_handle); horizon_function(ps, cgi, gii); + CCTK_TimerStopI(timer_handle); + CCTK_TimerPrintDataI(timer_handle, -1); ps.print_gridfn(nominal_gfns::gfn__H, "H.dat"); } else CCTK_VWarn(-1, __LINE__, __FILE__, CCTK_THORNSTRING, diff --git a/src/gr/horizon_function.cc b/src/gr/horizon_function.cc index 04bb186..42e8f1d 100644 --- a/src/gr/horizon_function.cc +++ b/src/gr/horizon_function.cc @@ -95,6 +95,9 @@ void horizon_function(patch_system& ps, { CCTK_VInfo(CCTK_THORNSTRING, " horizon function"); +// fill in values of ghosted gridfns in ghost zones +ps.synchronize_ghost_zones(ghosted_gfns::gfn__h, ghosted_gfns::gfn__h); + // compute xyz positions of grid points and xyz derivative coeffs xyz_posns_and_deriv_coeffs(ps); |