1 function rb_sim_data = stokes_rb_reconstruction(model, detailed_data, rb_sim_data)
2 %
function rb_sim_data = stokes_rb_reconstruction(model, detailed_data, rb_sim_data)
4 % detailed reconstruction
for stokes & lin_stat problems.
7 if ~isfield(rb_sim_data,
'uh')
8 rb_sim_data.uh =
Fem.DiscFunc([], detailed_data.df_info);
11 if get_rb_size(model, detailed_data) == 0
15 uh = detailed_data.RB(:,1:length(rb_sim_data.uN)) * rb_sim_data.uN;
17 if model.has_dirichlet_values
19 dir_comp = detailed_data.bc_info.dirichlet_dof_vector_components;
22 old_decomp_mode = model.decomp_mode;
23 model.decomp_mode = 2;
24 dir_coeff = model.dirichlet_values([], [], [], [], model);
25 model.decomp_mode = old_decomp_mode;
27 r_dir = lincomb_sequence(dir_comp, dir_coeff);
32 uh = uh + repmat(r_dir, 1, size(uh, 2));
35 rb_sim_data.uh.dofs = uh;