diff options
Diffstat (limited to 'src/ricci.F90')
| -rw-r--r-- | src/ricci.F90 | 10 |
1 files changed, 9 insertions, 1 deletions
diff --git a/src/ricci.F90 b/src/ricci.F90 index 83b144a..6575994 100644 --- a/src/ricci.F90 +++ b/src/ricci.F90 @@ -2,6 +2,8 @@ #include "cctk.h" +#undef DEBUG + module ricci implicit none private @@ -82,7 +84,9 @@ contains subroutine calc_ricci (gamma, dgamma, ri) CCTK_REAL, intent(in) :: gamma(3,3,3), dgamma(3,3,3,3) CCTK_REAL, intent(out) :: ri(3,3) - CCTK_REAL :: sum, cnt +#ifdef DEBUG + CCTK_REAL :: nrm, cnt +#endif integer :: i,j,k,l ! R_ij = Gamma^k_ij,k - Gamma^k_ik,j ! + Gamma^k_lk Gamma^l_ij - Gamma^k_lj Gamma^l_ki @@ -98,6 +102,7 @@ contains end do end do end do +#ifdef DEBUG ! check symmetries sum = 0 cnt = 0 @@ -108,12 +113,15 @@ contains end do end do if (sum > 1.0e-12 * cnt) call CCTK_WARN (0, "Ricci tensor is not symmetric") +#endif end subroutine calc_ricci subroutine calc_ricciderivs (gamma, dgamma, ddgamma, dri) CCTK_REAL, intent(in) :: gamma(3,3,3), dgamma(3,3,3,3), ddgamma(3,3,3,3,3) CCTK_REAL, intent(out) :: dri(3,3,3) +#if 0 CCTK_REAL :: nrm, cnt +#endif integer :: i,j,k,l,m ! R_ij = Gamma^k_ij,k - Gamma^k_ik,j ! + Gamma^k_lk Gamma^l_ij - Gamma^k_lj Gamma^l_ki |