elliptic_discrete_model.m

Go to the documentation of this file.

function local_model = elliptic_discrete_model(model)
% function local_model = elliptic_discrete_model(model)
% function creating a model with local functions out of a model
% with global functions. See detailed description for explanation.
% 
% All volume-data functions can be evaluated in many elements
% simultanously by e.g.
%
% local_model.source(grid,eindices,loc,params)
%
% where eindices is a n times 1 vector of element-indices and loc
% is a local coordinate in the reference triangle. the result is a
% n times 1 vector of resulting values of the source.
% and boundary-data functions can be evaluated by e.g.
%
% local_model.dirichlet_values(grid,eindices,face_index,lloc,params)
%
% where eindices is a n times 1 vector of element-indices and lloc
% is a 1d coordinate in [0,1]. the result is a n times 1 vector of
% resulting dirichlet-values on the edge given by face_index.
%
% to avoid conflicts in the "global" model, all data functions can
% still be evaluated with (glob,params), where glob are global
% coordinates.

% Immanuel Maier, 25.03.2011

local_model = model;

% convert all data functions of model:

if model.has_reaction
  local_model.reaction = @(varargin) ...
      discrete_volume_values(model.reaction,varargin{:});
end;

if model.has_advection
  local_model.velocity = @(varargin) ...
      discrete_volume_values(model.velocity,varargin{:});
end;

if model.has_diffusivity
  local_model.diffusivity_tensor = @(varargin) ...
      discrete_volume_values(model.diffusivity_tensor,varargin{:});
end;

if model.has_source
 local_model.source = @(varargin) ...
     discrete_volume_values(model.source,varargin{:});
end;

if model.has_dirichlet_values
  local_model.dirichlet_values = @(varargin) ...
      discrete_boundary_values(model.dirichlet_values,varargin{:});
end;

if model.has_neumann_values
  local_model.neumann_values = @(varargin) ...
      discrete_boundary_values(model.neumann_values,varargin{:});
end;

if model.has_robin_values
  local_model.robin_alpha = @(varargin) ...
      discrete_boundary_values(model.robin_alpha,varargin{:});
  local_model.robin_beta = @(varargin) ...
      discrete_boundary_values(model.robin_beta,varargin{:});
  local_model.robin_values = @(varargin) ...
      discrete_boundary_values(model.robin_values,varargin{:});
end

end