lin_evol_opt_reduced_data_subset_separate_bases.m

Go to the documentation of this file.

function reduced_data_subset = lin_evol_opt_reduced_data_subset_separate_bases(model,reduced_data)
%function reduced_data_subset = lin_evol_opt_reduced_data_subset_separate_bases(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 18.11.2011

reduced_data_subset = reduced_data;

incl_derivative = 0;

if isfield(reduced_data,'L_E_dd')
    incl_derivative=1;%with derivative caclulation
end

if incl_derivative
    
    nr_der = length(reduced_data.c0);
    
    if model.N > length(reduced_data.a0{1})
        error('N too large for current size of reduced basis!');
    end;
    N = model.N;
    
    if isfield(model,'N_der')
        N_der = model.N_der;
    else
        N_der = N.*ones(1,nr_der);
    end
    for i=1:nr_der
        if N_der(i) >length(reduced_data.c0{i}{1})
                error('N_der too large for current size of derivative reduced basis')
        end
    end
    
    
    %normal fields
      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
      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

    %derivative fields
    
    if isfield(reduced_data,'c0')
        reduced_data_subset.c0 = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.c0{p} = subblock_sequence(reduced_data.c0{p},1:N_der(p));
        end
    end
    if isfield(reduced_data,'L_E_dd')
        reduced_data_subset.L_E_dd = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.L_E_dd{p} = subblock_sequence(reduced_data.L_E_dd{p},1:N_der(p),1:N_der(p));
        end
    end
    if isfield(reduced_data,'L_I_dd')
        reduced_data_subset.L_I_dd = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.L_I_dd{p} = subblock_sequence(reduced_data.L_I_dd{p},1:N_der(p),1:N_der(p));
        end
    end
    if isfield(reduced_data,'dL_E_sd')
        reduced_data_subset.dL_E_sd = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.dL_E_sd{p} = subblock_sequence(reduced_data.dL_E_sd{p},1:N_der(p),1:N);
        end
    end
    if isfield(reduced_data,'dL_I_sd')
        reduced_data_subset.dL_I_sd = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.dL_I_sd{p} = subblock_sequence(reduced_data.dL_I_sd{p},1:N_der(p),1:N);
        end
    end
    if isfield(reduced_data,'db')
        reduced_data_subset.db = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.db{p} = subblock_sequence(reduced_data.db{p},1:N_der(p));
        end
    end
    if isfield(reduced_data,'K_E')
        reduced_data_subset.K_E = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_E{p} = subblock_sequence(reduced_data.K_E{p},1:N_der(p),1:N_der(p));
        end
    end
    if isfield(reduced_data,'K_I')
        reduced_data_subset.K_I = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_I{p} = subblock_sequence(reduced_data.K_I{p},1:N_der(p),1:N_der(p));
        end
    end
    if isfield(reduced_data,'K_EE')
        reduced_data_subset.K_EE = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_EE{p} = subblock_sequence(reduced_data.K_EE{p},1:N_der(p),1:N_der(p));
        end
    end
    if isfield(reduced_data,'K_II')
        reduced_data_subset.K_II = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_II{p} = subblock_sequence(reduced_data.K_II{p},1:N_der(p),1:N_der(p));
        end
    end
    if isfield(reduced_data,'K_IE')
        reduced_data_subset.K_IE = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_IE{p} = subblock_sequence(reduced_data.K_IE{p},1:N_der(p),1:N_der(p));
        end
    end
    if isfield(reduced_data,'K_EdEd')
        reduced_data_subset.K_EdEd = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_EdEd{p} = subblock_sequence(reduced_data.K_EdEd{p},1:N,1:N);
        end
    end
    if isfield(reduced_data,'K_IdId')
        reduced_data_subset.K_IdId = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_IdId{p} = subblock_sequence(reduced_data.K_IdId{p},1:N,1:N);
        end
    end
    if isfield(reduced_data,'K_Ed')
        reduced_data_subset.K_Ed = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_Ed{p} = subblock_sequence(reduced_data.K_Ed{p},1:N,1:N_der(p));
        end
    end
    if isfield(reduced_data,'K_Id')
        reduced_data_subset.K_Id = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_Id{p} = subblock_sequence(reduced_data.K_Id{p},1:N_der(p),1:N);
        end
    end
    if isfield(reduced_data,'K_bd')
        reduced_data_subset.K_bd = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_bd{p} = subblock_sequence(reduced_data.K_bd{p},1:N_der(p));
        end
    end
    if isfield(reduced_data,'K_IEd')
        reduced_data_subset.K_IEd = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_IEd{p} = subblock_sequence(reduced_data.K_IEd{p},1:N_der(p),1:N);
        end
    end
    if isfield(reduced_data,'K_IId')
        reduced_data_subset.K_IId = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_IId{p} = subblock_sequence(reduced_data.K_IId{p},1:N_der(p),1:N);
        end
    end
    if isfield(reduced_data,'K_Ibd')
        reduced_data_subset.K_Ibd = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_Ibd{p} = subblock_sequence(reduced_data.K_Ibd{p},1:N_der(p));
        end
    end
    if isfield(reduced_data,'K_EEd')
        reduced_data_subset.K_EEd = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_EEd{p} = subblock_sequence(reduced_data.K_EEd{p},1:N_der(p),1:N);
        end
    end
    if isfield(reduced_data,'K_EId')
        reduced_data_subset.K_EId = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_EId{p} = subblock_sequence(reduced_data.K_EId{p},1:N_der(p),1:N);
        end
    end
    if isfield(reduced_data,'K_Ebd')
        reduced_data_subset.K_Ebd = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_Ebd{p} = subblock_sequence(reduced_data.K_Ebd{p},1:N_der(p));
        end
    end
    if isfield(reduced_data,'K_EdId')
        reduced_data_subset.K_EdId = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_EdId{p} = subblock_sequence(reduced_data.K_EdId{p},1:N,1:N);
        end
    end
    if isfield(reduced_data,'K_Edbd')
        reduced_data_subset.K_Edbd = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_Edbd{p} = subblock_sequence(reduced_data.K_Edbd{p},1:N);
        end
    end
    if isfield(reduced_data,'K_Idbd')
        reduced_data_subset.K_Idbd= cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.K_Idbd{p} = subblock_sequence(reduced_data.K_Idbd{p},1:N);
        end
    end
    
    
else  %Extraction without derivative calculations
  % 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
  
  
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
    if incl_derivative
        reduced_data_subset.s_RB_der = cell(1,nr_der);
        for p=1:nr_der
            reduced_data_subset.s_RB_der{p} = reduced_data.s_RB_der{p}(1:N_der(p));
        end
    end
    
  end;

  reduced_data_subset.N = N;
  if incl_derivative
    reduced_data_subset.N_der = N_der;
  end