rbmatlab  1.16.09
 All Classes Namespaces Files Functions Variables Modules Pages
dom_dec_compute_error.m
1 function sim_data = dom_dec_compute_error(model,model_data, ...
2  sim_data)
3 %function sim_data = dom_dec_compute_error(model,model_data,sim_data)
4 %
5 
6 % IM, 06.03.2012
7 
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);
11 end;
12 
13 X_err = zeros(sim_data.numiter,1);
14 
15 eh = cell(1,2);
16 for i = 1:2
17  eh{i} = femdiscfunc([],model_data.df_infos{i});
18 end;
19 
20 for j = 1:sim_data.numiter
21  for i = 1:2
22  eh{i}.dofs = sim_data.all_dofs{i}(:,j) - ...
23  sim_data.base_sim_data.uh.dofs(model_data.masks{i});
24 
25  X_err(j) = X_err(j) + fem_l2_norm(eh{i})^2 ...
26  + fem_h10_norm(eh{i})^2;
27  end;
28 
29  X_err(j) = sqrt(X_err(j));
30 end;
31 
32 sim_data.X_err = X_err;