t_part_model.m

function t_model = t_part_model(model, params)
%function t_model = t_part_model(model, params)
%
% model for t-partition inheriting infos from base_model
% In this model layer the repartition of the time domain is fixed.
% 
% The implementation of the evolution scheme must accept fields
% model.starting_time_step and model.stopping_time_step to enable
% simulations restrained to specigfic partitions.
%
% The field 'basis_vector_overlap' fixes how many
% POD-Vectors/overlap-vectors are used for the initial basis in the
% following partition.
% The field 't_part_rb_generation_cut_error' should be set to 1. It fixes
% that only the error produced in the actual partition is regarded during
% the basis generation process.
% 
%
% Optional fields of params:
%   transisition_model: if this flag is set, then the projection between
%                       the reduced spaces is replaced by a transition
%                       phase (including projeciton in evolution
%                       step)(automatic creation only implemented for first
%                       order time integration)
%   no_middle_point_bisection: if flag is set to one, the bisection is
%                               done, where half the error is reached.
%
%   t_part_rb_generation_mode: 'fixed' or 'adaptive_error_restriction'
%
% Generated fields of t_model:
%   t_part_range:   cell-array giving the first and the last time step index in
%                   the respective time domain
%   t_part_transition_range:    cell-array giving the first and the last
%                               time step index of the transisiotn range
%                               between t-domains. (In Euler-time
%                               integrations the length of this range is
%                               1.)
%   
%
% Markus Dihlmann 11.02.11
%

transition_model =0;

if nargin>1
    if isfield(params,'transition_model')
        if params.transition_model
            transition_model =1;
        end
    end
end

%copy basically all information from the base model in the t-model (for
%basis generation this info is neccessary.)
t_model = model;
   
t_model.t_part_middle_point_bisection=1;

if nargin ==1
    params=[];
end
    
if isfield(params,'no_middle_point_bisection')
    if params.no_middle_point_bisection
        t_model.t_part_middle_point_bisection=0;
    end
end
   
if isfield(params,'t_part_rb_generation_mode')
    t_model.t_part_rb_generation_mode = params.t_part_rb_generation_mode;
else
    t_model.t_part_rb_generation_mode='adaptive_error_restriction';
end

%copy basically all information from the base model in the t-model (for
%basis generation this info is neccessary.)
t_model = model;

t_model.base_model = model;
t_model.transition_model = transition_model;
%t_model.model_type = 't-part model';

%copy time information from model
%t_model.T = model.T;
%t_model.nt = model.nt;
%t_model.dt = model.dt;
%t_model.base_model.save_time_indices = 0:model.nt;
%t_model.mu_names = model.mu_names;
%t_model.mu_ranges = model.mu_ranges;

t_model.set_mu=@set_mu_in_model_and_base_model;

%Time domain information
%-----------------------
%
if (nargin==1)||(~isfield(params,'t_part_range'))
    %initial t_part_model with only one time domain
    t_model.t_part_range{1} = [0,t_model.nt];
else
    %taking repartition of time domain from params
    t_part_range = params.t_part_range;
    t_model.t_part_range = t_part_range;
    %fixing transisiton ranges (ranges between time domains)
    if ~isfield(params,'t_part_transition_range')&&(length(t_part_range)>1)
        if transition_model
            %t_model.rb_operators = @rb_operators_t_part;
            for i=2:(length(t_part_range))
                %t_model.t_part_range{i} = [t_part_range{i}(1)+1,t_part_range{i}(2)];
                t_model.t_part_transition_range{i-1} = [t_part_range{i-1}(2), t_part_range{i}(1)+1];
            end
        end
    else
        if transition_model
            t_model.t_part_transition_range=params.t_part_transition_range;
        end
    end
end
        



% main function pointers
t_model.gen_model_data = @t_part_gen_model_data;
t_model.detailed_simulation = @t_part_detailed_simulation;
t_model.gen_reduced_data = @t_part_gen_reduced_data;
t_model.rb_simulation = @t_part_rb_simulation;
t_model.rb_reconstruction = @t_part_rb_reconstruction;
t_model.plot_sim_data = model.plot_sim_data;
t_model.gen_detailed_data = @t_part_gen_detailed_data;


%RB_generation fields
%t_model.t_part_rb_generation_mode = 'fixed'; %'adaptive'
t_model.t_part_rb_generation_cut_error = 1; %if set, only the error produced in the actual partition is used in greedy algorithm
t_model.basis_vector_overlap = 5; %how many basis vectors are used from the foregoing t-part RB as initial basis in the next
t_model.reduced_data_subset = @t_part_reduced_data_subset;@(model, reduced_data)reduced_data; %evtl. replace?!?
t_model.get_estimator_from_sim_data = model.get_estimator_from_sim_data;
t_model.RB_stop_epsilon_global = model.RB_stop_epsilon;
t_model.t_part_middle_point_bisection=1;