vi_detailed_simulation.m

function sim_data = vi_detailed_simulation(model,model_data)
%function sim_data = vi_detailed_simulation(model,model_data)
%
% function performing a detailed simulation of a vi-model (see
% elastic_rope_model) via "quadprog".
%
% Generated fields of sim_data:
%  U: primal solution vector
%  L: Lagrange multiplier vector

% I.Maier 30.05.2011

sim_data = [];

model.decomp_mode = 0;
[A,B,f,g] = model.operators(model,model_data);

warning off
options = optimset;
%options.TolFun = 0;
options.MaxIter = 1000; % required for large bases...
options.Display = 'off';
sim_data.U = quadprog(A,-f,B',g,[],[],[],[],[],options);
%[sim_data.U,FVAL,EXITFLAG,OUTPUT] = quadprog(A,-f,B',g,[],[],[],[],[],options);
warning on

sim_data.L = B\(f - A*sim_data.U);
sim_data.L = max(sim_data.L,0);
sim_data.A = A;
sim_data.B = B;
sim_data.f = f;
sim_data.g = g;