From f9fe9736f1da6d9b5c52978d2df84a2cc1988901 Mon Sep 17 00:00:00 2001 From: eschnett Date: Sun, 20 Jan 2013 15:10:24 +0000 Subject: Regenerate code after Kranc change git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/WeylScal4/trunk@131 4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843 --- src/WeylScal4_invars_calc_2nd.cc | 4 ++-- src/WeylScal4_invars_calc_4th.cc | 4 ++-- src/WeylScal4_invars_calc_Nth.cc | 4 ++-- src/WeylScal4_psi4_calc_2nd.cc | 4 ++-- src/WeylScal4_psi4_calc_4th.cc | 4 ++-- src/WeylScal4_psi4_calc_Nth.cc | 4 ++-- src/WeylScal4_psis_calc_2nd.cc | 4 ++-- src/WeylScal4_psis_calc_4th.cc | 4 ++-- src/WeylScal4_psis_calc_Nth.cc | 4 ++-- 9 files changed, 18 insertions(+), 18 deletions(-) diff --git a/src/WeylScal4_invars_calc_2nd.cc b/src/WeylScal4_invars_calc_2nd.cc index e466c5b..65eb394 100644 --- a/src/WeylScal4_invars_calc_2nd.cc +++ b/src/WeylScal4_invars_calc_2nd.cc @@ -155,7 +155,7 @@ static void WeylScal4_invars_calc_2nd_Body(cGH const * restrict const cctkGH, in CCTK_LOOP3STR(WeylScal4_invars_calc_2nd, 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; @@ -219,7 +219,7 @@ static void WeylScal4_invars_calc_2nd_Body(cGH const * restrict const cctkGH, in kmul(ToReal(-640.),kmadd(Psi1iL,kmul(kmadd(kmadd(Psi3rL,Psi4iL,kmul(Psi3iL,Psi4rL)),kmul(Psi1rL,Psi1rL),kmul(Psi0rL,kmul(Psi3iL,kmul(Psi3rL,Psi3rL)))),ToReal(-12.)),kmadd(kmadd(Psi3iL,kmul(Psi4iL,kmul(Psi1rL,kmul(Psi1rL,Psi1rL))),kmul(Psi0iL,kmul(Psi1iL,kmul(Psi3rL,kmul(Psi3rL,Psi3rL))))),ToReal(-4.),kmadd(kpow(Psi2rL,5.),ToReal(-3.),kmadd(Psi4rL,kmadd(kmul(Psi0rL,Psi0rL),kmul(Psi3iL,Psi3iL),kmadd(kmul(Psi0iL,Psi0iL),kmul(Psi3rL,Psi3rL),kmul(kmadd(kmul(Psi0iL,Psi0iL),kmul(Psi3iL,Psi3iL),kmul(kmul(Psi0rL,Psi0rL),kmul(Psi3rL,Psi3rL))),ToReal(-1.)))),kmadd(Psi0rL,kmadd(kmul(Psi1rL,Psi1rL),kmul(Psi4iL,Psi4iL),kmadd(kmul(Psi1iL,Psi1iL),kmul(Psi4rL,Psi4rL),kmul(kmadd(kmul(Psi1iL,Psi1iL),kmul(Psi4iL,Psi4iL),kmul(kmul(Psi1rL,Psi1rL),kmul(Psi4rL,Psi4rL))),ToReal(-1.)))),kmadd(kmadd(kmadd(Psi3rL,Psi4iL,kmul(Psi3iL,Psi4rL)),kmul(Psi1iL,kmul(Psi1iL,Psi1iL)),kmadd(Psi3rL,kmul(Psi4rL,kmul(Psi1rL,kmul(Psi1rL,Psi1rL))),kmadd(Psi0rL,kmadd(kmadd(Psi0iL,kmul(Psi3iL,Psi3rL),kmul(Psi1iL,kmul(Psi1rL,Psi4iL))),Psi4rL,kmul(Psi1iL,kmul(Psi3iL,kmul(Psi3iL,Psi3iL)))),kmul(Psi1rL,kmadd(Psi0iL,kmul(Psi3iL,kmul(Psi3iL,Psi3iL)),kmul(Psi0rL,kmul(Psi3rL,kmul(Psi3rL,Psi3rL)))))))),ToReal(4.),kmadd(ToReal(-2.),kmadd(Psi4iL,kmadd(Psi3iL,kmul(Psi3rL,kmul(Psi0iL,Psi0iL)),kmul(Psi0iL,kmadd(Psi4rL,kmul(Psi1iL,Psi1iL),kmul(Psi0rL,kmul(Psi3iL,Psi3iL))))),kmul(kmul(Psi2rL,kmul(Psi2rL,Psi2rL)),kmadd(Psi0iL,Psi4iL,kmadd(Psi1rL,kmul(Psi3rL,ToReal(-5.)),kmadd(Psi0rL,kmul(Psi4rL,ToReal(-1.)),kmul(Psi1iL,kmul(Psi3iL,ToReal(5.)))))))),kmadd(ToReal(2.),kmadd(Psi4iL,kmadd(Psi3iL,kmul(Psi3rL,kmul(Psi0rL,Psi0rL)),kmul(Psi0iL,kmadd(Psi4rL,kmul(Psi1rL,Psi1rL),kmul(Psi0rL,kmul(Psi3rL,Psi3rL))))),kmadd(kmul(Psi2iL,kmul(Psi2iL,Psi2iL)),kmadd(Psi0rL,Psi4iL,kmadd(Psi0iL,Psi4rL,kmul(kmadd(Psi1rL,Psi3iL,kmul(Psi1iL,Psi3rL)),ToReal(5.)))),kmul(Psi2iL,kmadd(Psi3rL,kmadd(Psi1iL,kmul(kmul(Psi2rL,Psi2rL),ToReal(-15.)),kmul(Psi3iL,kmul(kmul(Psi1iL,Psi1iL),ToReal(-8.)))),kmadd(Psi4iL,kmadd(Psi2rL,kmul(kmul(Psi1iL,Psi1iL),ToReal(-3.)),kmul(Psi4rL,kmadd(Psi0iL,Psi0iL,kmul(kmul(Psi0rL,Psi0rL),ToReal(-1.))))),kmadd(Psi0iL,kmul(ToReal(-1.),kmadd(Psi1iL,kmadd(Psi3rL,Psi4iL,kmul(Psi3iL,Psi4rL)),kmadd(Psi2rL,kmul(kmul(Psi3rL,Psi3rL),ToReal(-3.)),kmadd(Psi0rL,kmadd(Psi4rL,Psi4rL,kmul(kmul(Psi4iL,Psi4iL),ToReal(-1.))),kmul(kmadd(Psi4rL,kmul(Psi2rL,Psi2rL),kmul(Psi2rL,kmul(Psi3iL,Psi3iL))),ToReal(3.)))))),kmadd(Psi0rL,kmadd(Psi1iL,kmul(Psi3rL,Psi4rL),kmadd(Psi4iL,kmul(kmul(Psi2rL,Psi2rL),ToReal(-3.)),kmul(Psi3iL,kmadd(Psi1iL,kmul(Psi4iL,ToReal(-1.)),kmul(Psi2rL,kmul(Psi3rL,ToReal(6.))))))),kmadd(kmul(Psi1rL,Psi1rL),kmadd(Psi2rL,kmul(Psi4iL,ToReal(3.)),kmul(Psi3iL,kmul(Psi3rL,ToReal(8.)))),kmul(Psi1rL,kmadd(Psi3iL,kmul(kmul(Psi2rL,Psi2rL),ToReal(-15.)),kmadd(Psi4iL,kmadd(Psi0rL,Psi3rL,kmul(Psi0iL,kmul(Psi3iL,ToReal(-1.)))),kmadd(Psi4rL,kmadd(Psi0rL,Psi3iL,kmadd(Psi0iL,Psi3rL,kmul(Psi1iL,kmul(Psi2rL,ToReal(6.))))),kmul(Psi1iL,kmadd(kmul(Psi3iL,Psi3iL),ToReal(-8.),kmul(kmul(Psi3rL,Psi3rL),ToReal(8.)))))))))))))))),kmadd(ToReal(3.),kmadd(kmul(Psi2rL,Psi2rL),kmadd(Psi4rL,kmadd(Psi1iL,Psi1iL,kmul(kmul(Psi1rL,Psi1rL),ToReal(-1.))),kmadd(Psi0rL,kmadd(Psi3iL,Psi3iL,kmul(kmul(Psi3rL,Psi3rL),ToReal(-1.))),kmul(kmadd(Psi0iL,kmul(Psi3iL,Psi3rL),kmul(Psi1iL,kmul(Psi1rL,Psi4iL))),ToReal(2.)))),kmul(kmul(Psi2iL,Psi2iL),kmadd(kmadd(Psi0iL,kmul(Psi3iL,Psi3rL),kmul(Psi1iL,kmul(Psi1rL,Psi4iL))),ToReal(-2.),kmadd(Psi4rL,kmadd(Psi1rL,Psi1rL,kmul(kmul(Psi1iL,Psi1iL),ToReal(-1.))),kmadd(Psi0rL,kmadd(Psi3rL,Psi3rL,kmul(kmul(Psi3iL,Psi3iL),ToReal(-1.))),kmadd(Psi2rL,kmul(ToReal(2.),kmadd(Psi0iL,Psi4iL,kmadd(Psi1rL,kmul(Psi3rL,ToReal(-5.)),kmadd(Psi0rL,kmul(Psi4rL,ToReal(-1.)),kmul(Psi1iL,kmul(Psi3iL,ToReal(5.))))))),kmul(kmul(Psi2rL,kmul(Psi2rL,Psi2rL)),ToReal(10.)))))))),kmadd(Psi1rL,kmadd(Psi0rL,kmul(Psi3rL,kmul(kmul(Psi3iL,Psi3iL),ToReal(-12.))),kmadd(Psi0iL,kmul(Psi1iL,kmul(kmul(Psi4iL,Psi4iL),ToReal(-2.))),kmul(kmul(Psi1iL,Psi1iL),kmadd(Psi3rL,kmul(Psi4rL,ToReal(-12.)),kmul(Psi3iL,kmul(Psi4iL,ToReal(12.))))))),kmadd(Psi0iL,kmadd(Psi1rL,kmul(Psi3iL,kmul(kmul(Psi3rL,Psi3rL),ToReal(-12.))),kmul(Psi1iL,kmadd(Psi1rL,kmul(kmul(Psi4rL,Psi4rL),ToReal(2.)),kmul(Psi3rL,kmul(kmul(Psi3iL,Psi3iL),ToReal(12.)))))),kmul(Psi2rL,kmadd(kmul(Psi0iL,Psi0iL),kmul(Psi4iL,Psi4iL),kmadd(kmul(Psi0rL,Psi0rL),kmul(Psi4rL,Psi4rL),kmadd(kmul(kmul(Psi2iL,Psi2iL),kmul(Psi2iL,Psi2iL)),ToReal(-15.),kmadd(Psi0iL,kmul(Psi0rL,kmul(Psi4iL,kmul(Psi4rL,ToReal(-4.)))),kmadd(kmadd(kmul(Psi0rL,Psi0rL),kmul(Psi4iL,Psi4iL),kmul(kmul(Psi0iL,Psi0iL),kmul(Psi4rL,Psi4rL))),ToReal(-1.),kmadd(kmadd(Psi3iL,Psi3iL,kmul(kmul(Psi3rL,Psi3rL),ToReal(-1.))),kmadd(kmul(Psi1iL,Psi1iL),ToReal(-8.),kmul(kmul(Psi1rL,Psi1rL),ToReal(8.))),kmul(ToReal(2.),kmadd(Psi1rL,kmadd(kmadd(Psi0rL,Psi3iL,kmul(Psi0iL,Psi3rL)),Psi4iL,kmul(Psi4rL,kmadd(Psi0iL,Psi3iL,kmul(Psi0rL,kmul(Psi3rL,ToReal(-1.)))))),kmul(Psi1iL,kmadd(Psi0iL,kmul(Psi3rL,Psi4rL),kmadd(Psi0rL,kmadd(Psi3rL,Psi4iL,kmul(Psi3iL,Psi4rL)),kmul(Psi3iL,kmadd(Psi0iL,kmul(Psi4iL,ToReal(-1.)),kmul(Psi1rL,kmul(Psi3rL,ToReal(16.))))))))))))))))))))))))))))); /* 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(curvIi[index],curvIiL); vec_store_nta_partial(curvIr[index],curvIrL); vec_store_nta_partial(curvJ1[index],curvJ1L); diff --git a/src/WeylScal4_invars_calc_4th.cc b/src/WeylScal4_invars_calc_4th.cc index b6f6114..6e1d761 100644 --- a/src/WeylScal4_invars_calc_4th.cc +++ b/src/WeylScal4_invars_calc_4th.cc @@ -155,7 +155,7 @@ static void WeylScal4_invars_calc_4th_Body(cGH const * restrict const cctkGH, in CCTK_LOOP3STR(WeylScal4_invars_calc_4th, 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; @@ -219,7 +219,7 @@ static void WeylScal4_invars_calc_4th_Body(cGH const * restrict const cctkGH, in kmul(ToReal(-640.),kmadd(Psi1iL,kmul(kmadd(kmadd(Psi3rL,Psi4iL,kmul(Psi3iL,Psi4rL)),kmul(Psi1rL,Psi1rL),kmul(Psi0rL,kmul(Psi3iL,kmul(Psi3rL,Psi3rL)))),ToReal(-12.)),kmadd(kmadd(Psi3iL,kmul(Psi4iL,kmul(Psi1rL,kmul(Psi1rL,Psi1rL))),kmul(Psi0iL,kmul(Psi1iL,kmul(Psi3rL,kmul(Psi3rL,Psi3rL))))),ToReal(-4.),kmadd(kpow(Psi2rL,5.),ToReal(-3.),kmadd(Psi4rL,kmadd(kmul(Psi0rL,Psi0rL),kmul(Psi3iL,Psi3iL),kmadd(kmul(Psi0iL,Psi0iL),kmul(Psi3rL,Psi3rL),kmul(kmadd(kmul(Psi0iL,Psi0iL),kmul(Psi3iL,Psi3iL),kmul(kmul(Psi0rL,Psi0rL),kmul(Psi3rL,Psi3rL))),ToReal(-1.)))),kmadd(Psi0rL,kmadd(kmul(Psi1rL,Psi1rL),kmul(Psi4iL,Psi4iL),kmadd(kmul(Psi1iL,Psi1iL),kmul(Psi4rL,Psi4rL),kmul(kmadd(kmul(Psi1iL,Psi1iL),kmul(Psi4iL,Psi4iL),kmul(kmul(Psi1rL,Psi1rL),kmul(Psi4rL,Psi4rL))),ToReal(-1.)))),kmadd(kmadd(kmadd(Psi3rL,Psi4iL,kmul(Psi3iL,Psi4rL)),kmul(Psi1iL,kmul(Psi1iL,Psi1iL)),kmadd(Psi3rL,kmul(Psi4rL,kmul(Psi1rL,kmul(Psi1rL,Psi1rL))),kmadd(Psi0rL,kmadd(kmadd(Psi0iL,kmul(Psi3iL,Psi3rL),kmul(Psi1iL,kmul(Psi1rL,Psi4iL))),Psi4rL,kmul(Psi1iL,kmul(Psi3iL,kmul(Psi3iL,Psi3iL)))),kmul(Psi1rL,kmadd(Psi0iL,kmul(Psi3iL,kmul(Psi3iL,Psi3iL)),kmul(Psi0rL,kmul(Psi3rL,kmul(Psi3rL,Psi3rL)))))))),ToReal(4.),kmadd(ToReal(-2.),kmadd(Psi4iL,kmadd(Psi3iL,kmul(Psi3rL,kmul(Psi0iL,Psi0iL)),kmul(Psi0iL,kmadd(Psi4rL,kmul(Psi1iL,Psi1iL),kmul(Psi0rL,kmul(Psi3iL,Psi3iL))))),kmul(kmul(Psi2rL,kmul(Psi2rL,Psi2rL)),kmadd(Psi0iL,Psi4iL,kmadd(Psi1rL,kmul(Psi3rL,ToReal(-5.)),kmadd(Psi0rL,kmul(Psi4rL,ToReal(-1.)),kmul(Psi1iL,kmul(Psi3iL,ToReal(5.)))))))),kmadd(ToReal(2.),kmadd(Psi4iL,kmadd(Psi3iL,kmul(Psi3rL,kmul(Psi0rL,Psi0rL)),kmul(Psi0iL,kmadd(Psi4rL,kmul(Psi1rL,Psi1rL),kmul(Psi0rL,kmul(Psi3rL,Psi3rL))))),kmadd(kmul(Psi2iL,kmul(Psi2iL,Psi2iL)),kmadd(Psi0rL,Psi4iL,kmadd(Psi0iL,Psi4rL,kmul(kmadd(Psi1rL,Psi3iL,kmul(Psi1iL,Psi3rL)),ToReal(5.)))),kmul(Psi2iL,kmadd(Psi3rL,kmadd(Psi1iL,kmul(kmul(Psi2rL,Psi2rL),ToReal(-15.)),kmul(Psi3iL,kmul(kmul(Psi1iL,Psi1iL),ToReal(-8.)))),kmadd(Psi4iL,kmadd(Psi2rL,kmul(kmul(Psi1iL,Psi1iL),ToReal(-3.)),kmul(Psi4rL,kmadd(Psi0iL,Psi0iL,kmul(kmul(Psi0rL,Psi0rL),ToReal(-1.))))),kmadd(Psi0iL,kmul(ToReal(-1.),kmadd(Psi1iL,kmadd(Psi3rL,Psi4iL,kmul(Psi3iL,Psi4rL)),kmadd(Psi2rL,kmul(kmul(Psi3rL,Psi3rL),ToReal(-3.)),kmadd(Psi0rL,kmadd(Psi4rL,Psi4rL,kmul(kmul(Psi4iL,Psi4iL),ToReal(-1.))),kmul(kmadd(Psi4rL,kmul(Psi2rL,Psi2rL),kmul(Psi2rL,kmul(Psi3iL,Psi3iL))),ToReal(3.)))))),kmadd(Psi0rL,kmadd(Psi1iL,kmul(Psi3rL,Psi4rL),kmadd(Psi4iL,kmul(kmul(Psi2rL,Psi2rL),ToReal(-3.)),kmul(Psi3iL,kmadd(Psi1iL,kmul(Psi4iL,ToReal(-1.)),kmul(Psi2rL,kmul(Psi3rL,ToReal(6.))))))),kmadd(kmul(Psi1rL,Psi1rL),kmadd(Psi2rL,kmul(Psi4iL,ToReal(3.)),kmul(Psi3iL,kmul(Psi3rL,ToReal(8.)))),kmul(Psi1rL,kmadd(Psi3iL,kmul(kmul(Psi2rL,Psi2rL),ToReal(-15.)),kmadd(Psi4iL,kmadd(Psi0rL,Psi3rL,kmul(Psi0iL,kmul(Psi3iL,ToReal(-1.)))),kmadd(Psi4rL,kmadd(Psi0rL,Psi3iL,kmadd(Psi0iL,Psi3rL,kmul(Psi1iL,kmul(Psi2rL,ToReal(6.))))),kmul(Psi1iL,kmadd(kmul(Psi3iL,Psi3iL),ToReal(-8.),kmul(kmul(Psi3rL,Psi3rL),ToReal(8.)))))))))))))))),kmadd(ToReal(3.),kmadd(kmul(Psi2rL,Psi2rL),kmadd(Psi4rL,kmadd(Psi1iL,Psi1iL,kmul(kmul(Psi1rL,Psi1rL),ToReal(-1.))),kmadd(Psi0rL,kmadd(Psi3iL,Psi3iL,kmul(kmul(Psi3rL,Psi3rL),ToReal(-1.))),kmul(kmadd(Psi0iL,kmul(Psi3iL,Psi3rL),kmul(Psi1iL,kmul(Psi1rL,Psi4iL))),ToReal(2.)))),kmul(kmul(Psi2iL,Psi2iL),kmadd(kmadd(Psi0iL,kmul(Psi3iL,Psi3rL),kmul(Psi1iL,kmul(Psi1rL,Psi4iL))),ToReal(-2.),kmadd(Psi4rL,kmadd(Psi1rL,Psi1rL,kmul(kmul(Psi1iL,Psi1iL),ToReal(-1.))),kmadd(Psi0rL,kmadd(Psi3rL,Psi3rL,kmul(kmul(Psi3iL,Psi3iL),ToReal(-1.))),kmadd(Psi2rL,kmul(ToReal(2.),kmadd(Psi0iL,Psi4iL,kmadd(Psi1rL,kmul(Psi3rL,ToReal(-5.)),kmadd(Psi0rL,kmul(Psi4rL,ToReal(-1.)),kmul(Psi1iL,kmul(Psi3iL,ToReal(5.))))))),kmul(kmul(Psi2rL,kmul(Psi2rL,Psi2rL)),ToReal(10.)))))))),kmadd(Psi1rL,kmadd(Psi0rL,kmul(Psi3rL,kmul(kmul(Psi3iL,Psi3iL),ToReal(-12.))),kmadd(Psi0iL,kmul(Psi1iL,kmul(kmul(Psi4iL,Psi4iL),ToReal(-2.))),kmul(kmul(Psi1iL,Psi1iL),kmadd(Psi3rL,kmul(Psi4rL,ToReal(-12.)),kmul(Psi3iL,kmul(Psi4iL,ToReal(12.))))))),kmadd(Psi0iL,kmadd(Psi1rL,kmul(Psi3iL,kmul(kmul(Psi3rL,Psi3rL),ToReal(-12.))),kmul(Psi1iL,kmadd(Psi1rL,kmul(kmul(Psi4rL,Psi4rL),ToReal(2.)),kmul(Psi3rL,kmul(kmul(Psi3iL,Psi3iL),ToReal(12.)))))),kmul(Psi2rL,kmadd(kmul(Psi0iL,Psi0iL),kmul(Psi4iL,Psi4iL),kmadd(kmul(Psi0rL,Psi0rL),kmul(Psi4rL,Psi4rL),kmadd(kmul(kmul(Psi2iL,Psi2iL),kmul(Psi2iL,Psi2iL)),ToReal(-15.),kmadd(Psi0iL,kmul(Psi0rL,kmul(Psi4iL,kmul(Psi4rL,ToReal(-4.)))),kmadd(kmadd(kmul(Psi0rL,Psi0rL),kmul(Psi4iL,Psi4iL),kmul(kmul(Psi0iL,Psi0iL),kmul(Psi4rL,Psi4rL))),ToReal(-1.),kmadd(kmadd(Psi3iL,Psi3iL,kmul(kmul(Psi3rL,Psi3rL),ToReal(-1.))),kmadd(kmul(Psi1iL,Psi1iL),ToReal(-8.),kmul(kmul(Psi1rL,Psi1rL),ToReal(8.))),kmul(ToReal(2.),kmadd(Psi1rL,kmadd(kmadd(Psi0rL,Psi3iL,kmul(Psi0iL,Psi3rL)),Psi4iL,kmul(Psi4rL,kmadd(Psi0iL,Psi3iL,kmul(Psi0rL,kmul(Psi3rL,ToReal(-1.)))))),kmul(Psi1iL,kmadd(Psi0iL,kmul(Psi3rL,Psi4rL),kmadd(Psi0rL,kmadd(Psi3rL,Psi4iL,kmul(Psi3iL,Psi4rL)),kmul(Psi3iL,kmadd(Psi0iL,kmul(Psi4iL,ToReal(-1.)),kmul(Psi1rL,kmul(Psi3rL,ToReal(16.))))))))))))))))))))))))))))); /* 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(curvIi[index],curvIiL); vec_store_nta_partial(curvIr[index],curvIrL); vec_store_nta_partial(curvJ1[index],curvJ1L); diff --git a/src/WeylScal4_invars_calc_Nth.cc b/src/WeylScal4_invars_calc_Nth.cc index d1a9994..0ccb0f1 100644 --- a/src/WeylScal4_invars_calc_Nth.cc +++ b/src/WeylScal4_invars_calc_Nth.cc @@ -155,7 +155,7 @@ static void WeylScal4_invars_calc_Nth_Body(cGH const * restrict const cctkGH, in CCTK_LOOP3STR(WeylScal4_invars_calc_Nth, 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; @@ -219,7 +219,7 @@ static void WeylScal4_invars_calc_Nth_Body(cGH const * restrict const cctkGH, in kmul(ToReal(-640.),kmadd(Psi1iL,kmul(kmadd(kmadd(Psi3rL,Psi4iL,kmul(Psi3iL,Psi4rL)),kmul(Psi1rL,Psi1rL),kmul(Psi0rL,kmul(Psi3iL,kmul(Psi3rL,Psi3rL)))),ToReal(-12.)),kmadd(kmadd(Psi3iL,kmul(Psi4iL,kmul(Psi1rL,kmul(Psi1rL,Psi1rL))),kmul(Psi0iL,kmul(Psi1iL,kmul(Psi3rL,kmul(Psi3rL,Psi3rL))))),ToReal(-4.),kmadd(kpow(Psi2rL,5.),ToReal(-3.),kmadd(Psi4rL,kmadd(kmul(Psi0rL,Psi0rL),kmul(Psi3iL,Psi3iL),kmadd(kmul(Psi0iL,Psi0iL),kmul(Psi3rL,Psi3rL),kmul(kmadd(kmul(Psi0iL,Psi0iL),kmul(Psi3iL,Psi3iL),kmul(kmul(Psi0rL,Psi0rL),kmul(Psi3rL,Psi3rL))),ToReal(-1.)))),kmadd(Psi0rL,kmadd(kmul(Psi1rL,Psi1rL),kmul(Psi4iL,Psi4iL),kmadd(kmul(Psi1iL,Psi1iL),kmul(Psi4rL,Psi4rL),kmul(kmadd(kmul(Psi1iL,Psi1iL),kmul(Psi4iL,Psi4iL),kmul(kmul(Psi1rL,Psi1rL),kmul(Psi4rL,Psi4rL))),ToReal(-1.)))),kmadd(kmadd(kmadd(Psi3rL,Psi4iL,kmul(Psi3iL,Psi4rL)),kmul(Psi1iL,kmul(Psi1iL,Psi1iL)),kmadd(Psi3rL,kmul(Psi4rL,kmul(Psi1rL,kmul(Psi1rL,Psi1rL))),kmadd(Psi0rL,kmadd(kmadd(Psi0iL,kmul(Psi3iL,Psi3rL),kmul(Psi1iL,kmul(Psi1rL,Psi4iL))),Psi4rL,kmul(Psi1iL,kmul(Psi3iL,kmul(Psi3iL,Psi3iL)))),kmul(Psi1rL,kmadd(Psi0iL,kmul(Psi3iL,kmul(Psi3iL,Psi3iL)),kmul(Psi0rL,kmul(Psi3rL,kmul(Psi3rL,Psi3rL)))))))),ToReal(4.),kmadd(ToReal(-2.),kmadd(Psi4iL,kmadd(Psi3iL,kmul(Psi3rL,kmul(Psi0iL,Psi0iL)),kmul(Psi0iL,kmadd(Psi4rL,kmul(Psi1iL,Psi1iL),kmul(Psi0rL,kmul(Psi3iL,Psi3iL))))),kmul(kmul(Psi2rL,kmul(Psi2rL,Psi2rL)),kmadd(Psi0iL,Psi4iL,kmadd(Psi1rL,kmul(Psi3rL,ToReal(-5.)),kmadd(Psi0rL,kmul(Psi4rL,ToReal(-1.)),kmul(Psi1iL,kmul(Psi3iL,ToReal(5.)))))))),kmadd(ToReal(2.),kmadd(Psi4iL,kmadd(Psi3iL,kmul(Psi3rL,kmul(Psi0rL,Psi0rL)),kmul(Psi0iL,kmadd(Psi4rL,kmul(Psi1rL,Psi1rL),kmul(Psi0rL,kmul(Psi3rL,Psi3rL))))),kmadd(kmul(Psi2iL,kmul(Psi2iL,Psi2iL)),kmadd(Psi0rL,Psi4iL,kmadd(Psi0iL,Psi4rL,kmul(kmadd(Psi1rL,Psi3iL,kmul(Psi1iL,Psi3rL)),ToReal(5.)))),kmul(Psi2iL,kmadd(Psi3rL,kmadd(Psi1iL,kmul(kmul(Psi2rL,Psi2rL),ToReal(-15.)),kmul(Psi3iL,kmul(kmul(Psi1iL,Psi1iL),ToReal(-8.)))),kmadd(Psi4iL,kmadd(Psi2rL,kmul(kmul(Psi1iL,Psi1iL),ToReal(-3.)),kmul(Psi4rL,kmadd(Psi0iL,Psi0iL,kmul(kmul(Psi0rL,Psi0rL),ToReal(-1.))))),kmadd(Psi0iL,kmul(ToReal(-1.),kmadd(Psi1iL,kmadd(Psi3rL,Psi4iL,kmul(Psi3iL,Psi4rL)),kmadd(Psi2rL,kmul(kmul(Psi3rL,Psi3rL),ToReal(-3.)),kmadd(Psi0rL,kmadd(Psi4rL,Psi4rL,kmul(kmul(Psi4iL,Psi4iL),ToReal(-1.))),kmul(kmadd(Psi4rL,kmul(Psi2rL,Psi2rL),kmul(Psi2rL,kmul(Psi3iL,Psi3iL))),ToReal(3.)))))),kmadd(Psi0rL,kmadd(Psi1iL,kmul(Psi3rL,Psi4rL),kmadd(Psi4iL,kmul(kmul(Psi2rL,Psi2rL),ToReal(-3.)),kmul(Psi3iL,kmadd(Psi1iL,kmul(Psi4iL,ToReal(-1.)),kmul(Psi2rL,kmul(Psi3rL,ToReal(6.))))))),kmadd(kmul(Psi1rL,Psi1rL),kmadd(Psi2rL,kmul(Psi4iL,ToReal(3.)),kmul(Psi3iL,kmul(Psi3rL,ToReal(8.)))),kmul(Psi1rL,kmadd(Psi3iL,kmul(kmul(Psi2rL,Psi2rL),ToReal(-15.)),kmadd(Psi4iL,kmadd(Psi0rL,Psi3rL,kmul(Psi0iL,kmul(Psi3iL,ToReal(-1.)))),kmadd(Psi4rL,kmadd(Psi0rL,Psi3iL,kmadd(Psi0iL,Psi3rL,kmul(Psi1iL,kmul(Psi2rL,ToReal(6.))))),kmul(Psi1iL,kmadd(kmul(Psi3iL,Psi3iL),ToReal(-8.),kmul(kmul(Psi3rL,Psi3rL),ToReal(8.)))))))))))))))),kmadd(ToReal(3.),kmadd(kmul(Psi2rL,Psi2rL),kmadd(Psi4rL,kmadd(Psi1iL,Psi1iL,kmul(kmul(Psi1rL,Psi1rL),ToReal(-1.))),kmadd(Psi0rL,kmadd(Psi3iL,Psi3iL,kmul(kmul(Psi3rL,Psi3rL),ToReal(-1.))),kmul(kmadd(Psi0iL,kmul(Psi3iL,Psi3rL),kmul(Psi1iL,kmul(Psi1rL,Psi4iL))),ToReal(2.)))),kmul(kmul(Psi2iL,Psi2iL),kmadd(kmadd(Psi0iL,kmul(Psi3iL,Psi3rL),kmul(Psi1iL,kmul(Psi1rL,Psi4iL))),ToReal(-2.),kmadd(Psi4rL,kmadd(Psi1rL,Psi1rL,kmul(kmul(Psi1iL,Psi1iL),ToReal(-1.))),kmadd(Psi0rL,kmadd(Psi3rL,Psi3rL,kmul(kmul(Psi3iL,Psi3iL),ToReal(-1.))),kmadd(Psi2rL,kmul(ToReal(2.),kmadd(Psi0iL,Psi4iL,kmadd(Psi1rL,kmul(Psi3rL,ToReal(-5.)),kmadd(Psi0rL,kmul(Psi4rL,ToReal(-1.)),kmul(Psi1iL,kmul(Psi3iL,ToReal(5.))))))),kmul(kmul(Psi2rL,kmul(Psi2rL,Psi2rL)),ToReal(10.)))))))),kmadd(Psi1rL,kmadd(Psi0rL,kmul(Psi3rL,kmul(kmul(Psi3iL,Psi3iL),ToReal(-12.))),kmadd(Psi0iL,kmul(Psi1iL,kmul(kmul(Psi4iL,Psi4iL),ToReal(-2.))),kmul(kmul(Psi1iL,Psi1iL),kmadd(Psi3rL,kmul(Psi4rL,ToReal(-12.)),kmul(Psi3iL,kmul(Psi4iL,ToReal(12.))))))),kmadd(Psi0iL,kmadd(Psi1rL,kmul(Psi3iL,kmul(kmul(Psi3rL,Psi3rL),ToReal(-12.))),kmul(Psi1iL,kmadd(Psi1rL,kmul(kmul(Psi4rL,Psi4rL),ToReal(2.)),kmul(Psi3rL,kmul(kmul(Psi3iL,Psi3iL),ToReal(12.)))))),kmul(Psi2rL,kmadd(kmul(Psi0iL,Psi0iL),kmul(Psi4iL,Psi4iL),kmadd(kmul(Psi0rL,Psi0rL),kmul(Psi4rL,Psi4rL),kmadd(kmul(kmul(Psi2iL,Psi2iL),kmul(Psi2iL,Psi2iL)),ToReal(-15.),kmadd(Psi0iL,kmul(Psi0rL,kmul(Psi4iL,kmul(Psi4rL,ToReal(-4.)))),kmadd(kmadd(kmul(Psi0rL,Psi0rL),kmul(Psi4iL,Psi4iL),kmul(kmul(Psi0iL,Psi0iL),kmul(Psi4rL,Psi4rL))),ToReal(-1.),kmadd(kmadd(Psi3iL,Psi3iL,kmul(kmul(Psi3rL,Psi3rL),ToReal(-1.))),kmadd(kmul(Psi1iL,Psi1iL),ToReal(-8.),kmul(kmul(Psi1rL,Psi1rL),ToReal(8.))),kmul(ToReal(2.),kmadd(Psi1rL,kmadd(kmadd(Psi0rL,Psi3iL,kmul(Psi0iL,Psi3rL)),Psi4iL,kmul(Psi4rL,kmadd(Psi0iL,Psi3iL,kmul(Psi0rL,kmul(Psi3rL,ToReal(-1.)))))),kmul(Psi1iL,kmadd(Psi0iL,kmul(Psi3rL,Psi4rL),kmadd(Psi0rL,kmadd(Psi3rL,Psi4iL,kmul(Psi3iL,Psi4rL)),kmul(Psi3iL,kmadd(Psi0iL,kmul(Psi4iL,ToReal(-1.)),kmul(Psi1rL,kmul(Psi3rL,ToReal(16.))))))))))))))))))))))))))))); /* 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(curvIi[index],curvIiL); vec_store_nta_partial(curvIr[index],curvIrL); vec_store_nta_partial(curvJ1[index],curvJ1L); diff --git a/src/WeylScal4_psi4_calc_2nd.cc b/src/WeylScal4_psi4_calc_2nd.cc index 1a0084b..ddf5a89 100644 --- a/src/WeylScal4_psi4_calc_2nd.cc +++ b/src/WeylScal4_psi4_calc_2nd.cc @@ -170,7 +170,7 @@ static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int CCTK_LOOP3STR(WeylScal4_psi4_calc_2nd, 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; @@ -1313,7 +1313,7 @@ static void WeylScal4_psi4_calc_2nd_Body(cGH const * restrict const cctkGH, int kmadd(kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(n1,kmadd(n2,R4p1213,kmul(n3,R4p1313)),kmadd(n2,kmul(n3,R4p1323),kmul(R4p1223,kmul(n2,n2)))),kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(n1,kmadd(n2,R4p1212,kmul(n3,R4p1213)),kmul(n3,kmul(kmadd(n2,R4p1223,kmul(n3,R4p1323)),ToReal(-1.)))),kmul(nn,kmadd(kmadd(im1,kmul(rm1,kmadd(n1,Ro111,kmadd(n2,Ro121,kmul(n3,Ro131)))),kmadd(im2,kmul(rm2,kmadd(n1,Ro212,kmadd(n2,Ro222,kmul(n3,Ro232)))),kmul(im3,kmul(rm3,kmadd(n1,Ro313,kmadd(n2,Ro323,kmul(n3,Ro333))))))),ToReal(-2.),kmul(kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(n1,kadd(Ro112,Ro211),kmadd(n2,kadd(Ro122,Ro221),kmul(n3,kadd(Ro132,Ro231)))),kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(n1,kadd(Ro113,Ro311),kmadd(n2,kadd(Ro123,Ro321),kmul(n3,kadd(Ro133,Ro331)))),kmul(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(n1,kadd(Ro213,Ro312),kmadd(n2,kadd(Ro223,Ro322),kmul(n3,kadd(Ro233,Ro332))))))),ToReal(-1.)))))),ToReal(2.),kmadd(ToReal(-2.),kmadd(im2,kmul(rm2,kmadd(R4p1212,kmul(n1,n1),kmadd(R4p2323,kmul(n3,n3),kmul(n1,kmul(n3,kmul(R4p1223,ToReal(-2.))))))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(R4p1213,kmul(n1,n1),kmadd(n2,kmul(n3,kmul(R4p2323,ToReal(-1.))),kmul(n1,kmadd(n2,R4p1223,kmul(n3,kmul(R4p1323,ToReal(-1.))))))),kmadd(im1,kmul(rm1,kmadd(R4p1212,kmul(n2,n2),kmadd(R4p1313,kmul(n3,n3),kmul(n2,kmul(n3,kmul(R4p1213,ToReal(2.))))))),kmul(im3,kmul(rm3,kmadd(R4p1313,kmul(n1,n1),kmadd(R4p2323,kmul(n2,n2),kmul(n1,kmul(n2,kmul(R4p1323,ToReal(2.))))))))))),kmul(kmul(nn,nn),kmul(ToReal(-1.),kmadd(im1,kmadd(rm2,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo13,Rojo31),kmul(rm1,kmul(Rojo11,ToReal(2.))))),kmadd(im2,kmadd(rm1,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo23,Rojo32),kmul(rm2,kmul(Rojo22,ToReal(2.))))),kmul(im3,kmadd(rm1,kadd(Rojo13,Rojo31),kmadd(rm2,kadd(Rojo23,Rojo32),kmul(rm3,kmul(Rojo33,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(Psi4i[index],Psi4iL); vec_store_nta_partial(Psi4r[index],Psi4rL); } diff --git a/src/WeylScal4_psi4_calc_4th.cc b/src/WeylScal4_psi4_calc_4th.cc index 6d667ec..f4e387b 100644 --- a/src/WeylScal4_psi4_calc_4th.cc +++ b/src/WeylScal4_psi4_calc_4th.cc @@ -170,7 +170,7 @@ static void WeylScal4_psi4_calc_4th_Body(cGH const * restrict const cctkGH, int CCTK_LOOP3STR(WeylScal4_psi4_calc_4th, 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; @@ -1313,7 +1313,7 @@ static void WeylScal4_psi4_calc_4th_Body(cGH const * restrict const cctkGH, int kmadd(kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(n1,kmadd(n2,R4p1213,kmul(n3,R4p1313)),kmadd(n2,kmul(n3,R4p1323),kmul(R4p1223,kmul(n2,n2)))),kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(n1,kmadd(n2,R4p1212,kmul(n3,R4p1213)),kmul(n3,kmul(kmadd(n2,R4p1223,kmul(n3,R4p1323)),ToReal(-1.)))),kmul(nn,kmadd(kmadd(im1,kmul(rm1,kmadd(n1,Ro111,kmadd(n2,Ro121,kmul(n3,Ro131)))),kmadd(im2,kmul(rm2,kmadd(n1,Ro212,kmadd(n2,Ro222,kmul(n3,Ro232)))),kmul(im3,kmul(rm3,kmadd(n1,Ro313,kmadd(n2,Ro323,kmul(n3,Ro333))))))),ToReal(-2.),kmul(kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(n1,kadd(Ro112,Ro211),kmadd(n2,kadd(Ro122,Ro221),kmul(n3,kadd(Ro132,Ro231)))),kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(n1,kadd(Ro113,Ro311),kmadd(n2,kadd(Ro123,Ro321),kmul(n3,kadd(Ro133,Ro331)))),kmul(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(n1,kadd(Ro213,Ro312),kmadd(n2,kadd(Ro223,Ro322),kmul(n3,kadd(Ro233,Ro332))))))),ToReal(-1.)))))),ToReal(2.),kmadd(ToReal(-2.),kmadd(im2,kmul(rm2,kmadd(R4p1212,kmul(n1,n1),kmadd(R4p2323,kmul(n3,n3),kmul(n1,kmul(n3,kmul(R4p1223,ToReal(-2.))))))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(R4p1213,kmul(n1,n1),kmadd(n2,kmul(n3,kmul(R4p2323,ToReal(-1.))),kmul(n1,kmadd(n2,R4p1223,kmul(n3,kmul(R4p1323,ToReal(-1.))))))),kmadd(im1,kmul(rm1,kmadd(R4p1212,kmul(n2,n2),kmadd(R4p1313,kmul(n3,n3),kmul(n2,kmul(n3,kmul(R4p1213,ToReal(2.))))))),kmul(im3,kmul(rm3,kmadd(R4p1313,kmul(n1,n1),kmadd(R4p2323,kmul(n2,n2),kmul(n1,kmul(n2,kmul(R4p1323,ToReal(2.))))))))))),kmul(kmul(nn,nn),kmul(ToReal(-1.),kmadd(im1,kmadd(rm2,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo13,Rojo31),kmul(rm1,kmul(Rojo11,ToReal(2.))))),kmadd(im2,kmadd(rm1,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo23,Rojo32),kmul(rm2,kmul(Rojo22,ToReal(2.))))),kmul(im3,kmadd(rm1,kadd(Rojo13,Rojo31),kmadd(rm2,kadd(Rojo23,Rojo32),kmul(rm3,kmul(Rojo33,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(Psi4i[index],Psi4iL); vec_store_nta_partial(Psi4r[index],Psi4rL); } diff --git a/src/WeylScal4_psi4_calc_Nth.cc b/src/WeylScal4_psi4_calc_Nth.cc index 9d98af2..cdbd85c 100644 --- a/src/WeylScal4_psi4_calc_Nth.cc +++ b/src/WeylScal4_psi4_calc_Nth.cc @@ -170,7 +170,7 @@ static void WeylScal4_psi4_calc_Nth_Body(cGH const * restrict const cctkGH, int CCTK_LOOP3STR(WeylScal4_psi4_calc_Nth, 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; @@ -1313,7 +1313,7 @@ static void WeylScal4_psi4_calc_Nth_Body(cGH const * restrict const cctkGH, int kmadd(kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(n1,kmadd(n2,R4p1213,kmul(n3,R4p1313)),kmadd(n2,kmul(n3,R4p1323),kmul(R4p1223,kmul(n2,n2)))),kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(n1,kmadd(n2,R4p1212,kmul(n3,R4p1213)),kmul(n3,kmul(kmadd(n2,R4p1223,kmul(n3,R4p1323)),ToReal(-1.)))),kmul(nn,kmadd(kmadd(im1,kmul(rm1,kmadd(n1,Ro111,kmadd(n2,Ro121,kmul(n3,Ro131)))),kmadd(im2,kmul(rm2,kmadd(n1,Ro212,kmadd(n2,Ro222,kmul(n3,Ro232)))),kmul(im3,kmul(rm3,kmadd(n1,Ro313,kmadd(n2,Ro323,kmul(n3,Ro333))))))),ToReal(-2.),kmul(kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(n1,kadd(Ro112,Ro211),kmadd(n2,kadd(Ro122,Ro221),kmul(n3,kadd(Ro132,Ro231)))),kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(n1,kadd(Ro113,Ro311),kmadd(n2,kadd(Ro123,Ro321),kmul(n3,kadd(Ro133,Ro331)))),kmul(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(n1,kadd(Ro213,Ro312),kmadd(n2,kadd(Ro223,Ro322),kmul(n3,kadd(Ro233,Ro332))))))),ToReal(-1.)))))),ToReal(2.),kmadd(ToReal(-2.),kmadd(im2,kmul(rm2,kmadd(R4p1212,kmul(n1,n1),kmadd(R4p2323,kmul(n3,n3),kmul(n1,kmul(n3,kmul(R4p1223,ToReal(-2.))))))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(R4p1213,kmul(n1,n1),kmadd(n2,kmul(n3,kmul(R4p2323,ToReal(-1.))),kmul(n1,kmadd(n2,R4p1223,kmul(n3,kmul(R4p1323,ToReal(-1.))))))),kmadd(im1,kmul(rm1,kmadd(R4p1212,kmul(n2,n2),kmadd(R4p1313,kmul(n3,n3),kmul(n2,kmul(n3,kmul(R4p1213,ToReal(2.))))))),kmul(im3,kmul(rm3,kmadd(R4p1313,kmul(n1,n1),kmadd(R4p2323,kmul(n2,n2),kmul(n1,kmul(n2,kmul(R4p1323,ToReal(2.))))))))))),kmul(kmul(nn,nn),kmul(ToReal(-1.),kmadd(im1,kmadd(rm2,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo13,Rojo31),kmul(rm1,kmul(Rojo11,ToReal(2.))))),kmadd(im2,kmadd(rm1,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo23,Rojo32),kmul(rm2,kmul(Rojo22,ToReal(2.))))),kmul(im3,kmadd(rm1,kadd(Rojo13,Rojo31),kmadd(rm2,kadd(Rojo23,Rojo32),kmul(rm3,kmul(Rojo33,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(Psi4i[index],Psi4iL); vec_store_nta_partial(Psi4r[index],Psi4rL); } diff --git a/src/WeylScal4_psis_calc_2nd.cc b/src/WeylScal4_psis_calc_2nd.cc index 4fc2cff..a15e554 100644 --- a/src/WeylScal4_psis_calc_2nd.cc +++ b/src/WeylScal4_psis_calc_2nd.cc @@ -194,7 +194,7 @@ static void WeylScal4_psis_calc_2nd_Body(cGH const * restrict const cctkGH, int CCTK_LOOP3STR(WeylScal4_psis_calc_2nd, 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; @@ -1367,7 +1367,7 @@ static void WeylScal4_psis_calc_2nd_Body(cGH const * restrict const cctkGH, int kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmul(kmadd(ltet1,kmadd(ltet2,R4p1213,kmul(ltet3,R4p1313)),kmadd(ltet2,kmul(ltet3,R4p1323),kmul(R4p1223,kmul(ltet2,ltet2)))),ToReal(-2.)),kmadd(ToReal(2.),kmadd(im2,kmul(rm2,kmadd(R4p1212,kmul(ltet1,ltet1),kmadd(R4p2323,kmul(ltet3,ltet3),kmul(ltet1,kmul(ltet3,kmul(R4p1223,ToReal(-2.))))))),kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(ltet3,kmadd(ltet2,R4p1223,kmul(ltet3,R4p1323)),kmul(ltet1,kmul(kmadd(ltet2,R4p1212,kmul(ltet3,R4p1213)),ToReal(-1.)))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(R4p1213,kmul(ltet1,ltet1),kmadd(ltet2,kmul(ltet3,kmul(R4p2323,ToReal(-1.))),kmul(ltet1,kmadd(ltet2,R4p1223,kmul(ltet3,kmul(R4p1323,ToReal(-1.))))))),kmadd(im1,kmul(rm1,kmadd(R4p1212,kmul(ltet2,ltet2),kmadd(R4p1313,kmul(ltet3,ltet3),kmul(ltet2,kmul(ltet3,kmul(R4p1213,ToReal(2.))))))),kmadd(im3,kmul(rm3,kmadd(R4p1313,kmul(ltet1,ltet1),kmadd(R4p2323,kmul(ltet2,ltet2),kmul(ltet1,kmul(ltet2,kmul(R4p1323,ToReal(2.))))))),kmul(nn,kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(ltet1,kadd(Ro112,Ro211),kmadd(ltet2,kadd(Ro122,Ro221),kmul(ltet3,kadd(Ro132,Ro231)))),kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(ltet1,kadd(Ro113,Ro311),kmadd(ltet2,kadd(Ro123,Ro321),kmul(ltet3,kadd(Ro133,Ro331)))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(ltet1,kadd(Ro213,Ro312),kmadd(ltet2,kadd(Ro223,Ro322),kmul(ltet3,kadd(Ro233,Ro332)))),kmul(kmadd(im1,kmul(rm1,kmadd(ltet1,Ro111,kmadd(ltet2,Ro121,kmul(ltet3,Ro131)))),kmadd(im2,kmul(rm2,kmadd(ltet1,Ro212,kmadd(ltet2,Ro222,kmul(ltet3,Ro232)))),kmul(im3,kmul(rm3,kmadd(ltet1,Ro313,kmadd(ltet2,Ro323,kmul(ltet3,Ro333))))))),ToReal(2.))))))))))),kmul(kmul(nn,nn),kmadd(im1,kmadd(rm2,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo13,Rojo31),kmul(rm1,kmul(Rojo11,ToReal(2.))))),kmadd(im2,kmadd(rm1,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo23,Rojo32),kmul(rm2,kmul(Rojo22,ToReal(2.))))),kmul(im3,kmadd(rm1,kadd(Rojo13,Rojo31),kmadd(rm2,kadd(Rojo23,Rojo32),kmul(rm3,kmul(Rojo33,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(Psi0i[index],Psi0iL); vec_store_nta_partial(Psi0r[index],Psi0rL); vec_store_nta_partial(Psi1i[index],Psi1iL); diff --git a/src/WeylScal4_psis_calc_4th.cc b/src/WeylScal4_psis_calc_4th.cc index a2b0209..642b79c 100644 --- a/src/WeylScal4_psis_calc_4th.cc +++ b/src/WeylScal4_psis_calc_4th.cc @@ -194,7 +194,7 @@ static void WeylScal4_psis_calc_4th_Body(cGH const * restrict const cctkGH, int CCTK_LOOP3STR(WeylScal4_psis_calc_4th, 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; @@ -1367,7 +1367,7 @@ static void WeylScal4_psis_calc_4th_Body(cGH const * restrict const cctkGH, int kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmul(kmadd(ltet1,kmadd(ltet2,R4p1213,kmul(ltet3,R4p1313)),kmadd(ltet2,kmul(ltet3,R4p1323),kmul(R4p1223,kmul(ltet2,ltet2)))),ToReal(-2.)),kmadd(ToReal(2.),kmadd(im2,kmul(rm2,kmadd(R4p1212,kmul(ltet1,ltet1),kmadd(R4p2323,kmul(ltet3,ltet3),kmul(ltet1,kmul(ltet3,kmul(R4p1223,ToReal(-2.))))))),kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(ltet3,kmadd(ltet2,R4p1223,kmul(ltet3,R4p1323)),kmul(ltet1,kmul(kmadd(ltet2,R4p1212,kmul(ltet3,R4p1213)),ToReal(-1.)))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(R4p1213,kmul(ltet1,ltet1),kmadd(ltet2,kmul(ltet3,kmul(R4p2323,ToReal(-1.))),kmul(ltet1,kmadd(ltet2,R4p1223,kmul(ltet3,kmul(R4p1323,ToReal(-1.))))))),kmadd(im1,kmul(rm1,kmadd(R4p1212,kmul(ltet2,ltet2),kmadd(R4p1313,kmul(ltet3,ltet3),kmul(ltet2,kmul(ltet3,kmul(R4p1213,ToReal(2.))))))),kmadd(im3,kmul(rm3,kmadd(R4p1313,kmul(ltet1,ltet1),kmadd(R4p2323,kmul(ltet2,ltet2),kmul(ltet1,kmul(ltet2,kmul(R4p1323,ToReal(2.))))))),kmul(nn,kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(ltet1,kadd(Ro112,Ro211),kmadd(ltet2,kadd(Ro122,Ro221),kmul(ltet3,kadd(Ro132,Ro231)))),kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(ltet1,kadd(Ro113,Ro311),kmadd(ltet2,kadd(Ro123,Ro321),kmul(ltet3,kadd(Ro133,Ro331)))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(ltet1,kadd(Ro213,Ro312),kmadd(ltet2,kadd(Ro223,Ro322),kmul(ltet3,kadd(Ro233,Ro332)))),kmul(kmadd(im1,kmul(rm1,kmadd(ltet1,Ro111,kmadd(ltet2,Ro121,kmul(ltet3,Ro131)))),kmadd(im2,kmul(rm2,kmadd(ltet1,Ro212,kmadd(ltet2,Ro222,kmul(ltet3,Ro232)))),kmul(im3,kmul(rm3,kmadd(ltet1,Ro313,kmadd(ltet2,Ro323,kmul(ltet3,Ro333))))))),ToReal(2.))))))))))),kmul(kmul(nn,nn),kmadd(im1,kmadd(rm2,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo13,Rojo31),kmul(rm1,kmul(Rojo11,ToReal(2.))))),kmadd(im2,kmadd(rm1,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo23,Rojo32),kmul(rm2,kmul(Rojo22,ToReal(2.))))),kmul(im3,kmadd(rm1,kadd(Rojo13,Rojo31),kmadd(rm2,kadd(Rojo23,Rojo32),kmul(rm3,kmul(Rojo33,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(Psi0i[index],Psi0iL); vec_store_nta_partial(Psi0r[index],Psi0rL); vec_store_nta_partial(Psi1i[index],Psi1iL); diff --git a/src/WeylScal4_psis_calc_Nth.cc b/src/WeylScal4_psis_calc_Nth.cc index a136d5e..bd68c9d 100644 --- a/src/WeylScal4_psis_calc_Nth.cc +++ b/src/WeylScal4_psis_calc_Nth.cc @@ -194,7 +194,7 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int CCTK_LOOP3STR(WeylScal4_psis_calc_Nth, 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; @@ -1367,7 +1367,7 @@ static void WeylScal4_psis_calc_Nth_Body(cGH const * restrict const cctkGH, int kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmul(kmadd(ltet1,kmadd(ltet2,R4p1213,kmul(ltet3,R4p1313)),kmadd(ltet2,kmul(ltet3,R4p1323),kmul(R4p1223,kmul(ltet2,ltet2)))),ToReal(-2.)),kmadd(ToReal(2.),kmadd(im2,kmul(rm2,kmadd(R4p1212,kmul(ltet1,ltet1),kmadd(R4p2323,kmul(ltet3,ltet3),kmul(ltet1,kmul(ltet3,kmul(R4p1223,ToReal(-2.))))))),kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(ltet3,kmadd(ltet2,R4p1223,kmul(ltet3,R4p1323)),kmul(ltet1,kmul(kmadd(ltet2,R4p1212,kmul(ltet3,R4p1213)),ToReal(-1.)))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(R4p1213,kmul(ltet1,ltet1),kmadd(ltet2,kmul(ltet3,kmul(R4p2323,ToReal(-1.))),kmul(ltet1,kmadd(ltet2,R4p1223,kmul(ltet3,kmul(R4p1323,ToReal(-1.))))))),kmadd(im1,kmul(rm1,kmadd(R4p1212,kmul(ltet2,ltet2),kmadd(R4p1313,kmul(ltet3,ltet3),kmul(ltet2,kmul(ltet3,kmul(R4p1213,ToReal(2.))))))),kmadd(im3,kmul(rm3,kmadd(R4p1313,kmul(ltet1,ltet1),kmadd(R4p2323,kmul(ltet2,ltet2),kmul(ltet1,kmul(ltet2,kmul(R4p1323,ToReal(2.))))))),kmul(nn,kmadd(kmadd(im2,rm1,kmul(im1,rm2)),kmadd(ltet1,kadd(Ro112,Ro211),kmadd(ltet2,kadd(Ro122,Ro221),kmul(ltet3,kadd(Ro132,Ro231)))),kmadd(kmadd(im3,rm1,kmul(im1,rm3)),kmadd(ltet1,kadd(Ro113,Ro311),kmadd(ltet2,kadd(Ro123,Ro321),kmul(ltet3,kadd(Ro133,Ro331)))),kmadd(kmadd(im3,rm2,kmul(im2,rm3)),kmadd(ltet1,kadd(Ro213,Ro312),kmadd(ltet2,kadd(Ro223,Ro322),kmul(ltet3,kadd(Ro233,Ro332)))),kmul(kmadd(im1,kmul(rm1,kmadd(ltet1,Ro111,kmadd(ltet2,Ro121,kmul(ltet3,Ro131)))),kmadd(im2,kmul(rm2,kmadd(ltet1,Ro212,kmadd(ltet2,Ro222,kmul(ltet3,Ro232)))),kmul(im3,kmul(rm3,kmadd(ltet1,Ro313,kmadd(ltet2,Ro323,kmul(ltet3,Ro333))))))),ToReal(2.))))))))))),kmul(kmul(nn,nn),kmadd(im1,kmadd(rm2,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo13,Rojo31),kmul(rm1,kmul(Rojo11,ToReal(2.))))),kmadd(im2,kmadd(rm1,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo23,Rojo32),kmul(rm2,kmul(Rojo22,ToReal(2.))))),kmul(im3,kmadd(rm1,kadd(Rojo13,Rojo31),kmadd(rm2,kadd(Rojo23,Rojo32),kmul(rm3,kmul(Rojo33,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(Psi0i[index],Psi0iL); vec_store_nta_partial(Psi0r[index],Psi0rL); vec_store_nta_partial(Psi1i[index],Psi1iL); -- cgit v1.2.3