lin_evol_primal_dual_gen_detailed_data.m

Go to the documentation of this file.

function detailed_data = lin_evol_primal_dual_gen_detailed_data(model, model_data)
%detailed_data = lin_evol_primal_dual_gen_detailed_data(model, model_data)
%
% function generating both primal and dual detailed_data. Further writes
% both detailed_data_primal and detailed_data_dual into one structure:
% detailed_data. Used for the improved output and output estimator using
% the dual problem.
%
% special Note: in the european_option_pricing model a functional named
% 'theta' includes a partial timederivative and therefore produces a
% slightly different dual problem. Therefore the 'theta'-case is sometimes
% treated differently. Don't care about this case if you don't use the
% european_option_pricing model.

% Domink Garmatter 06.09 2012


% always use the energy error estimator for basis generation (remeber your
% final aim is the improved output and improved output estimator, which
% also uses the energy error estimator)
%model.RB_error_indicator = 'estimator';
%model.error_norm = 'energy';
% generate primal detailed_data
model.want_dual = 0;
detailed_data_primal = gen_detailed_data(model, model_data);
% generate dual detailed_data
model.want_dual = 1;
% if there are model fiels specifying the dual RB generation use them. 
% Otherwise the primal fields are automatically used. 
if isfield(model, 'dual_RB_stop_epsilon')
    model.RB_stop_epsilon = model.dual_RB_stop_epsilon;
end
if isfield(model, 'dual_RB_stop_Nmax')
    model.RB_stop_Nmax = model.dual_RB_stop_Nmax;
end
if isfield(model, 'dual_RB_stop_timeout')
    model.RB_stop_timeout = model.dual_RB_stop_timeout;
end
detailed_data_dual = gen_detailed_data(model, model_data);
% write both primal and dual detailed_data in one combined detailed_data
detailed_data = [];
detailed_data = structcpy(detailed_data, detailed_data_primal);
detailed_data.dual_RB = detailed_data_dual.RB;
detailed_data.dual_RB_info = detailed_data_dual.RB_info;
detailed_data.dual_N = detailed_data_dual.N;
if isfield(detailed_data_dual, 'L_I_comp')
    detailed_data.dual_L_I_comp = detailed_data_dual.L_I_comp;
end
if isfield(detailed_data_dual, 'L_E_comp')
    detailed_data.dual_L_E_comp = detailed_data_dual.L_E_comp;
end
if isfield(detailed_data_dual, 'b_comp')
    detailed_data.dual_b_comp = detailed_data_dual.b_comp;
end