comsol_lin_stat_gen_reduced_data.m

function reduced_data = comsol_lin_stat_gen_reduced_data(model,detailed_data)
%function reduced_data = comsol_lin_stat_gen_reduced_data(model,detailed_data)
%
% This function uses lin_stat_gen_reduced_data to generate the reduced data
% and then adds some comsol_specific fields to the reduced _data, the
% indices-vectors
% In order for the error estimators to work properly, it is needed, to cut
% off the comsol_DBC_data and only pass the eliminated data to
% lin_stat_gen_reduced_data! That is the main purpose of this function!
% 
%
% Oliver Zeeb, 19.09.2012



%Preparation: Cut the data in the comsol_matrices that belong to DBC-Dofs
% --> needed for correct error estimators!
if size(detailed_data.RB,1) == detailed_data.grid_info.ndofs %full vectors? --> eliminate DBC_data!
    detailed_data.RB=detailed_data.RB(detailed_data.ind_vectors.com_DOF_ind,:);% DBC eliminated from reduced Basis!
    in_prod_mat_string = char(model.get_inner_product_matrix);
    eval(['model.get_inner_product_matrix =',in_prod_mat_string,'_eliminated;']);
    model.get_dofs_from_sim_data = @(sim_data) sim_data.U(sim_data.ind_vectors.com_DOF_ind);
    model.comsol_get_eliminated_data = 1;
end

%use rbmatlab function to compute reduced_data
reduced_data = lin_stat_gen_reduced_data(model,detailed_data);

reduced_data.ind_vectors = detailed_data.ind_vectors;