Indices of points in canonical ordering for rectangular grid
[n1,n2,h1,h2,i,c,ii,cc] = poiindex(p,e,t,sd)
[n1,n2,h1,h2,i,c,ii,cc] = poiindex(p,e,t,sd) identifies
a given grid
sd as an evenly spaced rectangular
grid. If the grid is not rectangular,
n1 is 0 on
the number of points in the first and second directions,
ii are of
(n1-2)*(n2-2) and contain indices of interior
i contains indices of the original mesh,
ii contains indices of the canonical ordering.
n1*n2-(n1-2)*(n2-2) and contain indices
of border points.
In the canonical ordering, points are numbered from left to
right and then from bottom to top. Thus if
n1 = 3 and
= 5, then
ii = [5 8 11] and
= [1 2 3 4 6 7 9 10 12 13 14 15].