IntK Function

private function IntK(v, x)

Arguments

Type IntentOptional AttributesName
real(kind=rp), intent(in) :: v
real(kind=rp), intent(in) :: x

Return Value real(kind=rp)


Called by

proc~~intk~2~~CalledByGraph proc~intk~2 IntK proc~p_integral~2 P_integral proc~p_integral~2->proc~intk~2 proc~pr~2 PR proc~pr~2->proc~p_integral~2 proc~frexpr fRExPR proc~frexpr->proc~pr~2 proc~fre_hxpr fRE_HxPR proc~fre_hxpr->proc~pr~2

Contents

Source Code


Source Code

  FUNCTION IntK(v,x)
    REAL(rp) :: IntK
    REAL(rp), INTENT(IN) :: v
    REAL(rp), INTENT(IN) :: x

    IntK = (C_PI/SQRT(2.0_rp))*(1.0_rp - 0.25_rp*(4.0_rp*v**2 - 1.0_rp))*ERFC(SQRT(x))&
         + 0.25_rp*(4.0_rp*v**2 - 1.0_rp)*SQRT(0.5_rp*C_PI/x)*EXP(-x)
  END FUNCTION IntK