Right now the only way to fix user-specified internal edges in the resulting triangulation is to create multiple meshes and join them as explained in the README.
However, a more flexible interface would allow the user to specify a list of additional edges to include inside the domain. This requires changes to
the mesh smoothing logic (some nodes must stay fixed inside the domain; right now only boundary nodes are fixed)
the adaptive refinement logic (we must be able to identify possible new nodes popping up between the existing fixed nodes)
The latter would be fixed by the first approach described in #2.
Right now the only way to fix user-specified internal edges in the resulting triangulation is to create multiple meshes and join them as explained in the README.
However, a more flexible interface would allow the user to specify a list of additional edges to include inside the domain. This requires changes to
The latter would be fixed by the first approach described in #2.