diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-07-06 18:11:12 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-07-06 18:11:12 +0000 |
commit | 669d5b37c19d16380ba6c63a7195cbcea8d08f35 (patch) | |
tree | 43c35e9b941be21d80ba34711d3bbfe7825b3b6d /src/GRHydro_Source.F90 | |
parent | fcbaf698714e9594ef35fbba18530a4a1ab164ba (diff) |
GRHydro: add grid function for sqrt(detg)
* add new 1 tl grid function sdetg that stores the sqrt of the determinent
of the 3-metric.
* replace lots of re-computation of det by use of this grid function
From: Christian Ott <cott@bethe.tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@555 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Source.F90')
-rw-r--r-- | src/GRHydro_Source.F90 | 18 |
1 files changed, 8 insertions, 10 deletions
diff --git a/src/GRHydro_Source.F90 b/src/GRHydro_Source.F90 index 553c916..1923c2d 100644 --- a/src/GRHydro_Source.F90 +++ b/src/GRHydro_Source.F90 @@ -71,7 +71,7 @@ subroutine SourceTerms(CCTK_ARGUMENTS) CCTK_INT :: i, j, k, na, nb, nc CCTK_REAL :: one, two, half CCTK_REAL :: t00, t0a, t0b, t0c, taa, tab, tac, tbb, tbc, tcc - CCTK_REAL :: sqrtdet, det, uaa, uab, uac, ubb, ubc, ucc, rhoenthalpyW2 + CCTK_REAL :: uaa, uab, uac, ubb, ubc, ucc, rhoenthalpyW2 CCTK_REAL :: shifta, shiftb, shiftc, velashift, velbshift, velcshift CCTK_REAL :: vlowa, vlowb, vlowc CCTK_REAL :: da_betaa, da_betab, da_betac, db_betaa, db_betab,& @@ -191,7 +191,7 @@ subroutine SourceTerms(CCTK_ARGUMENTS) !$OMP PARALLEL DO PRIVATE(i, j, k, local_spatial_order,& !$OMP localgaa,localgab,localgac,localgbb,localgbc,localgcc,& - !$OMP det,sqrtdet,rhoenthalpyW2,shifta,shiftb,shiftc,& + !$OMP rhoenthalpyW2,shifta,shiftb,shiftc,& !$OMP da_betaa,da_betab,da_betac,db_betaa,db_betab,db_betac,& !$OMP dc_betaa,dc_betab,dc_betac,velashift,velbshift,velcshift,& !$OMP vlowa,vlowb,vlowc,t00,t0a,t0b,t0c,taa,tab,tac,tbb,tbc,tcc,& @@ -221,10 +221,8 @@ subroutine SourceTerms(CCTK_ARGUMENTS) localgbc = g23(i,j,k) localgcc = g33(i,j,k) - det = SPATIAL_DETERMINANT(localgaa, localgab, localgac,\ - localgbb, localgbc, localgcc) - sqrtdet = sqrt(det) - call UpperMetric(uaa, uab, uac, ubb, ubc, ucc, det, localgaa,& + call UpperMetric(uaa, uab, uac, ubb, ubc, ucc, & + sdetg(i,j,k)*sdetg(i,j,k), localgaa,& localgab, localgac, localgbb, localgbc, localgcc) !!$ All the matter variables (except velocity) always appear @@ -442,10 +440,10 @@ subroutine SourceTerms(CCTK_ARGUMENTS) invalp densrhs(i,j,k) = 0.d0 - srhs(i,j,k,1) = alp(i,j,k)*sqrtdet*sa_source - srhs(i,j,k,2) = alp(i,j,k)*sqrtdet*sb_source - srhs(i,j,k,3) = alp(i,j,k)*sqrtdet*sc_source - taurhs(i,j,k) = alp(i,j,k)*sqrtdet*tau_source + srhs(i,j,k,1) = alp(i,j,k)*sdetg(i,j,k)*sa_source + srhs(i,j,k,2) = alp(i,j,k)*sdetg(i,j,k)*sb_source + srhs(i,j,k,3) = alp(i,j,k)*sdetg(i,j,k)*sc_source + taurhs(i,j,k) = alp(i,j,k)*sdetg(i,j,k)*tau_source enddo enddo |