As discussed in this StackExchange question, there is a function available in SFCGAL that performs a Constrained Delaunay Triangulation. Some explanation and examples are available here.
There is a secret function available in SFCGAL that just needs to be exposed to use it.
Make sure to set the correct PostGIS version in '$libdir/postgis-2.4
.