LocalizedOperatorDefault.m

classdef LocalizedOperatorDefault < ILocalizedOperator

  properties (Constant)
    id = 'L_localized';
  end

  properties (Access = private)
    fptr;
  end

  methods
    function lod = LocalizedOperatorDefault(fptr)
      lod.fptr = fptr;
    end

    function [INC, INCoff] = apply(this, model, model_data, arg, NU_ind)
      [INC, INCoff] = this.fptr(model, model_data, arg, NU_ind);
    end

    function ipm = inner_product_matrix(this, model_data)
      % function ipm = inner_product_matrix(this, model_data)

      ipm = model_data.W;
    end
  end

  methods(Static)

    function ret_size = ret_size(model_data, NU_ind)
      if nargin < 2 || isempty(NU_ind)
        NU_ind = 1:model_data.grid.nelements;
      end
      ret_size = length(NU_ind);
    end

    function arg_size = arg_size(model_data, NU_ind)
      if nargin < 2 || isempty(NU_ind)
        NU_ind = 1:model_data.grid.nelements;
      end
      arg_size = length(NU_ind);
    end

    function eind = compute_TM_global(model_data, TM, rmodel)
      if isequal(rmodel.stencil_mode, 'vertex')
        eind = compute_TM_global_vertex(model_data, TM);
      else
        eind = compute_TM_global_edge(model_data, TM, rmodel.local_stencil_size);
      end
    end

  end
end