res = minimize(ConstrainedREPS._dual_function, eta_start,
                       jac=ConstrainedREPS._dual_function_diff,
                       bounds=((np.finfo(np.float32).eps, np.inf),),
                       args=(self._eps(), Jep, theta))