aboutsummaryrefslogtreecommitdiff
path: root/src/AHFinder_gau.F
diff options
context:
space:
mode:
authordiener <diener@89daf98e-ef62-4674-b946-b8ff9de2216c>2001-08-31 08:46:09 +0000
committerdiener <diener@89daf98e-ef62-4674-b946-b8ff9de2216c>2001-08-31 08:46:09 +0000
commit47ca138d1f1d463b5ecaaa56c4dfbd8b99691e19 (patch)
tree52c12572f8aba527184e4d5649c55fd90a31c95e /src/AHFinder_gau.F
parent72b8133852a9b7b039e58966d1935b0b5c62bf0f (diff)
Fixed several bugs in the call to CCTK_ReduceLocalScalar where the drift correction is communicated
to all processors from processor 0. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@233 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src/AHFinder_gau.F')
-rw-r--r--src/AHFinder_gau.F11
1 files changed, 8 insertions, 3 deletions
diff --git a/src/AHFinder_gau.F b/src/AHFinder_gau.F
index e73d92e..3c1f75c 100644
--- a/src/AHFinder_gau.F
+++ b/src/AHFinder_gau.F
@@ -46,6 +46,7 @@
CCTK_REAL trxi,xi2
CCTK_REAL zero,half,one,two,three,four,pi
CCTK_REAL aux,sina,cosa
+ CCTK_REAL omega_corr_all
CCTK_REAL, dimension(3,3) :: ug,xi
CCTK_REAL, dimension(2,2) :: ga,ua
@@ -604,7 +605,7 @@
write(0,*) 'AHFinder_gau: xi_drift = ',xi_drift
- omega_corr = - xi_drift
+ omega_corr = - 1.1D0 * xi_drift
. /(dble(ahf_findevery)*cctk_delta_time)
omega_cum = omega_cum + omega_corr
@@ -649,8 +650,12 @@
if (drift_correct_on.eq.1) then
- call CCTK_ReduceLocalScalar(omega_corr,cctkGH,-1,sum_handle,
- . error1,rerror,CCTK_VARIABLE_INT)
+ call CCTK_ReduceLocalScalar(ierror,cctkGH,-1,sum_handle,
+ . omega_corr,omega_corr_all,CCTK_VARIABLE_REAL)
+ if (ierror.ne.0) then
+ call CCTK_WARN(1,"Reduction failed!")
+ end if
+ omega_corr = omega_corr_all
ahf_omega_corr = omega_corr
ahf_omega_cum = omega_cum