diff options
Diffstat (limited to 'src/WeylScal4_psi4_calc_2nd.cc')
| -rw-r--r-- | src/WeylScal4_psi4_calc_2nd.cc | 4 |
1 files changed, 2 insertions, 2 deletions
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); } |