diff options
Diffstat (limited to 'src/macro/DCGDT_guts.h')
-rw-r--r-- | src/macro/DCGDT_guts.h | 44 |
1 files changed, 22 insertions, 22 deletions
diff --git a/src/macro/DCGDT_guts.h b/src/macro/DCGDT_guts.h index a88a831..e950f3e 100644 --- a/src/macro/DCGDT_guts.h +++ b/src/macro/DCGDT_guts.h @@ -55,38 +55,38 @@ #ifdef CCODE - DCGDT_DCGXXDT = - 2*DCGDT_A*DCGDT_KXX - DCGDT_DCGXYDT = - 2*DCGDT_A*DCGDT_KXY - DCGDT_DCGXZDT = - 2*DCGDT_A*DCGDT_KXZ - DCGDT_DCGYYDT = - 2*DCGDT_A*DCGDT_KYY - DCGDT_DCGYZDT = - 2*DCGDT_A*DCGDT_KYZ - DCGDT_DCGZZDT = - 2*DCGDT_A*DCGDT_KZZ - - if (conformal_state == CONFORMAL_METRIC) + DCGDT_DCGXXDT = - 2*DCGDT_A*DCGDT_KXX; + DCGDT_DCGXYDT = - 2*DCGDT_A*DCGDT_KXY; + DCGDT_DCGXZDT = - 2*DCGDT_A*DCGDT_KXZ; + DCGDT_DCGYYDT = - 2*DCGDT_A*DCGDT_KYY; + DCGDT_DCGYZDT = - 2*DCGDT_A*DCGDT_KYZ; + DCGDT_DCGZZDT = - 2*DCGDT_A*DCGDT_KZZ; + + if (*conformal_state == CONFORMAL_METRIC) { - DCGDT_IPSI4 = 1D0/(DCGDT_PSI**4) + DCGDT_IPSI4 = 1/(DCGDT_PSI*DCGDT_PSI*DCGDT_PSI*DCGDT_PSI); - DCGDT_DCGXXDT = DCGDT_DCGXXDT*DCGDT_IPSI4 - DCGDT_DCGXYDT = DCGDT_DCGXYDT*DCGDT_IPSI4 - DCGDT_DCGXZDT = DCGDT_DCGXZDT*DCGDT_IPSI4 - DCGDT_DCGYYDT = DCGDT_DCGYYDT*DCGDT_IPSI4 - DCGDT_DCGYZDT = DCGDT_DCGYZDT*DCGDT_IPSI4 - DCGDT_DCGZZDT = DCGDT_DCGZZDT*DCGDT_IPSI4 + DCGDT_DCGXXDT = DCGDT_DCGXXDT*DCGDT_IPSI4; + DCGDT_DCGXYDT = DCGDT_DCGXYDT*DCGDT_IPSI4; + DCGDT_DCGXZDT = DCGDT_DCGXZDT*DCGDT_IPSI4; + DCGDT_DCGYYDT = DCGDT_DCGYYDT*DCGDT_IPSI4; + DCGDT_DCGYZDT = DCGDT_DCGYZDT*DCGDT_IPSI4; + DCGDT_DCGZZDT = DCGDT_DCGZZDT*DCGDT_IPSI4; } - if (shift) + if (*shift_state != SHIFT_INACTIVE) { #include "CactusEinstein/Einstein/src/macro/LIEG_guts.h" - DCGDT_DCGXXDT = DCGDT_DCGXXDT + LIEG_LGXX - DCGDT_DCGXYDT = DCGDT_DCGXYDT + LIEG_LGXY - DCGDT_DCGXZDT = DCGDT_DCGXZDT + LIEG_LGXZ - DCGDT_DCGYYDT = DCGDT_DCGYYDT + LIEG_LGYY - DCGDT_DCGYZDT = DCGDT_DCGYZDT + LIEG_LGYZ - DCGDT_DCGZZDT = DCGDT_DCGZZDT + LIEG_LGZZ + DCGDT_DCGXXDT = DCGDT_DCGXXDT + LIEG_LGXX; + DCGDT_DCGXYDT = DCGDT_DCGXYDT + LIEG_LGXY; + DCGDT_DCGXZDT = DCGDT_DCGXZDT + LIEG_LGXZ; + DCGDT_DCGYYDT = DCGDT_DCGYYDT + LIEG_LGYY; + DCGDT_DCGYZDT = DCGDT_DCGYZDT + LIEG_LGYZ; + DCGDT_DCGZZDT = DCGDT_DCGZZDT + LIEG_LGZZ; } |