if not converged_pmax is True:
            pmax = np.minimum(pmax, b)