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/DETG_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/DETG_guts.h')
-rw-r--r-- | src/macro/DETG_guts.h | 57 |
1 files changed, 57 insertions, 0 deletions
diff --git a/src/macro/DETG_guts.h b/src/macro/DETG_guts.h new file mode 100644 index 0000000..a9025ec --- /dev/null +++ b/src/macro/DETG_guts.h @@ -0,0 +1,57 @@ +/* +DETG_guts.h + +Macro to calculate the determinants of the conformal and +physical 3-metric + +Gabrielle Allen, 5th November 1998 + +*/ + +#ifndef DETG_GUTS +#define DETG_GUTS + +#ifdef FCODE + + DETG_PSI4 = DETG_PSI**4 + + DETG_TEMPXX = DETG_GYY*DETG_GZZ - DETG_GYZ*DETG_GYZ + DETG_TEMPXY = DETG_GXZ*DETG_GYZ - DETG_GXY*DETG_GZZ + DETG_TEMPXZ = -DETG_GXZ*DETG_GYY + DETG_GXY*DETG_GYZ + DETG_TEMPYY = DETG_GXX*DETG_GZZ - DETG_GXZ*DETG_GXZ + DETG_TEMPYZ = DETG_GXY*DETG_GXZ - DETG_GXX*DETG_GYZ + DETG_TEMPZZ = DETG_GXX*DETG_GYY - DETG_GXY*DETG_GXY + +/* This is the conformal determinant */ + DETG_DETCG = (DETG_TEMPXX*DETG_GXX + + & DETG_TEMPXY*DETG_GXY + DETG_TEMPXZ*DETG_GXZ) + +/* This is the physical determinat */ + DETG_DETG = DETG_PSI4**3*DETG_DETCG + +#endif + +#ifdef CCODE + +#define Cal(x,y,z) ((x)?(y):(z)) +#define Power(x,y) (pow((double) (x), (double) (y))) + +DET_PSI4 = ((conformal_state != CONFORMAL_METRIC)?1:Power(DETG_PSI,4)); + +DETG_TEMPXX = DETG_GYY*DETG_GZZ - DETG_GYZ*DETG_GYZ; +DETG_TEMPXY = DETG_GXZ*DETG_GYZ - DETG_GXY*DETG_GZZ; +DETG_TEMPXZ = -DETG_GXZ*DETG_GYY + DETG_GXY*DETG_GYZ; +DETG_TEMPYY = DETG_GXX*DETG_GZZ - DETG_GXZ*DETG_GXZ; +DETG_TEMPYZ = DETG_GXY*DETG_GXZ - DETG_GXX*DETG_GYZ; +DETG_TEMPZZ = DETG_GXX*DETG_GYY - DETG_GXY*DETG_GXY; + +/* This is the conformal determinant */ +DETG_DETCG = (DETG_TEMPXX*DETG_GXX + + DETG_TEMPXY*DETG_GXY + DETG_TEMPXZ*DETG_GXZ); + +/* This is the physical determinat */ +DETG_DETG = DET_PSI4*DET_PSI4*DET_PSI4*DETG_DETCG + +#endif + +#endif |