diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-07-06 18:10:37 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-07-06 18:10:37 +0000 |
| commit | 31f6bcaeb06c94ae34c307158d15de09f90e5215 (patch) | |
| tree | 5d77d2ef8368db53a57c2f8cc460eb5cc57cca68 | |
| parent | 25ec5f611d758944670af91340fd85f3fc858cd3 (diff) | |
GRHydro: Don't compute Bcom twice.
From: Christian Reisswig <reisswig@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@547 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
| -rw-r--r-- | src/GRHydro_CalcBcom.F90 | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/src/GRHydro_CalcBcom.F90 b/src/GRHydro_CalcBcom.F90 index fc46eea..234343d 100644 --- a/src/GRHydro_CalcBcom.F90 +++ b/src/GRHydro_CalcBcom.F90 @@ -40,6 +40,11 @@ call CCTK_WARN(0,"No storage for GRHydro::bcom") end if +! Only compute if Tmunu is not computed! (Bcom is also computed with Tmunu) + if (stress_energy_state .ne. 0) then + return + endif + !$OMP PARALLEL DO PRIVATE(i,j,k,velxlow, velylow, velzlow,& !$OMP Bvecxlow,Bvecylow,Bveczlow, bdotv,dum1,dum2,b2,bxlow,bylow,bzlow) |