Triangulation routine, part of pdetools.jl