gen_reduced_model.m

Go to the documentation of this file.

function rmodel = gen_reduced_model(dmodel, bg_descr)
% function rmodel = gen_reduced_model(dmodel, bg_descr)
% generates an IReducedModel instance from a description structure.
%
% This function constructs the correct IReducedModel instance for a BasisGenDescr
% structure. The field 'rb_problem_type' should contain the package name of the
% 'ReducedModel' class. If e.g. 'descr.rb_problem_type' equals the string
% 'LinEvol', an object of type LinEvol.ReducedModel is returned.
%
% The class type of the returned object can also be chosen explictly by
% specifying the field 'descr.rmodel_constructor'.
%
% Parameters:
%       bg_descr: a Matlab structure of type BasisGenDescr.

if nargin == 1 || isempty(bg_descr);
  descr = dmodel.descr;
  fns    = fieldnames(descr);
  rbfnsI = cellfun(@(X) length(X)>3 && isequal(X(1:3), 'RB_'), fns);
  rbfns  = fns(rbfnsI);
  for i = 1:length(rbfns)
    rbfn = rbfns{i};
    bg_descr.(rbfn(4:end)) = descr.(rbfn);
    descr = rmfield(descr, rbfn);
  end
  bg_descr.rb_problem_type = descr.rb_problem_type;
end

if ~isfield(bg_descr, 'rb_problem_type')
  error('BasisGenDescr does not have the required field ''rb_problem_type''');
end

if ~isfield(bg_descr, 'rmodel_constructor') || isempty(bg_descr.rmodel_constructor)
  bg_descr.rmodel_constructor = eval(['@', bg_descr.rb_problem_type, '.ReducedModel']);
end

if ~isfield(bg_descr, 'detailed_data_constructor') || isempty(bg_descr.detailed_data_constructor)
  bg_descr.detailed_data_constructor = eval(['@', bg_descr.rb_problem_type, '.DetailedData']);
end

if ~isfield(bg_descr, 'reduced_data_constructor') || isempty(bg_descr.reduced_data_constructor)
  bg_descr.reduced_data_constructor = eval(['@', bg_descr.rb_problem_type, '.ReducedData']);
end

rmodel = bg_descr.rmodel_constructor(dmodel, bg_descr);

end