DetailedModel.m

classdef DetailedModel < Greedy.User.FVDetailedModelDefault
  % a test model doing nothing
  %
  % In this detailed model only the name of the called member function is
  % printed. Such it enables to set-up Greedy algorithms and check wether the
  % functions are called in the correct order.

  methods
    function dmei = DetailedModel(descr)
      % function dm = DetailedModel(descr)
      % constructor based on problem description
      %
      % Parameters:
      %   descr: structure of type ModelDescr describing the problem and the
      %          discretization
      dmei = dmei@Greedy.User.FVDetailedModelDefault(descr);
    end
  end

  methods
    function model_data   = gen_model_data(this)
      disp('calling model_data');
      gp.xrange = [0 1];
      gp.yrange = [0 1];
      if isfield(this.descr, 'gsize')
        gsize = this.descr.gsize;
      else
        gsize = [2 2];
      end
      gp.xnumintervals = gsize(1);
      gp.ynumintervals = gsize(2);
      model_data.grid = rectgrid(gp);
      model_data.W    = fv_mass_matrix([], model_data.grid, []);
    end

    function sim_data     = detailed_simulation(this, model_data)
      % function sim_data     = detailed_simulation(this, model_data)
      % @copybrief IDetailedModeldetailed_simulation()
      %
      % @copydetails IDetailedModeldetailed_simulation()
      disp('calling detailed_simulation');
      g = model_data.grid;
      funmode = false;
      miny = min(g.CY);
      xvals = g.CX(g.CY==miny);
      if isfield(this.descr, 'xfun')
        xdata = this.descr.xfun(this.descr, xvals);
        funmode = true;
      else
        xdata = ones(length(xvals), 1);
      end
      minx = min(g.CX);
      yvals = g.CY(g.CX==minx);
      if isfield(this.descr, 'yfun')
        ydata = this.descr.yfun(this.descr, yvals);
        funmode = true;
      else
        ydata = ones(length(yvals), 1);
      end

      if funmode
        sim_data.U = kron(ydata, xdata);
      else
        sim_data.U = sum(get_mu(this)) * ones(model_data.grid.nelements,this.descr.nt);
      end
    end

    function p = plot_sim_data(this, model_data, sim_data, plot_params)
      % function sim_data     = detailed_simulation(this, model_data[, plot_params])
      % @copybrief IDetailedModelplot_sim_data()
      %
      % @copydetails IDetailedModelplot_sim_data()

      if nargin <= 3
        plot_params = [];
      end
      p = fv_plot(model_data.grid, this.get_dofs_from_sim_data(sim_data), plot_params);
    end
  end

  methods(Static)
    function U            = get_dofs_from_sim_data(sim_data)
      disp('calling get_dofs_from_sim_data');
      U = sim_data.U;
    end
    
    function snapshot = get_dofs_at_time(sim_data, time_index)
      snapshot = sim_data(:,time_index);
    end
  end
end