POINT centroid(const POINT& p1, const POINT& p2, const POINT& p3) { NT x = p1.xcoord() + p2.xcoord() + p3.xcoord(); NT y = p1.ycoord() + p2.ycoord() + p3.ycoord(); NT z = p1.zcoord() + p2.zcoord() + p3.zcoord(); return POINT(x,y,z); }
POINT centroid(const list<POINT>& L) { POINT pi; NT x(0), y(0), z(0); forall(pi,L){ x = x + pi.xcoord(); y = y + pi.ycoord(); z = z + pi.zcoord(); }