gradmat.m

来自「解决麦克斯韦的matlab源文件」· M 代码 · 共 16 行

function G = Gradmat(mesh)
% Computes the topological gradient matrix of the mesh
%
% mesh -> Data structure for 2D unstructured mesh
% 
% the gradient matrix is the incidence matrix of oriented edges and
% vertices. 
% 

G = sparse(mesh.Ne,mesh.Nv);

for i=1:mesh.Ne
  G(i,mesh.ep(i,1)) = -1;
  G(i,mesh.ep(i,2)) = 1;
end