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/DKDT_declare.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/DKDT_declare.h')
-rw-r--r-- | src/macro/DKDT_declare.h | 111 |
1 files changed, 111 insertions, 0 deletions
diff --git a/src/macro/DKDT_declare.h b/src/macro/DKDT_declare.h new file mode 100644 index 0000000..b03ad4c --- /dev/null +++ b/src/macro/DKDT_declare.h @@ -0,0 +1,111 @@ +/*@@ + @header DKDT_declare.h + @date Jul 98 + @author Gabrielle Allen + @desc + Declarations for macro @seefile DKDT_guts.h to calculate the source term + in the evolution equation for the extrinsic curvature. That is + + d K_ij/dt = alpha*(R_ij-2 K_ik K^k_j + K_ij trK] - nabla_i nabla_j alpha + + + Lie_beta K_ij + + @enddesc +@@*/ + +#ifndef DKDT_DECLARE +#define DKDT_DECLARE + +#include "../../packages/CactusEinstein/Einstein/src/macro/RICCI_declare.h" +#include "../../packages/CactusEinstein/Einstein/src/macro/KK_declare.h" +#include "../../packages/CactusEinstein/Einstein/src/macro/TRK_declare.h" +#include "../../packages/CactusEinstein/Einstein/src/macro/CDCDA_declare.h" +#include "../../packages/CactusEinstein/Einstein/src/macro/LIEK_declare.h" + + +#ifdef FCODE + +/* Input variables */ +#undef DKDT_A +#define DKDT_A alp(i,j,k) +#undef DKDT_KXX +#define DKDT_KXX hxx(i,j,k) +#undef DKDT_KXY +#define DKDT_KXY hxy(i,j,k) +#undef DKDT_KXZ +#define DKDT_KXZ hxz(i,j,k) +#undef DKDT_KYY +#define DKDT_KYY hyy(i,j,k) +#undef DKDT_KYZ +#define DKDT_KYZ hyz(i,j,k) +#undef DKDT_KZZ +#define DKDT_KZZ hzz(i,j,k) + +/* Output variables */ +#undef DKDT_DKXXDT +#define DKDT_DKXXDT dkdt_dkxxdt +#undef DKDT_DKXYDT +#define DKDT_DKXYDT dkdt_dkxydt +#undef DKDT_DKXZDT +#define DKDT_DKXZDT dkdt_dkxzdt +#undef DKDT_DKYYDT +#define DKDT_DKYYDT dkdt_dkyydt +#undef DKDT_DKYZDT +#define DKDT_DKYZDT dkdt_dkyzdt +#undef DKDT_DKZZDT +#define DKDT_DKZZDT dkdt_dkzzdt + +/* Declare output variables */ + REAL DKDT_DKXXDT + REAL DKDT_DKXYDT + REAL DKDT_DKXZDT + REAL DKDT_DKYYDT + REAL DKDT_DKYZDT + REAL DKDT_DKZZDT + +#endif + + +#ifdef CCODE + +/* Input variables */ +#undef DKDT_A +#define DKDT_A alp[ijk] +#undef DKDT_KXX +#define DKDT_KXX hxx[ijk] +#undef DKDT_KXY +#define DKDT_KXY hxy[ijk] +#undef DKDT_KXZ +#define DKDT_KXZ hxz[ijk] +#undef DKDT_KYY +#define DKDT_KYY hyy[ijk] +#undef DKDT_KYZ +#define DKDT_KYZ hyz[ijk] +#undef DKDT_KZZ +#define DKDT_KZZ hzz[ijk] + +/* Output variables */ +#undef DKDT_DKXXDT +#define DKDT_DKXXDT dkdt_dkxxdt +#undef DKDT_DKXYDT +#define DKDT_DKXYDT dkdt_dkxydt +#undef DKDT_DKXZDT +#define DKDT_DKXZDT dkdt_dkxzdt +#undef DKDT_DKYYDT +#define DKDT_DKYYDT dkdt_dkyydt +#undef DKDT_DKYZDT +#define DKDT_DKYZDT dkdt_dkyzdt +#undef DKDT_DKZZDT +#define DKDT_DKZZDT dkdt_dkzzdt + +/* Declare output variables */ + double DKDT_DKXXDT; + double DKDT_DKXYDT; + double DKDT_DKXZDT; + double DKDT_DKYYDT; + double DKDT_DKYZDT; + double DKDT_DKZZDT; + +#endif + +#endif |