ENH: CONESTA use derived_eps = max(eps, self.eps) - mu * gM instead of derived_eps = eps - mu * gM