lin_evol_primal_dual_gen_reduced_data.m

Go to the documentation of this file.

function reduced_data = lin_evol_primal_dual_gen_reduced_data(model, detailed_data)
%reduced_data = lin_evol_primal_dual_gen_reduced_data(model, detailed_data)
%
% function generating both primal improved and dual reduced_data. Further
% writes both reduced_data_primal_improved and reduced_data_dual into one 
% structure: reduced_data.
% Used for the improved output and output estimator using the dual problem.

% Domink Garmatter 06.09 2012


% extract the dual detailed_data out of the connected detailed_data
[~, detailed_data_dual] = lin_evol_split_detailed_data(detailed_data);
% generate dual reduced_data
model.want_dual = 1;
model.want_improved_output = 0;
reduced_data_dual = lin_evol_gen_reduced_data_primal_dual(model, detailed_data_dual);
% generated primal improved reduced_data using the connected detailed_data
model.want_dual = 0;
model.want_improved_output = 1;
reduced_data_primal_improved = lin_evol_gen_reduced_data_primal_dual(model, detailed_data);
% connect both reduced data
reduced_data = [];
reduced_data = structcpy(reduced_data, reduced_data_dual); % simply copy the dual reduced_data
% and write alle the primal improved reduced_data - these fields are always like this
reduced_data.primal_improved_LL_I = reduced_data_primal_improved.LL_I;
reduced_data.primal_improved_LL_E = reduced_data_primal_improved.LL_E;
reduced_data.primal_improved_bb = reduced_data_primal_improved.bb;
reduced_data.primal_improved_K_II = reduced_data_primal_improved.K_II;
reduced_data.primal_improved_K_IE = reduced_data_primal_improved.K_IE;
reduced_data.primal_improved_K_EE = reduced_data_primal_improved.K_EE;
reduced_data.primal_improved_m_I = reduced_data_primal_improved.m_I;
reduced_data.primal_improved_m_E = reduced_data_primal_improved.m_E;
reduced_data.primal_improved_m = reduced_data_primal_improved.m;
reduced_data.primal_improved_LL_I_correct = reduced_data_primal_improved.LL_I_correct;
reduced_data.primal_improved_LL_E_correct = reduced_data_primal_improved.LL_E_correct;
reduced_data.primal_improved_bb_correct = reduced_data_primal_improved.bb_correct;
reduced_data.primal_improved_N = reduced_data_primal_improved.N;
reduced_data.primal_improved_a0 = reduced_data_primal_improved.a0;
reduced_data.primal_improved_s_RB = reduced_data_primal_improved.s_RB;
if isfield(reduced_data_primal_improved, 'scm_offline_data')
    reduced_data.primal_improved_scm_offline_data = reduced_data_primal_improved.scm_offline_data;
end