diff options
Diffstat (limited to 'src/macro/DYYDG_guts.h')
| -rw-r--r-- | src/macro/DYYDG_guts.h | 49 |
1 files changed, 35 insertions, 14 deletions
diff --git a/src/macro/DYYDG_guts.h b/src/macro/DYYDG_guts.h index 1896229..2dc495e 100644 --- a/src/macro/DYYDG_guts.h +++ b/src/macro/DYYDG_guts.h @@ -20,33 +20,54 @@ #ifdef FCODE - DYYDG_OODY2 = 1D0/(dy*dy) +#include "ADM_Derivative.h" /* Factor involving 2nd derivative of conformal factor */ IF (conformal_state .eq. 0) THEN DYYDG_FAC = 0 ELSE - DYYDG_FAC = DYDG_PSI4*(4*DYYDG_DYYDPSI_O_PSI + 12*DYDG_DYDPSI_O_PSI*DYDG_DYDPSI_O_PSI) + DYYDG_FAC = DYDG_PSI4*(4*DYYDG_DYYDPSI_O_PSI \ + + 12*DYDG_DYDPSI_O_PSI*DYDG_DYDPSI_O_PSI) ENDIF /* Now calculate the second deriatives */ - DYYDG_DYYDGXX = 2*DYDCG_DYDCGXX*DYDG_FAC+DYYDG_FAC*DYDG_GXX+DYDG_PSI4 \ - *DYYDG_OODY2*(DYDCG_GXX_JP-2*DYDG_GXX+DYDCG_GXX_JM) + if (spatial_order.eq.2) then + DYYDG_DYYDGXX = 2*DYDCG_DYDCGXX*DYDG_FAC + DYYDG_FAC*DYDG_GXX \ + + DYDG_PSI4*ADM_DYY_2(gxx,i,j,k) - DYYDG_DYYDGXY = 2*DYDCG_DYDCGXY*DYDG_FAC+DYYDG_FAC*DYDG_GXY+DYDG_PSI4 \ - *DYYDG_OODY2*(DYDCG_GXY_JP-2*DYDG_GXY+DYDCG_GXY_JM) + DYYDG_DYYDGXY = 2*DYDCG_DYDCGXY*DYDG_FAC + DYYDG_FAC*DYDG_GXY \ + + DYDG_PSI4*ADM_DYY_2(gxy,i,j,k) - DYYDG_DYYDGXZ = 2*DYDCG_DYDCGXZ*DYDG_FAC+DYYDG_FAC*DYDG_GXZ+DYDG_PSI4 \ - *DYYDG_OODY2*(DYDCG_GXZ_JP-2*DYDG_GXZ+DYDCG_GXZ_JM) + DYYDG_DYYDGXZ = 2*DYDCG_DYDCGXZ*DYDG_FAC + DYYDG_FAC*DYDG_GXZ \ + + DYDG_PSI4*ADM_DYY_2(gxz,i,j,k) - DYYDG_DYYDGYY = 2*DYDCG_DYDCGYY*DYDG_FAC+DYYDG_FAC*DYDG_GYY+DYDG_PSI4 \ - *DYYDG_OODY2*(DYDCG_GYY_JP-2*DYDG_GYY+DYDCG_GYY_JM) + DYYDG_DYYDGYY = 2*DYDCG_DYDCGYY*DYDG_FAC + DYYDG_FAC*DYDG_GYY \ + + DYDG_PSI4*ADM_DYY_2(gyy,i,j,k) - DYYDG_DYYDGYZ = 2*DYDCG_DYDCGYZ*DYDG_FAC+DYYDG_FAC*DYDG_GYZ+DYDG_PSI4 \ - *DYYDG_OODY2*(DYDCG_GYZ_JP-2*DYDG_GYZ+DYDCG_GYZ_JM) + DYYDG_DYYDGYZ = 2*DYDCG_DYDCGYZ*DYDG_FAC + DYYDG_FAC*DYDG_GYZ \ + + DYDG_PSI4*ADM_DYY_2(gyz,i,j,k) - DYYDG_DYYDGZZ = 2*DYDCG_DYDCGZZ*DYDG_FAC+DYYDG_FAC*DYDG_GZZ+DYDG_PSI4 \ - *DYYDG_OODY2*(DYDCG_GZZ_JP-2*DYDG_GZZ+DYDCG_GZZ_JM) + DYYDG_DYYDGZZ = 2*DYDCG_DYDCGZZ*DYDG_FAC + DYYDG_FAC*DYDG_GZZ \ + + DYDG_PSI4*ADM_DYY_2(gzz,i,j,k) + else + DYYDG_DYYDGXX = 2*DYDCG_DYDCGXX*DYDG_FAC + DYYDG_FAC*DYDG_GXX \ + + DYDG_PSI4*ADM_DYY_4(gxx,i,j,k) + + DYYDG_DYYDGXY = 2*DYDCG_DYDCGXY*DYDG_FAC + DYYDG_FAC*DYDG_GXY \ + + DYDG_PSI4*ADM_DYY_4(gxy,i,j,k) + + DYYDG_DYYDGXZ = 2*DYDCG_DYDCGXZ*DYDG_FAC + DYYDG_FAC*DYDG_GXZ \ + + DYDG_PSI4*ADM_DYY_4(gxz,i,j,k) + + DYYDG_DYYDGYY = 2*DYDCG_DYDCGYY*DYDG_FAC + DYYDG_FAC*DYDG_GYY \ + + DYDG_PSI4*ADM_DYY_4(gyy,i,j,k) + + DYYDG_DYYDGYZ = 2*DYDCG_DYDCGYZ*DYDG_FAC + DYYDG_FAC*DYDG_GYZ \ + + DYDG_PSI4*ADM_DYY_4(gyz,i,j,k) + + DYYDG_DYYDGZZ = 2*DYDCG_DYDCGZZ*DYDG_FAC + DYYDG_FAC*DYDG_GZZ \ + + DYDG_PSI4*ADM_DYY_4(gzz,i,j,k) + end if #endif |