diff options
Diffstat (limited to 'ML_ADM/src/ML_ADM_constraints_boundary.cc')
-rw-r--r-- | ML_ADM/src/ML_ADM_constraints_boundary.cc | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/ML_ADM/src/ML_ADM_constraints_boundary.cc b/ML_ADM/src/ML_ADM_constraints_boundary.cc index 01d49ca..e9cb98d 100644 --- a/ML_ADM/src/ML_ADM_constraints_boundary.cc +++ b/ML_ADM/src/ML_ADM_constraints_boundary.cc @@ -170,7 +170,7 @@ static void ML_ADM_constraints_boundary_Body(cGH const * restrict const cctkGH, CCTK_LOOP3STR(ML_ADM_constraints_boundary, 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; @@ -208,7 +208,7 @@ static void ML_ADM_constraints_boundary_Body(cGH const * restrict const cctkGH, CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED M3L = ToReal(0.); /* 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); vec_store_nta_partial(M1[index],M1L); vec_store_nta_partial(M2[index],M2L); |