diff options
Diffstat (limited to 'src/psi4_calc_4th.cc')
| -rw-r--r-- | src/psi4_calc_4th.cc | 586 |
1 files changed, 11 insertions, 575 deletions
diff --git a/src/psi4_calc_4th.cc b/src/psi4_calc_4th.cc index 48b61d6..8f307ec 100644 --- a/src/psi4_calc_4th.cc +++ b/src/psi4_calc_4th.cc @@ -606,23 +606,14 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC JacPDstandard4th11gxy; - CCTK_REAL_VEC JacPDstandard4th11gxz; CCTK_REAL_VEC JacPDstandard4th11gyy; CCTK_REAL_VEC JacPDstandard4th11gyz; CCTK_REAL_VEC JacPDstandard4th11gzz; - CCTK_REAL_VEC JacPDstandard4th12gxx; CCTK_REAL_VEC JacPDstandard4th12gxy; CCTK_REAL_VEC JacPDstandard4th12gxz; - CCTK_REAL_VEC JacPDstandard4th12gyy; CCTK_REAL_VEC JacPDstandard4th12gyz; CCTK_REAL_VEC JacPDstandard4th12gzz; - CCTK_REAL_VEC JacPDstandard4th13gxx; - CCTK_REAL_VEC JacPDstandard4th13gxy; CCTK_REAL_VEC JacPDstandard4th13gxz; - CCTK_REAL_VEC JacPDstandard4th13gyy; - CCTK_REAL_VEC JacPDstandard4th13gyz; - CCTK_REAL_VEC JacPDstandard4th13gzz; CCTK_REAL_VEC JacPDstandard4th1gxx; CCTK_REAL_VEC JacPDstandard4th1gxy; CCTK_REAL_VEC JacPDstandard4th1gxz; @@ -634,23 +625,14 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC JacPDstandard4th1kyy; CCTK_REAL_VEC JacPDstandard4th1kyz; CCTK_REAL_VEC JacPDstandard4th1kzz; - CCTK_REAL_VEC JacPDstandard4th21gxx; CCTK_REAL_VEC JacPDstandard4th21gxy; - CCTK_REAL_VEC JacPDstandard4th21gxz; - CCTK_REAL_VEC JacPDstandard4th21gyy; - CCTK_REAL_VEC JacPDstandard4th21gyz; - CCTK_REAL_VEC JacPDstandard4th21gzz; CCTK_REAL_VEC JacPDstandard4th22gxx; - CCTK_REAL_VEC JacPDstandard4th22gxy; CCTK_REAL_VEC JacPDstandard4th22gxz; - CCTK_REAL_VEC JacPDstandard4th22gyz; CCTK_REAL_VEC JacPDstandard4th22gzz; CCTK_REAL_VEC JacPDstandard4th23gxx; CCTK_REAL_VEC JacPDstandard4th23gxy; CCTK_REAL_VEC JacPDstandard4th23gxz; - CCTK_REAL_VEC JacPDstandard4th23gyy; CCTK_REAL_VEC JacPDstandard4th23gyz; - CCTK_REAL_VEC JacPDstandard4th23gzz; CCTK_REAL_VEC JacPDstandard4th2gxx; CCTK_REAL_VEC JacPDstandard4th2gxy; CCTK_REAL_VEC JacPDstandard4th2gxz; @@ -662,23 +644,14 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC JacPDstandard4th2kxz; CCTK_REAL_VEC JacPDstandard4th2kyz; CCTK_REAL_VEC JacPDstandard4th2kzz; - CCTK_REAL_VEC JacPDstandard4th31gxx; CCTK_REAL_VEC JacPDstandard4th31gxy; CCTK_REAL_VEC JacPDstandard4th31gxz; CCTK_REAL_VEC JacPDstandard4th31gyy; CCTK_REAL_VEC JacPDstandard4th31gyz; - CCTK_REAL_VEC JacPDstandard4th31gzz; - CCTK_REAL_VEC JacPDstandard4th32gxx; - CCTK_REAL_VEC JacPDstandard4th32gxy; - CCTK_REAL_VEC JacPDstandard4th32gxz; - CCTK_REAL_VEC JacPDstandard4th32gyy; CCTK_REAL_VEC JacPDstandard4th32gyz; - CCTK_REAL_VEC JacPDstandard4th32gzz; CCTK_REAL_VEC JacPDstandard4th33gxx; CCTK_REAL_VEC JacPDstandard4th33gxy; - CCTK_REAL_VEC JacPDstandard4th33gxz; CCTK_REAL_VEC JacPDstandard4th33gyy; - CCTK_REAL_VEC JacPDstandard4th33gyz; CCTK_REAL_VEC JacPDstandard4th3gxx; CCTK_REAL_VEC JacPDstandard4th3gxy; CCTK_REAL_VEC JacPDstandard4th3gxz; @@ -792,12 +765,6 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, JacPDstandard4th3kyz = kmadd(J13L,PDstandard4th1kyz,kmadd(J23L,PDstandard4th2kyz,kmul(J33L,PDstandard4th3kyz))); - JacPDstandard4th11gxy = - kmadd(dJ111L,PDstandard4th1gxy,kmadd(dJ211L,PDstandard4th2gxy,kmadd(dJ311L,PDstandard4th3gxy,kmadd(PDstandard4th11gxy,kmul(J11L,J11L),kmadd(PDstandard4th22gxy,kmul(J21L,J21L),kmadd(PDstandard4th33gxy,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gxy,kmul(J31L,PDstandard4th13gxy)),kmul(J21L,kmul(J31L,PDstandard4th23gxy))),ToReal(2.)))))))); - - JacPDstandard4th11gxz = - kmadd(dJ111L,PDstandard4th1gxz,kmadd(dJ211L,PDstandard4th2gxz,kmadd(dJ311L,PDstandard4th3gxz,kmadd(PDstandard4th11gxz,kmul(J11L,J11L),kmadd(PDstandard4th22gxz,kmul(J21L,J21L),kmadd(PDstandard4th33gxz,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gxz,kmul(J31L,PDstandard4th13gxz)),kmul(J21L,kmul(J31L,PDstandard4th23gxz))),ToReal(2.)))))))); - JacPDstandard4th11gyy = kmadd(dJ111L,PDstandard4th1gyy,kmadd(dJ211L,PDstandard4th2gyy,kmadd(dJ311L,PDstandard4th3gyy,kmadd(PDstandard4th11gyy,kmul(J11L,J11L),kmadd(PDstandard4th22gyy,kmul(J21L,J21L),kmadd(PDstandard4th33gyy,kmul(J31L,J31L),kmul(kmadd(J11L,kmadd(J21L,PDstandard4th12gyy,kmul(J31L,PDstandard4th13gyy)),kmul(J21L,kmul(J31L,PDstandard4th23gyy))),ToReal(2.)))))))); @@ -810,15 +777,9 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, JacPDstandard4th22gxx = kmadd(dJ122L,PDstandard4th1gxx,kmadd(dJ222L,PDstandard4th2gxx,kmadd(dJ322L,PDstandard4th3gxx,kmadd(PDstandard4th11gxx,kmul(J12L,J12L),kmadd(PDstandard4th22gxx,kmul(J22L,J22L),kmadd(PDstandard4th33gxx,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gxx,kmul(J32L,PDstandard4th13gxx)),kmul(J22L,kmul(J32L,PDstandard4th23gxx))),ToReal(2.)))))))); - JacPDstandard4th22gxy = - kmadd(dJ122L,PDstandard4th1gxy,kmadd(dJ222L,PDstandard4th2gxy,kmadd(dJ322L,PDstandard4th3gxy,kmadd(PDstandard4th11gxy,kmul(J12L,J12L),kmadd(PDstandard4th22gxy,kmul(J22L,J22L),kmadd(PDstandard4th33gxy,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gxy,kmul(J32L,PDstandard4th13gxy)),kmul(J22L,kmul(J32L,PDstandard4th23gxy))),ToReal(2.)))))))); - JacPDstandard4th22gxz = kmadd(dJ122L,PDstandard4th1gxz,kmadd(dJ222L,PDstandard4th2gxz,kmadd(dJ322L,PDstandard4th3gxz,kmadd(PDstandard4th11gxz,kmul(J12L,J12L),kmadd(PDstandard4th22gxz,kmul(J22L,J22L),kmadd(PDstandard4th33gxz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gxz,kmul(J32L,PDstandard4th13gxz)),kmul(J22L,kmul(J32L,PDstandard4th23gxz))),ToReal(2.)))))))); - JacPDstandard4th22gyz = - kmadd(dJ122L,PDstandard4th1gyz,kmadd(dJ222L,PDstandard4th2gyz,kmadd(dJ322L,PDstandard4th3gyz,kmadd(PDstandard4th11gyz,kmul(J12L,J12L),kmadd(PDstandard4th22gyz,kmul(J22L,J22L),kmadd(PDstandard4th33gyz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gyz,kmul(J32L,PDstandard4th13gyz)),kmul(J22L,kmul(J32L,PDstandard4th23gyz))),ToReal(2.)))))))); - JacPDstandard4th22gzz = kmadd(dJ122L,PDstandard4th1gzz,kmadd(dJ222L,PDstandard4th2gzz,kmadd(dJ322L,PDstandard4th3gzz,kmadd(PDstandard4th11gzz,kmul(J12L,J12L),kmadd(PDstandard4th22gzz,kmul(J22L,J22L),kmadd(PDstandard4th33gzz,kmul(J32L,J32L),kmul(kmadd(J12L,kmadd(J22L,PDstandard4th12gzz,kmul(J32L,PDstandard4th13gzz)),kmul(J22L,kmul(J32L,PDstandard4th23gzz))),ToReal(2.)))))))); @@ -828,69 +789,27 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, JacPDstandard4th33gxy = kmadd(dJ133L,PDstandard4th1gxy,kmadd(dJ233L,PDstandard4th2gxy,kmadd(dJ333L,PDstandard4th3gxy,kmadd(PDstandard4th11gxy,kmul(J13L,J13L),kmadd(PDstandard4th22gxy,kmul(J23L,J23L),kmadd(PDstandard4th33gxy,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gxy,kmul(J33L,PDstandard4th13gxy)),kmul(J23L,kmul(J33L,PDstandard4th23gxy))),ToReal(2.)))))))); - JacPDstandard4th33gxz = - kmadd(dJ133L,PDstandard4th1gxz,kmadd(dJ233L,PDstandard4th2gxz,kmadd(dJ333L,PDstandard4th3gxz,kmadd(PDstandard4th11gxz,kmul(J13L,J13L),kmadd(PDstandard4th22gxz,kmul(J23L,J23L),kmadd(PDstandard4th33gxz,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gxz,kmul(J33L,PDstandard4th13gxz)),kmul(J23L,kmul(J33L,PDstandard4th23gxz))),ToReal(2.)))))))); - JacPDstandard4th33gyy = kmadd(dJ133L,PDstandard4th1gyy,kmadd(dJ233L,PDstandard4th2gyy,kmadd(dJ333L,PDstandard4th3gyy,kmadd(PDstandard4th11gyy,kmul(J13L,J13L),kmadd(PDstandard4th22gyy,kmul(J23L,J23L),kmadd(PDstandard4th33gyy,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gyy,kmul(J33L,PDstandard4th13gyy)),kmul(J23L,kmul(J33L,PDstandard4th23gyy))),ToReal(2.)))))))); - JacPDstandard4th33gyz = - kmadd(dJ133L,PDstandard4th1gyz,kmadd(dJ233L,PDstandard4th2gyz,kmadd(dJ333L,PDstandard4th3gyz,kmadd(PDstandard4th11gyz,kmul(J13L,J13L),kmadd(PDstandard4th22gyz,kmul(J23L,J23L),kmadd(PDstandard4th33gyz,kmul(J33L,J33L),kmul(kmadd(J13L,kmadd(J23L,PDstandard4th12gyz,kmul(J33L,PDstandard4th13gyz)),kmul(J23L,kmul(J33L,PDstandard4th23gyz))),ToReal(2.)))))))); - - JacPDstandard4th12gxx = - kmadd(J12L,kmadd(J11L,PDstandard4th11gxx,kmadd(J21L,PDstandard4th12gxx,kmul(J31L,PDstandard4th13gxx))),kmadd(J11L,kmadd(J22L,PDstandard4th12gxx,kmul(J32L,PDstandard4th13gxx)),kmadd(dJ112L,PDstandard4th1gxx,kmadd(J22L,kmadd(J21L,PDstandard4th22gxx,kmul(J31L,PDstandard4th23gxx)),kmadd(dJ212L,PDstandard4th2gxx,kmadd(J32L,kmadd(J21L,PDstandard4th23gxx,kmul(J31L,PDstandard4th33gxx)),kmul(dJ312L,PDstandard4th3gxx))))))); - JacPDstandard4th12gxy = kmadd(J12L,kmadd(J11L,PDstandard4th11gxy,kmadd(J21L,PDstandard4th12gxy,kmul(J31L,PDstandard4th13gxy))),kmadd(J11L,kmadd(J22L,PDstandard4th12gxy,kmul(J32L,PDstandard4th13gxy)),kmadd(dJ112L,PDstandard4th1gxy,kmadd(J22L,kmadd(J21L,PDstandard4th22gxy,kmul(J31L,PDstandard4th23gxy)),kmadd(dJ212L,PDstandard4th2gxy,kmadd(J32L,kmadd(J21L,PDstandard4th23gxy,kmul(J31L,PDstandard4th33gxy)),kmul(dJ312L,PDstandard4th3gxy))))))); JacPDstandard4th12gxz = kmadd(J12L,kmadd(J11L,PDstandard4th11gxz,kmadd(J21L,PDstandard4th12gxz,kmul(J31L,PDstandard4th13gxz))),kmadd(J11L,kmadd(J22L,PDstandard4th12gxz,kmul(J32L,PDstandard4th13gxz)),kmadd(dJ112L,PDstandard4th1gxz,kmadd(J22L,kmadd(J21L,PDstandard4th22gxz,kmul(J31L,PDstandard4th23gxz)),kmadd(dJ212L,PDstandard4th2gxz,kmadd(J32L,kmadd(J21L,PDstandard4th23gxz,kmul(J31L,PDstandard4th33gxz)),kmul(dJ312L,PDstandard4th3gxz))))))); - JacPDstandard4th12gyy = - kmadd(J12L,kmadd(J11L,PDstandard4th11gyy,kmadd(J21L,PDstandard4th12gyy,kmul(J31L,PDstandard4th13gyy))),kmadd(J11L,kmadd(J22L,PDstandard4th12gyy,kmul(J32L,PDstandard4th13gyy)),kmadd(dJ112L,PDstandard4th1gyy,kmadd(J22L,kmadd(J21L,PDstandard4th22gyy,kmul(J31L,PDstandard4th23gyy)),kmadd(dJ212L,PDstandard4th2gyy,kmadd(J32L,kmadd(J21L,PDstandard4th23gyy,kmul(J31L,PDstandard4th33gyy)),kmul(dJ312L,PDstandard4th3gyy))))))); - JacPDstandard4th12gyz = kmadd(J12L,kmadd(J11L,PDstandard4th11gyz,kmadd(J21L,PDstandard4th12gyz,kmul(J31L,PDstandard4th13gyz))),kmadd(J11L,kmadd(J22L,PDstandard4th12gyz,kmul(J32L,PDstandard4th13gyz)),kmadd(dJ112L,PDstandard4th1gyz,kmadd(J22L,kmadd(J21L,PDstandard4th22gyz,kmul(J31L,PDstandard4th23gyz)),kmadd(dJ212L,PDstandard4th2gyz,kmadd(J32L,kmadd(J21L,PDstandard4th23gyz,kmul(J31L,PDstandard4th33gyz)),kmul(dJ312L,PDstandard4th3gyz))))))); JacPDstandard4th12gzz = kmadd(J12L,kmadd(J11L,PDstandard4th11gzz,kmadd(J21L,PDstandard4th12gzz,kmul(J31L,PDstandard4th13gzz))),kmadd(J11L,kmadd(J22L,PDstandard4th12gzz,kmul(J32L,PDstandard4th13gzz)),kmadd(dJ112L,PDstandard4th1gzz,kmadd(J22L,kmadd(J21L,PDstandard4th22gzz,kmul(J31L,PDstandard4th23gzz)),kmadd(dJ212L,PDstandard4th2gzz,kmadd(J32L,kmadd(J21L,PDstandard4th23gzz,kmul(J31L,PDstandard4th33gzz)),kmul(dJ312L,PDstandard4th3gzz))))))); - JacPDstandard4th13gxx = - kmadd(J13L,kmadd(J11L,PDstandard4th11gxx,kmadd(J21L,PDstandard4th12gxx,kmul(J31L,PDstandard4th13gxx))),kmadd(J11L,kmadd(J23L,PDstandard4th12gxx,kmul(J33L,PDstandard4th13gxx)),kmadd(dJ113L,PDstandard4th1gxx,kmadd(J23L,kmadd(J21L,PDstandard4th22gxx,kmul(J31L,PDstandard4th23gxx)),kmadd(dJ213L,PDstandard4th2gxx,kmadd(J33L,kmadd(J21L,PDstandard4th23gxx,kmul(J31L,PDstandard4th33gxx)),kmul(dJ313L,PDstandard4th3gxx))))))); - - JacPDstandard4th13gxy = - kmadd(J13L,kmadd(J11L,PDstandard4th11gxy,kmadd(J21L,PDstandard4th12gxy,kmul(J31L,PDstandard4th13gxy))),kmadd(J11L,kmadd(J23L,PDstandard4th12gxy,kmul(J33L,PDstandard4th13gxy)),kmadd(dJ113L,PDstandard4th1gxy,kmadd(J23L,kmadd(J21L,PDstandard4th22gxy,kmul(J31L,PDstandard4th23gxy)),kmadd(dJ213L,PDstandard4th2gxy,kmadd(J33L,kmadd(J21L,PDstandard4th23gxy,kmul(J31L,PDstandard4th33gxy)),kmul(dJ313L,PDstandard4th3gxy))))))); - JacPDstandard4th13gxz = kmadd(J13L,kmadd(J11L,PDstandard4th11gxz,kmadd(J21L,PDstandard4th12gxz,kmul(J31L,PDstandard4th13gxz))),kmadd(J11L,kmadd(J23L,PDstandard4th12gxz,kmul(J33L,PDstandard4th13gxz)),kmadd(dJ113L,PDstandard4th1gxz,kmadd(J23L,kmadd(J21L,PDstandard4th22gxz,kmul(J31L,PDstandard4th23gxz)),kmadd(dJ213L,PDstandard4th2gxz,kmadd(J33L,kmadd(J21L,PDstandard4th23gxz,kmul(J31L,PDstandard4th33gxz)),kmul(dJ313L,PDstandard4th3gxz))))))); - JacPDstandard4th13gyy = - kmadd(J13L,kmadd(J11L,PDstandard4th11gyy,kmadd(J21L,PDstandard4th12gyy,kmul(J31L,PDstandard4th13gyy))),kmadd(J11L,kmadd(J23L,PDstandard4th12gyy,kmul(J33L,PDstandard4th13gyy)),kmadd(dJ113L,PDstandard4th1gyy,kmadd(J23L,kmadd(J21L,PDstandard4th22gyy,kmul(J31L,PDstandard4th23gyy)),kmadd(dJ213L,PDstandard4th2gyy,kmadd(J33L,kmadd(J21L,PDstandard4th23gyy,kmul(J31L,PDstandard4th33gyy)),kmul(dJ313L,PDstandard4th3gyy))))))); - - JacPDstandard4th13gyz = - kmadd(J13L,kmadd(J11L,PDstandard4th11gyz,kmadd(J21L,PDstandard4th12gyz,kmul(J31L,PDstandard4th13gyz))),kmadd(J11L,kmadd(J23L,PDstandard4th12gyz,kmul(J33L,PDstandard4th13gyz)),kmadd(dJ113L,PDstandard4th1gyz,kmadd(J23L,kmadd(J21L,PDstandard4th22gyz,kmul(J31L,PDstandard4th23gyz)),kmadd(dJ213L,PDstandard4th2gyz,kmadd(J33L,kmadd(J21L,PDstandard4th23gyz,kmul(J31L,PDstandard4th33gyz)),kmul(dJ313L,PDstandard4th3gyz))))))); - - JacPDstandard4th13gzz = - kmadd(J13L,kmadd(J11L,PDstandard4th11gzz,kmadd(J21L,PDstandard4th12gzz,kmul(J31L,PDstandard4th13gzz))),kmadd(J11L,kmadd(J23L,PDstandard4th12gzz,kmul(J33L,PDstandard4th13gzz)),kmadd(dJ113L,PDstandard4th1gzz,kmadd(J23L,kmadd(J21L,PDstandard4th22gzz,kmul(J31L,PDstandard4th23gzz)),kmadd(dJ213L,PDstandard4th2gzz,kmadd(J33L,kmadd(J21L,PDstandard4th23gzz,kmul(J31L,PDstandard4th33gzz)),kmul(dJ313L,PDstandard4th3gzz))))))); - - JacPDstandard4th21gxx = - kmadd(J12L,kmadd(J11L,PDstandard4th11gxx,kmadd(J21L,PDstandard4th12gxx,kmul(J31L,PDstandard4th13gxx))),kmadd(J11L,kmadd(J22L,PDstandard4th12gxx,kmul(J32L,PDstandard4th13gxx)),kmadd(dJ112L,PDstandard4th1gxx,kmadd(J22L,kmadd(J21L,PDstandard4th22gxx,kmul(J31L,PDstandard4th23gxx)),kmadd(dJ212L,PDstandard4th2gxx,kmadd(J32L,kmadd(J21L,PDstandard4th23gxx,kmul(J31L,PDstandard4th33gxx)),kmul(dJ312L,PDstandard4th3gxx))))))); - JacPDstandard4th21gxy = kmadd(J12L,kmadd(J11L,PDstandard4th11gxy,kmadd(J21L,PDstandard4th12gxy,kmul(J31L,PDstandard4th13gxy))),kmadd(J11L,kmadd(J22L,PDstandard4th12gxy,kmul(J32L,PDstandard4th13gxy)),kmadd(dJ112L,PDstandard4th1gxy,kmadd(J22L,kmadd(J21L,PDstandard4th22gxy,kmul(J31L,PDstandard4th23gxy)),kmadd(dJ212L,PDstandard4th2gxy,kmadd(J32L,kmadd(J21L,PDstandard4th23gxy,kmul(J31L,PDstandard4th33gxy)),kmul(dJ312L,PDstandard4th3gxy))))))); - JacPDstandard4th21gxz = - kmadd(J12L,kmadd(J11L,PDstandard4th11gxz,kmadd(J21L,PDstandard4th12gxz,kmul(J31L,PDstandard4th13gxz))),kmadd(J11L,kmadd(J22L,PDstandard4th12gxz,kmul(J32L,PDstandard4th13gxz)),kmadd(dJ112L,PDstandard4th1gxz,kmadd(J22L,kmadd(J21L,PDstandard4th22gxz,kmul(J31L,PDstandard4th23gxz)),kmadd(dJ212L,PDstandard4th2gxz,kmadd(J32L,kmadd(J21L,PDstandard4th23gxz,kmul(J31L,PDstandard4th33gxz)),kmul(dJ312L,PDstandard4th3gxz))))))); - - JacPDstandard4th21gyy = - kmadd(J12L,kmadd(J11L,PDstandard4th11gyy,kmadd(J21L,PDstandard4th12gyy,kmul(J31L,PDstandard4th13gyy))),kmadd(J11L,kmadd(J22L,PDstandard4th12gyy,kmul(J32L,PDstandard4th13gyy)),kmadd(dJ112L,PDstandard4th1gyy,kmadd(J22L,kmadd(J21L,PDstandard4th22gyy,kmul(J31L,PDstandard4th23gyy)),kmadd(dJ212L,PDstandard4th2gyy,kmadd(J32L,kmadd(J21L,PDstandard4th23gyy,kmul(J31L,PDstandard4th33gyy)),kmul(dJ312L,PDstandard4th3gyy))))))); - - JacPDstandard4th21gyz = - kmadd(J12L,kmadd(J11L,PDstandard4th11gyz,kmadd(J21L,PDstandard4th12gyz,kmul(J31L,PDstandard4th13gyz))),kmadd(J11L,kmadd(J22L,PDstandard4th12gyz,kmul(J32L,PDstandard4th13gyz)),kmadd(dJ112L,PDstandard4th1gyz,kmadd(J22L,kmadd(J21L,PDstandard4th22gyz,kmul(J31L,PDstandard4th23gyz)),kmadd(dJ212L,PDstandard4th2gyz,kmadd(J32L,kmadd(J21L,PDstandard4th23gyz,kmul(J31L,PDstandard4th33gyz)),kmul(dJ312L,PDstandard4th3gyz))))))); - - JacPDstandard4th21gzz = - kmadd(J12L,kmadd(J11L,PDstandard4th11gzz,kmadd(J21L,PDstandard4th12gzz,kmul(J31L,PDstandard4th13gzz))),kmadd(J11L,kmadd(J22L,PDstandard4th12gzz,kmul(J32L,PDstandard4th13gzz)),kmadd(dJ112L,PDstandard4th1gzz,kmadd(J22L,kmadd(J21L,PDstandard4th22gzz,kmul(J31L,PDstandard4th23gzz)),kmadd(dJ212L,PDstandard4th2gzz,kmadd(J32L,kmadd(J21L,PDstandard4th23gzz,kmul(J31L,PDstandard4th33gzz)),kmul(dJ312L,PDstandard4th3gzz))))))); - JacPDstandard4th23gxx = kmadd(J13L,kmadd(J12L,PDstandard4th11gxx,kmadd(J22L,PDstandard4th12gxx,kmul(J32L,PDstandard4th13gxx))),kmadd(J12L,kmadd(J23L,PDstandard4th12gxx,kmul(J33L,PDstandard4th13gxx)),kmadd(dJ123L,PDstandard4th1gxx,kmadd(J23L,kmadd(J22L,PDstandard4th22gxx,kmul(J32L,PDstandard4th23gxx)),kmadd(dJ223L,PDstandard4th2gxx,kmadd(J33L,kmadd(J22L,PDstandard4th23gxx,kmul(J32L,PDstandard4th33gxx)),kmul(dJ323L,PDstandard4th3gxx))))))); @@ -900,18 +819,9 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, JacPDstandard4th23gxz = kmadd(J13L,kmadd(J12L,PDstandard4th11gxz,kmadd(J22L,PDstandard4th12gxz,kmul(J32L,PDstandard4th13gxz))),kmadd(J12L,kmadd(J23L,PDstandard4th12gxz,kmul(J33L,PDstandard4th13gxz)),kmadd(dJ123L,PDstandard4th1gxz,kmadd(J23L,kmadd(J22L,PDstandard4th22gxz,kmul(J32L,PDstandard4th23gxz)),kmadd(dJ223L,PDstandard4th2gxz,kmadd(J33L,kmadd(J22L,PDstandard4th23gxz,kmul(J32L,PDstandard4th33gxz)),kmul(dJ323L,PDstandard4th3gxz))))))); - JacPDstandard4th23gyy = - kmadd(J13L,kmadd(J12L,PDstandard4th11gyy,kmadd(J22L,PDstandard4th12gyy,kmul(J32L,PDstandard4th13gyy))),kmadd(J12L,kmadd(J23L,PDstandard4th12gyy,kmul(J33L,PDstandard4th13gyy)),kmadd(dJ123L,PDstandard4th1gyy,kmadd(J23L,kmadd(J22L,PDstandard4th22gyy,kmul(J32L,PDstandard4th23gyy)),kmadd(dJ223L,PDstandard4th2gyy,kmadd(J33L,kmadd(J22L,PDstandard4th23gyy,kmul(J32L,PDstandard4th33gyy)),kmul(dJ323L,PDstandard4th3gyy))))))); - JacPDstandard4th23gyz = kmadd(J13L,kmadd(J12L,PDstandard4th11gyz,kmadd(J22L,PDstandard4th12gyz,kmul(J32L,PDstandard4th13gyz))),kmadd(J12L,kmadd(J23L,PDstandard4th12gyz,kmul(J33L,PDstandard4th13gyz)),kmadd(dJ123L,PDstandard4th1gyz,kmadd(J23L,kmadd(J22L,PDstandard4th22gyz,kmul(J32L,PDstandard4th23gyz)),kmadd(dJ223L,PDstandard4th2gyz,kmadd(J33L,kmadd(J22L,PDstandard4th23gyz,kmul(J32L,PDstandard4th33gyz)),kmul(dJ323L,PDstandard4th3gyz))))))); - JacPDstandard4th23gzz = - kmadd(J13L,kmadd(J12L,PDstandard4th11gzz,kmadd(J22L,PDstandard4th12gzz,kmul(J32L,PDstandard4th13gzz))),kmadd(J12L,kmadd(J23L,PDstandard4th12gzz,kmul(J33L,PDstandard4th13gzz)),kmadd(dJ123L,PDstandard4th1gzz,kmadd(J23L,kmadd(J22L,PDstandard4th22gzz,kmul(J32L,PDstandard4th23gzz)),kmadd(dJ223L,PDstandard4th2gzz,kmadd(J33L,kmadd(J22L,PDstandard4th23gzz,kmul(J32L,PDstandard4th33gzz)),kmul(dJ323L,PDstandard4th3gzz))))))); - - JacPDstandard4th31gxx = - kmadd(J13L,kmadd(J11L,PDstandard4th11gxx,kmadd(J21L,PDstandard4th12gxx,kmul(J31L,PDstandard4th13gxx))),kmadd(J11L,kmadd(J23L,PDstandard4th12gxx,kmul(J33L,PDstandard4th13gxx)),kmadd(dJ113L,PDstandard4th1gxx,kmadd(J23L,kmadd(J21L,PDstandard4th22gxx,kmul(J31L,PDstandard4th23gxx)),kmadd(dJ213L,PDstandard4th2gxx,kmadd(J33L,kmadd(J21L,PDstandard4th23gxx,kmul(J31L,PDstandard4th33gxx)),kmul(dJ313L,PDstandard4th3gxx))))))); - JacPDstandard4th31gxy = kmadd(J13L,kmadd(J11L,PDstandard4th11gxy,kmadd(J21L,PDstandard4th12gxy,kmul(J31L,PDstandard4th13gxy))),kmadd(J11L,kmadd(J23L,PDstandard4th12gxy,kmul(J33L,PDstandard4th13gxy)),kmadd(dJ113L,PDstandard4th1gxy,kmadd(J23L,kmadd(J21L,PDstandard4th22gxy,kmul(J31L,PDstandard4th23gxy)),kmadd(dJ213L,PDstandard4th2gxy,kmadd(J33L,kmadd(J21L,PDstandard4th23gxy,kmul(J31L,PDstandard4th33gxy)),kmul(dJ313L,PDstandard4th3gxy))))))); @@ -924,26 +834,8 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, JacPDstandard4th31gyz = kmadd(J13L,kmadd(J11L,PDstandard4th11gyz,kmadd(J21L,PDstandard4th12gyz,kmul(J31L,PDstandard4th13gyz))),kmadd(J11L,kmadd(J23L,PDstandard4th12gyz,kmul(J33L,PDstandard4th13gyz)),kmadd(dJ113L,PDstandard4th1gyz,kmadd(J23L,kmadd(J21L,PDstandard4th22gyz,kmul(J31L,PDstandard4th23gyz)),kmadd(dJ213L,PDstandard4th2gyz,kmadd(J33L,kmadd(J21L,PDstandard4th23gyz,kmul(J31L,PDstandard4th33gyz)),kmul(dJ313L,PDstandard4th3gyz))))))); - JacPDstandard4th31gzz = - kmadd(J13L,kmadd(J11L,PDstandard4th11gzz,kmadd(J21L,PDstandard4th12gzz,kmul(J31L,PDstandard4th13gzz))),kmadd(J11L,kmadd(J23L,PDstandard4th12gzz,kmul(J33L,PDstandard4th13gzz)),kmadd(dJ113L,PDstandard4th1gzz,kmadd(J23L,kmadd(J21L,PDstandard4th22gzz,kmul(J31L,PDstandard4th23gzz)),kmadd(dJ213L,PDstandard4th2gzz,kmadd(J33L,kmadd(J21L,PDstandard4th23gzz,kmul(J31L,PDstandard4th33gzz)),kmul(dJ313L,PDstandard4th3gzz))))))); - - JacPDstandard4th32gxx = - kmadd(J13L,kmadd(J12L,PDstandard4th11gxx,kmadd(J22L,PDstandard4th12gxx,kmul(J32L,PDstandard4th13gxx))),kmadd(J12L,kmadd(J23L,PDstandard4th12gxx,kmul(J33L,PDstandard4th13gxx)),kmadd(dJ123L,PDstandard4th1gxx,kmadd(J23L,kmadd(J22L,PDstandard4th22gxx,kmul(J32L,PDstandard4th23gxx)),kmadd(dJ223L,PDstandard4th2gxx,kmadd(J33L,kmadd(J22L,PDstandard4th23gxx,kmul(J32L,PDstandard4th33gxx)),kmul(dJ323L,PDstandard4th3gxx))))))); - - JacPDstandard4th32gxy = - kmadd(J13L,kmadd(J12L,PDstandard4th11gxy,kmadd(J22L,PDstandard4th12gxy,kmul(J32L,PDstandard4th13gxy))),kmadd(J12L,kmadd(J23L,PDstandard4th12gxy,kmul(J33L,PDstandard4th13gxy)),kmadd(dJ123L,PDstandard4th1gxy,kmadd(J23L,kmadd(J22L,PDstandard4th22gxy,kmul(J32L,PDstandard4th23gxy)),kmadd(dJ223L,PDstandard4th2gxy,kmadd(J33L,kmadd(J22L,PDstandard4th23gxy,kmul(J32L,PDstandard4th33gxy)),kmul(dJ323L,PDstandard4th3gxy))))))); - - JacPDstandard4th32gxz = - kmadd(J13L,kmadd(J12L,PDstandard4th11gxz,kmadd(J22L,PDstandard4th12gxz,kmul(J32L,PDstandard4th13gxz))),kmadd(J12L,kmadd(J23L,PDstandard4th12gxz,kmul(J33L,PDstandard4th13gxz)),kmadd(dJ123L,PDstandard4th1gxz,kmadd(J23L,kmadd(J22L,PDstandard4th22gxz,kmul(J32L,PDstandard4th23gxz)),kmadd(dJ223L,PDstandard4th2gxz,kmadd(J33L,kmadd(J22L,PDstandard4th23gxz,kmul(J32L,PDstandard4th33gxz)),kmul(dJ323L,PDstandard4th3gxz))))))); - - JacPDstandard4th32gyy = - kmadd(J13L,kmadd(J12L,PDstandard4th11gyy,kmadd(J22L,PDstandard4th12gyy,kmul(J32L,PDstandard4th13gyy))),kmadd(J12L,kmadd(J23L,PDstandard4th12gyy,kmul(J33L,PDstandard4th13gyy)),kmadd(dJ123L,PDstandard4th1gyy,kmadd(J23L,kmadd(J22L,PDstandard4th22gyy,kmul(J32L,PDstandard4th23gyy)),kmadd(dJ223L,PDstandard4th2gyy,kmadd(J33L,kmadd(J22L,PDstandard4th23gyy,kmul(J32L,PDstandard4th33gyy)),kmul(dJ323L,PDstandard4th3gyy))))))); - JacPDstandard4th32gyz = kmadd(J13L,kmadd(J12L,PDstandard4th11gyz,kmadd(J22L,PDstandard4th12gyz,kmul(J32L,PDstandard4th13gyz))),kmadd(J12L,kmadd(J23L,PDstandard4th12gyz,kmul(J33L,PDstandard4th13gyz)),kmadd(dJ123L,PDstandard4th1gyz,kmadd(J23L,kmadd(J22L,PDstandard4th22gyz,kmul(J32L,PDstandard4th23gyz)),kmadd(dJ223L,PDstandard4th2gyz,kmadd(J33L,kmadd(J22L,PDstandard4th23gyz,kmul(J32L,PDstandard4th33gyz)),kmul(dJ323L,PDstandard4th3gyz))))))); - - JacPDstandard4th32gzz = - kmadd(J13L,kmadd(J12L,PDstandard4th11gzz,kmadd(J22L,PDstandard4th12gzz,kmul(J32L,PDstandard4th13gzz))),kmadd(J12L,kmadd(J23L,PDstandard4th12gzz,kmul(J33L,PDstandard4th13gzz)),kmadd(dJ123L,PDstandard4th1gzz,kmadd(J23L,kmadd(J22L,PDstandard4th22gzz,kmul(J32L,PDstandard4th23gzz)),kmadd(dJ223L,PDstandard4th2gzz,kmadd(J33L,kmadd(J22L,PDstandard4th23gzz,kmul(J32L,PDstandard4th33gzz)),kmul(dJ323L,PDstandard4th3gzz))))))); } else { @@ -1013,10 +905,6 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, JacPDstandard4th3kyz = PDstandard4th3kyz; - JacPDstandard4th11gxy = PDstandard4th11gxy; - - JacPDstandard4th11gxz = PDstandard4th11gxz; - JacPDstandard4th11gyy = PDstandard4th11gyy; JacPDstandard4th11gyz = PDstandard4th11gyz; @@ -1025,74 +913,36 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, JacPDstandard4th22gxx = PDstandard4th22gxx; - JacPDstandard4th22gxy = PDstandard4th22gxy; - JacPDstandard4th22gxz = PDstandard4th22gxz; - JacPDstandard4th22gyz = PDstandard4th22gyz; - JacPDstandard4th22gzz = PDstandard4th22gzz; JacPDstandard4th33gxx = PDstandard4th33gxx; JacPDstandard4th33gxy = PDstandard4th33gxy; - JacPDstandard4th33gxz = PDstandard4th33gxz; - JacPDstandard4th33gyy = PDstandard4th33gyy; - JacPDstandard4th33gyz = PDstandard4th33gyz; - - JacPDstandard4th12gxx = PDstandard4th12gxx; - JacPDstandard4th12gxy = PDstandard4th12gxy; JacPDstandard4th12gxz = PDstandard4th12gxz; - JacPDstandard4th12gyy = PDstandard4th12gyy; - JacPDstandard4th12gyz = PDstandard4th12gyz; JacPDstandard4th12gzz = PDstandard4th12gzz; - JacPDstandard4th13gxx = PDstandard4th13gxx; - - JacPDstandard4th13gxy = PDstandard4th13gxy; - JacPDstandard4th13gxz = PDstandard4th13gxz; - JacPDstandard4th13gyy = PDstandard4th13gyy; - - JacPDstandard4th13gyz = PDstandard4th13gyz; - - JacPDstandard4th13gzz = PDstandard4th13gzz; - - JacPDstandard4th21gxx = PDstandard4th12gxx; - JacPDstandard4th21gxy = PDstandard4th12gxy; - JacPDstandard4th21gxz = PDstandard4th12gxz; - - JacPDstandard4th21gyy = PDstandard4th12gyy; - - JacPDstandard4th21gyz = PDstandard4th12gyz; - - JacPDstandard4th21gzz = PDstandard4th12gzz; - JacPDstandard4th23gxx = PDstandard4th23gxx; JacPDstandard4th23gxy = PDstandard4th23gxy; JacPDstandard4th23gxz = PDstandard4th23gxz; - JacPDstandard4th23gyy = PDstandard4th23gyy; - JacPDstandard4th23gyz = PDstandard4th23gyz; - JacPDstandard4th23gzz = PDstandard4th23gzz; - - JacPDstandard4th31gxx = PDstandard4th13gxx; - JacPDstandard4th31gxy = PDstandard4th13gxy; JacPDstandard4th31gxz = PDstandard4th13gxz; @@ -1101,19 +951,7 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, JacPDstandard4th31gyz = PDstandard4th13gyz; - JacPDstandard4th31gzz = PDstandard4th13gzz; - - JacPDstandard4th32gxx = PDstandard4th23gxx; - - JacPDstandard4th32gxy = PDstandard4th23gxy; - - JacPDstandard4th32gxz = PDstandard4th23gxz; - - JacPDstandard4th32gyy = PDstandard4th23gyy; - JacPDstandard4th32gyz = PDstandard4th23gyz; - - JacPDstandard4th32gzz = PDstandard4th23gzz; } CCTK_REAL_VEC detg = @@ -1317,444 +1155,42 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC nn = isqrt2; - CCTK_REAL_VEC R1111 = ToReal(0.); - - CCTK_REAL_VEC R1112 = - kmul(ksub(JacPDstandard4th21gxx,JacPDstandard4th12gxx),ToReal(0.5)); - - CCTK_REAL_VEC R1113 = - kmul(ksub(JacPDstandard4th31gxx,JacPDstandard4th13gxx),ToReal(0.5)); - - CCTK_REAL_VEC R1121 = - kmul(ksub(JacPDstandard4th21gxx,JacPDstandard4th12gxx),ToReal(0.5)); - - CCTK_REAL_VEC R1122 = - ksub(JacPDstandard4th21gxy,JacPDstandard4th12gxy); - - CCTK_REAL_VEC R1123 = - kmul(kadd(JacPDstandard4th21gxz,ksub(JacPDstandard4th31gxy,kadd(JacPDstandard4th13gxy,JacPDstandard4th12gxz))),ToReal(0.5)); - - CCTK_REAL_VEC R1131 = - kmul(ksub(JacPDstandard4th31gxx,JacPDstandard4th13gxx),ToReal(0.5)); - - CCTK_REAL_VEC R1132 = - kmul(kadd(JacPDstandard4th21gxz,ksub(JacPDstandard4th31gxy,kadd(JacPDstandard4th13gxy,JacPDstandard4th12gxz))),ToReal(0.5)); - - CCTK_REAL_VEC R1133 = - ksub(JacPDstandard4th31gxz,JacPDstandard4th13gxz); - - CCTK_REAL_VEC R1211 = - kmul(ksub(JacPDstandard4th12gxx,JacPDstandard4th21gxx),ToReal(0.5)); - CCTK_REAL_VEC R1212 = kmul(ToReal(0.5),kadd(JacPDstandard4th12gxy,kadd(JacPDstandard4th21gxy,kmadd(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmadd(gzzL,kmul(gamma311,gamma322),kmul(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)))))))))),ToReal(-2.),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma121,gamma121),ToReal(2.)),kmadd(gyyL,kmul(kmul(gamma221,gamma221),ToReal(2.)),kmadd(gzzL,kmul(kmul(gamma321,gamma321),ToReal(2.)),kmul(kmadd(gyzL,kmul(gamma221,gamma321),kmul(gamma121,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)))),ToReal(4.))))),JacPDstandard4th22gxx),JacPDstandard4th11gyy))))); CCTK_REAL_VEC R1213 = kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard4th12gxz,kadd(JacPDstandard4th31gxy,kmadd(kmadd(gyyL,kmul(gamma211,gamma232),kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmadd(gzzL,kmul(gamma311,gamma332),kmul(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2.),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma221),kmadd(gyyL,kmul(gamma221,gamma231),kmadd(gxzL,kmul(gamma131,gamma321),kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))))),ToReal(2.),JacPDstandard4th23gxx),JacPDstandard4th11gyz)))))); - CCTK_REAL_VEC R1221 = - kmadd(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322))),kmadd(gamma121,kmul(kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),ToReal(-2.)),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyy,kadd(JacPDstandard4th22gxx,kmadd(gyzL,kmul(gamma221,kmul(gamma321,ToReal(-4.))),kmadd(gyyL,kmul(kmul(gamma221,gamma221),ToReal(-2.)),kmadd(gzzL,kmul(kmul(gamma321,gamma321),ToReal(-2.)),ksub(kmsub(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmul(gzzL,kmul(gamma311,gamma322))))))),ToReal(2.),JacPDstandard4th21gxy),JacPDstandard4th12gxy)))))),kmul(gxxL,kmul(gamma121,gamma121))))); - - CCTK_REAL_VEC R1222 = - kmul(ksub(JacPDstandard4th21gyy,JacPDstandard4th12gyy),ToReal(0.5)); - CCTK_REAL_VEC R1223 = kmadd(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard4th22gxz,kadd(JacPDstandard4th31gyy,kmadd(kmadd(gyyL,kmul(gamma221,gamma232),kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmadd(gzzL,kmul(gamma321,gamma332),kmul(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2.),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma222),kmadd(gyyL,kmul(gamma222,gamma231),kmadd(gxzL,kmul(gamma131,gamma322),kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))))),ToReal(2.),JacPDstandard4th23gxy),JacPDstandard4th12gyz)))))); - CCTK_REAL_VEC R1231 = - kmadd(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyz,kadd(JacPDstandard4th32gxx,kmadd(kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gamma131,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))),ToReal(-2.),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmul(gzzL,kmul(gamma311,gamma332))))),ToReal(2.),kmul(gyyL,kmadd(gamma221,kmul(gamma231,ToReal(-2.)),kmul(gamma211,kmul(gamma232,ToReal(2.)))))),JacPDstandard4th21gxz),JacPDstandard4th13gxy)))),kmul(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); - - CCTK_REAL_VEC R1232 = - kmadd(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard4th21gyz,kadd(JacPDstandard4th32gxy,kmadd(kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gamma131,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))),ToReal(-2.),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmul(gzzL,kmul(gamma321,gamma332))))),ToReal(2.),kmul(gyyL,kmadd(gamma222,kmul(gamma231,ToReal(-2.)),kmul(gamma221,kmul(gamma232,ToReal(2.)))))),JacPDstandard4th22gxz),JacPDstandard4th13gyy)))),kmul(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); - - CCTK_REAL_VEC R1233 = - kmul(kadd(JacPDstandard4th31gyz,ksub(JacPDstandard4th32gxz,kadd(JacPDstandard4th23gxz,JacPDstandard4th13gyz))),ToReal(0.5)); - - CCTK_REAL_VEC R1311 = - kmul(ksub(JacPDstandard4th13gxx,JacPDstandard4th31gxx),ToReal(0.5)); - - CCTK_REAL_VEC R1312 = - kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard4th13gxy,kadd(JacPDstandard4th21gxz,kmadd(kmadd(gyyL,kmul(gamma211,gamma232),kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmadd(gzzL,kmul(gamma311,gamma332),kmul(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2.),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma221),kmadd(gyyL,kmul(gamma221,gamma231),kmadd(gxzL,kmul(gamma131,gamma321),kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))))),ToReal(2.),JacPDstandard4th32gxx),JacPDstandard4th11gyz)))))); - CCTK_REAL_VEC R1313 = kmul(ToReal(0.5),kadd(JacPDstandard4th13gxz,kadd(JacPDstandard4th31gxz,kmadd(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmadd(gzzL,kmul(gamma311,gamma333),kmul(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2.),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma131,gamma131),ToReal(2.)),kmadd(gyyL,kmul(kmul(gamma231,gamma231),ToReal(2.)),kmadd(gzzL,kmul(kmul(gamma331,gamma331),ToReal(2.)),kmul(kmadd(gyzL,kmul(gamma231,gamma331),kmul(gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))),ToReal(4.))))),JacPDstandard4th33gxx),JacPDstandard4th11gzz))))); - CCTK_REAL_VEC R1321 = - kmadd(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyz,kadd(JacPDstandard4th23gxx,kmadd(kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gamma131,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))),ToReal(-2.),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmul(gzzL,kmul(gamma311,gamma332))))),ToReal(2.),kmul(gyyL,kmadd(gamma221,kmul(gamma231,ToReal(-2.)),kmul(gamma211,kmul(gamma232,ToReal(2.)))))),JacPDstandard4th31gxy),JacPDstandard4th12gxz)))),kmul(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); - - CCTK_REAL_VEC R1322 = - kmul(kadd(JacPDstandard4th21gyz,ksub(JacPDstandard4th23gxy,kadd(JacPDstandard4th32gxy,JacPDstandard4th12gyz))),ToReal(0.5)); - CCTK_REAL_VEC R1323 = kmadd(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmul(ToReal(0.5),kadd(JacPDstandard4th23gxz,kadd(JacPDstandard4th31gyz,kmadd(kmadd(gyyL,kmul(gamma221,gamma233),kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmadd(gzzL,kmul(gamma321,gamma333),kmul(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))))))))),ToReal(-2.),ksub(kmsub(kmadd(gxyL,kmul(gamma132,gamma231),kmadd(gyyL,kmul(gamma231,gamma232),kmadd(gxzL,kmul(gamma132,gamma331),kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))))),ToReal(2.),JacPDstandard4th33gxy),JacPDstandard4th12gzz)))))); - CCTK_REAL_VEC R1331 = - kmadd(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmadd(gamma131,kmul(kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),ToReal(-2.)),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gzz,kadd(JacPDstandard4th33gxx,kmadd(gyzL,kmul(gamma231,kmul(gamma331,ToReal(-4.))),kmadd(gyyL,kmul(kmul(gamma231,gamma231),ToReal(-2.)),kmadd(gzzL,kmul(kmul(gamma331,gamma331),ToReal(-2.)),ksub(kmsub(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmul(gzzL,kmul(gamma311,gamma333))))))),ToReal(2.),JacPDstandard4th31gxz),JacPDstandard4th13gxz)))))),kmul(gxxL,kmul(gamma131,gamma131))))); - - CCTK_REAL_VEC R1332 = - kmadd(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmsub(ToReal(0.5),kadd(JacPDstandard4th21gzz,kadd(JacPDstandard4th33gxy,kmadd(kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gamma132,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))),ToReal(-2.),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmul(gzzL,kmul(gamma321,gamma333))))),ToReal(2.),kmul(gyyL,kmadd(gamma231,kmul(gamma232,ToReal(-2.)),kmul(gamma221,kmul(gamma233,ToReal(2.)))))),JacPDstandard4th32gxz),JacPDstandard4th13gyz)))),kmul(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))))); - - CCTK_REAL_VEC R1333 = - kmul(ksub(JacPDstandard4th31gzz,JacPDstandard4th13gzz),ToReal(0.5)); - - CCTK_REAL_VEC R2111 = - kmul(ksub(JacPDstandard4th12gxx,JacPDstandard4th21gxx),ToReal(0.5)); - - CCTK_REAL_VEC R2112 = - kmadd(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322))),kmadd(gamma121,kmul(kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),ToReal(-2.)),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyy,kadd(JacPDstandard4th22gxx,kmadd(gyzL,kmul(gamma221,kmul(gamma321,ToReal(-4.))),kmadd(gyyL,kmul(kmul(gamma221,gamma221),ToReal(-2.)),kmadd(gzzL,kmul(kmul(gamma321,gamma321),ToReal(-2.)),ksub(kmsub(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmul(gzzL,kmul(gamma311,gamma322))))))),ToReal(2.),JacPDstandard4th21gxy),JacPDstandard4th12gxy)))))),kmul(gxxL,kmul(gamma121,gamma121))))); - - CCTK_REAL_VEC R2113 = - kmadd(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyz,kadd(JacPDstandard4th32gxx,kmadd(kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gamma131,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))),ToReal(-2.),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmul(gzzL,kmul(gamma311,gamma332))))),ToReal(2.),kmul(gyyL,kmadd(gamma221,kmul(gamma231,ToReal(-2.)),kmul(gamma211,kmul(gamma232,ToReal(2.)))))),JacPDstandard4th21gxz),JacPDstandard4th13gxy)))),kmul(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); - - CCTK_REAL_VEC R2121 = - kmul(ToReal(0.5),kadd(JacPDstandard4th12gxy,kadd(JacPDstandard4th21gxy,kmadd(kmadd(gxyL,kmul(gamma122,gamma211),kmadd(gyyL,kmul(gamma211,gamma222),kmadd(gxzL,kmul(gamma122,gamma311),kmadd(gyzL,kmul(gamma222,gamma311),kmadd(gyzL,kmul(gamma211,gamma322),kmadd(gzzL,kmul(gamma311,gamma322),kmul(gamma111,kmadd(gxxL,gamma122,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)))))))))),ToReal(-2.),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma121,gamma121),ToReal(2.)),kmadd(gyyL,kmul(kmul(gamma221,gamma221),ToReal(2.)),kmadd(gzzL,kmul(kmul(gamma321,gamma321),ToReal(2.)),kmul(kmadd(gyzL,kmul(gamma221,gamma321),kmul(gamma121,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)))),ToReal(4.))))),JacPDstandard4th22gxx),JacPDstandard4th11gyy))))); - - CCTK_REAL_VEC R2122 = - kmul(ksub(JacPDstandard4th21gyy,JacPDstandard4th12gyy),ToReal(0.5)); - - CCTK_REAL_VEC R2123 = - kmadd(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard4th21gyz,kadd(JacPDstandard4th32gxy,kmadd(kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gamma131,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))),ToReal(-2.),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmul(gzzL,kmul(gamma321,gamma332))))),ToReal(2.),kmul(gyyL,kmadd(gamma222,kmul(gamma231,ToReal(-2.)),kmul(gamma221,kmul(gamma232,ToReal(2.)))))),JacPDstandard4th22gxz),JacPDstandard4th13gyy)))),kmul(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); - - CCTK_REAL_VEC R2131 = - kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard4th12gxz,kadd(JacPDstandard4th31gxy,kmadd(kmadd(gyyL,kmul(gamma211,gamma232),kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmadd(gzzL,kmul(gamma311,gamma332),kmul(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2.),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma221),kmadd(gyyL,kmul(gamma221,gamma231),kmadd(gxzL,kmul(gamma131,gamma321),kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))))),ToReal(2.),JacPDstandard4th23gxx),JacPDstandard4th11gyz)))))); - - CCTK_REAL_VEC R2132 = - kmadd(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard4th22gxz,kadd(JacPDstandard4th31gyy,kmadd(kmadd(gyyL,kmul(gamma221,gamma232),kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmadd(gzzL,kmul(gamma321,gamma332),kmul(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2.),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma222),kmadd(gyyL,kmul(gamma222,gamma231),kmadd(gxzL,kmul(gamma131,gamma322),kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))))),ToReal(2.),JacPDstandard4th23gxy),JacPDstandard4th12gyz)))))); - - CCTK_REAL_VEC R2133 = - kmul(kadd(JacPDstandard4th31gyz,ksub(JacPDstandard4th32gxz,kadd(JacPDstandard4th23gxz,JacPDstandard4th13gyz))),ToReal(0.5)); - - CCTK_REAL_VEC R2211 = - ksub(JacPDstandard4th12gxy,JacPDstandard4th21gxy); - - CCTK_REAL_VEC R2212 = - kmul(ksub(JacPDstandard4th12gyy,JacPDstandard4th21gyy),ToReal(0.5)); - - CCTK_REAL_VEC R2213 = - kmul(kadd(JacPDstandard4th12gyz,ksub(JacPDstandard4th32gxy,kadd(JacPDstandard4th23gxy,JacPDstandard4th21gyz))),ToReal(0.5)); - - CCTK_REAL_VEC R2221 = - kmul(ksub(JacPDstandard4th12gyy,JacPDstandard4th21gyy),ToReal(0.5)); - - CCTK_REAL_VEC R2222 = ToReal(0.); - - CCTK_REAL_VEC R2223 = - kmul(ksub(JacPDstandard4th32gyy,JacPDstandard4th23gyy),ToReal(0.5)); - - CCTK_REAL_VEC R2231 = - kmul(kadd(JacPDstandard4th12gyz,ksub(JacPDstandard4th32gxy,kadd(JacPDstandard4th23gxy,JacPDstandard4th21gyz))),ToReal(0.5)); - - CCTK_REAL_VEC R2232 = - kmul(ksub(JacPDstandard4th32gyy,JacPDstandard4th23gyy),ToReal(0.5)); - - CCTK_REAL_VEC R2233 = - ksub(JacPDstandard4th32gyz,JacPDstandard4th23gyz); - - CCTK_REAL_VEC R2311 = - kmul(kadd(JacPDstandard4th12gxz,ksub(JacPDstandard4th13gxy,kadd(JacPDstandard4th31gxy,JacPDstandard4th21gxz))),ToReal(0.5)); - - CCTK_REAL_VEC R2312 = - kmadd(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard4th13gyy,kadd(JacPDstandard4th22gxz,kmadd(kmadd(gyyL,kmul(gamma221,gamma232),kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmadd(gzzL,kmul(gamma321,gamma332),kmul(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2.),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma222),kmadd(gyyL,kmul(gamma222,gamma231),kmadd(gxzL,kmul(gamma131,gamma322),kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))))),ToReal(2.),JacPDstandard4th32gxy),JacPDstandard4th21gyz)))))); - - CCTK_REAL_VEC R2313 = - kmadd(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmul(ToReal(0.5),kadd(JacPDstandard4th13gyz,kadd(JacPDstandard4th32gxz,kmadd(kmadd(gyyL,kmul(gamma221,gamma233),kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmadd(gzzL,kmul(gamma321,gamma333),kmul(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))))))))),ToReal(-2.),ksub(kmsub(kmadd(gxyL,kmul(gamma132,gamma231),kmadd(gyyL,kmul(gamma231,gamma232),kmadd(gxzL,kmul(gamma132,gamma331),kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))))),ToReal(2.),JacPDstandard4th33gxy),JacPDstandard4th21gzz)))))); - - CCTK_REAL_VEC R2321 = - kmadd(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard4th12gyz,kadd(JacPDstandard4th23gxy,kmadd(kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gamma131,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))),ToReal(-2.),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmul(gzzL,kmul(gamma321,gamma332))))),ToReal(2.),kmul(gyyL,kmadd(gamma222,kmul(gamma231,ToReal(-2.)),kmul(gamma221,kmul(gamma232,ToReal(2.)))))),JacPDstandard4th31gyy),JacPDstandard4th22gxz)))),kmul(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); - - CCTK_REAL_VEC R2322 = - kmul(ksub(JacPDstandard4th23gyy,JacPDstandard4th32gyy),ToReal(0.5)); - CCTK_REAL_VEC R2323 = kmul(ToReal(0.5),kadd(JacPDstandard4th23gyz,kadd(JacPDstandard4th32gyz,kmadd(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmadd(gzzL,kmul(gamma322,gamma333),kmul(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2.),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma132,gamma132),ToReal(2.)),kmadd(gyyL,kmul(kmul(gamma232,gamma232),ToReal(2.)),kmadd(gzzL,kmul(kmul(gamma332,gamma332),ToReal(2.)),kmul(kmadd(gyzL,kmul(gamma232,gamma332),kmul(gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))),ToReal(4.))))),JacPDstandard4th33gyy),JacPDstandard4th22gzz))))); - CCTK_REAL_VEC R2331 = - kmadd(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmsub(ToReal(0.5),kadd(JacPDstandard4th12gzz,kadd(JacPDstandard4th33gxy,kmadd(kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gamma132,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))),ToReal(-2.),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmul(gzzL,kmul(gamma321,gamma333))))),ToReal(2.),kmul(gyyL,kmadd(gamma231,kmul(gamma232,ToReal(-2.)),kmul(gamma221,kmul(gamma233,ToReal(2.)))))),JacPDstandard4th31gyz),JacPDstandard4th23gxz)))),kmul(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))))); - - CCTK_REAL_VEC R2332 = - kmadd(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmadd(gamma132,kmul(kmadd(gxyL,gamma232,kmul(gxzL,gamma332)),ToReal(-2.)),kmsub(ToReal(0.5),kadd(JacPDstandard4th22gzz,kadd(JacPDstandard4th33gyy,kmadd(gyzL,kmul(gamma232,kmul(gamma332,ToReal(-4.))),kmadd(gyyL,kmul(kmul(gamma232,gamma232),ToReal(-2.)),kmadd(gzzL,kmul(kmul(gamma332,gamma332),ToReal(-2.)),ksub(kmsub(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmul(gzzL,kmul(gamma322,gamma333))))))),ToReal(2.),JacPDstandard4th32gyz),JacPDstandard4th23gyz)))))),kmul(gxxL,kmul(gamma132,gamma132))))); - - CCTK_REAL_VEC R2333 = - kmul(ksub(JacPDstandard4th32gzz,JacPDstandard4th23gzz),ToReal(0.5)); - - CCTK_REAL_VEC R3111 = - kmul(ksub(JacPDstandard4th13gxx,JacPDstandard4th31gxx),ToReal(0.5)); - - CCTK_REAL_VEC R3112 = - kmadd(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gyz,kadd(JacPDstandard4th23gxx,kmadd(kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gamma131,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))),ToReal(-2.),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmul(gzzL,kmul(gamma311,gamma332))))),ToReal(2.),kmul(gyyL,kmadd(gamma221,kmul(gamma231,ToReal(-2.)),kmul(gamma211,kmul(gamma232,ToReal(2.)))))),JacPDstandard4th31gxy),JacPDstandard4th12gxz)))),kmul(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); - - CCTK_REAL_VEC R3113 = - kmadd(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmadd(gamma131,kmul(kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),ToReal(-2.)),kmsub(ToReal(0.5),kadd(JacPDstandard4th11gzz,kadd(JacPDstandard4th33gxx,kmadd(gyzL,kmul(gamma231,kmul(gamma331,ToReal(-4.))),kmadd(gyyL,kmul(kmul(gamma231,gamma231),ToReal(-2.)),kmadd(gzzL,kmul(kmul(gamma331,gamma331),ToReal(-2.)),ksub(kmsub(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmul(gzzL,kmul(gamma311,gamma333))))))),ToReal(2.),JacPDstandard4th31gxz),JacPDstandard4th13gxz)))))),kmul(gxxL,kmul(gamma131,gamma131))))); - - CCTK_REAL_VEC R3121 = - kmadd(gamma121,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard4th13gxy,kadd(JacPDstandard4th21gxz,kmadd(kmadd(gyyL,kmul(gamma211,gamma232),kmadd(gyzL,kmul(gamma232,gamma311),kmadd(gamma132,kmadd(gxyL,gamma211,kmul(gxzL,gamma311)),kmadd(gyzL,kmul(gamma211,gamma332),kmadd(gzzL,kmul(gamma311,gamma332),kmul(gamma111,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2.),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma221),kmadd(gyyL,kmul(gamma221,gamma231),kmadd(gxzL,kmul(gamma131,gamma321),kmadd(gyzL,kmul(gamma231,gamma321),kmadd(gyzL,kmul(gamma221,gamma331),kmul(gzzL,kmul(gamma321,gamma331))))))),ToReal(2.),JacPDstandard4th32gxx),JacPDstandard4th11gyz)))))); - - CCTK_REAL_VEC R3122 = - kmul(kadd(JacPDstandard4th21gyz,ksub(JacPDstandard4th23gxy,kadd(JacPDstandard4th32gxy,JacPDstandard4th12gyz))),ToReal(0.5)); - - CCTK_REAL_VEC R3123 = - kmadd(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmsub(ToReal(0.5),kadd(JacPDstandard4th21gzz,kadd(JacPDstandard4th33gxy,kmadd(kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gamma132,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))),ToReal(-2.),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmul(gzzL,kmul(gamma321,gamma333))))),ToReal(2.),kmul(gyyL,kmadd(gamma231,kmul(gamma232,ToReal(-2.)),kmul(gamma221,kmul(gamma233,ToReal(2.)))))),JacPDstandard4th32gxz),JacPDstandard4th13gyz)))),kmul(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))))); - - CCTK_REAL_VEC R3131 = - kmul(ToReal(0.5),kadd(JacPDstandard4th13gxz,kadd(JacPDstandard4th31gxz,kmadd(kmadd(gxyL,kmul(gamma133,gamma211),kmadd(gyyL,kmul(gamma211,gamma233),kmadd(gxzL,kmul(gamma133,gamma311),kmadd(gyzL,kmul(gamma233,gamma311),kmadd(gyzL,kmul(gamma211,gamma333),kmadd(gzzL,kmul(gamma311,gamma333),kmul(gamma111,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2.),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma131,gamma131),ToReal(2.)),kmadd(gyyL,kmul(kmul(gamma231,gamma231),ToReal(2.)),kmadd(gzzL,kmul(kmul(gamma331,gamma331),ToReal(2.)),kmul(kmadd(gyzL,kmul(gamma231,gamma331),kmul(gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))),ToReal(4.))))),JacPDstandard4th33gxx),JacPDstandard4th11gzz))))); - - CCTK_REAL_VEC R3132 = - kmadd(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmul(ToReal(0.5),kadd(JacPDstandard4th23gxz,kadd(JacPDstandard4th31gyz,kmadd(kmadd(gyyL,kmul(gamma221,gamma233),kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmadd(gzzL,kmul(gamma321,gamma333),kmul(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))))))))),ToReal(-2.),ksub(kmsub(kmadd(gxyL,kmul(gamma132,gamma231),kmadd(gyyL,kmul(gamma231,gamma232),kmadd(gxzL,kmul(gamma132,gamma331),kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))))),ToReal(2.),JacPDstandard4th33gxy),JacPDstandard4th12gzz)))))); - - CCTK_REAL_VEC R3133 = - kmul(ksub(JacPDstandard4th31gzz,JacPDstandard4th13gzz),ToReal(0.5)); - - CCTK_REAL_VEC R3211 = - kmul(kadd(JacPDstandard4th12gxz,ksub(JacPDstandard4th13gxy,kadd(JacPDstandard4th31gxy,JacPDstandard4th21gxz))),ToReal(0.5)); - - CCTK_REAL_VEC R3212 = - kmadd(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmsub(ToReal(0.5),kadd(JacPDstandard4th12gyz,kadd(JacPDstandard4th23gxy,kmadd(kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gamma131,kmadd(gxyL,gamma222,kmul(gxzL,gamma322)),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))),ToReal(-2.),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmul(gzzL,kmul(gamma321,gamma332))))),ToReal(2.),kmul(gyyL,kmadd(gamma222,kmul(gamma231,ToReal(-2.)),kmul(gamma221,kmul(gamma232,ToReal(2.)))))),JacPDstandard4th31gyy),JacPDstandard4th22gxz)))),kmul(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)))))); - - CCTK_REAL_VEC R3213 = - kmadd(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmsub(ToReal(0.5),kadd(JacPDstandard4th12gzz,kadd(JacPDstandard4th33gxy,kmadd(kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gamma132,kmadd(gxyL,gamma231,kmul(gxzL,gamma331)),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))),ToReal(-2.),ksub(ksub(kmadd(kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmul(gzzL,kmul(gamma321,gamma333))))),ToReal(2.),kmul(gyyL,kmadd(gamma231,kmul(gamma232,ToReal(-2.)),kmul(gamma221,kmul(gamma233,ToReal(2.)))))),JacPDstandard4th31gyz),JacPDstandard4th23gxz)))),kmul(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))))); - - CCTK_REAL_VEC R3221 = - kmadd(gamma122,kmadd(gxxL,gamma131,kmadd(gxyL,gamma231,kmul(gxzL,gamma331))),kmul(ToReal(0.5),kadd(JacPDstandard4th13gyy,kadd(JacPDstandard4th22gxz,kmadd(kmadd(gyyL,kmul(gamma221,gamma232),kmadd(gyzL,kmul(gamma232,gamma321),kmadd(gamma132,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma332),kmadd(gzzL,kmul(gamma321,gamma332),kmul(gamma121,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))))))))),ToReal(-2.),ksub(kmsub(kmadd(gxyL,kmul(gamma131,gamma222),kmadd(gyyL,kmul(gamma222,gamma231),kmadd(gxzL,kmul(gamma131,gamma322),kmadd(gyzL,kmul(gamma231,gamma322),kmadd(gyzL,kmul(gamma222,gamma331),kmul(gzzL,kmul(gamma322,gamma331))))))),ToReal(2.),JacPDstandard4th32gxy),JacPDstandard4th21gyz)))))); - - CCTK_REAL_VEC R3222 = - kmul(ksub(JacPDstandard4th23gyy,JacPDstandard4th32gyy),ToReal(0.5)); - - CCTK_REAL_VEC R3223 = - kmadd(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))),kmadd(gamma132,kmul(kmadd(gxyL,gamma232,kmul(gxzL,gamma332)),ToReal(-2.)),kmsub(ToReal(0.5),kadd(JacPDstandard4th22gzz,kadd(JacPDstandard4th33gyy,kmadd(gyzL,kmul(gamma232,kmul(gamma332,ToReal(-4.))),kmadd(gyyL,kmul(kmul(gamma232,gamma232),ToReal(-2.)),kmadd(gzzL,kmul(kmul(gamma332,gamma332),ToReal(-2.)),ksub(kmsub(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmul(gzzL,kmul(gamma322,gamma333))))))),ToReal(2.),JacPDstandard4th32gyz),JacPDstandard4th23gyz)))))),kmul(gxxL,kmul(gamma132,gamma132))))); - - CCTK_REAL_VEC R3231 = - kmadd(gamma131,kmadd(gxxL,gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332))),kmul(ToReal(0.5),kadd(JacPDstandard4th13gyz,kadd(JacPDstandard4th32gxz,kmadd(kmadd(gyyL,kmul(gamma221,gamma233),kmadd(gyzL,kmul(gamma233,gamma321),kmadd(gamma133,kmadd(gxyL,gamma221,kmul(gxzL,gamma321)),kmadd(gyzL,kmul(gamma221,gamma333),kmadd(gzzL,kmul(gamma321,gamma333),kmul(gamma121,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333))))))))),ToReal(-2.),ksub(kmsub(kmadd(gxyL,kmul(gamma132,gamma231),kmadd(gyyL,kmul(gamma231,gamma232),kmadd(gxzL,kmul(gamma132,gamma331),kmadd(gyzL,kmul(gamma232,gamma331),kmadd(gyzL,kmul(gamma231,gamma332),kmul(gzzL,kmul(gamma331,gamma332))))))),ToReal(2.),JacPDstandard4th33gxy),JacPDstandard4th21gzz)))))); - - CCTK_REAL_VEC R3232 = - kmul(ToReal(0.5),kadd(JacPDstandard4th23gyz,kadd(JacPDstandard4th32gyz,kmadd(kmadd(gxyL,kmul(gamma133,gamma222),kmadd(gyyL,kmul(gamma222,gamma233),kmadd(gxzL,kmul(gamma133,gamma322),kmadd(gyzL,kmul(gamma233,gamma322),kmadd(gyzL,kmul(gamma222,gamma333),kmadd(gzzL,kmul(gamma322,gamma333),kmul(gamma122,kmadd(gxxL,gamma133,kmadd(gxyL,gamma233,kmul(gxzL,gamma333)))))))))),ToReal(-2.),ksub(ksub(kmadd(gxxL,kmul(kmul(gamma132,gamma132),ToReal(2.)),kmadd(gyyL,kmul(kmul(gamma232,gamma232),ToReal(2.)),kmadd(gzzL,kmul(kmul(gamma332,gamma332),ToReal(2.)),kmul(kmadd(gyzL,kmul(gamma232,gamma332),kmul(gamma132,kmadd(gxyL,gamma232,kmul(gxzL,gamma332)))),ToReal(4.))))),JacPDstandard4th33gyy),JacPDstandard4th22gzz))))); - - CCTK_REAL_VEC R3233 = - kmul(ksub(JacPDstandard4th32gzz,JacPDstandard4th23gzz),ToReal(0.5)); - - CCTK_REAL_VEC R3311 = - ksub(JacPDstandard4th13gxz,JacPDstandard4th31gxz); - - CCTK_REAL_VEC R3312 = - kmul(kadd(JacPDstandard4th13gyz,ksub(JacPDstandard4th23gxz,kadd(JacPDstandard4th32gxz,JacPDstandard4th31gyz))),ToReal(0.5)); - - CCTK_REAL_VEC R3313 = - kmul(ksub(JacPDstandard4th13gzz,JacPDstandard4th31gzz),ToReal(0.5)); - - CCTK_REAL_VEC R3321 = - kmul(kadd(JacPDstandard4th13gyz,ksub(JacPDstandard4th23gxz,kadd(JacPDstandard4th32gxz,JacPDstandard4th31gyz))),ToReal(0.5)); - - CCTK_REAL_VEC R3322 = - ksub(JacPDstandard4th23gyz,JacPDstandard4th32gyz); - - CCTK_REAL_VEC R3323 = - kmul(ksub(JacPDstandard4th23gzz,JacPDstandard4th32gzz),ToReal(0.5)); - - CCTK_REAL_VEC R3331 = - kmul(ksub(JacPDstandard4th13gzz,JacPDstandard4th31gzz),ToReal(0.5)); - - CCTK_REAL_VEC R3332 = - kmul(ksub(JacPDstandard4th23gzz,JacPDstandard4th32gzz),ToReal(0.5)); - - CCTK_REAL_VEC R3333 = ToReal(0.); - - CCTK_REAL_VEC R4p1111 = R1111; - - CCTK_REAL_VEC R4p1112 = R1112; - - CCTK_REAL_VEC R4p1113 = R1113; - - CCTK_REAL_VEC R4p1121 = R1121; - - CCTK_REAL_VEC R4p1122 = R1122; - - CCTK_REAL_VEC R4p1123 = R1123; - - CCTK_REAL_VEC R4p1131 = R1131; - - CCTK_REAL_VEC R4p1132 = R1132; - - CCTK_REAL_VEC R4p1133 = R1133; - - CCTK_REAL_VEC R4p1211 = R1211; - CCTK_REAL_VEC R4p1212 = kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R1212)); CCTK_REAL_VEC R4p1213 = kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)); - CCTK_REAL_VEC R4p1221 = - kadd(R1221,kmsub(kxyL,kxyL,kmul(kxxL,kyyL))); - - CCTK_REAL_VEC R4p1222 = R1222; - CCTK_REAL_VEC R4p1223 = kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R1223)); - CCTK_REAL_VEC R4p1231 = - kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R1231)); - - CCTK_REAL_VEC R4p1232 = - kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R1232)); - - CCTK_REAL_VEC R4p1233 = R1233; - - CCTK_REAL_VEC R4p1311 = R1311; - - CCTK_REAL_VEC R4p1312 = - kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1312)); - CCTK_REAL_VEC R4p1313 = kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R1313)); - CCTK_REAL_VEC R4p1321 = - kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R1321)); - - CCTK_REAL_VEC R4p1322 = R1322; - CCTK_REAL_VEC R4p1323 = kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323)); - CCTK_REAL_VEC R4p1331 = - kadd(R1331,kmsub(kxzL,kxzL,kmul(kxxL,kzzL))); - - CCTK_REAL_VEC R4p1332 = - kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R1332)); - - CCTK_REAL_VEC R4p1333 = R1333; - - CCTK_REAL_VEC R4p2111 = R2111; - - CCTK_REAL_VEC R4p2112 = - kadd(R2112,kmsub(kxyL,kxyL,kmul(kxxL,kyyL))); - - CCTK_REAL_VEC R4p2113 = - kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R2113)); - - CCTK_REAL_VEC R4p2121 = - kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R2121)); - - CCTK_REAL_VEC R4p2122 = R2122; - - CCTK_REAL_VEC R4p2123 = - kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R2123)); - - CCTK_REAL_VEC R4p2131 = - kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R2131)); - - CCTK_REAL_VEC R4p2132 = - kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R2132)); - - CCTK_REAL_VEC R4p2133 = R2133; - - CCTK_REAL_VEC R4p2211 = R2211; - - CCTK_REAL_VEC R4p2212 = R2212; - - CCTK_REAL_VEC R4p2213 = R2213; - - CCTK_REAL_VEC R4p2221 = R2221; - - CCTK_REAL_VEC R4p2222 = R2222; - - CCTK_REAL_VEC R4p2223 = R2223; - - CCTK_REAL_VEC R4p2231 = R2231; - - CCTK_REAL_VEC R4p2232 = R2232; - - CCTK_REAL_VEC R4p2233 = R2233; - - CCTK_REAL_VEC R4p2311 = R2311; - - CCTK_REAL_VEC R4p2312 = - kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R2312)); - - CCTK_REAL_VEC R4p2313 = - kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R2313)); - - CCTK_REAL_VEC R4p2321 = - kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R2321)); - - CCTK_REAL_VEC R4p2322 = R2322; - CCTK_REAL_VEC R4p2323 = kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R2323)); - CCTK_REAL_VEC R4p2331 = - kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R2331)); - - CCTK_REAL_VEC R4p2332 = - kadd(R2332,kmsub(kyzL,kyzL,kmul(kyyL,kzzL))); - - CCTK_REAL_VEC R4p2333 = R2333; - - CCTK_REAL_VEC R4p3111 = R3111; - - CCTK_REAL_VEC R4p3112 = - kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R3112)); - - CCTK_REAL_VEC R4p3113 = - kadd(R3113,kmsub(kxzL,kxzL,kmul(kxxL,kzzL))); - - CCTK_REAL_VEC R4p3121 = - kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R3121)); - - CCTK_REAL_VEC R4p3122 = R3122; - - CCTK_REAL_VEC R4p3123 = - kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R3123)); - - CCTK_REAL_VEC R4p3131 = - kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R3131)); - - CCTK_REAL_VEC R4p3132 = - kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R3132)); - - CCTK_REAL_VEC R4p3133 = R3133; - - CCTK_REAL_VEC R4p3211 = R3211; - - CCTK_REAL_VEC R4p3212 = - kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R3212)); - - CCTK_REAL_VEC R4p3213 = - kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R3213)); - - CCTK_REAL_VEC R4p3221 = - kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R3221)); - - CCTK_REAL_VEC R4p3222 = R3222; - - CCTK_REAL_VEC R4p3223 = - kadd(R3223,kmsub(kyzL,kyzL,kmul(kyyL,kzzL))); - - CCTK_REAL_VEC R4p3231 = - kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R3231)); - - CCTK_REAL_VEC R4p3232 = - kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R3232)); - - CCTK_REAL_VEC R4p3233 = R3233; - - CCTK_REAL_VEC R4p3311 = R3311; - - CCTK_REAL_VEC R4p3312 = R3312; - - CCTK_REAL_VEC R4p3313 = R3313; - - CCTK_REAL_VEC R4p3321 = R3321; - - CCTK_REAL_VEC R4p3322 = R3322; - - CCTK_REAL_VEC R4p3323 = R3323; - - CCTK_REAL_VEC R4p3331 = R3331; - - CCTK_REAL_VEC R4p3332 = R3332; - - CCTK_REAL_VEC R4p3333 = R3333; - CCTK_REAL_VEC Ro111 = ToReal(0.); CCTK_REAL_VEC Ro112 = @@ -1828,37 +1264,37 @@ static void psi4_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL_VEC Ro333 = ToReal(0.); CCTK_REAL_VEC Rojo11 = - kmadd(gInv11,R1111,kmadd(gInv12,R1112,kmadd(gInv13,R1113,kmadd(gInv21,R1211,kmadd(gInv31,R1311,kmadd(gInv23,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)),kmadd(gInv32,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1312)),kmadd(gInv22,kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R1212)),kmul(gInv33,kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R1313))))))))))); + kmadd(kadd(gInv23,gInv32),kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)),kmadd(gInv22,kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R1212)),kmul(gInv33,kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R1313))))); CCTK_REAL_VEC Rojo12 = - kmadd(gInv11,R1121,kmadd(gInv21,R1221,kmadd(gInv22,R1222,kmadd(gInv23,R1223,kmadd(gInv31,R1321,kmadd(gInv12,kadd(R1122,kmsub(kxyL,kxyL,kmul(kxxL,kyyL))),kmadd(gInv32,kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R1322)),kmadd(gInv13,kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R1123)),kmul(gInv33,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323))))))))))); + kmadd(gInv23,R1223,knmsub(gInv21,R1212,knmsub(gInv31,R1213,kmadd(gInv12,kmsub(kxyL,kxyL,kmul(kxxL,kyyL)),kmadd(gInv32,kmsub(kxyL,kyzL,kmul(kxzL,kyyL)),kmadd(gInv13,kmsub(kxyL,kxzL,kmul(kxxL,kyzL)),kmul(gInv33,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323))))))))); CCTK_REAL_VEC Rojo13 = - kmadd(gInv11,R1131,kmadd(gInv21,R1231,kmadd(gInv31,R1331,kmadd(gInv32,R1332,kmadd(gInv33,R1333,kmadd(gInv12,kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R1132)),kmadd(gInv22,kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R1232)),kmadd(gInv13,kadd(R1133,kmsub(kxzL,kxzL,kmul(kxxL,kzzL))),kmul(gInv23,kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R1233))))))))))); + kmadd(kxzL,kmul(kyzL,gInv23),kmadd(gInv13,kmul(kxzL,kxzL),knmsub(gInv21,R1213,knmsub(gInv31,R1313,knmsub(gInv32,R1323,kmadd(gInv12,kmsub(kxyL,kxzL,kmul(kxxL,kyzL)),kmsub(gInv22,kmsub(kxzL,kyyL,kmadd(kxyL,kyzL,R1223)),kmul(kzzL,kmadd(kxyL,gInv23,kmul(kxxL,gInv13)))))))))); CCTK_REAL_VEC Rojo21 = - kmadd(gInv11,R2111,kmadd(gInv12,R2112,kmadd(gInv13,R2113,kmadd(gInv22,R2212,kmadd(gInv32,R2312,kmadd(gInv21,kadd(R2211,kmsub(kxyL,kxyL,kmul(kxxL,kyyL))),kmadd(gInv23,kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R2213)),kmadd(gInv31,kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R2311)),kmul(gInv33,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R2313))))))))))); + kmadd(gInv32,R1223,knmsub(gInv12,R1212,knmsub(gInv13,R1213,kmadd(gInv21,kmsub(kxyL,kxyL,kmul(kxxL,kyyL)),kmadd(gInv23,kmsub(kxyL,kyzL,kmul(kxzL,kyyL)),kmadd(gInv31,kmsub(kxyL,kxzL,kmul(kxxL,kyzL)),kmul(gInv33,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323))))))))); CCTK_REAL_VEC Rojo22 = - kmadd(gInv12,R2122,kmadd(gInv21,R2221,kmadd(gInv22,R2222,kmadd(gInv23,R2223,kmadd(gInv32,R2322,kmadd(gInv13,kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R2123)),kmadd(gInv31,kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R2321)),kmadd(gInv11,kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R2121)),kmul(gInv33,kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R2323))))))))))); + kmadd(kadd(gInv13,gInv31),kmsub(kxzL,kyyL,kmadd(kxyL,kyzL,R1223)),kmadd(gInv11,kmadd(kxxL,kyyL,knmsub(kxyL,kxyL,R1212)),kmul(gInv33,kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R2323))))); CCTK_REAL_VEC Rojo23 = - kmadd(gInv12,R2132,kmadd(gInv22,R2232,kmadd(gInv31,R2331,kmadd(gInv32,R2332,kmadd(gInv33,R2333,kmadd(gInv11,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R2131)),kmadd(gInv21,kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R2231)),kmadd(gInv13,kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R2133)),kmul(gInv23,kadd(R2233,kmsub(kyzL,kyzL,kmul(kyyL,kzzL)))))))))))); + kmadd(gInv12,R1223,knmsub(gInv31,R1323,knmsub(gInv32,R2323,kmadd(gInv11,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)),kmadd(gInv21,kmsub(kxyL,kyzL,kmul(kxzL,kyyL)),kmadd(gInv13,kmsub(kxzL,kyzL,kmul(kxyL,kzzL)),kmul(gInv23,kmsub(kyzL,kyzL,kmul(kyyL,kzzL))))))))); CCTK_REAL_VEC Rojo31 = - kmadd(gInv11,R3111,kmadd(gInv12,R3112,kmadd(gInv13,R3113,kmadd(gInv23,R3213,kmadd(gInv33,R3313,kmadd(gInv21,kmadd(kxyL,kxzL,knmsub(kxxL,kyzL,R3211)),kmadd(gInv22,kmadd(kxzL,kyyL,knmsub(kxyL,kyzL,R3212)),kmadd(gInv31,kadd(R3311,kmsub(kxzL,kxzL,kmul(kxxL,kzzL))),kmul(gInv32,kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R3312))))))))))); + kmadd(kxzL,kmul(kyzL,gInv32),kmadd(gInv31,kmul(kxzL,kxzL),knmsub(gInv12,R1213,knmsub(gInv13,R1313,knmsub(gInv23,R1323,kmadd(gInv21,kmsub(kxyL,kxzL,kmul(kxxL,kyzL)),kmsub(gInv22,kmsub(kxzL,kyyL,kmadd(kxyL,kyzL,R1223)),kmul(kzzL,kmadd(kxyL,gInv32,kmul(kxxL,gInv31)))))))))); CCTK_REAL_VEC Rojo32 = - kmadd(gInv13,R3123,kmadd(gInv21,R3221,kmadd(gInv22,R3222,kmadd(gInv23,R3223,kmadd(gInv33,R3323,kmadd(gInv11,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R3121)),kmadd(gInv12,kmadd(kxyL,kyzL,knmsub(kxzL,kyyL,R3122)),kmadd(gInv31,kmadd(kxzL,kyzL,knmsub(kxyL,kzzL,R3321)),kmul(gInv32,kadd(R3322,kmsub(kyzL,kyzL,kmul(kyyL,kzzL)))))))))))); + kmadd(gInv21,R1223,knmsub(gInv13,R1323,knmsub(gInv23,R2323,kmadd(gInv11,kmadd(kxxL,kyzL,knmsub(kxyL,kxzL,R1213)),kmadd(gInv12,kmsub(kxyL,kyzL,kmul(kxzL,kyyL)),kmadd(gInv31,kmsub(kxzL,kyzL,kmul(kxyL,kzzL)),kmul(gInv32,kmsub(kyzL,kyzL,kmul(kyyL,kzzL))))))))); CCTK_REAL_VEC Rojo33 = - kmadd(gInv13,R3133,kmadd(gInv23,R3233,kmadd(gInv31,R3331,kmadd(gInv32,R3332,kmadd(gInv33,R3333,kmadd(gInv12,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R3132)),kmadd(gInv21,kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R3231)),kmadd(gInv11,kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R3131)),kmul(gInv22,kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R3232))))))))))); + kmadd(kadd(gInv12,gInv21),kmadd(kxyL,kzzL,knmsub(kxzL,kyzL,R1323)),kmadd(gInv11,kmadd(kxxL,kzzL,knmsub(kxzL,kxzL,R1313)),kmul(gInv22,kmadd(kyyL,kzzL,knmsub(kyzL,kyzL,R2323))))); CCTK_REAL_VEC Psi4rL = - kmadd(R4p1112,kmul(rmbar1,kmul(rmbar2,kmul(n1,n1))),kmadd(R4p1211,kmul(rmbar1,kmul(rmbar2,kmul(n1,n1))),kmadd(R4p1113,kmul(rmbar1,kmul(rmbar3,kmul(n1,n1))),kmadd(R4p1311,kmul(rmbar1,kmul(rmbar3,kmul(n1,n1))),kmadd(R4p1213,kmul(rmbar2,kmul(rmbar3,kmul(n1,n1))),kmadd(R4p1312,kmul(rmbar2,kmul(rmbar3,kmul(n1,n1))),kmadd(R4p2122,kmul(rmbar1,kmul(rmbar2,kmul(n2,n2))),kmadd(R4p2221,kmul(rmbar1,kmul(rmbar2,kmul(n2,n2))),kmadd(R4p2123,kmul(rmbar1,kmul(rmbar3,kmul(n2,n2))),kmadd(R4p2321,kmul(rmbar1,kmul(rmbar3,kmul(n2,n2))),kmadd(R4p2223,kmul(rmbar2,kmul(rmbar3,kmul(n2,n2))),kmadd(R4p2322,kmul(rmbar2,kmul(rmbar3,kmul(n2,n2))),kmadd(R4p3132,kmul(rmbar1,kmul(rmbar2,kmul(n3,n3))),kmadd(R4p3231,kmul(rmbar1,kmul(rmbar2,kmul(n3,n3))),kmadd(R4p3133,kmul(rmbar1,kmul(rmbar3,kmul(n3,n3))),kmadd(R4p3331,kmul(rmbar1,kmul(rmbar3,kmul(n3,n3))),kmadd(R4p3233,kmul(rmbar2,kmul(rmbar3,kmul(n3,n3))),kmadd(R4p3332,kmul(rmbar2,kmul(rmbar3,kmul(n3,n3))),kmadd(rmbar1,kmul(rmbar2,kmul(Rojo12,kmul(nn,nn))),kmadd(rmbar1,kmul(rmbar3,kmul(Rojo13,kmul(nn,nn))),kmadd(rmbar1,kmul(rmbar2,kmul(Rojo21,kmul(nn,nn))),kmadd(rmbar2,kmul(rmbar3,kmul(Rojo23,kmul(nn,nn))),kmadd(rmbar1,kmul(rmbar3,kmul(Rojo31,kmul(nn,nn))),kmadd(rmbar2,kmul(rmbar3,kmul(Rojo32,kmul(nn,nn))),kmadd(n1,kmul(n2,kmul(R4p1121,kmul(rmbar1,rmbar1))),kmadd(n1,kmul(n3,kmul(R4p1131,kmul(rmbar1,rmbar1))),kmadd(n1,kmul(n2,kmul(R4p2111,kmul(rmbar1,rmbar1))),kmadd(n2,kmul(n3,kmul(R4p2131,kmul(rmbar1,rmbar1))),kmadd(n1,kmul(n3,kmul(R4p3111,kmul(rmbar1,rmbar1))),kmadd(n2,kmul(n3,kmul(R4p3121,kmul(rmbar1,rmbar1))),kmadd(R4p1111,kmul(kmul(n1,n1),kmul(rmbar1,rmbar1)),kmadd(R4p2121,kmul(kmul(n2,n2),kmul(rmbar1,rmbar1)),kmadd(R4p3131,kmul(kmul(n3,n3),kmul(rmbar1,rmbar1)),kmadd(Rojo11,kmul(kmul(nn,nn),kmul(rmbar1,rmbar1)),kmadd(n1,kmul(n2,kmul(R4p1222,kmul(rmbar2,rmbar2))),kmadd(n1,kmul(n3,kmul(R4p1232,kmul(rmbar2,rmbar2))),kmadd(n1,kmul(n2,kmul(R4p2212,kmul(rmbar2,rmbar2))),kmadd(n2,kmul(n3,kmul(R4p2232,kmul(rmbar2,rmbar2))),kmadd(n1,kmul(n3,kmul(R4p3212,kmul(rmbar2,rmbar2))),kmadd(n2,kmul(n3,kmul(R4p3222,kmul(rmbar2,rmbar2))),kmadd(R4p1212,kmul(kmul(n1,n1),kmul(rmbar2,rmbar2)),kmadd(R4p2222,kmul(kmul(n2,n2),kmul(rmbar2,rmbar2)),kmadd(R4p3232,kmul(kmul(n3,n3),kmul(rmbar2,rmbar2)),kmadd(Rojo22,kmul(kmul(nn,nn),kmul(rmbar2,rmbar2)),kmadd(n1,kmul(n2,kmul(R4p1323,kmul(rmbar3,rmbar3))),kmadd(n1,kmul(n3,kmul(R4p1333,kmul(rmbar3,rmbar3))),kmadd(n1,kmul(n2,kmul(R4p2313,kmul(rmbar3,rmbar3))),kmadd(n2,kmul(n3,kmul(R4p2333,kmul(rmbar3,rmbar3))),kmadd(n1,kmul(n3,kmul(R4p3313,kmul(rmbar3,rmbar3))),kmadd(n2,kmul(n3,kmul(R4p3323,kmul(rmbar3,rmbar3))),kmadd(n2,kmul(nn,kmul(Ro323,kmul(kmul(imbar3,imbar3),ToReal(-2.)))),kmadd(n3,kmul(nn,kmul(Ro333,kmul(kmul(imbar3,imbar3),ToReal(-2.)))),knmsub(n3,kmul(kmul(imbar3,imbar3),kmadd(n1,R4p3313,kmul(n2,R4p2333))),kmadd(kmul(imbar3,imbar3),kmsub(n1,kmul(nn,kmul(Ro313,ToReal(-2.))),kmul(n2,kmul(n3,R4p3323))),kmadd(R4p1313,kmsub(kmul(n1,n1),kmul(rmbar3,rmbar3),kmul(kmul(imbar3,imbar3),kmul(n1,n1))),kmadd(R4p2323,kmsub(kmul(n2,n2),kmul(rmbar3,rmbar3),kmul(kmul(imbar3,imbar3),kmul(n2,n2))),kmadd(R4p3333,kmsub(kmul(n3,n3),kmul(rmbar3,rmbar3),kmul(kmul(imbar3,imbar3),kmul(n3,n3))),kmadd(Rojo33,kmsub(kmul(nn,nn),kmul(rmbar3,rmbar3),kmul(kmul(imbar3,imbar3),kmul(nn,nn))),kmadd(n1,kmul(nn,kmul(Ro111,kmul(kmul(rmbar1,rmbar1),ToReal(2.)))),kmadd(n2,kmul(nn,kmul(Ro121,kmul(kmul(rmbar1,rmbar1),ToReal(2.)))),kmadd(n3,kmul(nn,kmul(Ro131,kmul(kmul(rmbar1,rmbar1),ToReal(2.)))),kmadd(n1,kmul(nn,kmul(Ro212,kmul(kmul(rmbar2,rmbar2),ToReal(2.)))),kmadd(n2,kmul(nn,kmul(Ro222,kmul(kmul(rmbar2,rmbar2),ToReal(2.)))),kmadd(n3,kmul(nn,kmul(Ro232,kmul(kmul(rmbar2,rmbar2),ToReal(2.)))),kmadd(n1,kmul(nn,kmul(Ro313,kmul(kmul(rmbar3,rmbar3),ToReal(2.)))),kmadd(n2,kmul(nn,kmul(Ro323,kmul(kmul(rmbar3,rmbar3),ToReal(2.)))),kmadd(n3,kmul(nn,kmul(Ro333,kmul(kmul(rmbar3,rmbar3),ToReal(2.)))),kmadd(n1,kmadd(n3,kmul(R4p1132,kmul(rmbar1,rmbar2)),kmadd(n2,kmul(R4p1221,kmul(rmbar1,rmbar2)),kmadd(n3,kmul(R4p1231,kmul(rmbar1,rmbar2)),kmadd(n2,kmul(R4p1123,kmul(rmbar1,rmbar3)),kmadd(n3,kmul(R4p1331,kmul(rmbar1,rmbar3)),kmadd(n2,kmul(R4p2113,kmul(rmbar1,rmbar3)),kmadd(n3,kmul(R4p1233,kmul(rmbar2,rmbar3)),kmadd(n2,kmul(R4p1322,kmul(rmbar2,rmbar3)),kmadd(n3,kmul(R4p1332,kmul(rmbar2,rmbar3)),kmadd(n2,kmul(R4p2213,kmul(rmbar2,rmbar3)),kmadd(n2,kmsub(R4p1122,kmul(rmbar1,rmbar2),kmul(R4p1323,kmul(imbar3,imbar3))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro112,ToReal(2.)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro113,ToReal(2.)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro211,ToReal(2.)))),kmadd(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro213,ToReal(2.)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro311,ToReal(2.)))),kmul(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro312,ToReal(2.))))))))))))))))))))),kmadd(n2,kmadd(n3,kmul(R4p2132,kmul(rmbar1,rmbar2)),kmadd(n1,kmul(R4p2211,kmul(rmbar1,rmbar2)),kmadd(n3,kmul(R4p2231,kmul(rmbar1,rmbar2)),kmadd(n1,kmul(R4p1321,kmul(rmbar1,rmbar3)),kmadd(n3,kmul(R4p2133,kmul(rmbar1,rmbar3)),kmadd(n1,kmul(R4p2311,kmul(rmbar1,rmbar3)),kmadd(n3,kmul(R4p3321,kmul(rmbar1,rmbar3)),kmadd(n1,kmul(R4p1223,kmul(rmbar2,rmbar3)),kmadd(n3,kmul(R4p2233,kmul(rmbar2,rmbar3)),kmadd(n1,kmul(R4p2312,kmul(rmbar2,rmbar3)),kmadd(n1,kmsub(R4p2112,kmul(rmbar1,rmbar2),kmul(R4p2313,kmul(imbar3,imbar3))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro122,ToReal(2.)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro123,ToReal(2.)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro221,ToReal(2.)))),kmadd(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro223,ToReal(2.)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro321,ToReal(2.)))),kmul(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro322,ToReal(2.))))))))))))))))))))),kmsub(n3,kmadd(n2,kmul(R4p3122,kmul(rmbar1,rmbar2)),kmadd(n1,kmul(R4p3211,kmul(rmbar1,rmbar2)),kmadd(n2,kmul(R4p3221,kmul(rmbar1,rmbar2)),kmadd(n1,kmul(R4p1133,kmul(rmbar1,rmbar3)),kmadd(n2,kmul(R4p2331,kmul(rmbar1,rmbar3)),kmadd(n1,kmul(R4p3113,kmul(rmbar1,rmbar3)),kmadd(n2,kmul(R4p3123,kmul(rmbar1,rmbar3)),kmadd(n1,kmul(R4p3311,kmul(rmbar1,rmbar3)),kmadd(n2,kmul(R4p2332,kmul(rmbar2,rmbar3)),kmadd(n1,kmul(R4p3213,kmul(rmbar2,rmbar3)),kmadd(n2,kmul(R4p3223,kmul(rmbar2,rmbar3)),kmadd(n1,kmul(R4p3312,kmul(rmbar2,rmbar3)),kmadd(n2,kmul(R4p3322,kmul(rmbar2,rmbar3)),kmadd(n1,kmsub(R4p3112,kmul(rmbar1,rmbar2),kmul(R4p1333,kmul(imbar3,imbar3))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro132,ToReal(2.)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro133,ToReal(2.)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro231,ToReal(2.)))),kmadd(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro233,ToReal(2.)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro331,ToReal(2.)))),kmul(nn,kmul(rmbar2,kmul(rmbar3,kmul(Ro332,ToReal(2.)))))))))))))))))))))))),kmadd(kmul(imbar1,imbar1),kmadd(n2,kmul(n3,R4p2131),kmadd(n2,kmul(n3,R4p3121),kmadd(R4p1111,kmul(n1,n1),kmadd(R4p2121,kmul(n2,n2),kmadd(R4p3131,kmul(n3,n3),kmadd(Rojo11,kmul(nn,nn),kmadd(n2,kmul(nn,kmul(Ro121,ToReal(2.))),kmadd(n3,kmul(nn,kmul(Ro131,ToReal(2.))),kmul(n1,kmadd(n2,kadd(R4p1121,R4p2111),kmadd(n3,kadd(R4p1131,R4p3111),kmul(nn,kmul(Ro111,ToReal(2.)))))))))))))),kmadd(kmul(imbar2,imbar2),kmadd(n2,kmul(n3,R4p2232),kmadd(n2,kmul(n3,R4p3222),kmadd(R4p1212,kmul(n1,n1),kmadd(R4p2222,kmul(n2,n2),kmadd(R4p3232,kmul(n3,n3),kmadd(Rojo22,kmul(nn,nn),kmadd(n2,kmul(nn,kmul(Ro222,ToReal(2.))),kmadd(n3,kmul(nn,kmul(Ro232,ToReal(2.))),kmul(n1,kmadd(n2,kadd(R4p1222,R4p2212),kmadd(n3,kadd(R4p1232,R4p3212),kmul(nn,kmul(Ro212,ToReal(2.)))))))))))))),kmadd(imbar2,kmul(imbar3,kmadd(kadd(R4p1213,R4p1312),kmul(n1,n1),kmadd(kadd(R4p2223,R4p2322),kmul(n2,n2),kmadd(R4p3233,kmul(n3,n3),kmadd(R4p3332,kmul(n3,n3),kmadd(Rojo23,kmul(nn,nn),kmadd(Rojo32,kmul(nn,nn),kmadd(n3,kmul(nn,kmul(Ro233,ToReal(2.))),kmadd(n3,kmul(nn,kmul(Ro332,ToReal(2.))),kmadd(n1,kmadd(n2,kadd(R4p1223,kadd(R4p1322,kadd(R4p2213,R4p2312))),kmadd(n3,kadd(R4p1233,kadd(R4p1332,kadd(R4p3213,R4p3312))),kmul(nn,kmul(kadd(Ro213,Ro312),ToReal(2.))))),kmul(n2,kmadd(n3,kadd(R4p2233,kadd(R4p2332,kadd(R4p3223,R4p3322))),kmul(nn,kmul(kadd(Ro223,Ro322),ToReal(2.))))))))))))))),kmul(imbar1,kmadd(imbar2,kmadd(kadd(R4p1112,R4p1211),kmul(n1,n1),kmadd(kadd(R4p2122,R4p2221),kmul(n2,n2),kmadd(R4p3132,kmul(n3,n3),kmadd(R4p3231,kmul(n3,n3),kmadd(Rojo12,kmul(nn,nn),kmadd(Rojo21,kmul(nn,nn),kmadd(n3,kmul(nn,kmul(Ro132,ToReal(2.))),kmadd(n3,kmul(nn,kmul(Ro231,ToReal(2.))),kmadd(n1,kmadd(n2,kadd(R4p1122,kadd(R4p1221,kadd(R4p2112,R4p2211))),kmadd(n3,kadd(R4p1132,kadd(R4p1231,kadd(R4p3112,R4p3211))),kmul(nn,kmul(kadd(Ro112,Ro211),ToReal(2.))))),kmul(n2,kmadd(n3,kadd(R4p2132,kadd(R4p2231,kadd(R4p3122,R4p3221))),kmul(nn,kmul(kadd(Ro122,Ro221),ToReal(2.)))))))))))))),kmul(imbar3,kmadd(kadd(R4p1113,R4p1311),kmul(n1,n1),kmadd(kadd(R4p2123,R4p2321),kmul(n2,n2),kmadd(R4p3133,kmul(n3,n3),kmadd(R4p3331,kmul(n3,n3),kmadd(Rojo13,kmul(nn,nn),kmadd(Rojo31,kmul(nn,nn),kmadd(n3,kmul(nn,kmul(Ro133,ToReal(2.))),kmadd(n3,kmul(nn,kmul(Ro331,ToReal(2.))),kmadd(n1,kmadd(n2,kadd(R4p1123,kadd(R4p1321,kadd(R4p2113,R4p2311))),kmadd(n3,kadd(R4p1133,kadd(R4p1331,kadd(R4p3113,R4p3311))),kmul(nn,kmul(kadd(Ro113,Ro311),ToReal(2.))))),kmul(n2,kmadd(n3,kadd(R4p2133,kadd(R4p2331,kadd(R4p3123,R4p3321))),kmul(nn,kmul(kadd(Ro123,Ro321),ToReal(2.)))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))); + kmadd(rmbar1,kmul(rmbar2,kmul(Rojo12,kmul(nn,nn))),kmadd(rmbar1,kmul(rmbar3,kmul(Rojo13,kmul(nn,nn))),kmadd(rmbar1,kmul(rmbar2,kmul(Rojo21,kmul(nn,nn))),kmadd(rmbar2,kmul(rmbar3,kmul(Rojo23,kmul(nn,nn))),kmadd(rmbar1,kmul(rmbar3,kmul(Rojo31,kmul(nn,nn))),kmadd(rmbar2,kmul(rmbar3,kmul(Rojo32,kmul(nn,nn))),kmadd(R4p1212,kmul(kmul(n2,n2),kmul(rmbar1,rmbar1)),kmadd(Rojo11,kmul(kmul(nn,nn),kmul(rmbar1,rmbar1)),kmadd(R4p1212,kmul(kmul(n1,n1),kmul(rmbar2,rmbar2)),kmadd(Rojo22,kmul(kmul(nn,nn),kmul(rmbar2,rmbar2)),kmadd(R4p1313,kmul(kmul(n1,n1),kmul(rmbar3,rmbar3)),kmadd(R4p2323,kmul(kmul(n2,n2),kmul(rmbar3,rmbar3)),kmadd(kmadd(n1,kmul(n2,kmul(R4p1212,kmul(rmbar1,rmbar2))),kmadd(n1,kmul(n3,kmul(R4p1213,kmul(rmbar1,rmbar2))),kmadd(n1,kmul(n3,kmul(R4p1313,kmul(rmbar1,rmbar3))),kmadd(n1,kmul(n3,kmul(R4p1323,kmul(rmbar2,rmbar3))),kmadd(n2,kmul(n3,kmul(R4p2323,kmul(rmbar2,rmbar3))),kmadd(n1,kmul(n2,kmul(R4p1323,kmul(imbar3,imbar3))),kmadd(n1,kmul(nn,kmul(Ro313,kmul(imbar3,imbar3))),kmadd(n3,kmul(nn,kmul(Ro333,kmul(imbar3,imbar3))),kmul(R4p1223,kmul(rmbar1,kmul(rmbar3,kmul(n2,n2)))))))))))),ToReal(-2.),kmadd(n1,kmul(n3,kmul(R4p1223,kmul(kmul(rmbar2,rmbar2),ToReal(-2.)))),kmadd(R4p1313,kmsub(kmul(n3,n3),kmul(rmbar1,rmbar1),kmul(kmul(imbar3,imbar3),kmul(n1,n1))),kmadd(R4p2323,kmsub(kmul(n3,n3),kmul(rmbar2,rmbar2),kmul(kmul(imbar3,imbar3),kmul(n2,n2))),kmadd(Rojo33,kmsub(kmul(nn,nn),kmul(rmbar3,rmbar3),kmul(kmul(imbar3,imbar3),kmul(nn,nn))),kmadd(kmadd(n1,kmul(nn,kmul(rmbar1,kmul(rmbar2,Ro112))),kmadd(n1,kmul(nn,kmul(rmbar1,kmul(rmbar3,Ro113))),kmadd(n2,kmul(nn,kmul(rmbar1,kmul(rmbar3,Ro123))),kmadd(n3,kmul(nn,kmul(rmbar1,kmul(rmbar2,Ro132))),kmadd(n3,kmul(nn,kmul(rmbar1,kmul(rmbar3,Ro133))),kmadd(n1,kmul(nn,kmul(rmbar1,kmul(rmbar2,Ro211))),kmadd(n1,kmul(nn,kmul(rmbar2,kmul(rmbar3,Ro213))),kmadd(n2,kmul(nn,kmul(rmbar2,kmul(rmbar3,Ro223))),kmadd(n3,kmul(nn,kmul(rmbar1,kmul(rmbar2,Ro231))),kmadd(n3,kmul(nn,kmul(rmbar2,kmul(rmbar3,Ro233))),kmadd(n1,kmul(nn,kmul(rmbar1,kmul(rmbar3,Ro311))),kmadd(n1,kmul(nn,kmul(rmbar2,kmul(rmbar3,Ro312))),kmadd(n2,kmul(nn,kmul(rmbar2,kmul(rmbar3,Ro322))),kmadd(n3,kmul(nn,kmul(rmbar1,kmul(rmbar3,Ro331))),kmadd(n3,kmul(nn,kmul(rmbar2,kmul(rmbar3,Ro332))),kmadd(R4p1213,kmul(rmbar2,kmul(rmbar3,kmul(n1,n1))),kmul(R4p1323,kmul(rmbar1,kmul(rmbar2,kmul(n3,n3)))))))))))))))))))),ToReal(2.),kmadd(n1,kmul(nn,kmul(Ro111,kmul(kmul(rmbar1,rmbar1),ToReal(2.)))),kmadd(n2,kmul(nn,kmul(Ro121,kmul(kmul(rmbar1,rmbar1),ToReal(2.)))),kmadd(n3,kmul(nn,kmul(Ro131,kmul(kmul(rmbar1,rmbar1),ToReal(2.)))),kmadd(n1,kmul(nn,kmul(Ro212,kmul(kmul(rmbar2,rmbar2),ToReal(2.)))),kmadd(n2,kmul(nn,kmul(Ro222,kmul(kmul(rmbar2,rmbar2),ToReal(2.)))),kmadd(n3,kmul(nn,kmul(Ro232,kmul(kmul(rmbar2,rmbar2),ToReal(2.)))),kmadd(n1,kmul(n2,kmul(R4p1323,kmul(kmul(rmbar3,rmbar3),ToReal(2.)))),kmadd(n1,kmul(nn,kmul(Ro313,kmul(kmul(rmbar3,rmbar3),ToReal(2.)))),kmadd(n2,kmul(nn,kmul(Ro323,kmul(kmul(rmbar3,rmbar3),ToReal(2.)))),kmadd(n3,kmul(nn,kmul(Ro333,kmul(kmul(rmbar3,rmbar3),ToReal(2.)))),knmsub(kmul(imbar2,imbar2),kmadd(R4p1212,kmul(n1,n1),kmadd(R4p2323,kmul(n3,n3),kmadd(Rojo22,kmul(nn,nn),kmadd(n1,kmul(n3,kmul(R4p1223,ToReal(-2.))),kmadd(n1,kmul(nn,kmul(Ro212,ToReal(2.))),kmadd(n2,kmul(nn,kmul(Ro222,ToReal(2.))),kmul(n3,kmul(nn,kmul(Ro232,ToReal(2.)))))))))),kmadd(n2,kmadd(n1,kmul(R4p1213,kmul(rmbar1,kmul(rmbar3,ToReal(-2.)))),kmadd(n3,kmul(R4p1323,kmul(rmbar1,kmul(rmbar3,ToReal(-2.)))),kmadd(nn,kmul(Ro323,kmul(kmul(imbar3,imbar3),ToReal(-2.))),kmadd(n3,kmul(R4p1223,kmul(rmbar1,kmul(rmbar2,ToReal(2.)))),kmadd(n1,kmul(R4p1223,kmul(rmbar2,kmul(rmbar3,ToReal(2.)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro122,ToReal(2.)))),kmadd(nn,kmul(rmbar1,kmul(rmbar2,kmul(Ro221,ToReal(2.)))),kmadd(nn,kmul(rmbar1,kmul(rmbar3,kmul(Ro321,ToReal(2.)))),kmul(n3,kmul(R4p1213,kmul(kmul(rmbar1,rmbar1),ToReal(2.)))))))))))),knmsub(kmul(imbar1,imbar1),kmadd(R4p1212,kmul(n2,n2),kmadd(R4p1313,kmul(n3,n3),kmadd(n2,kmul(kmadd(n3,R4p1213,kmul(nn,Ro121)),ToReal(2.)),kmadd(n3,kmul(nn,kmul(Ro131,ToReal(2.))),kmul(nn,kmadd(nn,Rojo11,kmul(n1,kmul(Ro111,ToReal(2.))))))))),kmsub(imbar1,kmul(imbar3,kmadd(R4p1223,kmul(kmul(n2,n2),ToReal(2.)),kmadd(n1,kmul(kmadd(n2,R4p1213,kmsub(n3,R4p1313,kmul(nn,kadd(Ro113,Ro311)))),ToReal(2.)),kmsub(n2,kmul(kmsub(n3,R4p1323,kmul(nn,kadd(Ro123,Ro321))),ToReal(2.)),kmul(nn,kmadd(nn,kadd(Rojo13,Rojo31),kmul(n3,kmul(kadd(Ro133,Ro331),ToReal(2.))))))))),kmul(imbar2,kmadd(imbar3,kmadd(R4p1213,kmul(kmul(n1,n1),ToReal(2.)),kmadd(n1,kmul(kmadd(n2,R4p1223,kmsub(nn,kadd(Ro213,Ro312),kmul(n3,R4p1323))),ToReal(2.)),kmadd(n2,kmadd(n3,kmul(R4p2323,ToReal(-2.)),kmul(nn,kmul(kadd(Ro223,Ro322),ToReal(2.)))),kmul(nn,kmadd(nn,kadd(Rojo23,Rojo32),kmul(n3,kmul(kadd(Ro233,Ro332),ToReal(2.)))))))),kmul(imbar1,kmadd(Rojo12,kmul(nn,nn),kmadd(Rojo21,kmul(nn,nn),kmadd(n1,kmul(ToReal(-2.),kmadd(n2,R4p1212,kmsub(n3,R4p1213,kmul(nn,kadd(Ro112,Ro211))))),kmadd(n3,kmul(nn,kmul(Ro132,ToReal(2.))),kmadd(n2,kmul(kmadd(n3,R4p1223,kmul(nn,kadd(Ro122,Ro221))),ToReal(2.)),kmadd(n3,kmul(nn,kmul(Ro231,ToReal(2.))),kmul(R4p1323,kmul(kmul(n3,n3),ToReal(2.)))))))))))))))))))))))))))))))))))))))))))); CCTK_REAL_VEC Psi4iL = - knmadd(im1,kmadd(R4p3132,kmul(rm2,kmul(n3,n3)),kmadd(R4p3231,kmul(rm2,kmul(n3,n3)),kmadd(R4p3133,kmul(rm3,kmul(n3,n3)),kmadd(R4p3331,kmul(rm3,kmul(n3,n3)),kmadd(rm2,kmul(Rojo12,kmul(nn,nn)),kmadd(rm3,kmul(Rojo13,kmul(nn,nn)),kmadd(rm2,kmul(Rojo21,kmul(nn,nn)),kmadd(rm3,kmul(Rojo31,kmul(nn,nn)),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro132,ToReal(2.)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro133,ToReal(2.)))),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro231,ToReal(2.)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro331,ToReal(2.)))),kmadd(R4p3131,kmul(rm1,kmul(kmul(n3,n3),ToReal(2.))),kmadd(rm1,kmul(Rojo11,kmul(kmul(nn,nn),ToReal(2.))),kmadd(kmul(n1,n1),kmadd(R4p1112,rm2,kmadd(R4p1211,rm2,kmadd(R4p1113,rm3,kmadd(R4p1311,rm3,kmul(R4p1111,kmul(rm1,ToReal(2.))))))),kmadd(kmul(n2,n2),kmadd(R4p2122,rm2,kmadd(R4p2221,rm2,kmadd(R4p2123,rm3,kmadd(R4p2321,rm3,kmul(R4p2121,kmul(rm1,ToReal(2.))))))),kmadd(n1,kmadd(n2,kmadd(R4p1122,rm2,kmadd(R4p1221,rm2,kmadd(R4p2112,rm2,kmadd(R4p2211,rm2,kmadd(R4p1123,rm3,kmadd(R4p1321,rm3,kmadd(R4p2113,rm3,kmadd(R4p2311,rm3,kmadd(R4p1121,kmul(rm1,ToReal(2.)),kmul(R4p2111,kmul(rm1,ToReal(2.)))))))))))),kmadd(n3,kmadd(R4p1132,rm2,kmadd(R4p1231,rm2,kmadd(R4p3112,rm2,kmadd(R4p3211,rm2,kmadd(R4p1133,rm3,kmadd(R4p1331,rm3,kmadd(R4p3113,rm3,kmadd(R4p3311,rm3,kmadd(R4p1131,kmul(rm1,ToReal(2.)),kmul(R4p3111,kmul(rm1,ToReal(2.)))))))))))),kmul(nn,kmul(ToReal(2.),kmadd(rm2,kadd(Ro112,Ro211),kmadd(rm3,kadd(Ro113,Ro311),kmul(rm1,kmul(Ro111,ToReal(2.))))))))),kmadd(n2,kmadd(n3,kmadd(R4p2132,rm2,kmadd(R4p2231,rm2,kmadd(R4p3122,rm2,kmadd(R4p3221,rm2,kmadd(R4p2133,rm3,kmadd(R4p2331,rm3,kmadd(R4p3123,rm3,kmadd(R4p3321,rm3,kmadd(R4p2131,kmul(rm1,ToReal(2.)),kmul(R4p3121,kmul(rm1,ToReal(2.)))))))))))),kmul(nn,kmul(ToReal(2.),kmadd(rm2,kadd(Ro122,Ro221),kmadd(rm3,kadd(Ro123,Ro321),kmul(rm1,kmul(Ro121,ToReal(2.)))))))),kmul(n3,kmul(nn,kmul(rm1,kmul(Ro131,ToReal(4.))))))))))))))))))))))),kmadd(im3,kmadd(R4p3133,kmul(rm1,kmul(n3,n3)),kmadd(R4p3331,kmul(rm1,kmul(n3,n3)),kmadd(R4p3233,kmul(rm2,kmul(n3,n3)),kmadd(R4p3332,kmul(rm2,kmul(n3,n3)),kmadd(rm1,kmul(Rojo13,kmul(nn,nn)),kmadd(rm2,kmul(Rojo23,kmul(nn,nn)),kmadd(rm1,kmul(Rojo31,kmul(nn,nn)),kmadd(rm2,kmul(Rojo32,kmul(nn,nn)),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro133,ToReal(2.)))),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro233,ToReal(2.)))),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro331,ToReal(2.)))),kmadd(n3,kmul(nn,kmul(rm2,kmul(Ro332,ToReal(2.)))),kmadd(R4p3333,kmul(rm3,kmul(kmul(n3,n3),ToReal(2.))),kmadd(rm3,kmul(Rojo33,kmul(kmul(nn,nn),ToReal(2.))),kmadd(kmul(n1,n1),kmadd(R4p1113,rm1,kmadd(R4p1311,rm1,kmadd(R4p1213,rm2,kmadd(R4p1312,rm2,kmul(R4p1313,kmul(rm3,ToReal(2.))))))),kmadd(kmul(n2,n2),kmadd(R4p2123,rm1,kmadd(R4p2321,rm1,kmadd(R4p2223,rm2,kmadd(R4p2322,rm2,kmul(R4p2323,kmul(rm3,ToReal(2.))))))),kmadd(n1,kmadd(n2,kmadd(R4p1123,rm1,kmadd(R4p1321,rm1,kmadd(R4p2113,rm1,kmadd(R4p2311,rm1,kmadd(R4p1223,rm2,kmadd(R4p1322,rm2,kmadd(R4p2213,rm2,kmadd(R4p2312,rm2,kmadd(R4p1323,kmul(rm3,ToReal(2.)),kmul(R4p2313,kmul(rm3,ToReal(2.)))))))))))),kmadd(n3,kmadd(R4p1133,rm1,kmadd(R4p1331,rm1,kmadd(R4p3113,rm1,kmadd(R4p3311,rm1,kmadd(R4p1233,rm2,kmadd(R4p1332,rm2,kmadd(R4p3213,rm2,kmadd(R4p3312,rm2,kmadd(R4p1333,kmul(rm3,ToReal(2.)),kmul(R4p3313,kmul(rm3,ToReal(2.)))))))))))),kmul(nn,kmul(ToReal(2.),kmadd(rm1,kadd(Ro113,Ro311),kmadd(rm2,kadd(Ro213,Ro312),kmul(rm3,kmul(Ro313,ToReal(2.))))))))),kmadd(n2,kmadd(n3,kmadd(R4p2133,rm1,kmadd(R4p2331,rm1,kmadd(R4p3123,rm1,kmadd(R4p3321,rm1,kmadd(R4p2233,rm2,kmadd(R4p2332,rm2,kmadd(R4p3223,rm2,kmadd(R4p3322,rm2,kmadd(R4p2333,kmul(rm3,ToReal(2.)),kmul(R4p3323,kmul(rm3,ToReal(2.)))))))))))),kmul(nn,kmul(ToReal(2.),kmadd(rm1,kadd(Ro123,Ro321),kmadd(rm2,kadd(Ro223,Ro322),kmul(rm3,kmul(Ro323,ToReal(2.)))))))),kmul(n3,kmul(nn,kmul(rm3,kmul(Ro333,ToReal(4.))))))))))))))))))))))),kmul(im2,kmadd(R4p3132,kmul(rm1,kmul(n3,n3)),kmadd(R4p3231,kmul(rm1,kmul(n3,n3)),kmadd(R4p3233,kmul(rm3,kmul(n3,n3)),kmadd(R4p3332,kmul(rm3,kmul(n3,n3)),kmadd(rm1,kmul(Rojo12,kmul(nn,nn)),kmadd(rm1,kmul(Rojo21,kmul(nn,nn)),kmadd(rm3,kmul(Rojo23,kmul(nn,nn)),kmadd(rm3,kmul(Rojo32,kmul(nn,nn)),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro132,ToReal(2.)))),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro231,ToReal(2.)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro233,ToReal(2.)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro332,ToReal(2.)))),kmadd(R4p3232,kmul(rm2,kmul(kmul(n3,n3),ToReal(2.))),kmadd(rm2,kmul(Rojo22,kmul(kmul(nn,nn),ToReal(2.))),kmadd(kmul(n1,n1),kmadd(R4p1112,rm1,kmadd(R4p1211,rm1,kmadd(R4p1213,rm3,kmadd(R4p1312,rm3,kmul(R4p1212,kmul(rm2,ToReal(2.))))))),kmadd(kmul(n2,n2),kmadd(R4p2122,rm1,kmadd(R4p2221,rm1,kmadd(R4p2223,rm3,kmadd(R4p2322,rm3,kmul(R4p2222,kmul(rm2,ToReal(2.))))))),kmadd(n1,kmadd(n2,kmadd(R4p1122,rm1,kmadd(R4p1221,rm1,kmadd(R4p2112,rm1,kmadd(R4p2211,rm1,kmadd(R4p1223,rm3,kmadd(R4p1322,rm3,kmadd(R4p2213,rm3,kmadd(R4p2312,rm3,kmadd(R4p1222,kmul(rm2,ToReal(2.)),kmul(R4p2212,kmul(rm2,ToReal(2.)))))))))))),kmadd(n3,kmadd(R4p1132,rm1,kmadd(R4p1231,rm1,kmadd(R4p3112,rm1,kmadd(R4p3211,rm1,kmadd(R4p1233,rm3,kmadd(R4p1332,rm3,kmadd(R4p3213,rm3,kmadd(R4p3312,rm3,kmadd(R4p1232,kmul(rm2,ToReal(2.)),kmul(R4p3212,kmul(rm2,ToReal(2.)))))))))))),kmul(nn,kmul(ToReal(2.),kmadd(rm1,kadd(Ro112,Ro211),kmadd(rm3,kadd(Ro213,Ro312),kmul(rm2,kmul(Ro212,ToReal(2.))))))))),kmadd(n2,kmadd(n3,kmadd(R4p2132,rm1,kmadd(R4p2231,rm1,kmadd(R4p3122,rm1,kmadd(R4p3221,rm1,kmadd(R4p2233,rm3,kmadd(R4p2332,rm3,kmadd(R4p3223,rm3,kmadd(R4p3322,rm3,kmadd(R4p2232,kmul(rm2,ToReal(2.)),kmul(R4p3222,kmul(rm2,ToReal(2.)))))))))))),kmul(nn,kmul(ToReal(2.),kmadd(rm1,kadd(Ro122,Ro221),kmadd(rm3,kadd(Ro223,Ro322),kmul(rm2,kmul(Ro222,ToReal(2.)))))))),kmul(n3,kmul(nn,kmul(rm2,kmul(Ro232,ToReal(4.)))))))))))))))))))))))))); + knmadd(im1,kmadd(kmadd(R4p1313,rm1,kmul(R4p1323,rm2)),kmul(kmul(n3,n3),ToReal(2.)),kmadd(kmul(n2,n2),kmul(kmsub(R4p1212,rm1,kmul(R4p1223,rm3)),ToReal(2.)),kmadd(n2,kmul(ToReal(2.),knmsub(n1,kmadd(R4p1212,rm2,kmul(R4p1213,rm3)),kmadd(n3,kmadd(R4p1223,rm2,kmsub(R4p1213,kmul(rm1,ToReal(2.)),kmul(R4p1323,rm3))),kmul(nn,kmadd(rm2,kadd(Ro122,Ro221),kmadd(rm3,kadd(Ro123,Ro321),kmul(rm1,kmul(Ro121,ToReal(2.))))))))),kmadd(n3,kmul(ToReal(2.),kmsub(nn,kmadd(rm2,kadd(Ro132,Ro231),kmadd(rm3,kadd(Ro133,Ro331),kmul(rm1,kmul(Ro131,ToReal(2.))))),kmul(n1,kmadd(R4p1213,rm2,kmul(R4p1313,rm3))))),kmul(nn,kmadd(n1,kmul(ToReal(2.),kmadd(rm2,kadd(Ro112,Ro211),kmadd(rm3,kadd(Ro113,Ro311),kmul(rm1,kmul(Ro111,ToReal(2.)))))),kmul(nn,kmadd(rm2,kadd(Rojo12,Rojo21),kmadd(rm3,kadd(Rojo13,Rojo31),kmul(rm1,kmul(Rojo11,ToReal(2.)))))))))))),kmadd(im2,kmadd(rm1,kmul(Rojo12,kmul(nn,nn)),kmadd(rm1,kmul(Rojo21,kmul(nn,nn)),kmadd(rm3,kmul(Rojo23,kmul(nn,nn)),kmadd(rm3,kmul(Rojo32,kmul(nn,nn)),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro132,ToReal(2.)))),kmadd(n3,kmul(nn,kmul(rm1,kmul(Ro231,ToReal(2.)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro233,ToReal(2.)))),kmadd(n3,kmul(nn,kmul(rm3,kmul(Ro332,ToReal(2.)))),kmadd(kmadd(R4p1212,rm2,kmul(R4p1213,rm3)),kmul(kmul(n1,n1),ToReal(2.)),kmadd(R4p1323,kmul(rm1,kmul(kmul(n3,n3),ToReal(2.))),kmadd(R4p2323,kmul(rm2,kmul(kmul(n3,n3),ToReal(2.))),kmadd(rm2,kmul(Rojo22,kmul(kmul(nn,nn),ToReal(2.))),kmadd(n1,kmul(ToReal(2.),kmadd(n2,kmsub(R4p1223,rm3,kmul(R4p1212,rm1)),kmsub(nn,kmadd(rm1,kadd(Ro112,Ro211),kmadd(rm3,kadd(Ro213,Ro312),kmul(rm2,kmul(Ro212,ToReal(2.))))),kmul(n3,kmadd(R4p1213,rm1,kmadd(R4p1323,rm3,kmul(R4p1223,kmul(rm2,ToReal(2.))))))))),kmadd(n2,kmul(ToReal(2.),kmadd(n3,kmsub(R4p1223,rm1,kmul(R4p2323,rm3)),kmul(nn,kmadd(rm1,kadd(Ro122,Ro221),kmadd(rm3,kadd(Ro223,Ro322),kmul(rm2,kmul(Ro222,ToReal(2.)))))))),kmul(n3,kmul(nn,kmul(rm2,kmul(Ro232,ToReal(4.))))))))))))))))))),kmul(im3,kmadd(kmadd(R4p1213,rm2,kmul(R4p1313,rm3)),kmul(kmul(n1,n1),ToReal(2.)),kmadd(kmul(n2,n2),kmadd(R4p1223,kmul(rm1,ToReal(-2.)),kmul(R4p2323,kmul(rm3,ToReal(2.)))),kmadd(n1,kmul(ToReal(2.),knmsub(n3,kmadd(R4p1313,rm1,kmul(R4p1323,rm2)),kmadd(n2,kmadd(R4p1223,rm2,kmsub(R4p1323,kmul(rm3,ToReal(2.)),kmul(R4p1213,rm1))),kmul(nn,kmadd(rm1,kadd(Ro113,Ro311),kmadd(rm2,kadd(Ro213,Ro312),kmul(rm3,kmul(Ro313,ToReal(2.))))))))),kmadd(n2,kmul(ToReal(2.),kmsub(nn,kmadd(rm1,kadd(Ro123,Ro321),kmadd(rm2,kadd(Ro223,Ro322),kmul(rm3,kmul(Ro323,ToReal(2.))))),kmul(n3,kmadd(R4p1323,rm1,kmul(R4p2323,rm2))))),kmul(nn,kmadd(n3,kmul(ToReal(2.),kmadd(rm1,kadd(Ro133,Ro331),kmadd(rm2,kadd(Ro233,Ro332),kmul(rm3,kmul(Ro333,ToReal(2.)))))),kmul(nn,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); |