lin_evol_gen_detailed_data.m

function detailed_data=lin_evol_gen_detailed_data(model,model_data)
%function detailed_data=lin_evol_gen_detailed_data(model,model_data)
%
% method, which prepares reduced_data, which is meant as 
% data, that may be dependent on H only, and is not required during the 
% online-simulation, but it may be used for reconstruction purposes or
% computation of online_quantities.  i.e. reduced basis vectors, 
% colateral reduced basis spaces, grid, etc. can be stored here.
%
% So detailed_data is produced and used only in offline-phase algorithms.
%
% allowed dependency of data: H
% allowed dependency of computation: H
% Unknown at this stage: mu, Nmax, N
%
% Required fields of model:
%  gridtype : rectgrid or triagrid to be generated
%  mu_names : cell array of parameter-names
%  mu_ranges : cell array of parameter-intervals
%
% Further fields are required for rb generation,
% see rb_basis_generation
%
% generated fields of detailed_data:
%
%  RB    : reduced basis columns 1,...,Nmax
%  grid  : grid to be used in the subsequent stages
%  RB_info : depending on generation method some detailed information

% Bernard Haasdonk 15.5.2007
% edited by Dominik Garmatter 11.09 2012

if(isempty(model_data))
  model_data = gen_model_data(model);
end
detailed_data = [];

%% begin of edit by Dominik Garmatter 11.09 2012
% components are now computet and stored in the detailed_data.

if ~isfield(model, 'want_dual')
    model.want_dual = 0;
end
% computing the operators dependant on the controlfield 'model.want_dual'
% 1 = dual Operators are computet, 0 = primal ones are.
% 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.
model.decomp_mode = 1;
if model.want_dual % dual case
  [~, ~, ~, detailed_data.L_I_comp, detailed_data.L_E_comp] = model.operators_ptr(model, model_data);
  if strcmp(model.name_output_functional, 'theta') == 1
    detailed_data.b_comp = model.operators_output(model, model_data); % v_l the Riesz_representant of the theta-functional
  end
else % primal case
  [detailed_data.L_I_comp, detailed_data.L_E_comp, detailed_data.b_comp] = model.operators_ptr(model, model_data);
end
%% end of edit

skip_rb_generation = false;
detailed_data = structcpy(detailed_data, model_data);
if isfield(model_data, 'mexptr') && isfield(model, 'detailedfname')
  if model_data.mexptr('is_valid_rb')
    if exist(model.detailedfname, 'file')
      tmp = load(model.detailedfname);
      if length(tmp.detailed_data.RB) == model_data.mexptr('get_rb_size')
        disp('loading already generated rb space!');
        detailed_data = tmp.detailed_data;
        skip_rb_generation = true;
      else
        model.mexptr('reset_rb');
      end
    else
      model.mexptr('reset_rb');
    end
  end
end
if ~skip_rb_generation
  detailed_data = rb_basis_generation(model,detailed_data);
  if isfield(model, 'detailedfname')
    save(model.detailedfname, 'model', 'detailed_data');
  end
end