vi_reduced_data_subset.m

function reduced_data_subset = vi_reduced_data_subset(model,reduced_data)
%function reduced_data_subset = vi_reduced_data_subset(model,reduced_data)
%
% function generating all parameter independent reduced data for
% the RB-Method for Parametrized Variational Inequalities.
%
% Generated fields of reduced_data:
%  AN_comp: components of matrix AN
%  fN_comp: components of vector fN
%  gN_comp: components of vector gN
%  BN: this matrix needs no decomposition

% B. Haasdonk 29.09.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 > length(reduced_data.AN_comp{1})
    error('N too large for current size of reduced basis!');
  end;
  N = model.N;
  if isfield(reduced_data, 'AN_comp')
    reduced_data_subset.AN_comp   = subblock_sequence(reduced_data.AN_comp,1:N,1:N);
  end
  if isfield(reduced_data, 'fN_comp')
    reduced_data_subset.fN_comp   = subblock_sequence(reduced_data.fN_comp,1:N);
  end
  if isfield(reduced_data, 'gN_comp')
%    reduced_data_subset.gN_comp   = subblock_sequence(reduced_data.gN_comp,1:N);
    reduced_data_subset.gN_comp   = reduced_data.gN_comp;
  end
  if isfield(reduced_data, 'BN')
    reduced_data_subset.BN    = reduced_data.BN(1:N,:);
  end
  
  reduced_data_subset.N = N;

end