From 6720d80eab38f9f24584bcf657bc85c1708ab701 Mon Sep 17 00:00:00 2001 From: Erik Schnetter Date: Sat, 19 Jan 2013 22:45:29 -0500 Subject: Regenerate code after Kranc change --- ML_BSSN/src/ML_BSSN_Advect.cc | 4 ++-- ML_BSSN/src/ML_BSSN_Dissipation.cc | 4 ++-- ML_BSSN/src/ML_BSSN_InitGamma.cc | 4 ++-- ML_BSSN/src/ML_BSSN_InitRHS.cc | 4 ++-- ML_BSSN/src/ML_BSSN_Minkowski.cc | 4 ++-- ML_BSSN/src/ML_BSSN_RHS1.cc | 4 ++-- ML_BSSN/src/ML_BSSN_RHS2.cc | 4 ++-- ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc | 4 ++-- ML_BSSN/src/ML_BSSN_boundary.cc | 4 ++-- ML_BSSN/src/ML_BSSN_constraints1.cc | 4 ++-- ML_BSSN/src/ML_BSSN_constraints2.cc | 4 ++-- ML_BSSN/src/ML_BSSN_convertFromADMBase.cc | 4 ++-- ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc | 4 ++-- ML_BSSN/src/ML_BSSN_convertToADMBase.cc | 4 ++-- ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc | 4 ++-- ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc | 4 ++-- ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc | 4 ++-- ML_BSSN/src/ML_BSSN_enforce.cc | 4 ++-- 18 files changed, 36 insertions(+), 36 deletions(-) (limited to 'ML_BSSN') diff --git a/ML_BSSN/src/ML_BSSN_Advect.cc b/ML_BSSN/src/ML_BSSN_Advect.cc index 94ea8e6..d74a205 100644 --- a/ML_BSSN/src/ML_BSSN_Advect.cc +++ b/ML_BSSN/src/ML_BSSN_Advect.cc @@ -233,7 +233,7 @@ static void ML_BSSN_Advect_Body(cGH const * restrict const cctkGH, int const dir CCTK_LOOP3STR(ML_BSSN_Advect, 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; @@ -2068,7 +2068,7 @@ static void ML_BSSN_Advect_Body(cGH const * restrict const cctkGH, int const dir + ShiftAdvectionCoeff)))),ToReal(ShiftBCoeff),B3rhsL); /* 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(alpharhs[index],alpharhsL); vec_store_nta_partial(Arhs[index],ArhsL); vec_store_nta_partial(At11rhs[index],At11rhsL); diff --git a/ML_BSSN/src/ML_BSSN_Dissipation.cc b/ML_BSSN/src/ML_BSSN_Dissipation.cc index a2645cb..d2a0715 100644 --- a/ML_BSSN/src/ML_BSSN_Dissipation.cc +++ b/ML_BSSN/src/ML_BSSN_Dissipation.cc @@ -233,7 +233,7 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons CCTK_LOOP3STR(ML_BSSN_Dissipation, 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; @@ -1239,7 +1239,7 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons kadd(B3rhsL,kmadd(epsdiss1,JacPDdissipationNth1B3,kmadd(epsdiss2,JacPDdissipationNth2B3,kmul(epsdiss3,JacPDdissipationNth3B3)))); /* 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(alpharhs[index],alpharhsL); vec_store_nta_partial(Arhs[index],ArhsL); vec_store_nta_partial(At11rhs[index],At11rhsL); diff --git a/ML_BSSN/src/ML_BSSN_InitGamma.cc b/ML_BSSN/src/ML_BSSN_InitGamma.cc index 968ac6b..c37ebc8 100644 --- a/ML_BSSN/src/ML_BSSN_InitGamma.cc +++ b/ML_BSSN/src/ML_BSSN_InitGamma.cc @@ -197,7 +197,7 @@ static void ML_BSSN_InitGamma_Body(cGH const * restrict const cctkGH, int const CCTK_LOOP3STR(ML_BSSN_InitGamma, 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; @@ -241,7 +241,7 @@ static void ML_BSSN_InitGamma_Body(cGH const * restrict const cctkGH, int const CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED B3L = 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(A[index],AL); vec_store_nta_partial(B1[index],B1L); vec_store_nta_partial(B2[index],B2L); diff --git a/ML_BSSN/src/ML_BSSN_InitRHS.cc b/ML_BSSN/src/ML_BSSN_InitRHS.cc index dde5f60..90a93bc 100644 --- a/ML_BSSN/src/ML_BSSN_InitRHS.cc +++ b/ML_BSSN/src/ML_BSSN_InitRHS.cc @@ -197,7 +197,7 @@ static void ML_BSSN_InitRHS_Body(cGH const * restrict const cctkGH, int const di CCTK_LOOP3STR(ML_BSSN_InitRHS, 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; @@ -277,7 +277,7 @@ static void ML_BSSN_InitRHS_Body(cGH const * restrict const cctkGH, int const di CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED B3rhsL = 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(alpharhs[index],alpharhsL); vec_store_nta_partial(Arhs[index],ArhsL); vec_store_nta_partial(At11rhs[index],At11rhsL); diff --git a/ML_BSSN/src/ML_BSSN_Minkowski.cc b/ML_BSSN/src/ML_BSSN_Minkowski.cc index d5ac436..812f9e1 100644 --- a/ML_BSSN/src/ML_BSSN_Minkowski.cc +++ b/ML_BSSN/src/ML_BSSN_Minkowski.cc @@ -197,7 +197,7 @@ static void ML_BSSN_Minkowski_Body(cGH const * restrict const cctkGH, int const CCTK_LOOP3STR(ML_BSSN_Minkowski, 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; @@ -278,7 +278,7 @@ static void ML_BSSN_Minkowski_Body(cGH const * restrict const cctkGH, int const CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED B3L = 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(A[index],AL); vec_store_nta_partial(alpha[index],alphaL); vec_store_nta_partial(At11[index],At11L); diff --git a/ML_BSSN/src/ML_BSSN_RHS1.cc b/ML_BSSN/src/ML_BSSN_RHS1.cc index 31b0601..66a02ac 100644 --- a/ML_BSSN/src/ML_BSSN_RHS1.cc +++ b/ML_BSSN/src/ML_BSSN_RHS1.cc @@ -230,7 +230,7 @@ static void ML_BSSN_RHS1_Body(cGH const * restrict const cctkGH, int const dir, CCTK_LOOP3STR(ML_BSSN_RHS1, 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; @@ -1398,7 +1398,7 @@ static void ML_BSSN_RHS1_Body(cGH const * restrict const cctkGH, int const dir, kmul(kmadd(B3L,kmul(eta,ToReal(-1.*BetaDriver)),dotXt3),ToReal(ShiftBCoeff)); /* 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(alpharhs[index],alpharhsL); vec_store_nta_partial(Arhs[index],ArhsL); vec_store_nta_partial(B1rhs[index],B1rhsL); diff --git a/ML_BSSN/src/ML_BSSN_RHS2.cc b/ML_BSSN/src/ML_BSSN_RHS2.cc index 02b35fe..d27c0b7 100644 --- a/ML_BSSN/src/ML_BSSN_RHS2.cc +++ b/ML_BSSN/src/ML_BSSN_RHS2.cc @@ -209,7 +209,7 @@ static void ML_BSSN_RHS2_Body(cGH const * restrict const cctkGH, int const dir, CCTK_LOOP3STR(ML_BSSN_RHS2, 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; @@ -1819,7 +1819,7 @@ static void ML_BSSN_RHS2_Body(cGH const * restrict const cctkGH, int const dir, kmadd(At33L,kmul(kadd(JacPDstandardNth1beta1,kadd(JacPDstandardNth2beta2,JacPDstandardNth3beta3)),ToReal(-0.666666666666666666666666666667)),kmadd(em4phi,kmadd(g33,kmul(trAts,ToReal(-0.333333333333333333333333333333)),Ats33),kmadd(alphaL,kmadd(At33L,trKL,kmadd(kmadd(At13L,Atm13,kmadd(At23L,Atm23,kmul(At33L,Atm33))),ToReal(-2.),kmul(em4phi,kmul(ToReal(-25.132741228718345),kmadd(g33,kmul(trS,ToReal(-0.333333333333333333333333333333)),eTzzL))))),kmul(kmadd(At13L,JacPDstandardNth3beta1,kmadd(At23L,JacPDstandardNth3beta2,kmul(At33L,JacPDstandardNth3beta3))),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(At11rhs[index],At11rhsL); vec_store_nta_partial(At12rhs[index],At12rhsL); vec_store_nta_partial(At13rhs[index],At13rhsL); diff --git a/ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc b/ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc index c69d486..08fc3dd 100644 --- a/ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc +++ b/ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc @@ -233,7 +233,7 @@ static void ML_BSSN_RHSStaticBoundary_Body(cGH const * restrict const cctkGH, in CCTK_LOOP3STR(ML_BSSN_RHSStaticBoundary, 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; @@ -313,7 +313,7 @@ static void ML_BSSN_RHSStaticBoundary_Body(cGH const * restrict const cctkGH, in CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED B3rhsL = 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(alpharhs[index],alpharhsL); vec_store_nta_partial(Arhs[index],ArhsL); vec_store_nta_partial(At11rhs[index],At11rhsL); diff --git a/ML_BSSN/src/ML_BSSN_boundary.cc b/ML_BSSN/src/ML_BSSN_boundary.cc index e838f14..e8e87b5 100644 --- a/ML_BSSN/src/ML_BSSN_boundary.cc +++ b/ML_BSSN/src/ML_BSSN_boundary.cc @@ -233,7 +233,7 @@ static void ML_BSSN_boundary_Body(cGH const * restrict const cctkGH, int const d CCTK_LOOP3STR(ML_BSSN_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; @@ -314,7 +314,7 @@ static void ML_BSSN_boundary_Body(cGH const * restrict const cctkGH, int const d CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED B3L = 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(A[index],AL); vec_store_nta_partial(alpha[index],alphaL); vec_store_nta_partial(At11[index],At11L); diff --git a/ML_BSSN/src/ML_BSSN_constraints1.cc b/ML_BSSN/src/ML_BSSN_constraints1.cc index 096e4d4..5c880fd 100644 --- a/ML_BSSN/src/ML_BSSN_constraints1.cc +++ b/ML_BSSN/src/ML_BSSN_constraints1.cc @@ -209,7 +209,7 @@ static void ML_BSSN_constraints1_Body(cGH const * restrict const cctkGH, int con CCTK_LOOP3STR(ML_BSSN_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; @@ -1578,7 +1578,7 @@ static void ML_BSSN_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_BSSN_constraints1); 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); diff --git a/ML_BSSN/src/ML_BSSN_convertFromADMBase.cc b/ML_BSSN/src/ML_BSSN_convertFromADMBase.cc index 49f4382..8b5be6d 100644 --- a/ML_BSSN/src/ML_BSSN_convertFromADMBase.cc +++ b/ML_BSSN/src/ML_BSSN_convertFromADMBase.cc @@ -197,7 +197,7 @@ static void ML_BSSN_convertFromADMBase_Body(cGH const * restrict const cctkGH, i CCTK_LOOP3STR(ML_BSSN_convertFromADMBase, 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; @@ -334,7 +334,7 @@ static void ML_BSSN_convertFromADMBase_Body(cGH const * restrict const cctkGH, i CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED beta3L = betazL; /* 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(alpha[index],alphaL); vec_store_nta_partial(At11[index],At11L); vec_store_nta_partial(At12[index],At12L); diff --git a/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc b/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc index 253b6a0..e9db745 100644 --- a/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc +++ b/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc @@ -215,7 +215,7 @@ static void ML_BSSN_convertFromADMBaseGamma_Body(cGH const * restrict const cctk CCTK_LOOP3STR(ML_BSSN_convertFromADMBaseGamma, 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; @@ -861,7 +861,7 @@ static void ML_BSSN_convertFromADMBaseGamma_Body(cGH const * restrict const cctk } /* 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(A[index],AL); vec_store_nta_partial(B1[index],B1L); vec_store_nta_partial(B2[index],B2L); diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBase.cc b/ML_BSSN/src/ML_BSSN_convertToADMBase.cc index 24c249e..36d0c7d 100644 --- a/ML_BSSN/src/ML_BSSN_convertToADMBase.cc +++ b/ML_BSSN/src/ML_BSSN_convertToADMBase.cc @@ -197,7 +197,7 @@ static void ML_BSSN_convertToADMBase_Body(cGH const * restrict const cctkGH, int CCTK_LOOP3STR(ML_BSSN_convertToADMBase, 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; @@ -292,7 +292,7 @@ static void ML_BSSN_convertToADMBase_Body(cGH const * restrict const cctkGH, int CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED betazL = beta3L; /* 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(alp[index],alpL); vec_store_nta_partial(betax[index],betaxL); vec_store_nta_partial(betay[index],betayL); diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc index 3f243e2..d10b5d6 100644 --- a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc +++ b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc @@ -212,7 +212,7 @@ static void ML_BSSN_convertToADMBaseDtLapseShift_Body(cGH const * restrict const CCTK_LOOP3STR(ML_BSSN_convertToADMBaseDtLapseShift, 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; @@ -857,7 +857,7 @@ static void ML_BSSN_convertToADMBaseDtLapseShift_Body(cGH const * restrict const + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff))))); /* 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(dtalp[index],dtalpL); vec_store_nta_partial(dtbetax[index],dtbetaxL); vec_store_nta_partial(dtbetay[index],dtbetayL); diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc index 2f00e65..f50e8c2 100644 --- a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc +++ b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc @@ -212,7 +212,7 @@ static void ML_BSSN_convertToADMBaseDtLapseShiftBoundary_Body(cGH const * restri CCTK_LOOP3STR(ML_BSSN_convertToADMBaseDtLapseShiftBoundary, 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; @@ -291,7 +291,7 @@ static void ML_BSSN_convertToADMBaseDtLapseShiftBoundary_Body(cGH const * restri } /* 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(dtalp[index],dtalpL); vec_store_nta_partial(dtbetax[index],dtbetaxL); vec_store_nta_partial(dtbetay[index],dtbetayL); diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc b/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc index 6de84e5..dd877ee 100644 --- a/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc +++ b/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc @@ -197,7 +197,7 @@ static void ML_BSSN_convertToADMBaseFakeDtLapseShift_Body(cGH const * restrict c CCTK_LOOP3STR(ML_BSSN_convertToADMBaseFakeDtLapseShift, 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; @@ -276,7 +276,7 @@ static void ML_BSSN_convertToADMBaseFakeDtLapseShift_Body(cGH const * restrict c } /* 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(dtalp[index],dtalpL); vec_store_nta_partial(dtbetax[index],dtbetaxL); vec_store_nta_partial(dtbetay[index],dtbetayL); diff --git a/ML_BSSN/src/ML_BSSN_enforce.cc b/ML_BSSN/src/ML_BSSN_enforce.cc index 71b2032..62f0a87 100644 --- a/ML_BSSN/src/ML_BSSN_enforce.cc +++ b/ML_BSSN/src/ML_BSSN_enforce.cc @@ -197,7 +197,7 @@ static void ML_BSSN_enforce_Body(cGH const * restrict const cctkGH, int const di CCTK_LOOP3STR(ML_BSSN_enforce, 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; @@ -283,7 +283,7 @@ static void ML_BSSN_enforce_Body(cGH const * restrict const cctkGH, int const di alphaL = kfmax(alphaL,ToReal(MinimumLapse)); /* 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(alpha[index],alphaL); vec_store_nta_partial(At11[index],At11L); vec_store_nta_partial(At12[index],At12L); -- cgit v1.2.3