diff options
| author | evans <evans@b1d164ef-f17a-46e7-89d4-021c7118ef4e> | 1999-03-09 21:04:02 +0000 |
|---|---|---|
| committer | evans <evans@b1d164ef-f17a-46e7-89d4-021c7118ef4e> | 1999-03-09 21:04:02 +0000 |
| commit | 0a13cbb3d0f5ebf208d962d374a197c66bd13b1b (patch) | |
| tree | 19582fa856a1193d17170ad971353f58e7be36c9 /src/macro/DZZDG_guts.h | |
| parent | c58897b8676212b9aee16d6d2f80ca2adeb50c6c (diff) | |
added the Einstien macros.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinBase/ADMMacros/trunk@2 b1d164ef-f17a-46e7-89d4-021c7118ef4e
Diffstat (limited to 'src/macro/DZZDG_guts.h')
| -rw-r--r-- | src/macro/DZZDG_guts.h | 79 |
1 files changed, 79 insertions, 0 deletions
diff --git a/src/macro/DZZDG_guts.h b/src/macro/DZZDG_guts.h new file mode 100644 index 0000000..ef0a946 --- /dev/null +++ b/src/macro/DZZDG_guts.h @@ -0,0 +1,79 @@ +/*@@ + @header DZZDG_guts.h + @date Jun 98 + @author Gabrielle Allen + @desc + Macro to calculate the (first and) second derivatives of the + physical metric with respect to z + + The macro is defined in terms of standard variables in + @seefile DZZDG_declare.h + + The macro uses @seefile DZDG_guts.h and @seefile DZDG_declare.h + @enddesc +@@*/ + +#ifndef DZZDG_GUTS +#define DZZDG_GUTS + +#include "../../packages/CactusEinstein/Einstein/src/macro/DZDG_guts.h" + +#ifdef FCODE + + DZZDG_OODZ2 = 1D0/(dz*dz) + + /* Factor involving 2nd derivative of conformal factor */ + IF (conformal_state /= CONFORMAL_METRIC) THEN + DZZDG_FAC = 0 + ELSE + DZZDG_FAC = DZDG_PSI4*(4*DZZDG_DZZDPSI_O_PSI + 12*DZDG_DZDPSI_O_PSI*DZDG_DZDPSI_O_PSI) + ENDIF + + DZZDG_DZZDGXX = 2*DZDCG_DZDCGXX*DZDG_FAC+DZZDG_FAC*DZDG_GXX+DZDG_PSI4 + &*DZZDG_OODZ2*(DZDCG_GXX_KP-2*DZDG_GXX+DZDCG_GXX_KM) + + DZZDG_DZZDGXY = 2*DZDCG_DZDCGXY*DZDG_FAC+DZZDG_FAC*DZDG_GXY+DZDG_PSI4 + &*DZZDG_OODZ2*(DZDCG_GXY_KP-2*DZDG_GXY+DZDCG_GXY_KM) + + DZZDG_DZZDGXZ = 2*DZDCG_DZDCGXZ*DZDG_FAC+DZZDG_FAC*DZDG_GXZ+DZDG_PSI4 + &*DZZDG_OODZ2*(DZDCG_GXZ_KP-2*DZDG_GXZ+DZDCG_GXZ_KM) + + DZZDG_DZZDGYY = 2*DZDCG_DZDCGYY*DZDG_FAC+DZZDG_FAC*DZDG_GYY+DZDG_PSI4 + &*DZZDG_OODZ2*(DZDCG_GYY_KP-2*DZDG_GYY+DZDCG_GYY_KM) + + DZZDG_DZZDGYZ = 2*DZDCG_DZDCGYZ*DZDG_FAC+DZZDG_FAC*DZDG_GYZ+DZDG_PSI4 + &*DZZDG_OODZ2*(DZDCG_GYZ_KP-2*DZDG_GYZ+DZDCG_GYZ_KM) + + DZZDG_DZZDGZZ = 2*DZDCG_DZDCGZZ*DZDG_FAC+DZZDG_FAC*DZDG_GZZ+DZDG_PSI4 + &*DZZDG_OODZ2*(DZDCG_GZZ_KP-2*DZDG_GZZ+DZDCG_GZZ_KM) + +#endif + +#ifdef CCODE + + /* Factor involving 2nd derivative of conformal factor */ + DZZDG_FAC = ((conformal_state != CONFORMAL_METRIC)?0: + DZDG_PSI4*(4*DZZDG_DZZDPSI_O_PSI + 12*DZDG_DZDPSI_O_PSI*DZDG_DZDPSI_O_PSI)); + + /* Now calculate the second deriatives */ + DZZDG_DZZDGXX = 2*DZDCG_DZDCGXX*DZDG_FAC+DZZDG_FAC*DZDG_GXX+DZDG_PSI4 + *DZZDG_OODZ2*(DZDCG_GXX_KP-2*DZDG_GXX+DZDCG_GXX_KM); + + DZZDG_DZZDGXY = 2*DZDCG_DZDCGXY*DZDG_FAC+DZZDG_FAC*DZDG_GXY+DZDG_PSI4 + *DZZDG_OODZ2*(DZDCG_GXY_KP-2*DZDG_GXY+DZDCG_GXY_KM); + + DZZDG_DZZDGXZ = 2*DZDCG_DZDCGXZ*DZDG_FAC+DZZDG_FAC*DZDG_GXZ+DZDG_PSI4 + *DZZDG_OODZ2*(DZDCG_GXZ_KP-2*DZDG_GXZ+DZDCG_GXZ_KM); + + DZZDG_DZZDGYY = 2*DZDCG_DZDCGYY*DZDG_FAC+DZZDG_FAC*DZDG_GYY+DZDG_PSI4 + *DZZDG_OODZ2*(DZDCG_GYY_KP-2*DZDG_GYY+DZDCG_GYY_KM); + + DZZDG_DZZDGYZ = 2*DZDCG_DZDCGYZ*DZDG_FAC+DZZDG_FAC*DZDG_GYZ+DZDG_PSI4 + *DZZDG_OODZ2*(DZDCG_GYZ_KP-2*DZDG_GYZ+DZDCG_GYZ_KM); + + DZZDG_DZZDGZZ = 2*DZDCG_DZDCGZZ*DZDG_FAC+DZZDG_FAC*DZDG_GZZ+DZDG_PSI4 + *DZZDG_OODZ2*(DZDCG_GZZ_KP-2*DZDG_GZZ+DZDCG_GZZ_KM); + +#endif + +#endif |