pure function funB_m1_0(b) result(ret) real(dp), intent(in) :: b real(dp) :: ret ret = funB(b, -1.0_dp, 0.0_dp) end function funB_m1_0