comsol_ThermalBlock3D_model.m

function model = comsol_ThermalBlock3D_model(varargin)
%function model = comsol_ThermalBlock3D_model(varargin)
%
% This function defines the Comsol-3D-ThermalBlock-Model and adds the
% necessary fields to the model!
% 
% Oliver Zeeb, 2012



%%%%%%%%%%%%%%%%%%%%%%%%%%
% CHOOSE THE MODEL HERE!!!
%%%%%%%%%%%%%%%%%%%%%%%%%%

if isempty(varargin)
    model=gen_rbmatlab_model_from_comsol_model('ThermalBlock_3D_5x5x5_LinLagr_Einfluss_unten_3Params');
    %model=gen_rbmatlab_model_from_comsol_model('ThermalBlock_3D_5x5x5_LinLagr_Einfluss_unten_8Params');
    %model=gen_rbmatlab_model_from_comsol_model('ThermalBlock_3D_5x5x5_LinLagr_Einfluss_unten_125Params');
else
    model=gen_rbmatlab_model_from_comsol_model(varargin{1});
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% RB_GENERATION_MODE
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
model.RB_generation_mode = 'greedy_uniform_fixed';
model.RB_numintervals = 3 * ones(size(model.RB_numintervals));



%%%%%%%%%%%%%%%%%%%
%BASIS GENERATION
%%%%%%%%%%%%%%%%%%%

%set initial mu and mu_ranges
length_mu = length(get_mu(model));
model.mu_ranges = cell(length_mu,1);
mu_min = 1;
mu_max = 100;
for k=1:length_mu
    model.mu_ranges{k} = [mu_min, mu_max];
end
model=set_mu(model, ones(1,length_mu)*(mu_min+mu_max)/2);

model.get_inner_product_matrix = @(detailed_data)detailed_data.inner_product_matrices.L2;
model.RB_stop_epsilon = 1e-12;
model.RB_stop_Nmax = 20;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% MODEL OPERATORS
%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% set model operator
model.operators = @comsol_ThermalBlock3D_operators;
model.comsol_calculate_inner_product_matrices = @comsol_ThermalBlock3D_calculate_inner_product_matrices;
model.operators_output = @comsol_ThermalBlock3D_operators_output;
model.compute_output_functional = 1;

%error estimation quantities
model.coercivity_alpha = @(model) min(get_mu(model));
model.coercivity_alpha_LB = min(min(cell2mat(model.mu_ranges)));