comsol_operators_output_default.m

function v = comsol_operators_output_default(model,model_data)
%function v = comsol_operators_output_default(model,model_data)
%
% This function returns the output defined in the COMSOL model as for
% example a line integral (in comsol_model.result.numerical. ...).
%
% PLEASE WRITE AN OWN OUTPUT FUNCTION FOR EVERY MODEL THAT SUPPORTS THE
% THE FIELD model.use_comsol=0/1;
% See for example the functions used for the ThermalBlock models!
%
%This function is for linear output functionals with no explicit mu-dependency 
%and therefor no affine decomposition!!!
%
% Oliver Zeeb, 2012.

switch model.decomp_mode
    case 0
        error('not yet implemented')
        % is this case needed???
        
    case 1 %model_data is detailed_data
        %store original solution:
        u_save = model_data.comsol_model.sol(model.comsol_tags.sol).getU;
        %calculate the output for every basis function:
        v_comp = zeros(size(model_data.RB,2),1);
        for k = 1:size(model_data.RB,2)
            model_data.comsol_model.sol(model.comsol_tags.sol).setU(model_data.RB(:,k));
            model_data.comsol_model.sol(model.comsol_tags.sol).createSolution;
            v_comp(k) = model_data.comsol_model.result.numerical(model.comsol_tags.output).getReal();
        end
        v{1} = v_comp;
        %reset to original solution:
        model_data.comsol_model.sol(model.comsol_tags.sol).setU(u_save);
        model_data.comsol_model.sol(model.comsol_tags.sol).updateSolution;
        
    case 2
        v(1) = 1;
end