function full_H_combined = get_full_H_with_diagonal(obj)
% Get the full approximate Hessian, including diagonal terms that are included in obj.full_H
    full_H_combined = obj.full_H + eye(obj.K_max).*sum(obj.min_eig_sub(obj.active));
end
