Utilities

StaticConformal provides aliased functions to convert between physical and conformal 3-metric values. It is very important to understand that these functions apply the conversion in place. That is, if gxx contains the conformal metric value, when the routine is exited it will now contain the physical metric value. These functions do not change the value of conformal_state and should be used with due care. (These functions are for example used by some analysis thorns who work only with the physical metric, they apply the transformation on entry to the analysis routine and switch it back on exit).

Convert from conformal to physical:

subroutine ConfToPhysInPlace (nx, ny, nz,
                              psi,
                              gxx, gxy, gxz, gyy, gyz, gzz)
   implicit none
   CCTK_INT,                         intent(in)    :: nx, ny, nz
   CCTK_REAL, dimension(nx, ny, nz), intent(in)    :: psi
   CCTK_REAL, dimension(nx, ny, nz), intent(inout) :: gxx, gxy, gxz, gyy, gyz, gzz
end subroutine ConfToPhysInPlace

Convert from physical to conformal:

subroutine PhysToConfInPlace (nx, ny, nz,
                              psi,
                              gxx, gxy, gxz, gyy, gyz, gzz)
   implicit none
   CCTK_INT,                         intent(in)    :: nx, ny, nz
   CCTK_REAL, dimension(nx, ny, nz), intent(in)    :: psi
   CCTK_REAL, dimension(nx, ny, nz), intent(inout) :: gxx, gxy, gxz, gyy, gyz, gzz
end subroutine ConfToPhysInPlace