diff options
Diffstat (limited to 'src/macro/DZZDG_guts.h')
| -rw-r--r-- | src/macro/DZZDG_guts.h | 14 |
1 files changed, 4 insertions, 10 deletions
diff --git a/src/macro/DZZDG_guts.h b/src/macro/DZZDG_guts.h index 1d2242e..a0b53b4 100644 --- a/src/macro/DZZDG_guts.h +++ b/src/macro/DZZDG_guts.h @@ -22,13 +22,8 @@ #include "ADM_Derivative.h" - /* Factor involving 2nd derivative of conformal factor */ - IF (conformal_state .eq. 0) THEN - DZZDG_FAC = 0 - ELSE - DZZDG_FAC = DZDG_PSI4*(4*DZZDG_DZZDPSI_O_PSI \ - + 12*DZDG_DZDPSI_O_PSI*DZDG_DZDPSI_O_PSI) - ENDIF + /* Factor involving 2nd derivative of conformal factor, nowadays zero */ + DZZDG_FAC = 0 if (local_spatial_order.eq.2) then DZZDG_DZZDGXX = 2*DZDCG_DZDCGXX*DZDG_FAC + DZZDG_FAC*DZDG_GXX \ @@ -71,9 +66,8 @@ #ifdef CCODE - /* Factor involving 2nd derivative of conformal factor */ - DZZDG_FAC = ((*conformal_state == 0)?0: - DZDG_PSI4*(4*DZZDG_DZZDPSI_O_PSI + 12*DZDG_DZDPSI_O_PSI*DZDG_DZDPSI_O_PSI)); + /* Factor involving 2nd derivative of conformal factor, nowadays zero */ + DZZDG_FAC = 0; /* Now calculate the second deriatives */ DZZDG_DZZDGXX = 2*DZDCG_DZDCGXX*DZDG_FAC+DZZDG_FAC*DZDG_GXX+DZDG_PSI4 |