1
function index = meshindex(p, t)
3
% MESHINDEX - COMPUTE APPROXIMATE MULTI-ADAPTIVE EFFICIENCY INDEX FOR MESH
5
% Usage: index = meshindex(p, t)
7
% p - points (exported from PDE Toolbox)
8
% t - triangles (exported from PDE Toolbox)
10
% Copyright (C) 2005 Anders Logg.
11
% Licensed under the GNU LGPL Version 2.1.
14
hlist = hmin*ones(size(p, 2), 1);
30
h = min([h0, h1, h2]);
32
hlist(i0) = min(hlist(i0), h);
33
hlist(i1) = min(hlist(i1), h);
34
hlist(i2) = min(hlist(i2), h);
40
index = (length(hlist) / hmin) / sum(1.0 ./ hlist);