diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_constraints2.cc')
-rw-r--r-- | ML_BSSN/src/ML_BSSN_constraints2.cc | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/ML_BSSN/src/ML_BSSN_constraints2.cc b/ML_BSSN/src/ML_BSSN_constraints2.cc index a132d2c..362f533 100644 --- a/ML_BSSN/src/ML_BSSN_constraints2.cc +++ b/ML_BSSN/src/ML_BSSN_constraints2.cc @@ -218,7 +218,7 @@ static void ML_BSSN_constraints2_Body(cGH const * restrict const cctkGH, int con CCTK_LOOP3STR(ML_BSSN_constraints2, 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; @@ -949,7 +949,7 @@ static void ML_BSSN_constraints2_Body(cGH const * restrict const cctkGH, int con kmadd(At11L,gtu11,kmadd(At22L,gtu22,kmadd(At33L,gtu33,kmul(kmadd(At12L,gtu12,kmadd(At13L,gtu13,kmul(At23L,gtu23))),ToReal(2.))))); /* 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(cA[index],cAL); vec_store_nta_partial(cS[index],cSL); vec_store_nta_partial(cXt1[index],cXt1L); |