1. A method for generating a constrained Delaunay triangulation for a planar domain with boundaries and internal features, comprising: approximating the boundaries and internal features of the domain with polylines; constructing unconstrained Delaunay triangulation for the domain; modifying the unconstrained Delaunay triangulation to conform triangle sides to the polylines; correcting the modified constrained triangulation to make it a constrained Delaunay triangulation.