SeparableFunctionNode.m

classdef SeparableFunctionNode < Greedy.User.IReducedDataNode

  properties
    bb;

    sfhandle;
  end

  properties (SetAccess = private, Dependent)
    id;
    N;
  end

  methods
    function sfn = SeparableFunctionNode(rmodel, detailed_data, sfhandle)
      if nargin == 0
        error('Greedy.DataTree.Reduced.SeparableFunctionNode needs three constructor arguments');
      elseif nargin == 2 && isa(rmodel, 'IReducedModel') ...
          && isa(detailed_data, 'Greedy.DataTree.Reduced.SeparableFunctionNode')
        copy = detailed_data;
        copy_extract(sfn, copy, rmodel);
      elseif nargin == 3
        sfn.sfhandle = sfhandle;
        fill_fields(sfn, rmodel, detailed_data);
      else
        error('Wrong constructor arguments');
      end
    end

    function yesno = needs_subset_copy(this, rmodel)
      % function yesno = needs_subset_copy(this, rmodel)
      % @copybrief ::GreedyUser.IReducedDataNode.needs_subset_copy()
      %
      % @copydetails ::GreedyUser.IReducedDataNode.needs_subset_copy()
      %
      % Parameters:
      %   rmodel: of type TwoPhaseFlow.ReducedModel
      rn = get_by_description(rmodel.N, this.id);
      yesno = this.N ~= rn;

    end

    function N = get.N(this)
      N = size(this.bb,1);
    end

    function id = get.id(this)
      id = this.sfhandle.id;
    end

    function conds = get_conds(this)
      conds = [];
    end

    function lc = lincomb_sequence(this, descr)

      coeffs = this.sfhandle.coefficients(descr);
      lc = this.bb * coeffs;
    end
  end

  methods(Access=private)

    function fill_fields(this, rmodel, detailed_data)
      this.bb = this.sfhandle.project_components(rmodel, detailed_data);
    end

    function copy_extract(this, reduced_data, rmodel)
      this.sfhandle = reduced_data.sfhandle;
      rn = get_by_description(rmodel.N, this.id);

      if rn > this.N
        error('N too large for current size of reduced data!');
      end

      this.bb = reduced_data.bb(1:rn);
    end
  end
end