1 function sim_data = dom_dec_compute_error(model,model_data, ...
3 %
function sim_data = dom_dec_compute_error(model,model_data,sim_data)
8 if ~isfield(sim_data,
'base_sim_data')
9 sim_data.base_sim_data = detailed_simulation(model.base_model, ...
10 model_data.base_model_data);
13 X_err = zeros(sim_data.numiter,1);
17 eh{i} = femdiscfunc([],model_data.df_infos{i});
20 for j = 1:sim_data.numiter
22 eh{i}.dofs = sim_data.all_dofs{i}(:,j) - ...
23 sim_data.base_sim_data.uh.dofs(model_data.masks{i});
25 X_err(j) = X_err(j) + fem_l2_norm(eh{i})^2 ...
26 + fem_h10_norm(eh{i})^2;
29 X_err(j) = sqrt(X_err(j));
32 sim_data.X_err = X_err;