aboutsummaryrefslogtreecommitdiff
path: root/src/AHFinder.F
diff options
context:
space:
mode:
authordiener <diener@89daf98e-ef62-4674-b946-b8ff9de2216c>2002-04-30 12:50:27 +0000
committerdiener <diener@89daf98e-ef62-4674-b946-b8ff9de2216c>2002-04-30 12:50:27 +0000
commita5f32fbd5bdd6794a8eab5cd19b5a4ef953ab58f (patch)
treed8adfadeaa5bf4a38432ad18ffb901ebdb7cafc3 /src/AHFinder.F
parent28e8940b0905cc35c44fc40b93681e974e2897df (diff)
AHFinder converted to the changes in Einstein. Passes the testsuites on the
xeon. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@295 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src/AHFinder.F')
-rw-r--r--src/AHFinder.F15
1 files changed, 9 insertions, 6 deletions
diff --git a/src/AHFinder.F b/src/AHFinder.F
index 9ead22e..d01cd04 100644
--- a/src/AHFinder.F
+++ b/src/AHFinder.F
@@ -54,6 +54,7 @@
CCTK_REAL, parameter :: one = 1.0D0
CCTK_REAL, dimension (1:3) :: ahfmass
+ CCTK_INT :: save_conformal_state
character*200 logf
@@ -190,12 +191,13 @@
! I need to remember the original value of the flag to change
! it back at the end if necesary.
- if (conformal_state == CONFORMAL_METRIC) then
-
+ if ( CCTK_Equals(metric_type,"static conformal") ) then
ahfconformal = .true.
- conformal_state = NOCONFORMAL_METRIC
+ save_conformal_state = conformal_state
+ conformal_state = 0
- call ConfToPhys(nx,ny,nz,psi,gxx,gxy,gxz,gyy,gyz,gzz)
+ call StaticConformal_ConfToPhysInPlace (nx,ny,nz,psi,
+ & gxx,gxy,gxz,gyy,gyz,gzz)
end if
@@ -1604,9 +1606,10 @@
if (ahfconformal) then
- conformal_state = CONFORMAL_METRIC
+ conformal_state = save_conformal_state
- call PhysToConf(nx,ny,nz,psi,gxx,gxy,gxz,gyy,gyz,gzz)
+ call StaticConformal_PhysToConfInPlace (nx,ny,nz,psi,
+ & gxx,gxy,gxz,gyy,gyz,gzz)
end if