diff options
Diffstat (limited to 'ML_CCZ4/src/ML_CCZ4_constraints1.cc')
| -rw-r--r-- | ML_CCZ4/src/ML_CCZ4_constraints1.cc | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/ML_CCZ4/src/ML_CCZ4_constraints1.cc b/ML_CCZ4/src/ML_CCZ4_constraints1.cc index 9e3c5ab..d79e218 100644 --- a/ML_CCZ4/src/ML_CCZ4_constraints1.cc +++ b/ML_CCZ4/src/ML_CCZ4_constraints1.cc @@ -209,7 +209,7 @@ static void ML_CCZ4_constraints1_Body(cGH const * restrict const cctkGH, int con CCTK_LOOP3STR(ML_CCZ4_constraints1, i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2], cctk_ash[0],cctk_ash[1],cctk_ash[2], - imin,imax, CCTK_REAL_VEC_SIZE) + kimin,kimax, CCTK_REAL_VEC_SIZE) { ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k; @@ -1617,7 +1617,7 @@ static void ML_CCZ4_constraints1_Body(cGH const * restrict const cctkGH, int con kadd(trR,kmadd(rho,ToReal(-50.26548245743669),kmadd(kmadd(Atm12,Atm21,kmadd(Atm13,Atm31,kmul(Atm23,Atm32))),ToReal(-2.),kmadd(kmadd(Atm11,Atm11,kmadd(Atm22,Atm22,kmul(Atm33,Atm33))),ToReal(-1.),kmul(kmul(trKL,trKL),ToReal(0.666666666666666666666666666667)))))); /* Copy local copies back to grid functions */ - vec_store_partial_prepare(i,lc_imin,lc_imax); + vec_store_partial_prepare(i,kimin,kimax); vec_store_nta_partial(H[index],HL); } CCTK_ENDLOOP3STR(ML_CCZ4_constraints1); |