getmesh.m

来自「五点差分型多重网格方法：各种插值算子的比较）」· M 代码 · 共 40 行

function mesh = getmesh(node,elem,Dirichlet,Neumann,type,solu)
% GETMESH generates initial mesh data structure.
% 
% USAGE
%    mesh = getmesh(node,elem,Dirichlet,Neumann,type,solu)
%
% INPUT 
%       node:  array of nodes (required)
%       elem:  array of element indices (required)
%  Dirichlet:  Dirichlet boundary edges (required)
%    Neumann:  Neumann boundary edges (required)
%       type:  node in use (optional)
%       solu:  initial solution (optional)
%
% OUTPUT
%    mesh:  current mesh
%

% L. Chen & C. Zhang 10-15-2006

%--------------------------------------------------------------------------
% Generate opitional data
%--------------------------------------------------------------------------
if nargin <= 5, solu = zeros(size(node,1),1); end
if nargin <= 4, type = uint8(ones(size(node,1),1)); end

%--------------------------------------------------------------------------
% Label the mesh by the longest edge rule
%--------------------------------------------------------------------------
elem = label(node,elem);

%--------------------------------------------------------------------------
% Generate mesh data structure
%--------------------------------------------------------------------------
mesh = struct('node',node, 'elem',elem, 'type',type, 'solu',solu, ...
              'Dirichlet',Dirichlet, 'Neumann',Neumann);
      
%--------------------------------------------------------------------------
% End of function GETMESH
%--------------------------------------------------------------------------