lin_evol_primal_dual_rb_simulation.m

Go to the documentation of this file.

function simulation_data = lin_evol_primal_dual_rb_simulation(model, reduced_data)
%simulation_data = lin_evol_primal_dual_rb_simulation(model, reduced_data)
%
% function which first computes a dual rb_simulation using dual
% reduced_data and afterwards computes a primal improved rb_simulation
% using primal improved reduced_data. This finally yields the improved
% output and improved output estimator!
%
% simulation_data yields the final results:
% - primal rb coefficients
% - primal energy error estimator
% - improved output
% - improved output estimator

% Dominik Garmatter 06.09 2012


% first split up the reduced_data in the dual and primal improved parts
[reduced_data_primal_improved, reduced_data_dual] = lin_evol_split_reduced_data(reduced_data);
% you want the energy error for the improved output estimator
if ~strcmp(model.error_norm ,'energy')
    warning('Primal dual error estimator only implemented for energy error norm');
    model.error_norm = 'energy';
end
% do a dual rb_simulation to recieve the dual rb coefficients and the dual
% energy error estimator - both required for the improved output and
% improved output estimator
model.want_dual = 1;
rb_sim_data_dual = lin_evol_rb_simulation_primal_dual(model, reduced_data_dual);
% now do an improved primal rb_simulation
model.want_dual = 0;
model.want_improved_output = 1;
% the dual rb coefficients and the dual energy error estimator are required
% for the improved simulation
reduced_data_primal_improved.dual_sim_data = rb_sim_data_dual.a;
reduced_data_primal_improved.dual_Delta = rb_sim_data_dual.Delta;
simulation_data = lin_evol_rb_simulation_primal_dual(model, reduced_data_primal_improved);