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

if model.decomp_mode == 2
    
    v = model.output_function([],[],[],model);
else
    
    dirichlet_gids = model_data.df_info.dirichlet_gids;
    model_data.df_info.dirichlet_gids = [];
    
    v = fem_rhs_volume_part_assembly(model.output_function,model,model_data.df_info);
    
    model_data.df_info.dirichlet_gids = dirichlet_gids;
end

% 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;