porsche_circle_x3opt_y3opt_r10_get_Hessian_eq_constr.m

function hessian_H_j = porsche_circle_x3opt_y3opt_r10_get_Hessian_eq_constr(model,data)
%function hessian_H_j = porsche_get_Hessian_eq_constr(model,data)
%
% This function calculates the matrices \nabla^2 h_j in my diploma-thesis,
% either detailed or reduced according to model.optimimzation.opt_mode.
% The corresponding data (model_data or reduced_data) must be passed to
% this function.
%
% 

nr_constr = model.optimization.nr_constraints;
hessian_H_j = cell(1,nr_constr);
l=1;
for k=1:nr_constr
    hessian_H_j{k} = zeros(length(model.mu_names));
    hessian_H_j{k}(l,l) = 2;
    hessian_H_j{k}(l+1,l+1) = 2;
    l=l+2;
end