Default: delaunay_reduction = 0 (off)
performs a reduction of the local Delaunay tetrahedralization: If a boundary element triangle intersects with a local Delaunay
tetrahedralization, the corresponding tet is going to deleted. Advise: This option should only be activated when necessary
because it is a very expensive operation