[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 547)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/07/06 01:10 PM

Modified:
 /trunk/src/
  GRHydro_CalcBcom.F90

Log:
 GRHydro: Don't compute Bcom twice.
 
 From: Christian Reisswig <reisswig at tapir.caltech.edu>

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_CalcBcom.F90
Delta lines: +5 -0
===================================================================
--- trunk/src/GRHydro_CalcBcom.F90	2013-07-06 18:10:32 UTC (rev 546)
+++ trunk/src/GRHydro_CalcBcom.F90	2013-07-06 18:10:37 UTC (rev 547)
@@ -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)