Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=dp), | intent(in) | :: | b | |||
real(kind=dp), | intent(in) | :: | x | |||
real(kind=dp), | intent(in) | :: | y |
pure function funB(b, x, y) result(ret)
real(dp), intent(in) :: b, x, y
real(dp) :: ret
!ASSERT( (b + x) > 0.0_dp)
!ASSERT( (b + y) > 0.0_dp)
ret = sqrt(2.0_dp / ((b + x) * (b + y)))
end function funB