get_triangle_midpoints.m

Go to the documentation of this file.

function midpoints = get_triangle_midpoints(p,t)
%computes triangle midpoints for given mesh data

ntria = size(t,2);              %nof triangles

tp=t(1:3,:);                    %indices of all triangles
tp(:);

%x
px=p(1,tp);                     %corresponding x values

px=reshape(px,3,ntria);         %corr. x values sorted

cx=mean(px)';                   

%y
py=p(2,tp);

py=reshape(py,3,ntria);

cy=mean(py)';

midpoints=[cx,cy];