lin_stat_reduced_data_subset.m

Go to the documentation of this file.

function reduced_data_subset = lin_stat_reduced_data_subset(model,reduced_data)
%function reduced_data_subset = lin_evol_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
%
% Required fields of reduced_data:
%  N    : number of reduced basis vectors in the reduced_data,
%         must be larger than model.N!!

% Bernard Haasdonk 27.2.2011

reduced_data_subset = reduced_data;

if reduced_data.N ~= model.N
  % extract correct N-sized submatrices and subvectors from reduced_data
  if model.N > reduced_data.N
    error('N too large for current size of reduced basis!');
  end;
  N = model.N;
  if isfield(reduced_data, 'AN_comp')
    dummy = zeros(reduced_data.N);
    dummy(1:N,1:N)= 1;
    ind = find(dummy);
    reduced_data_subset.AN_comp   = reduced_data.AN_comp(ind,:);
  end
  if isfield(reduced_data, 'fN_comp')
    reduced_data_subset.fN_comp = reduced_data.fN_comp(1:N,:);
  end
  if isfield(reduced_data, 'lN_comp')
    reduced_data_subset.lN_comp = reduced_data.lN_comp(1:N,:);
  end

  Q_f = size(reduced_data_subset.fN_comp,2);
  Q_a = size(reduced_data_subset.AN_comp,2);
  
  ind = reshape(1:reduced_data.N*Q_a,reduced_data.N,Q_a)'; 
  ind = ind(:,1:N)';
  ind = [(1:Q_f)'; ind(:) + Q_f];  
  reduced_data_subset.G = reduced_data.G(ind,ind);
  
  reduced_data_subset.N = N;

end