dom_dec_compute_error.m

Go to the documentation of this file.

function sim_data = dom_dec_compute_error(model,model_data, ...
                                          sim_data)
% function sim_data =
% dom_dec_compute_error(model,model_data,sim_data)
%

% IM, 06.03.2012

if ~isfield(sim_data, 'base_sim_data')
  sim_data.base_sim_data = detailed_simulation(model.base_model, ...
                                    model_data.base_model_data);
end;

X_err = zeros(sim_data.numiter,1);

eh = cell(1,2);
for i = 1:2
  eh{i} = femdiscfunc([],model_data.df_infos{i});
end;

for j = 1:sim_data.numiter
  for i = 1:2
    eh{i}.dofs = sim_data.all_dofs{i}(:,j) - ...
        sim_data.base_sim_data.uh.dofs(model_data.masks{i});

    X_err(j) = X_err(j) + fem_l2_norm(eh{i})^2 ...
        + fem_h10_norm(eh{i})^2;
  end;
  
  X_err(j) = sqrt(X_err(j));
end;  

sim_data.X_err = X_err;