Following work of Heindl and of Powell and Sabin, each triangle of an arbitrary (regular) triangulation $\Delta$ of a polygonal region $\Omega$ in $\Bbb R^2$ is subdivided into twelve triangles, using ...