ReducedModel.m

classdef ReducedModel < IReducedModel
  % reduced model for linear stationary problems as given by a
  % LinStat.DetailedModel.

  properties
    enable_error_estimator = false;
  end

  methods
    function rm = ReducedModel(dmodel, bg_descr)
      % function rm = ReducedModel(detailed_model, basis_generator)
      % Constructor for the reduced model.
      %
      % Parameters:
      %   dmodel: of type LinStat.DetailedModel

      % this implements a copy constructor if necessary...
      rm = rm@IReducedModel(dmodel, bg_descr);

    end

    %
    %
    % Parameters:
    %   reduced_data: of type LinStat.ReducedData
    rb_sim_data = rb_simulation(this, reduced_data);

    function rb_sim_data = rb_reconstruction(this, detailed_data, rb_sim_data)
    % function rb_sim_data = rb_reconstruction(this, detailed_data, rb_sim_data)
    % @copybrief IReducedModel.rb_reconstruction()
    %
    % Parameters:
    %   detailed_data: of type LinStat.DetailedData
    %   rb_sim_data:   struct holding reduced simulation data returned by IReducedModelrb_simulation().
    %
    % Return values:
    %   rb_sim_data: struct holding the reduced simulation results and their
    %                reconstructions.
    %
    % (trivial) function computing a detailed reconstruction by linear
    % combination of the coefficients in the simulation data with the
    % orthonormal reduced basis RB

    % Bernard Haasdonk 22.2.2011

      if ~isfield(rb_sim_data,'uh')
        rb_sim_data.uh = femdiscfunc([],detailed_data.model_data.df_info);
      end;
      rb_sim_data.uh.dofs = ...
        detailed_data.RB(:,1:length(rb_sim_data.uN)) * rb_sim_data.uN;
    end

    function c = copy(this)
      % function c = copy(this)
      % @copybrief IReducedModel.copy()
      %
      % Return values:
      %   c:   an object of type LinStat.ReducedModel which is a deep copy of
      %        this object.
      c = LinStat.ReducedModel(this);
    end

  end
  methods (Static)
    function Delta        = get_estimators_from_sim_data(rb_sim_data)
      % function Delta        = get_estimators_from_sim_data(sim_data)
      % @copybrief IReducedModel.get_estimators_from_sim_data()
      %
      % @copydetails IReducedModel.get_estimators_from_sim_data()
      Delta = rb_sim_data.Delta;
    end

    function Delta        = get_estimator_from_sim_data(rb_sim_data)
      % function Delta        = get_estimator_from_sim_data(sim_data)
      % @copybrief IReducedModel.get_estimator_from_sim_data()
      %
      % @copydetails IReducedModel.get_estimator_from_sim_data()
      Delta = rb_sim_data.Delta(end);
    end


  end
end