pure function minFunD_0(b) result(ret) real(dp), intent(in) :: b real(dp) :: ret ret = -funD_0(b) end function minFunD_0