pure function funOverB_0(b) result(ret) real(dp), intent(in) :: b real(dp) :: ret !ASSERT( b > 0.0_dp) ret = funOverB(b, 0.0_dp) end function funOverB_0