lin_evol_opt_reduced_data_subset.m

Go to the documentation of this file.

function reduced_data_subset = lin_evol_opt_reduced_data_subset(model,reduced_data)
%function reduced_data_subset = lin_evol_opt_reduced_data_subset(model,reduced_data)
% method which modifies reduced_data, which is the data, that will
% be passed to the online-simulation algorithm.
%
% Typically, this routine only does a submatrix extraction of the 'Nmax' sized
% offline-objects to produce 'N' sized objects for the real simulation.
%
% Required fields of model:
%  N       : number of reduced basis vectors to choose
%   N_der  : number of reduced basis vectors in derivative solution
%
% Optional fields of model:
%  name_output_functional : name of an output functional
%
% Required fields of reduced_data:
%  N    : number of reduced basis vectors in the reduced_data.
%
% optional fields of reduced_data:
%  a0   : a cell array sequence of N-vector components of initial data
%         projection
%  LL_E : a cell array sequence of N x N component-Matrices of explicit
%         operator evaluations.
%  LL_I : a cell array sequence of N x N component-Matrices of implicit
%         operator evaluations.
%  bb   : a cell array sequence of N-vector components of the offset
%  bb_I : a cell array sequence of N-vector components of the implicit offset
%  ...
%...
%
% generated fields of reduced_data_subset:
%  a0   : a cell array sequence of N-vector components of initial data
%         projection
%  LL_E : a cell array sequence of N x N component-Matrices of explicit
%         operator evaluations.
%  LL_I : a cell array sequence of N x N component-Matrices of implicit
%         operator evaluations.
%  bb   : a cell array sequence of N-vector components of the offset
%  bb_I : a cell array sequence of N-vector components of the implicit offset
% .,...
%...
%s_RB     : in case of a given output functional, this is a vector of output
%             functional values of the reduced basis
%  s_l2norm : in case of a given output functional, this is the `L^2`-norm of
%             the output functional
%  N    : number of reduced basis vectors in the reduced_data.
%
% \note The fields for 'reduced_data_subset' are only generated if they exist in
% 'reduced_data'.

% Markus Dihlmann 04.07.2010

reduced_data_subset = reduced_data;

if reduced_data.N ~= model.N
  % extract correct N-sized submatrices and subvectors from reduced_data
  if model.N > length(reduced_data.a0{1})
    error('N too large for current size of reduced basis!');
  end;
  N = model.N;
  if isfield(reduced_data, 'a0')
    reduced_data_subset.a0   = subblock_sequence(reduced_data.a0,1:N);
  end
  if isfield(reduced_data, 'LL_I')
    reduced_data_subset.LL_I = subblock_sequence(reduced_data.LL_I,1:N,1:N);
  end
  if isfield(reduced_data, 'LL_E')
    reduced_data_subset.LL_E = subblock_sequence(reduced_data.LL_E,1:N,1:N);
  end
  if isfield(reduced_data, 'bb')
    reduced_data_subset.bb   = subblock_sequence(reduced_data.bb,1:N);
  end
  if isfield(reduced_data, 'bb_I')
    reduced_data_subset.bb_I = subblock_sequence(reduced_data.bb_I,1:N);
  end
  if isfield(reduced_data, 'M_E')
    reduced_data_subset.M_E = subblock_sequence(reduced_data.M_E,1:N,1:N);
  end
  if isfield(reduced_data, 'M_b')
    reduced_data_subset.M_b = subblock_sequence(reduced_data.M_b,1:N);
  end
  if isfield(reduced_data, 'M_EE')
    reduced_data_subset.M_EE = subblock_sequence(reduced_data.M_EE,1:N,1:N);
  end
  if isfield(reduced_data, 'M_Eb')
    reduced_data_subset.M_Eb  = subblock_sequence(reduced_data.M_Eb,1:N);
  end
% NICH BENÖTIGT, DA EH SKALAR?  
%   if isfield(reduced_data, 'M_bb')
%     reduced_data_subset.M_bb  = subblock_sequence(reduced_data.M_bb,1);
%   end
  if isfield(reduced_data, 'M_I')
    reduced_data_subset.M_I    = subblock_sequence(reduced_data.M_I,1:N,1:N);
  end
  if isfield(reduced_data, 'M_II')
    reduced_data_subset.M_II    = subblock_sequence(reduced_data.M_II,1:N,1:N);
  end
  if isfield(reduced_data, 'M_IE')
    reduced_data_subset.M_IE    = subblock_sequence(reduced_data.M_IE,1:N,1:N);
  end
  if isfield(reduced_data, 'M_Ib')
    reduced_data_subset.M_Ib    = subblock_sequence(reduced_data.M_Ib,1:N);
  end
  if isfield(reduced_data, 'M_EdEd')
    reduced_data_subset.M_EdEd    = subblock_sequence(reduced_data.M_EdEd,1:N,1:N);
  end
  if isfield(reduced_data, 'M_IdId')
    reduced_data_subset.M_IdId    = subblock_sequence(reduced_data.M_IdId,1:N,1:N);
  end
% NICHT BENÖTIGT, DA EH SKALAR?  
%   if isfield(reduced_data, 'M_bdbd')
%     reduced_data_subset.M_bdbd    = subblock_sequence(reduced_data.M_bdbd,1);
%   end
  if isfield(reduced_data, 'M_Ed')
    reduced_data_subset.M_Ed    = subblock_sequence(reduced_data.M_Ed,1:N,1:N);
  end
  if isfield(reduced_data, 'M_Id')
    reduced_data_subset.M_Id    = subblock_sequence(reduced_data.M_Id,1:N,1:N);
  end
  if isfield(reduced_data, 'M_bd')
    reduced_data_subset.M_bd    = subblock_sequence(reduced_data.M_bd,1:N);
  end
  if isfield(reduced_data, 'M_IEd')
    reduced_data_subset.M_IEd    = subblock_sequence(reduced_data.M_IEd,1:N,1:N);
  end
  if isfield(reduced_data, 'M_IId')
    reduced_data_subset.M_IId    = subblock_sequence(reduced_data.M_IId,1:N,1:N);
  end
  if isfield(reduced_data, 'M_Ibd')
    reduced_data_subset.M_Ibd    = subblock_sequence(reduced_data.M_Ibd,1:N);
  end
  if isfield(reduced_data, 'M_EEd')
    reduced_data_subset.M_EEd    = subblock_sequence(reduced_data.M_EEd,1:N,1:N);
  end
  if isfield(reduced_data, 'M_EId')
    reduced_data_subset.M_EId    = subblock_sequence(reduced_data.M_EId,1:N,1:N);
  end
  if isfield(reduced_data, 'M_Ebd')
    reduced_data_subset.M_Ebd    = subblock_sequence(reduced_data.M_Ebd,1:N);
  end
  if isfield(reduced_data, 'M_EdId')
    reduced_data_subset.M_EdId    = subblock_sequence(reduced_data.M_EdId,1:N,1:N);
  end
  if isfield(reduced_data, 'M_Edbd')
    reduced_data_subset.M_Edbd    = subblock_sequence(reduced_data.M_Edbd,1:N);
  end
  if isfield(reduced_data, 'M_Idbd')
    reduced_data_subset.M_Idbd    = subblock_sequence(reduced_data.M_Idbd,1:N);
  end
  

  if isfield(model,'name_output_functional')
    reduced_data_subset.s_RB = (reduced_data.s_RB(1:N));
    if isfield(reduced_data,'s_l2norm')
        reduced_data_subset.s_l2norm = reduced_data.s_l2norm;
    end
  end;

  reduced_data_subset.N = N;

end