ISeparableFunction.m

classdef ISeparableFunction < handle

  properties (Constant, Abstract)
    id;
  end

  methods (Abstract)

    res = evaluate(this, descr, model_data);

    comps = components(this, descr, model_data);

    coeffs = coefficients(this, descr);

    A = inner_product_matrix(this, descr, model_data);
  end

  methods
    function fa = projection(this, model, detailed_data)
      model_data = detailed_data.model_data;
      A = inner_product_matrix(this, model, model_data);
      f = evaluate(this, model, model_data);

      fa = f' * A * detailed_data.RB;
    end

    function acomp = project_components(this, rmodel, detailed_data)
      model_data = detailed_data.model_data;
      A = inner_product_matrix(this, rmodel, model_data);

      Nmax = get_N(rmodel, this.id);

      comps = components(this, rmodel.descr, model_data);

      acomp = cellfun(@(X) X' * A * detailed_data.RB(:,1:Nmax), ...
                      comps, 'UniformOutput', false);
      acomp = cat(1, acomp{:});

      if rmodel.debug && length(acomp) > 1
        if max(abs(detailed_data.RB(:,Nmax) * acomp{1}'-acomp{1}))>0.001
          U0 = [acomp{1},acomp{2}]; 
          U0_RB = [detailed_data.RB(:,Nmax)*acomp{1}',detailed_data.RB*a0{2}']; 
          plot_params = [];
          plot_params.plot = model.plot;
          plot_params.title = [this.id, 'Components'];
          plot_sequence(U0,model_data.grid,plot_params);
          plot_params.title = [this.id, 'RB'];
          plot_sequence(U0_RB,model_data.grid,plot_params);
          error('rb projection badly approximating separable function !!!');
        end
      end
    end
  end
end