fem_operators_output_volume_integral.m

function v = fem_operators_output_volume_integral(model,model_data)
%function v = fem_operators_output_volume_integral(model,model_data)
%
% function computing the output vectors of an elliptic problem with
% finite element discretization and output functional consisting of
% a weighted volume integral
%
%    `s(u) = \int_{\Omega} f(x) u(x)`
%
% required fields of model:
%    output_function:  weight function f

% B. Haasdonk 22.2.2011

disp('output to be implemented, currently dummy implementation!');

% dummy output: mean over computational domain
vec = model_data.df_info.l2_inner_product_matrix * ...
      ones(model.df_info.ndofs,1);
if model.decomp_mode == 0 % == complete
  v = vec;
elseif model.decomp_mode == 1 % == components
  v ={vec};
else % == 2 coefficients
  v = 1;
end;