t_part_rb_reconstruction.m

Go to the documentation of this file.

function sim_data = t_part_rb_reconstruction(model, detailed_data, sim_data)
%function sim_data = t_part_rb_reconstruction(model, detailed_data,
%sim_data)
%
% function making an rb reconstruciton for t-partition models.
% In t-partition there is an overlap at the border of time-partitions. At
% the border the initial value of the next t-partition is discarded and the
% last timestep solution of the foregoing t-partition is kept.
%
% 

% Markus Dihlmann 15.02.2011

base_model = model.base_model;

nr_t_parts = length(sim_data.t_part_sim_data); %number of t-partitions

t_part_sim_data  = sim_data.t_part_sim_data;

%rb_reconstruction of the sim_data in every partition
for ind = 1:nr_t_parts
    t_part_sim_data{ind} = rb_reconstruction(base_model, detailed_data.t_part_detailed_data{ind},t_part_sim_data{ind});
end

%putting together global sim_data cutting off first step in each partition
%(except initial partition);
sim_data.U = t_part_sim_data{1}.U;

if nr_t_parts>1
    if ~model.transition_model
        for ind = 2:nr_t_parts
            len=size(t_part_sim_data{ind}.U,2);
            sim_data.U = [sim_data.U,t_part_sim_data{ind}.U(:,2:len)];
        end
    else
        for ind = 2:nr_t_parts
            sim_data.U = [sim_data.U,t_part_sim_data{ind}.U];
        end
    end
end


sim_data.t_part_sim_data = t_part_sim_data;