if pmax is None:
            pmax, converged_pmax = self._search_p_min_max(phatfree0, p_up, p_opt, 'max')