t_part_gen_reduced_data.m

function reduced_data = t_part_gen_reduced_data(model, detailed_data)
%function reduced_data = t_part_gen_reduced_data(model, detailed_data)
%
% function generating the reduced data for a t-part model.
%   To the normal reduced data of the base model a transfer matrix T is
%   generated, which projects the coefficients from the foregoing reduced
%   space to the next.
%
%   Required fields of detailed_data:
%       t_part_detailed_data:   the detailed data including the reduced
%                               bases from the t-partitions
%
%   Generated fields of reduced_data:
%       reduced_data_t_partitions:  Cell array containing the reduced data
%                                   for every t-partition.
%
%

% Markus Dihlmann 15.02.2010

if ~isfield(model,'t_part_in_basis_generation')
    model.t_part_in_basis_generation =0;
end

if ~isfield(model,'transition_model')
  model.transition_model =0;
end


if model.t_part_in_basis_generation
    %model is in basis generation mode
    %means that the detailed_data RB field has to be added temporarily as a
    %new field of tp_detailed_data
    if isfield(detailed_data,'t_part_detailed_data')
        %we are at least in the second domain in basis generation
        tp_detailed_data = detailed_data.t_part_detailed_data;
        len=length(tp_detailed_data);
        if isfield(model,'t_part_id_in_basisgen')
            %tp_detailed_data{model.t_part_id_in_basisgen} = [];
            %tp_detailed_data{model.t_part_id_in_basisgen} = structcpy(tp_detailed_data{model.t_part_id_in_basisgen}, detailed_data);
            if isfield(detailed_data,'RB')
                tp_detailed_data{model.t_part_id_in_basisgen}.RB = detailed_data.RB;
                tp_detailed_data{model.t_part_id_in_basisgen}.RB_info = detailed_data.RB_info;
            end
            tp_detailed_data{model.t_part_id_in_basisgen}.W = detailed_data.W;
            tp_detailed_data{model.t_part_id_in_basisgen}.grid = detailed_data.grid;
        else
            %append actual basis as "new t_partition_detailed_data"
            tp_detailed_data{len+1} = [];
            %tp_detailed_data{len+1} = structcpy(tp_detailed_data{len+1}, detailed_data);
            tp_detailed_data{len+1}.RB = detailed_data.RB;
            tp_detailed_data{len+1}.RB_info = detailed_data.RB_info;
            tp_detailed_data{len+1}.W = detailed_data.W;
            tp_detailed_data{len+1}.grid = detailed_data.grid;
        end
    else
        tp_detailed_data{1} = [];
        tp_detailed_data{1} = structcpy(tp_detailed_data{1}, detailed_data);
    end
    
else
    if ~isfield(detailed_data,'t_part_detailed_data')
        error('No t-partition model. Can not create reduced_data');
    else
        tp_detailed_data = detailed_data.t_part_detailed_data;
    end

end


reduced_data = [];

base_model = model.base_model;


nr_t_parts = length(tp_detailed_data);

if (isfield(model,'t_part_for_simulation'))&&~isempty(model.t_part_for_simulation)
    stop_part = model.t_part_for_simulation;
else
    stop_part = nr_t_parts;
end


t_part_reduced_data = cell(1,stop_part);
if model.transition_model
  if stop_part>1
    t_part_transition_reduced_data = cell(1,stop_part-1);
  end
end


%generate reduced_data for every partition without transition part
for t_part_indx = 1:stop_part
    
    %usual reduced_data
    t_part_reduced_data{t_part_indx} = gen_reduced_data(base_model, tp_detailed_data{t_part_indx});
    
    %Transfer matrix
    if t_part_indx>1
        T=tp_detailed_data{t_part_indx}.RB' * detailed_data.W * tp_detailed_data{t_part_indx-1}.RB;
        t_part_reduced_data{t_part_indx}.T = T;
    end
    
    if t_part_indx>1 &&model.transition_model
      %add the mixed explicit operator to reduced_data
      t_part_transition_reduced_data{t_part_indx-1}=t_part_reduced_data{t_part_indx};
      %decomp_old = model.decomp_mode;
      model.decomp_mode =1;
      [L_I, L_E, b] = model.operators_ptr(model, detailed_data);
      for i=1:length(L_E)
        t_part_transition_reduced_data{t_part_indx-1}.LL_E{i} = ...
            tp_detailed_data{t_part_indx}.RB'*detailed_data.W* ...
            L_E{i}*tp_detailed_data{t_part_indx-1}.RB;
      end
      %nr_comp = length(t_part_transition_reduced_data{t_part_indx-1}.LL_E);
      %  t_part_transition_reduced_data{t_part_indx-1}.LL_E{1}=T;
        
      %Fehlermatrizen?!?!?!?!?
    end
    
    
end

reduced_data.t_part_reduced_data = t_part_reduced_data;
if model.transition_model
    reduced_data.t_part_transition_reduced_data = t_part_transition_reduced_data;
end
if model.t_part_in_basis_generation
    reduced_data.N = t_part_reduced_data{stop_part}.N;
end