comsol_ThermalBlock3D_operators_output.m

function output = comsol_ThermalBlock3D_operators_output(model,model_data)
%function output = comsol_ThermalBlock3D_operators_output(model,model_data)
%
% This function returns the output data (components and coefficients) for
% the comsol 3D Thermal Block model.
% Since this is a compliant case model, the output functional is equal to
% the rhs (this has to be scaled with the flux vector of the bottom 
% boundary, in case it is not = 1 !
%
% Oliver Zeeb, 08/2013



%%% LINEAR COMPLIANT CASE --> output = f;
if model.use_comsol
    switch model.decomp_mode
        case 0
            disp(' in comsol_ThermalBlock3D_operators_output: decomp_mode 0 not yet implemented!')
            keyboard
            
        case 1
            flux_factor = str2double(model_data.comsol_model.physics('c').feature('flux1').getString('g'));
            if model.comsol_get_eliminated_data %cut the matrices --> return eliminated matrices!
                output = {model_data.operators.f_comp_eliminated{1} ./ flux_factor};
            else
                output = {model_data.operators.f_comp_full{1} ./ flux_factor};
            end
            
        case 2
            output = 1;
    end
    
else
    switch model.decomp_mode
        case 0
            disp(' in comsol_ThermalBlock3D_operators_output: decomp_mode 0 not yet implemented!')
            keyboard
            
        case 1
            if model.comsol_get_eliminated_data %cut the matrices --> return eliminated matrices!
                output = model_data.operators.output_comp_eliminated;
            else
                output = model_data.operators.output_comp_full;
            end
            
        case 2
            output = 1;
    end
end