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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2014/04/15 02:49 PM

Modified:
 /trunk/src/
  GRHydro_CalcUpdate.F90

Log:
 GRHydro: rewrite MHD rhs in symmetric way

File Changes:

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

File [modified]: GRHydro_CalcUpdate.F90
Delta lines: +12 -12
===================================================================
--- trunk/src/GRHydro_CalcUpdate.F90	2014-04-15 19:49:34 UTC (rev 627)
+++ trunk/src/GRHydro_CalcUpdate.F90	2014-04-15 19:49:37 UTC (rev 628)
@@ -303,18 +303,18 @@
     do k = GRHydro_stencil + 1, cctk_lsh(3) - GRHydro_stencil
       do j = GRHydro_stencil + 1, cctk_lsh(2) - GRHydro_stencil
         do i = GRHydro_stencil + 1, cctk_lsh(1) - GRHydro_stencil
-          Bconsrhs(i,j,k,1) = - 0.5d0 * (Evec(i-1,j  ,k-1,2)-Evec(i-1,j  ,k  ,2)) / CCTK_DELTA_SPACE(3) &
-                              - 0.5d0 * (Evec(i-1,j  ,k  ,3)-Evec(i-1,j-1,k  ,3)) / CCTK_DELTA_SPACE(2) &
-                              - 0.5d0 * (Evec(i  ,j  ,k-1,2)-Evec(i  ,j  ,k  ,2)) / CCTK_DELTA_SPACE(3) &
-                              - 0.5d0 * (Evec(i  ,j  ,k  ,3)-Evec(i  ,j-1,k  ,3)) / CCTK_DELTA_SPACE(2)
-          Bconsrhs(i,j,k,2) = - 0.5d0 * (Evec(i-1,j-1,k  ,3)-Evec(i  ,j-1,k  ,3)) / CCTK_DELTA_SPACE(1) &
-                              - 0.5d0 * (Evec(i  ,j-1,k  ,1)-Evec(i  ,j-1,k-1,1)) / CCTK_DELTA_SPACE(3) &
-                              - 0.5d0 * (Evec(i-1,j  ,k  ,3)-Evec(i  ,j  ,k  ,3)) / CCTK_DELTA_SPACE(1) &
-                              - 0.5d0 * (Evec(i  ,j  ,k  ,1)-Evec(i  ,j  ,k-1,1)) / CCTK_DELTA_SPACE(3)
-          Bconsrhs(i,j,k,3) = - 0.5d0 * (Evec(i  ,j-1,k-1,1)-Evec(i  ,j  ,k-1,1)) / CCTK_DELTA_SPACE(2) &
-                              - 0.5d0 * (Evec(i  ,j  ,k-1,2)-Evec(i-1,j  ,k-1,2)) / CCTK_DELTA_SPACE(1) &
-                              - 0.5d0 * (Evec(i  ,j-1,k  ,1)-Evec(i  ,j  ,k  ,1)) / CCTK_DELTA_SPACE(2) &
-                              - 0.5d0 * (Evec(i  ,j  ,k  ,2)-Evec(i-1,j  ,k  ,2)) / CCTK_DELTA_SPACE(1)
+          Bconsrhs(i,j,k,1) = - 0.5d0 * ((Evec(i-1,j  ,k-1,2)-Evec(i-1,j  ,k  ,2)) &
+                                       + (Evec(i  ,j  ,k-1,2)-Evec(i  ,j  ,k  ,2))) / CCTK_DELTA_SPACE(3) &
+                              - 0.5d0 * ((Evec(i-1,j  ,k  ,3)-Evec(i-1,j-1,k  ,3)) &
+                                       + (Evec(i  ,j  ,k  ,3)-Evec(i  ,j-1,k  ,3))) / CCTK_DELTA_SPACE(2)
+          Bconsrhs(i,j,k,2) = - 0.5d0 * ((Evec(i-1,j-1,k  ,3)-Evec(i  ,j-1,k  ,3)) &
+                                       + (Evec(i-1,j  ,k  ,3)-Evec(i  ,j  ,k  ,3))) / CCTK_DELTA_SPACE(1) &
+                              - 0.5d0 * ((Evec(i  ,j-1,k  ,1)-Evec(i  ,j-1,k-1,1)) &
+                                       + (Evec(i  ,j  ,k  ,1)-Evec(i  ,j  ,k-1,1))) / CCTK_DELTA_SPACE(3)
+          Bconsrhs(i,j,k,3) = - 0.5d0 * ((Evec(i  ,j-1,k-1,1)-Evec(i  ,j  ,k-1,1)) &
+                                       + (Evec(i  ,j-1,k  ,1)-Evec(i  ,j  ,k  ,1))) / CCTK_DELTA_SPACE(2) &
+                              - 0.5d0 * ((Evec(i  ,j  ,k-1,2)-Evec(i-1,j  ,k-1,2)) &
+                                       + (Evec(i  ,j  ,k  ,2)-Evec(i-1,j  ,k  ,2))) / CCTK_DELTA_SPACE(1)
         enddo
       enddo
     enddo