get_ham Subroutine

subroutine get_ham(tpof, i, j, k, l, hc)

Wrapper SUBROUTINE to get the value of the Hamiltonian constraint computed using ID read on the refined mesh, at a given mesh point (i,j,k), on the refinement level l, for the tpo object tpof

Arguments

Type IntentOptional Attributes Name
class(tpo), intent(inout) :: tpof

tpo object to use
integer, intent(in) :: i

index of the desired point
integer, intent(in) :: j

index of the desired point
integer, intent(in) :: k

index of the desired point
integer, intent(in) :: l

Index of the refinement level
double precision, intent(out) :: hc

Value of the Hamiltonian constraint at


Contents

Source Code

get_ham

Source Code

  SUBROUTINE get_ham( tpof, i, j, k, l, hc )
  !# Wrapper SUBROUTINE to get the value of the Hamiltonian constraint
  !  computed using ID read on the refined mesh,
  !  at a given mesh point `(i,j,k)`, on the refinement level `l`, for the
  !  [[tpo]] object `tpof`

    IMPLICIT NONE

    CLASS(tpo), INTENT(INOUT):: tpof
    !! [[tpo]] object to use
    INTEGER, INTENT(IN):: i
    !! \(x\) index of the desired point
    INTEGER, INTENT(IN):: j
    !! \(y\) index of the desired point
    INTEGER, INTENT(IN):: k
    !! \(z\) index of the desired point
    INTEGER, INTENT(IN):: l
    !! Index of the refinement level
    DOUBLE PRECISION, INTENT(OUT):: hc
    !! Value of the Hamiltonian constraint at \((i,j,k)\)

    hc= tpof% get_HC( i, j, k, l )

  END SUBROUTINE get_ham