diff options
Diffstat (limited to 'src/psis_calc_4th.cc')
| -rw-r--r-- | src/psis_calc_4th.cc | 271 |
1 files changed, 141 insertions, 130 deletions
diff --git a/src/psis_calc_4th.cc b/src/psis_calc_4th.cc index 6afc465..9aa697e 100644 --- a/src/psis_calc_4th.cc +++ b/src/psis_calc_4th.cc @@ -65,6 +65,9 @@ static void psis_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, ptrdiff_t const di = 1; ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0); ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0); + ptrdiff_t const cdi = sizeof(CCTK_REAL) * di; + ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj; + ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk; CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0)); CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1)); CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2)); @@ -1122,136 +1125,144 @@ static void psis_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, + gInv31*R3331 + gInv32*R3332 + gInv33*R3333 + gInv11*(kxxL*kzzL + R3131 - SQR(kxzL)) + gInv22*(kyyL*kzzL + R3232 - SQR(kyzL)); - CCTK_REAL Psi4rL = 2*nn*(-(n1*((imbar2*imbar3 - rmbar2*rmbar3)*Ro213 + - imbar1*(imbar2*(Ro112 + Ro211) + imbar3*(Ro113 + Ro311)) - - rmbar1*(rmbar2*(Ro112 + Ro211) + rmbar3*(Ro113 + Ro311)) + - (imbar2*imbar3 - rmbar2*rmbar3)*Ro312 + Ro111*(SQR(imbar1) - - SQR(rmbar1)) + Ro212*(SQR(imbar2) - SQR(rmbar2)) + Ro313*(SQR(imbar3) - - SQR(rmbar3)))) - n2*((imbar2*imbar3 - rmbar2*rmbar3)*Ro223 + - imbar1*(imbar2*(Ro122 + Ro221) + imbar3*(Ro123 + Ro321)) - - rmbar1*(rmbar2*(Ro122 + Ro221) + rmbar3*(Ro123 + Ro321)) + - (imbar2*imbar3 - rmbar2*rmbar3)*Ro322 + Ro121*(SQR(imbar1) - - SQR(rmbar1)) + Ro222*(SQR(imbar2) - SQR(rmbar2)) + Ro323*(SQR(imbar3) - - SQR(rmbar3))) - n3*((imbar2*imbar3 - rmbar2*rmbar3)*Ro233 + - imbar1*(imbar2*(Ro132 + Ro231) + imbar3*(Ro133 + Ro331)) - - rmbar1*(rmbar2*(Ro132 + Ro231) + rmbar3*(Ro133 + Ro331)) + - (imbar2*imbar3 - rmbar2*rmbar3)*Ro332 + Ro131*(SQR(imbar1) - - SQR(rmbar1)) + Ro232*(SQR(imbar2) - SQR(rmbar2)) + Ro333*(SQR(imbar3) - - SQR(rmbar3)))) - SQR(nn)*((imbar2*imbar3 - rmbar2*rmbar3)*Rojo23 + - imbar1*(imbar2*(Rojo12 + Rojo21) + imbar3*(Rojo13 + Rojo31)) - - rmbar1*(rmbar2*(Rojo12 + Rojo21) + rmbar3*(Rojo13 + Rojo31)) + - (imbar2*imbar3 - rmbar2*rmbar3)*Rojo32 + Rojo11*(SQR(imbar1) - - SQR(rmbar1)) + Rojo22*(SQR(imbar2) - SQR(rmbar2)) + Rojo33*(SQR(imbar3) - - SQR(rmbar3))) + n1*(-(imbar1*(imbar2*(n1*(R4p1112 + R4p1211) + - n2*(R4p1122 + R4p1221) + n3*(R4p1132 + R4p1231)) + imbar3*(n1*(R4p1113 - + R4p1311) + n2*(R4p1123 + R4p1321) + n3*(R4p1133 + R4p1331)))) - - imbar2*imbar3*(n1*(R4p1213 + R4p1312) + n2*(R4p1223 + R4p1322) + - n3*(R4p1233 + R4p1332)) + n1*R4p1211*rmbar1*rmbar2 + - n2*R4p1221*rmbar1*rmbar2 + n3*R4p1231*rmbar1*rmbar2 + - n1*R4p1113*rmbar1*rmbar3 + n2*R4p1123*rmbar1*rmbar3 + - n3*R4p1133*rmbar1*rmbar3 + n1*R4p1311*rmbar1*rmbar3 + - n2*R4p1321*rmbar1*rmbar3 + n3*R4p1331*rmbar1*rmbar3 + - n1*R4p1213*rmbar2*rmbar3 + n2*R4p1223*rmbar2*rmbar3 + - n3*R4p1233*rmbar2*rmbar3 + n1*R4p1312*rmbar2*rmbar3 + - n2*R4p1322*rmbar2*rmbar3 + n3*R4p1332*rmbar2*rmbar3 - (n1*R4p1111 + - n2*R4p1121 + n3*R4p1131)*SQR(imbar1) - (n1*R4p1212 + n2*R4p1222 + - n3*R4p1232)*SQR(imbar2) + n1*(R4p1112*rmbar1*rmbar2 - - R4p1313*SQR(imbar3)) + n2*(R4p1122*rmbar1*rmbar2 - R4p1323*SQR(imbar3)) - + n3*(R4p1132*rmbar1*rmbar2 - R4p1333*SQR(imbar3)) + - n1*R4p1111*SQR(rmbar1) + n2*R4p1121*SQR(rmbar1) + - n3*R4p1131*SQR(rmbar1) + n1*R4p1212*SQR(rmbar2) + - n2*R4p1222*SQR(rmbar2) + n3*R4p1232*SQR(rmbar2) + - n1*R4p1313*SQR(rmbar3) + n2*R4p1323*SQR(rmbar3) + - n3*R4p1333*SQR(rmbar3)) + n2*(-(imbar1*(imbar2*(n1*(R4p2112 + R4p2211) - + n2*(R4p2122 + R4p2221) + n3*(R4p2132 + R4p2231)) + - imbar3*(n1*(R4p2113 + R4p2311) + n2*(R4p2123 + R4p2321) + n3*(R4p2133 + - R4p2331)))) - imbar2*imbar3*(n1*(R4p2213 + R4p2312) + n2*(R4p2223 + - R4p2322) + n3*(R4p2233 + R4p2332)) + n1*R4p2211*rmbar1*rmbar2 + - n2*R4p2221*rmbar1*rmbar2 + n3*R4p2231*rmbar1*rmbar2 + - n1*R4p2113*rmbar1*rmbar3 + n2*R4p2123*rmbar1*rmbar3 + - n3*R4p2133*rmbar1*rmbar3 + n1*R4p2311*rmbar1*rmbar3 + - n2*R4p2321*rmbar1*rmbar3 + n3*R4p2331*rmbar1*rmbar3 + - n1*R4p2213*rmbar2*rmbar3 + n2*R4p2223*rmbar2*rmbar3 + - n3*R4p2233*rmbar2*rmbar3 + n1*R4p2312*rmbar2*rmbar3 + - n2*R4p2322*rmbar2*rmbar3 + n3*R4p2332*rmbar2*rmbar3 - (n1*R4p2111 + - n2*R4p2121 + n3*R4p2131)*SQR(imbar1) - (n1*R4p2212 + n2*R4p2222 + - n3*R4p2232)*SQR(imbar2) + n1*(R4p2112*rmbar1*rmbar2 - - R4p2313*SQR(imbar3)) + n2*(R4p2122*rmbar1*rmbar2 - R4p2323*SQR(imbar3)) - + n3*(R4p2132*rmbar1*rmbar2 - R4p2333*SQR(imbar3)) + - n1*R4p2111*SQR(rmbar1) + n2*R4p2121*SQR(rmbar1) + - n3*R4p2131*SQR(rmbar1) + n1*R4p2212*SQR(rmbar2) + - n2*R4p2222*SQR(rmbar2) + n3*R4p2232*SQR(rmbar2) + - n1*R4p2313*SQR(rmbar3) + n2*R4p2323*SQR(rmbar3) + - n3*R4p2333*SQR(rmbar3)) + n3*(-(imbar1*(imbar2*(n1*(R4p3112 + R4p3211) - + n2*(R4p3122 + R4p3221) + n3*(R4p3132 + R4p3231)) + - imbar3*(n1*(R4p3113 + R4p3311) + n2*(R4p3123 + R4p3321) + n3*(R4p3133 + - R4p3331)))) - imbar2*imbar3*(n1*(R4p3213 + R4p3312) + n2*(R4p3223 + - R4p3322) + n3*(R4p3233 + R4p3332)) + n1*R4p3211*rmbar1*rmbar2 + - n2*R4p3221*rmbar1*rmbar2 + n3*R4p3231*rmbar1*rmbar2 + - n1*R4p3113*rmbar1*rmbar3 + n2*R4p3123*rmbar1*rmbar3 + - n3*R4p3133*rmbar1*rmbar3 + n1*R4p3311*rmbar1*rmbar3 + - n2*R4p3321*rmbar1*rmbar3 + n3*R4p3331*rmbar1*rmbar3 + - n1*R4p3213*rmbar2*rmbar3 + n2*R4p3223*rmbar2*rmbar3 + - n3*R4p3233*rmbar2*rmbar3 + n1*R4p3312*rmbar2*rmbar3 + - n2*R4p3322*rmbar2*rmbar3 + n3*R4p3332*rmbar2*rmbar3 - (n1*R4p3111 + - n2*R4p3121 + n3*R4p3131)*SQR(imbar1) - (n1*R4p3212 + n2*R4p3222 + - n3*R4p3232)*SQR(imbar2) + n1*(R4p3112*rmbar1*rmbar2 - - R4p3313*SQR(imbar3)) + n2*(R4p3122*rmbar1*rmbar2 - R4p3323*SQR(imbar3)) - + n3*(R4p3132*rmbar1*rmbar2 - R4p3333*SQR(imbar3)) + - n1*R4p3111*SQR(rmbar1) + n2*R4p3121*SQR(rmbar1) + - n3*R4p3131*SQR(rmbar1) + n1*R4p3212*SQR(rmbar2) + - n2*R4p3222*SQR(rmbar2) + n3*R4p3232*SQR(rmbar2) + - n1*R4p3313*SQR(rmbar3) + n2*R4p3323*SQR(rmbar3) + - n3*R4p3333*SQR(rmbar3)); - - CCTK_REAL Psi4iL = -(n1*(im1*(n1*(2*R4p1111*rm1 + R4p1112*rm2 + - R4p1211*rm2 + R4p1113*rm3 + R4p1311*rm3) + n2*(2*R4p1121*rm1 + - R4p1122*rm2 + R4p1221*rm2 + R4p1123*rm3 + R4p1321*rm3) + - n3*(2*R4p1131*rm1 + R4p1132*rm2 + R4p1231*rm2 + R4p1133*rm3 + - R4p1331*rm3)) + im2*(n1*(R4p1112*rm1 + R4p1211*rm1 + 2*R4p1212*rm2 + - R4p1213*rm3 + R4p1312*rm3) + n2*(R4p1122*rm1 + R4p1221*rm1 + - 2*R4p1222*rm2 + R4p1223*rm3 + R4p1322*rm3) + n3*(R4p1132*rm1 + - R4p1231*rm1 + 2*R4p1232*rm2 + R4p1233*rm3 + R4p1332*rm3)) + - im3*(n1*(R4p1113*rm1 + R4p1311*rm1 + R4p1213*rm2 + R4p1312*rm2 + - 2*R4p1313*rm3) + n2*(R4p1123*rm1 + R4p1321*rm1 + R4p1223*rm2 + - R4p1322*rm2 + 2*R4p1323*rm3) + n3*(R4p1133*rm1 + R4p1331*rm1 + - R4p1233*rm2 + R4p1332*rm2 + 2*R4p1333*rm3)))) - - n2*(im1*(n1*(2*R4p2111*rm1 + R4p2112*rm2 + R4p2211*rm2 + R4p2113*rm3 + - R4p2311*rm3) + n2*(2*R4p2121*rm1 + R4p2122*rm2 + R4p2221*rm2 + - R4p2123*rm3 + R4p2321*rm3) + n3*(2*R4p2131*rm1 + R4p2132*rm2 + - R4p2231*rm2 + R4p2133*rm3 + R4p2331*rm3)) + im2*(n1*(R4p2112*rm1 + - R4p2211*rm1 + 2*R4p2212*rm2 + R4p2213*rm3 + R4p2312*rm3) + - n2*(R4p2122*rm1 + R4p2221*rm1 + 2*R4p2222*rm2 + R4p2223*rm3 + - R4p2322*rm3) + n3*(R4p2132*rm1 + R4p2231*rm1 + 2*R4p2232*rm2 + - R4p2233*rm3 + R4p2332*rm3)) + im3*(n1*(R4p2113*rm1 + R4p2311*rm1 + - R4p2213*rm2 + R4p2312*rm2 + 2*R4p2313*rm3) + n2*(R4p2123*rm1 + - R4p2321*rm1 + R4p2223*rm2 + R4p2322*rm2 + 2*R4p2323*rm3) + - n3*(R4p2133*rm1 + R4p2331*rm1 + R4p2233*rm2 + R4p2332*rm2 + - 2*R4p2333*rm3))) - n3*(im1*(n1*(2*R4p3111*rm1 + R4p3112*rm2 + - R4p3211*rm2 + R4p3113*rm3 + R4p3311*rm3) + n2*(2*R4p3121*rm1 + - R4p3122*rm2 + R4p3221*rm2 + R4p3123*rm3 + R4p3321*rm3) + - n3*(2*R4p3131*rm1 + R4p3132*rm2 + R4p3231*rm2 + R4p3133*rm3 + - R4p3331*rm3)) + im2*(n1*(R4p3112*rm1 + R4p3211*rm1 + 2*R4p3212*rm2 + - R4p3213*rm3 + R4p3312*rm3) + n2*(R4p3122*rm1 + R4p3221*rm1 + - 2*R4p3222*rm2 + R4p3223*rm3 + R4p3322*rm3) + n3*(R4p3132*rm1 + - R4p3231*rm1 + 2*R4p3232*rm2 + R4p3233*rm3 + R4p3332*rm3)) + - im3*(n1*(R4p3113*rm1 + R4p3311*rm1 + R4p3213*rm2 + R4p3312*rm2 + - 2*R4p3313*rm3) + n2*(R4p3123*rm1 + R4p3321*rm1 + R4p3223*rm2 + - R4p3322*rm2 + 2*R4p3323*rm3) + n3*(R4p3133*rm1 + R4p3331*rm1 + - R4p3233*rm2 + R4p3332*rm2 + 2*R4p3333*rm3))) - - 2*nn*(im1*(n1*(2*rm1*Ro111 + rm2*(Ro112 + Ro211) + rm3*(Ro113 + Ro311)) - + n2*(2*rm1*Ro121 + rm2*(Ro122 + Ro221) + rm3*(Ro123 + Ro321)) + - n3*(2*rm1*Ro131 + rm2*(Ro132 + Ro231) + rm3*(Ro133 + Ro331))) + - im2*(n1*(rm1*(Ro112 + Ro211) + 2*rm2*Ro212 + rm3*(Ro213 + Ro312)) + - n2*(rm1*(Ro122 + Ro221) + 2*rm2*Ro222 + rm3*(Ro223 + Ro322)) + - n3*(rm1*(Ro132 + Ro231) + 2*rm2*Ro232 + rm3*(Ro233 + Ro332))) + - im3*(n1*(rm1*(Ro113 + Ro311) + rm2*(Ro213 + Ro312) + 2*rm3*Ro313) + - n2*(rm1*(Ro123 + Ro321) + rm2*(Ro223 + Ro322) + 2*rm3*Ro323) + - n3*(rm1*(Ro133 + Ro331) + rm2*(Ro233 + Ro332) + 2*rm3*Ro333))) - - (im1*(2*rm1*Rojo11 + rm2*(Rojo12 + Rojo21) + rm3*(Rojo13 + Rojo31)) + - im2*(rm1*(Rojo12 + Rojo21) + 2*rm2*Rojo22 + rm3*(Rojo23 + Rojo32)) + - im3*(rm1*(Rojo13 + Rojo31) + rm2*(Rojo23 + Rojo32) + - 2*rm3*Rojo33))*SQR(nn); + CCTK_REAL Psi4rL = n3*(-(n2*R4p2333) - n1*R4p3313)*SQR(imbar3) + + n1*(n3*R4p1132*rmbar1*rmbar2 + n2*R4p1221*rmbar1*rmbar2 + + n3*R4p1231*rmbar1*rmbar2 + n2*R4p1123*rmbar1*rmbar3 + + n3*R4p1331*rmbar1*rmbar3 + n2*R4p2113*rmbar1*rmbar3 + + n3*R4p1233*rmbar2*rmbar3 + n2*R4p1322*rmbar2*rmbar3 + + n3*R4p1332*rmbar2*rmbar3 + n2*R4p2213*rmbar2*rmbar3 + + 2*nn*rmbar1*rmbar2*Ro112 + 2*nn*rmbar1*rmbar3*Ro113 + + 2*nn*rmbar1*rmbar2*Ro211 + 2*nn*rmbar2*rmbar3*Ro213 + + 2*nn*rmbar1*rmbar3*Ro311 + 2*nn*rmbar2*rmbar3*Ro312 + + n2*(R4p1122*rmbar1*rmbar2 - R4p1323*SQR(imbar3))) + + n3*(n2*R4p3122*rmbar1*rmbar2 + n1*R4p3211*rmbar1*rmbar2 + + n2*R4p3221*rmbar1*rmbar2 + n1*R4p1133*rmbar1*rmbar3 + + n2*R4p2331*rmbar1*rmbar3 + n1*R4p3113*rmbar1*rmbar3 + + n2*R4p3123*rmbar1*rmbar3 + n1*R4p3311*rmbar1*rmbar3 + + n2*R4p2332*rmbar2*rmbar3 + n1*R4p3213*rmbar2*rmbar3 + + n2*R4p3223*rmbar2*rmbar3 + n1*R4p3312*rmbar2*rmbar3 + + n2*R4p3322*rmbar2*rmbar3 + 2*nn*rmbar1*rmbar2*Ro132 + + 2*nn*rmbar1*rmbar3*Ro133 + 2*nn*rmbar1*rmbar2*Ro231 + + 2*nn*rmbar2*rmbar3*Ro233 + 2*nn*rmbar1*rmbar3*Ro331 + + 2*nn*rmbar2*rmbar3*Ro332 + n1*(R4p3112*rmbar1*rmbar2 - + R4p1333*SQR(imbar3))) + n2*(n3*R4p2132*rmbar1*rmbar2 + + n1*R4p2211*rmbar1*rmbar2 + n3*R4p2231*rmbar1*rmbar2 + + n1*R4p1321*rmbar1*rmbar3 + n3*R4p2133*rmbar1*rmbar3 + + n1*R4p2311*rmbar1*rmbar3 + n3*R4p3321*rmbar1*rmbar3 + + n1*R4p1223*rmbar2*rmbar3 + n3*R4p2233*rmbar2*rmbar3 + + n1*R4p2312*rmbar2*rmbar3 + 2*nn*rmbar1*rmbar2*Ro122 + + 2*nn*rmbar1*rmbar3*Ro123 + 2*nn*rmbar1*rmbar2*Ro221 + + 2*nn*rmbar2*rmbar3*Ro223 + 2*nn*rmbar1*rmbar3*Ro321 + + 2*nn*rmbar2*rmbar3*Ro322 + n1*(R4p2112*rmbar1*rmbar2 - + R4p2313*SQR(imbar3))) + R4p1112*rmbar1*rmbar2*SQR(n1) + + R4p1211*rmbar1*rmbar2*SQR(n1) + R4p1113*rmbar1*rmbar3*SQR(n1) + + R4p1311*rmbar1*rmbar3*SQR(n1) + R4p1213*rmbar2*rmbar3*SQR(n1) + + R4p1312*rmbar2*rmbar3*SQR(n1) + R4p2122*rmbar1*rmbar2*SQR(n2) + + R4p2221*rmbar1*rmbar2*SQR(n2) + R4p2123*rmbar1*rmbar3*SQR(n2) + + R4p2321*rmbar1*rmbar3*SQR(n2) + R4p2223*rmbar2*rmbar3*SQR(n2) + + R4p2322*rmbar2*rmbar3*SQR(n2) + SQR(imbar3)*(-2*n2*nn*Ro323 - + R4p2323*SQR(n2)) + R4p3132*rmbar1*rmbar2*SQR(n3) + + R4p3231*rmbar1*rmbar2*SQR(n3) + R4p3133*rmbar1*rmbar3*SQR(n3) + + R4p3331*rmbar1*rmbar3*SQR(n3) + R4p3233*rmbar2*rmbar3*SQR(n3) + + R4p3332*rmbar2*rmbar3*SQR(n3) + SQR(imbar3)*(-2*n3*nn*Ro333 - + R4p3333*SQR(n3)) + rmbar1*rmbar2*Rojo12*SQR(nn) + + rmbar1*rmbar3*Rojo13*SQR(nn) + rmbar1*rmbar2*Rojo21*SQR(nn) + + rmbar2*rmbar3*Rojo23*SQR(nn) + rmbar1*rmbar3*Rojo31*SQR(nn) + + rmbar2*rmbar3*Rojo32*SQR(nn) - SQR(imbar1)*(n2*n3*R4p2131 + + n2*n3*R4p3121 + n1*(n2*(R4p1121 + R4p2111) + n3*(R4p1131 + R4p3111) + + 2*nn*Ro111) + 2*n2*nn*Ro121 + 2*n3*nn*Ro131 + R4p1111*SQR(n1) + + R4p2121*SQR(n2) + R4p3131*SQR(n3) + Rojo11*SQR(nn)) - + SQR(imbar2)*(n2*n3*R4p2232 + n2*n3*R4p3222 + n1*(n2*(R4p1222 + R4p2212) + + n3*(R4p1232 + R4p3212) + 2*nn*Ro212) + 2*n2*nn*Ro222 + 2*n3*nn*Ro232 + + R4p1212*SQR(n1) + R4p2222*SQR(n2) + R4p3232*SQR(n3) + Rojo22*SQR(nn)) + - imbar2*imbar3*(2*n3*nn*Ro233 + n1*(n2*(R4p1223 + R4p1322 + R4p2213 + + R4p2312) + n3*(R4p1233 + R4p1332 + R4p3213 + R4p3312) + 2*nn*(Ro213 + + Ro312)) + n2*(n3*(R4p2233 + R4p2332 + R4p3223 + R4p3322) + 2*nn*(Ro223 + + Ro322)) + 2*n3*nn*Ro332 + (R4p1213 + R4p1312)*SQR(n1) + (R4p2223 + + R4p2322)*SQR(n2) + R4p3233*SQR(n3) + R4p3332*SQR(n3) + Rojo23*SQR(nn) + + Rojo32*SQR(nn)) + SQR(imbar3)*(-(n2*n3*R4p3323) - 2*n1*nn*Ro313 - + R4p1313*SQR(n1) - Rojo33*SQR(nn)) - imbar1*(imbar2*(2*n3*nn*Ro132 + + n1*(n2*(R4p1122 + R4p1221 + R4p2112 + R4p2211) + n3*(R4p1132 + R4p1231 + + R4p3112 + R4p3211) + 2*nn*(Ro112 + Ro211)) + n2*(n3*(R4p2132 + + R4p2231 + R4p3122 + R4p3221) + 2*nn*(Ro122 + Ro221)) + 2*n3*nn*Ro231 + + (R4p1112 + R4p1211)*SQR(n1) + (R4p2122 + R4p2221)*SQR(n2) + + R4p3132*SQR(n3) + R4p3231*SQR(n3) + Rojo12*SQR(nn) + Rojo21*SQR(nn)) + + imbar3*(2*n3*nn*Ro133 + n1*(n2*(R4p1123 + R4p1321 + R4p2113 + R4p2311) + + n3*(R4p1133 + R4p1331 + R4p3113 + R4p3311) + 2*nn*(Ro113 + Ro311)) + + n2*(n3*(R4p2133 + R4p2331 + R4p3123 + R4p3321) + 2*nn*(Ro123 + Ro321)) + + 2*n3*nn*Ro331 + (R4p1113 + R4p1311)*SQR(n1) + (R4p2123 + + R4p2321)*SQR(n2) + R4p3133*SQR(n3) + R4p3331*SQR(n3) + Rojo13*SQR(nn) + + Rojo31*SQR(nn))) + n1*n2*R4p1121*SQR(rmbar1) + + n1*n3*R4p1131*SQR(rmbar1) + n1*n2*R4p2111*SQR(rmbar1) + + n2*n3*R4p2131*SQR(rmbar1) + n1*n3*R4p3111*SQR(rmbar1) + + n2*n3*R4p3121*SQR(rmbar1) + 2*n1*nn*Ro111*SQR(rmbar1) + + 2*n2*nn*Ro121*SQR(rmbar1) + 2*n3*nn*Ro131*SQR(rmbar1) + + R4p1111*SQR(n1)*SQR(rmbar1) + R4p2121*SQR(n2)*SQR(rmbar1) + + R4p3131*SQR(n3)*SQR(rmbar1) + Rojo11*SQR(nn)*SQR(rmbar1) + + n1*n2*R4p1222*SQR(rmbar2) + n1*n3*R4p1232*SQR(rmbar2) + + n1*n2*R4p2212*SQR(rmbar2) + n2*n3*R4p2232*SQR(rmbar2) + + n1*n3*R4p3212*SQR(rmbar2) + n2*n3*R4p3222*SQR(rmbar2) + + 2*n1*nn*Ro212*SQR(rmbar2) + 2*n2*nn*Ro222*SQR(rmbar2) + + 2*n3*nn*Ro232*SQR(rmbar2) + R4p1212*SQR(n1)*SQR(rmbar2) + + R4p2222*SQR(n2)*SQR(rmbar2) + R4p3232*SQR(n3)*SQR(rmbar2) + + Rojo22*SQR(nn)*SQR(rmbar2) + n1*n2*R4p1323*SQR(rmbar3) + + n1*n3*R4p1333*SQR(rmbar3) + n1*n2*R4p2313*SQR(rmbar3) + + n2*n3*R4p2333*SQR(rmbar3) + n1*n3*R4p3313*SQR(rmbar3) + + n2*n3*R4p3323*SQR(rmbar3) + 2*n1*nn*Ro313*SQR(rmbar3) + + 2*n2*nn*Ro323*SQR(rmbar3) + 2*n3*nn*Ro333*SQR(rmbar3) + + R4p1313*SQR(n1)*SQR(rmbar3) + R4p2323*SQR(n2)*SQR(rmbar3) + + R4p3333*SQR(n3)*SQR(rmbar3) + Rojo33*SQR(nn)*SQR(rmbar3); + + CCTK_REAL Psi4iL = -(im1*(4*n3*nn*rm1*Ro131 + 2*n3*nn*rm2*Ro132 + + 2*n3*nn*rm3*Ro133 + 2*n3*nn*rm2*Ro231 + n1*(n2*(2*R4p1121*rm1 + + 2*R4p2111*rm1 + R4p1122*rm2 + R4p1221*rm2 + R4p2112*rm2 + R4p2211*rm2 + + R4p1123*rm3 + R4p1321*rm3 + R4p2113*rm3 + R4p2311*rm3) + + n3*(2*R4p1131*rm1 + 2*R4p3111*rm1 + R4p1132*rm2 + R4p1231*rm2 + + R4p3112*rm2 + R4p3211*rm2 + R4p1133*rm3 + R4p1331*rm3 + R4p3113*rm3 + + R4p3311*rm3) + 2*nn*(2*rm1*Ro111 + rm2*(Ro112 + Ro211) + rm3*(Ro113 + + Ro311))) + n2*(n3*(2*R4p2131*rm1 + 2*R4p3121*rm1 + R4p2132*rm2 + + R4p2231*rm2 + R4p3122*rm2 + R4p3221*rm2 + R4p2133*rm3 + R4p2331*rm3 + + R4p3123*rm3 + R4p3321*rm3) + 2*nn*(2*rm1*Ro121 + rm2*(Ro122 + Ro221) + + rm3*(Ro123 + Ro321))) + 2*n3*nn*rm3*Ro331 + (2*R4p1111*rm1 + + R4p1112*rm2 + R4p1211*rm2 + R4p1113*rm3 + R4p1311*rm3)*SQR(n1) + + (2*R4p2121*rm1 + R4p2122*rm2 + R4p2221*rm2 + R4p2123*rm3 + + R4p2321*rm3)*SQR(n2) + 2*R4p3131*rm1*SQR(n3) + R4p3132*rm2*SQR(n3) + + R4p3231*rm2*SQR(n3) + R4p3133*rm3*SQR(n3) + R4p3331*rm3*SQR(n3) + + 2*rm1*Rojo11*SQR(nn) + rm2*Rojo12*SQR(nn) + rm3*Rojo13*SQR(nn) + + rm2*Rojo21*SQR(nn) + rm3*Rojo31*SQR(nn))) - im2*(2*n3*nn*rm1*Ro132 + + 2*n3*nn*rm1*Ro231 + 4*n3*nn*rm2*Ro232 + 2*n3*nn*rm3*Ro233 + + n1*(n2*(R4p1122*rm1 + R4p1221*rm1 + R4p2112*rm1 + R4p2211*rm1 + + 2*R4p1222*rm2 + 2*R4p2212*rm2 + R4p1223*rm3 + R4p1322*rm3 + R4p2213*rm3 + + R4p2312*rm3) + n3*(R4p1132*rm1 + R4p1231*rm1 + R4p3112*rm1 + + R4p3211*rm1 + 2*R4p1232*rm2 + 2*R4p3212*rm2 + R4p1233*rm3 + R4p1332*rm3 + + R4p3213*rm3 + R4p3312*rm3) + 2*nn*(rm1*(Ro112 + Ro211) + 2*rm2*Ro212 + + rm3*(Ro213 + Ro312))) + n2*(n3*(R4p2132*rm1 + R4p2231*rm1 + + R4p3122*rm1 + R4p3221*rm1 + 2*R4p2232*rm2 + 2*R4p3222*rm2 + R4p2233*rm3 + + R4p2332*rm3 + R4p3223*rm3 + R4p3322*rm3) + 2*nn*(rm1*(Ro122 + Ro221) + + 2*rm2*Ro222 + rm3*(Ro223 + Ro322))) + 2*n3*nn*rm3*Ro332 + + (R4p1112*rm1 + R4p1211*rm1 + 2*R4p1212*rm2 + R4p1213*rm3 + + R4p1312*rm3)*SQR(n1) + (R4p2122*rm1 + R4p2221*rm1 + 2*R4p2222*rm2 + + R4p2223*rm3 + R4p2322*rm3)*SQR(n2) + R4p3132*rm1*SQR(n3) + + R4p3231*rm1*SQR(n3) + 2*R4p3232*rm2*SQR(n3) + R4p3233*rm3*SQR(n3) + + R4p3332*rm3*SQR(n3) + rm1*Rojo12*SQR(nn) + rm1*Rojo21*SQR(nn) + + 2*rm2*Rojo22*SQR(nn) + rm3*Rojo23*SQR(nn) + rm3*Rojo32*SQR(nn)) - + im3*(2*n3*nn*rm1*Ro133 + 2*n3*nn*rm2*Ro233 + n1*(n2*(R4p1123*rm1 + + R4p1321*rm1 + R4p2113*rm1 + R4p2311*rm1 + R4p1223*rm2 + R4p1322*rm2 + + R4p2213*rm2 + R4p2312*rm2 + 2*R4p1323*rm3 + 2*R4p2313*rm3) + + n3*(R4p1133*rm1 + R4p1331*rm1 + R4p3113*rm1 + R4p3311*rm1 + R4p1233*rm2 + + R4p1332*rm2 + R4p3213*rm2 + R4p3312*rm2 + 2*R4p1333*rm3 + + 2*R4p3313*rm3) + 2*nn*(rm1*(Ro113 + Ro311) + rm2*(Ro213 + Ro312) + + 2*rm3*Ro313)) + n2*(n3*(R4p2133*rm1 + R4p2331*rm1 + R4p3123*rm1 + + R4p3321*rm1 + R4p2233*rm2 + R4p2332*rm2 + R4p3223*rm2 + R4p3322*rm2 + + 2*R4p2333*rm3 + 2*R4p3323*rm3) + 2*nn*(rm1*(Ro123 + Ro321) + rm2*(Ro223 + + Ro322) + 2*rm3*Ro323)) + 2*n3*nn*rm1*Ro331 + 2*n3*nn*rm2*Ro332 + + 4*n3*nn*rm3*Ro333 + (R4p1113*rm1 + R4p1311*rm1 + R4p1213*rm2 + + R4p1312*rm2 + 2*R4p1313*rm3)*SQR(n1) + (R4p2123*rm1 + R4p2321*rm1 + + R4p2223*rm2 + R4p2322*rm2 + 2*R4p2323*rm3)*SQR(n2) + + R4p3133*rm1*SQR(n3) + R4p3331*rm1*SQR(n3) + R4p3233*rm2*SQR(n3) + + R4p3332*rm2*SQR(n3) + 2*R4p3333*rm3*SQR(n3) + rm1*Rojo13*SQR(nn) + + rm2*Rojo23*SQR(nn) + rm1*Rojo31*SQR(nn) + rm2*Rojo32*SQR(nn) + + 2*rm3*Rojo33*SQR(nn)); /* Copy local copies back to grid functions */ Psi4i[index] = Psi4iL; |