From 1fdb87dd0e04ad5cf89fab3a650446cfbdb2aff0 Mon Sep 17 00:00:00 2001 From: diener Date: Fri, 1 Oct 2004 19:40:44 +0000 Subject: Fixed a bug in the RK2 re-initialization routines that were introduced when converting to Carpet. These routines are never used (which is why this bug wasn't discovered before) and may be removed at a later time. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/EHFinder/trunk@188 2a26948c-0e4f-0410-aee8-f1d3e353619c --- src/EHFinder_ReInitialize.F90 | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/EHFinder_ReInitialize.F90 b/src/EHFinder_ReInitialize.F90 index 5daf5fd..337c3a1 100644 --- a/src/EHFinder_ReInitialize.F90 +++ b/src/EHFinder_ReInitialize.F90 @@ -175,7 +175,7 @@ subroutine EHFinder_ReInitializeRK2_1(CCTK_ARGUMENTS) ! Set the step size from the courant factor and the minimum of the ! grid spacing. - h = hfac*minval(CCTK_DELTA_SPACE) + h = hfac*min(CCTK_DELTA_SPACE(1),CCTK_DELTA_SPACE(2),CCTK_DELTA_SPACE(3)) ! If centered finite differences are required... if ( CCTK_EQUALS ( pde_differences, 'centered' ) ) then @@ -375,7 +375,8 @@ subroutine EHFinder_ReInitializeRK2_2(CCTK_ARGUMENTS) ! If the maximum change on all processors are small enough signal that ! we are done with re-initialization by setting re_init_control = 0. - if ( all ( maxdf < h*minval(CCTK_DELTA_SPACE)**2 ) ) then + if ( all ( maxdf < h*min(CCTK_DELTA_SPACE(1),CCTK_DELTA_SPACE(2), & + CCTK_DELTA_SPACE(3))**2 ) ) then ! pugh_re_init_control = 0 re_init_control = 0 end if -- cgit v1.2.3