From 31f6bcaeb06c94ae34c307158d15de09f90e5215 Mon Sep 17 00:00:00 2001 From: rhaas Date: Sat, 6 Jul 2013 18:10:37 +0000 Subject: GRHydro: Don't compute Bcom twice. From: Christian Reisswig git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@547 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45 --- src/GRHydro_CalcBcom.F90 | 5 +++++ 1 file changed, 5 insertions(+) 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) -- cgit v1.2.3