comsol_lin_evol_rb_simulation.m

function rb_sim_data = comsol_lin_evol_rb_simulation(model,reduced_data)
%function rb_sim_data = comsol_lin_evol_rb_simulation(model,reduced_data)
%
% This function does a reduced simulation using the rbmatlab function 
% lin_evol_rb_simulation. 
% Due to the way comsol assembles the matrices the coefficients calculated
% by rbmatlab must be "shifted" so thaht also the initial solution is
% correct. We cannot use 0 as an initial reduced basis, since that qould
% crash all the error estimators!
%
% Oliver Zeeb, 2013

rb_sim_data=lin_evol_rb_simulation(model,reduced_data);

% since COMSOL does a homogenization the initial values must be added to
% the solution calculated using rbmatlab.
% so the first "calculated" solution must be 0 so that after adding the
% initial values detailed and reduced simulation have the same values for
% t=0!
rb_sim_data.a(:,end) = [];
rb_sim_data.a=[zeros(size(rb_sim_data.a,1),1), rb_sim_data.a];