На самом деле я понял это сам. Это довольно просто:
def optim2 (obj_fun c, initial_theta, bounds): (1) optimResult = opt.minimize (obj_fun c, initial_theta, method = 'BFGS', jac = 'False') theta_opt = optimResult.x func_min = optimResult.fun возвращает theta_opt, func_min
OR
(2) theta_opt, func_min, Convergence_dict = opt.fmin_l_bfgs_b (obj_fun c, bounds), ) вернуть theta_opt, func_min
gp = GaussianProcessRegressor (ядро = ядро, альфа = 0,015, normalize_y = True, оптимизатор = optim2)