refine.m

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

function refmesh = refine(mesh)
% Regular refinement of unstructured triangular mesh
%
% mesh -> data structure for unstructured 2D triangulation
%

Nv_ref = mesh.Nv + mesh.Ne;
Nt_ref = 4*mesh.Nt;

x = zeros(Nv_ref,1);
y = zeros(Nv_ref,1);
tri = [];

x(1:mesh.Nv,1) = mesh.vt(:,1);
y(1:mesh.Nv,1) = mesh.vt(:,2);

for i=1:mesh.Ne
  mp = sum(mesh.vt(mesh.ep(i,:),:))/2;
  x(i+mesh.Nv,1) = mp(1);
  y(i+mesh.Nv,1) = mp(2);
end

for i=1:mesh.Nt
  v = mesh.trv(i,:);
  e = mesh.tre(i,:) + mesh.Nv;
  tri = [tri; v(1) e(3) e(2); e(3) v(2) e(1); e(2) e(1) v(3); e(3) e(1) e(2)];
end

refmesh = convert(x,y,tri);