diff options
author | diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c> | 2004-10-01 19:40:44 +0000 |
---|---|---|
committer | diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c> | 2004-10-01 19:40:44 +0000 |
commit | 1fdb87dd0e04ad5cf89fab3a650446cfbdb2aff0 (patch) | |
tree | 9136262c281a6e1c467f82a106a0f972a9027fbe | |
parent | 51643738d6f24f963950615d2321dca5ebc86a28 (diff) |
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
-rw-r--r-- | src/EHFinder_ReInitialize.F90 | 5 |
1 files 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 |