rb_reconstruction_derivative.m

function rb_sim_data = rb_reconstruction_derivative(model, detailed_data, rb_sim_data)
%function rb_sim_data = rb_reconstruction_default(model_data, detailed_data, rb_sim_data)
%
% function computing a detailed reconstruction by linear
% combination of the coefficients in the simulation data with the
% orthonormal reduced basis RB
% The sensitivity derivative solution is also reconstructed. Either using
% the common reduced basis or separte derivative information reduced bases.

% Markus Dihlmann 24.06.2010

rb_sim_data.U = detailed_data.RB(:,1:size(rb_sim_data.a,1)) * rb_sim_data.a;

if model.compute_derivative_info
    if ~isfield(detailed_data,'RB_der')
        for i = 1:size(rb_sim_data.c,3)
            rb_sim_data.U_der(:,:,i) = detailed_data.RB(:,1:size(rb_sim_data.a,1)) * rb_sim_data.c(:,:,i);
        end
    else
        indx_mu=find(model.optimization.params_to_optimize);
        for i = 1:sum(model.optimization.params_to_optimize)
            rb_sim_data.U_der(:,:,i) = detailed_data.RB_der{indx_mu(i)}(:,1:size(rb_sim_data.c{i},1)) * rb_sim_data.c{i}(:,:);
        end
    end
end
%calculate output
model.decomp_mode=0;
v = model.operators_output(model,detailed_data);
  rb_sim_data.output = (v(:)') * rb_sim_data.U;
  if model.compute_derivative_info
  for i=1:size(rb_sim_data.U_der,3)
    rb_sim_data.output_derivative(:,i) = (v(:)') * rb_sim_data.U_der(:,:,i);
  end
  end

%| \docupdate