[0] TetGenMeshDelaunizeVertices(): Constructing Delaunay tetrahedralization.
[0] TetGenMeshDelaunizeVertices():   Sorting vertices by a bsp-tree.
[0] TetGenMeshBTreeSort():   Depth 0, 8 verts.
[0] TetGenMeshBTreeSort():   Bbox (0., 0., 0.),(1., 1., 1.). x-axis
[0] TetGenMeshBTreeSort():     leftsize = 4, rightsize = 4
[0] TetGenMeshDelaunizeVertices():   Number of tree nodes: 2.
[0] TetGenMeshDelaunizeVertices():   Maximum tree node size: 4.
[0] TetGenMeshDelaunizeVertices():   Maximum tree depth: 0.
[0] TetGenMeshDelaunizeVertices():   Incrementally inserting vertices.
[0] TetGenMeshDelaunayIncrFlip():     Create the first tet (0, 3, 5, 1).
[0] TetGenMeshDelaunayIncrFlip():     Incrementally inserting points.
[0] TetGenMeshInsertVertexBW():     Insert point 7
[0] TetGenMeshInsertVertexBW():     Walk distance (# tets): 2
[0] TetGenMeshInsertVertexBW():     Located (4) tet (3, 0, 5, 1).
[0] TetGenMeshInsertVertexBW():     Insert a hull vertex.
[0] TetGenMeshInsertVertexBW():     Cavity formed: 2 tets, 2 faces.
[0] TetGenMeshInsertVertexBW():     Insert point 6
[0] TetGenMeshInsertVertexBW():     Walk distance (# tets): 2
[0] TetGenMeshInsertVertexBW():     Located (4) tet (0, 5, 7, 1).
[0] TetGenMeshInsertVertexBW():     Insert a hull vertex.
[0] TetGenMeshInsertVertexBW():     Cavity formed: 2 tets, 2 faces.
[0] TetGenMeshInsertVertexBW():     Insert point 2
[0] TetGenMeshBTreeSearch():     Get point 6 (cell size 1).
[0] TetGenMeshInsertVertexBW():     Walk distance (# tets): 2
[0] TetGenMeshInsertVertexBW():     Located (4) tet (3, 0, 6, 7).
[0] TetGenMeshInsertVertexBW():     Insert a hull vertex.
[0] TetGenMeshInsertVertexBW():     Cavity formed: 2 tets, 2 faces.
[0] TetGenMeshInsertVertexBW():     Insert point 4
[0] TetGenMeshBTreeSearch():     Get point 6 (cell size 2).
[0] TetGenMeshInsertVertexBW():     Walk distance (# tets): 3
[0] TetGenMeshInsertVertexBW():     Located (4) tet (6, 0, 5, 7).
[0] TetGenMeshInsertVertexBW():     Insert a hull vertex.
[0] TetGenMeshInsertVertexBW():     Cavity formed: 2 tets, 2 faces.
[0] TetGenMeshMeshSurface(): Creating surface mesh.
[0] TetGenMeshMakeIndex2PointMap():   Constructing mapping from indices to points.
[0] TetGenMeshTriangulate():     4 vertices, 4 segments, 0 holes, shmark: 1.
[0] TetGenMeshGetFacetAbovePoint():     Chosen abovepoint 0 for facet 1.
[0] TetGenMeshRecoverSegment():     Insert seg (4, 5).
[0] TetGenMeshRecoverSegment():     Insert seg (5, 7).
[0] TetGenMeshRecoverSegment():     Insert seg (7, 6).
[0] TetGenMeshRecoverSegment():     Insert seg (6, 4).
[0] TetGenMeshTriangulate():     4 vertices, 4 segments, 0 holes, shmark: 2.
[0] TetGenMeshGetFacetAbovePoint():     Chosen abovepoint 6 for facet 2.
[0] TetGenMeshRecoverSegment():     Insert seg (1, 0).
[0] TetGenMeshRecoverSegment():     Insert seg (0, 2).
[0] TetGenMeshRecoverSegment():     Insert seg (2, 3).
[0] TetGenMeshRecoverSegment():     Insert seg (3, 1).
[0] TetGenMeshTriangulate():     4 vertices, 4 segments, 0 holes, shmark: 3.
[0] TetGenMeshGetFacetAbovePoint():     Chosen abovepoint 6 for facet 3.
[0] TetGenMeshRecoverSegment():     Insert seg (0, 1).
[0] TetGenMeshRecoverSegment():     Insert seg (1, 5).
[0] TetGenMeshRecoverSegment():     Insert seg (5, 4).
[0] TetGenMeshRecoverSegment():     Insert seg (4, 0).
[0] TetGenMeshTriangulate():     4 vertices, 4 segments, 0 holes, shmark: 4.
[0] TetGenMeshGetFacetAbovePoint():     Chosen abovepoint 0 for facet 4.
[0] TetGenMeshRecoverSegment():     Insert seg (6, 7).
[0] TetGenMeshRecoverSegment():     Insert seg (7, 3).
[0] TetGenMeshRecoverSegment():     Insert seg (3, 2).
[0] TetGenMeshRecoverSegment():     Insert seg (2, 6).
[0] TetGenMeshTriangulate():     4 vertices, 4 segments, 0 holes, shmark: 5.
[0] TetGenMeshGetFacetAbovePoint():     Chosen abovepoint 7 for facet 5.
[0] TetGenMeshRecoverSegment():     Insert seg (0, 4).
[0] TetGenMeshRecoverSegment():     Insert seg (4, 6).
[0] TetGenMeshRecoverSegment():     Insert seg (6, 2).
[0] TetGenMeshRecoverSegment():     Insert seg (2, 0).
[0] TetGenMeshTriangulate():     4 vertices, 4 segments, 0 holes, shmark: 6.
[0] TetGenMeshGetFacetAbovePoint():     Chosen abovepoint 0 for facet 6.
[0] TetGenMeshRecoverSegment():     Insert seg (5, 1).
[0] TetGenMeshRecoverSegment():     Insert seg (1, 3).
[0] TetGenMeshRecoverSegment():     Insert seg (3, 7).
[0] TetGenMeshRecoverSegment():     Insert seg (7, 5).
[0] TetGenMeshUnifySegments():   Unifying segments.
[0] TetGenMeshMakeSubfaceMap():   Constructing mapping from points to subfaces.
[0] TetGenMeshUnifySegments():     Identifying 2 segments of (5  4).
[0] TetGenMeshUnifySegments():     Bond subfaces (5, 4, 7) and (5, 4, 0).
[0] TetGenMeshUnifySegments():     Bond subfaces (5, 4, 0) and (5, 4, 7).
[0] TetGenMeshUnifySegments():     Identifying 2 segments of (5  7).
[0] TetGenMeshUnifySegments():     Bond subfaces (5, 7, 4) and (5, 7, 3).
[0] TetGenMeshUnifySegments():     Bond subfaces (5, 7, 3) and (5, 7, 4).
[0] TetGenMeshUnifySegments():     Identifying 2 segments of (6  7).
[0] TetGenMeshUnifySegments():     Bond subfaces (6, 7, 4) and (6, 7, 3).
[0] TetGenMeshUnifySegments():     Bond subfaces (6, 7, 3) and (6, 7, 4).
[0] TetGenMeshUnifySegments():     Identifying 2 segments of (4  6).
[0] TetGenMeshUnifySegments():     Bond subfaces (4, 6, 7) and (4, 6, 0).
[0] TetGenMeshUnifySegments():     Bond subfaces (4, 6, 0) and (4, 6, 7).
[0] TetGenMeshUnifySegments():     Identifying 2 segments of (0  1).
[0] TetGenMeshUnifySegments():     Bond subfaces (0, 1, 2) and (0, 1, 5).
[0] TetGenMeshUnifySegments():     Bond subfaces (0, 1, 5) and (0, 1, 2).
[0] TetGenMeshUnifySegments():     Identifying 2 segments of (0  2).
[0] TetGenMeshUnifySegments():     Bond subfaces (0, 2, 1) and (0, 2, 6).
[0] TetGenMeshUnifySegments():     Bond subfaces (0, 2, 6) and (0, 2, 1).
[0] TetGenMeshUnifySegments():     Identifying 2 segments of (3  2).
[0] TetGenMeshUnifySegments():     Bond subfaces (3, 2, 1) and (3, 2, 6).
[0] TetGenMeshUnifySegments():     Bond subfaces (3, 2, 6) and (3, 2, 1).
[0] TetGenMeshUnifySegments():     Identifying 2 segments of (1  3).
[0] TetGenMeshUnifySegments():     Bond subfaces (1, 3, 2) and (1, 3, 5).
[0] TetGenMeshUnifySegments():     Bond subfaces (1, 3, 5) and (1, 3, 2).
[0] TetGenMeshUnifySegments():     Identifying 2 segments of (1  5).
[0] TetGenMeshUnifySegments():     Bond subfaces (1, 5, 0) and (1, 5, 3).
[0] TetGenMeshUnifySegments():     Bond subfaces (1, 5, 3) and (1, 5, 0).
[0] TetGenMeshUnifySegments():     Identifying 2 segments of (0  4).
[0] TetGenMeshUnifySegments():     Bond subfaces (0, 4, 5) and (0, 4, 6).
[0] TetGenMeshUnifySegments():     Bond subfaces (0, 4, 6) and (0, 4, 5).
[0] TetGenMeshUnifySegments():     Identifying 2 segments of (7  3).
[0] TetGenMeshUnifySegments():     Bond subfaces (7, 3, 6) and (7, 3, 5).
[0] TetGenMeshUnifySegments():     Bond subfaces (7, 3, 5) and (7, 3, 6).
[0] TetGenMeshUnifySegments():     Identifying 2 segments of (6  2).
[0] TetGenMeshUnifySegments():     Bond subfaces (6, 2, 3) and (6, 2, 0).
[0] TetGenMeshUnifySegments():     Bond subfaces (6, 2, 0) and (6, 2, 3).
[0] TetGenMeshMergeFacets():   Merging coplanar facets.
[0] TetGenMeshMakePoint2SegMap():   Constructing mapping from points to segments.
[0] TetGenMeshMarkAcuteVertices():   Marking acute vertices.
[0] TetGenMeshMakeSegmentMap():   Constructing mapping from points to segments.
[0] TetGenMeshMarkAcuteVertices():   0 acute vertices.
[0] TetGenMeshFormSkeleton(): Recovering boundaries.
[0] TetGenMeshDelaunizeSegments2():   Delaunizing segments.
[0] TetGenMeshScoutSegment2():     Scout seg (0, 1).
[0] TetGenMeshFindDirection2():       From tet (0, 7, 6, 4) to 1.
[0] TetGenMeshFindDirection2():       From tet (0, 7, 4, 5) to 1.
[0] TetGenMeshFindDirection2():       From tet (0, 7, 5, 1) to 1.
[0] TetGenMeshScoutSegment2():     Scout seg (0, 2).
[0] TetGenMeshFindDirection2():       From tet (0, 7, 6, 4) to 2.
[0] TetGenMeshFindDirection2():       From tet (0, 6, 7, 2) to 2.
[0] TetGenMeshScoutSegment2():     Scout seg (1, 5).
[0] TetGenMeshScoutSegment2():     Scout seg (3, 2).
[0] TetGenMeshFindDirection2():       From tet (3, 0, 7, 2) to 2.
[0] TetGenMeshScoutSegment2():     Scout seg (6, 7).
[0] TetGenMeshScoutSegment2():     Scout seg (6, 2).
[0] TetGenMeshFindDirection2():       From tet (6, 0, 7, 4) to 2.
[0] TetGenMeshFindDirection2():       From tet (6, 7, 0, 2) to 2.
[0] TetGenMeshScoutSegment2():     Scout seg (5, 4).
[0] TetGenMeshFindDirection2():       From tet (5, 7, 0, 4) to 4.
[0] TetGenMeshScoutSegment2():     Scout seg (5, 7).
[0] TetGenMeshScoutSegment2():     Scout seg (7, 3).
[0] TetGenMeshFindDirection2():       From tet (7, 6, 0, 4) to 3.
[0] TetGenMeshFindDirection2():       From tet (7, 0, 6, 2) to 3.
[0] TetGenMeshFindDirection2():       From tet (7, 0, 2, 3) to 3.
[0] TetGenMeshScoutSegment2():     Scout seg (1, 3).
[0] TetGenMeshFindDirection2():       From tet (1, 0, 5, 7) to 3.
[0] TetGenMeshFindDirection2():       From tet (1, 0, 7, 3) to 3.
[0] TetGenMeshScoutSegment2():     Scout seg (0, 4).
[0] TetGenMeshFindDirection2():       From tet (0, 7, 6, 4) to 4.
[0] TetGenMeshScoutSegment2():     Scout seg (4, 6).
[0] TetGenMeshDelaunizeSegments2():   0 protecting points.
[0] TetGenMeshConstrainedFacets2():   Constraining facets.
[0] TetGenMeshConstrainedFacets2():   Recover facet #1: 2 subfaces, 4 vertices.
[0] TetGenMeshFindDirection2():       From tet (3, 0, 7, 2) to 5.
[0] TetGenMeshFindDirection2():       From tet (3, 7, 0, 1) to 5.
[0] TetGenMeshScoutSubface():     Scout subface (5, 7, 3) (11).
[0] TetGenMeshScoutCrossTet():     Found a co-facet face (5, 7, 1) op (0).
[0] TetGenMeshFlip22Sub():     Flip subedge (5, 3) to (7, 1).
[0] TetGenMeshFindDirection2():       From tet (1, 0, 5, 7) to 3.
[0] TetGenMeshFindDirection2():       From tet (1, 0, 7, 3) to 3.
[0] TetGenMeshScoutSubface():     Scout subface (1, 3, 7) (11).
[0] TetGenMeshFindDirection2():       From tet (7, 6, 0, 4) to 1.
[0] TetGenMeshFindDirection2():       From tet (7, 4, 0, 5) to 1.
[0] TetGenMeshFindDirection2():       From tet (7, 5, 0, 1) to 1.
[0] TetGenMeshScoutSubface():     Scout subface (7, 1, 5) (11).
[0] TetGenMeshConstrainedFacets2():   Recover facet #2: 2 subfaces, 4 vertices.
[0] TetGenMeshFindDirection2():       From tet (4, 6, 7, 0) to 5.
[0] TetGenMeshFindDirection2():       From tet (4, 0, 7, 5) to 5.
[0] TetGenMeshScoutSubface():     Scout subface (4, 5, 7) (10).
[0] TetGenMeshFindDirection2():       From tet (7, 6, 0, 4) to 4.
[0] TetGenMeshScoutSubface():     Scout subface (7, 4, 6) (10).
[0] TetGenMeshConstrainedFacets2():   Recover facet #3: 2 subfaces, 4 vertices.
[0] TetGenMeshScoutSubface():     Scout subface (6, 0, 2) (9).
[0] TetGenMeshFindDirection2():       From tet (0, 7, 6, 4) to 4.
[0] TetGenMeshScoutSubface():     Scout subface (0, 4, 6) (9).
[0] TetGenMeshConstrainedFacets2():   Recover facet #4: 2 subfaces, 4 vertices.
[0] TetGenMeshScoutSubface():     Scout subface (5, 0, 4) (8).
[0] TetGenMeshFindDirection2():       From tet (0, 7, 6, 4) to 1.
[0] TetGenMeshFindDirection2():       From tet (0, 7, 4, 5) to 1.
[0] TetGenMeshFindDirection2():       From tet (0, 7, 5, 1) to 1.
[0] TetGenMeshScoutSubface():     Scout subface (0, 1, 5) (8).
[0] TetGenMeshConstrainedFacets2():   Recover facet #5: 2 subfaces, 4 vertices.
[0] TetGenMeshScoutSubface():     Scout subface (1, 0, 2) (7).
[0] TetGenMeshScoutCrossTet():     Found a co-facet face (1, 0, 3) op (7).
[0] TetGenMeshFlip22Sub():     Flip subedge (1, 2) to (0, 3).
[0] TetGenMeshFindDirection2():       From tet (0, 7, 6, 4) to 3.
[0] TetGenMeshFindDirection2():       From tet (0, 6, 7, 2) to 3.
[0] TetGenMeshFindDirection2():       From tet (0, 2, 7, 3) to 3.
[0] TetGenMeshScoutSubface():     Scout subface (0, 3, 2) (7).
[0] TetGenMeshFindDirection2():       From tet (3, 0, 7, 2) to 1.
[0] TetGenMeshFindDirection2():       From tet (3, 7, 0, 1) to 1.
[0] TetGenMeshScoutSubface():     Scout subface (3, 1, 0) (7).
[0] TetGenMeshConstrainedFacets2():   Recover facet #6: 2 subfaces, 4 vertices.
[0] TetGenMeshScoutSubface():     Scout subface (6, 7, 3) (2).
[0] TetGenMeshScoutCrossTet():     Found a co-facet face (7, 6, 2) op (0).
[0] TetGenMeshFlip22Sub():     Flip subedge (6, 3) to (7, 2).
[0] TetGenMeshFindDirection2():       From tet (7, 6, 0, 4) to 2.
[0] TetGenMeshFindDirection2():       From tet (7, 0, 6, 2) to 2.
[0] TetGenMeshScoutSubface():     Scout subface (7, 2, 3) (2).
[0] TetGenMeshScoutSubface():     Scout subface (2, 6, 7) (2).
[0] TetGenMeshConstrainedFacets2():   3 subedge flips  0 cavities remeshed.
[0] TetGenMeshCarveHoles(): Removing exterior tetrahedra.
[0] TetGenMeshInfectHull():   Marking concavities for elimination.
[0] TetGenMeshPlague():   Marking neighbors of marked tetrahedra.
[0] TetGenMeshRemoveHoleTets():   Deleting marked tetrahedra.
[0] TetGenMeshMakePoint2TetMap():   Constructing mapping from points to tetrahedra.
[0] TetGenMeshOptimize(): Repairing mesh.
[0] TetGenMeshEnforceQuality(): Adding Steiner points to enforce quality.
[0] TetGenMeshMarkSharpSegments():   Marking sharp segments.
[0] TetGenMeshDecideFeaturePointSizes():   Deciding feature-point sizes.
[0] TetGenMeshMakeSegmentMap():   Constructing mapping from points to segments.
[0] TetGenMeshDecideFeaturePointSizes():   0 feature points.
[0] TetGenMeshEnforceQuality():   0 split points.
[0] TetGenMeshEnforceQuality():   0 split points.
[0] TetGenMeshEnforceQuality():   0 refinement points.
[0] TetGenMeshEnforceQuality():   Totally added 0 points.
[0] TetGenMeshOptimize(): Optimizing mesh.
[0] TetGenMeshOptimize():   level = 3.
[0] TetGenMeshOptimize():   level = 3 pass 0.
[0] TetGenMeshOutputNodes(): Writing nodes.
[0] TetGenMeshOutputElements(): Writing elements.
[0] TetGenMeshOutputSubfaces(): Writing faces.
PetscSection Object: 1 MPI processes
  type not yet set
3 fields
  field 0 with 1 components
Process 0:
  (   0) dim  0 offset   0
  (   1) dim  0 offset   3
  (   2) dim  0 offset   6
  (   3) dim  0 offset   9
  (   4) dim  0 offset  12
  (   5) dim  0 offset  15
  (   6) dim  1 offset  18 constrained 0
  (   7) dim  1 offset  19 constrained 0
  (   8) dim  1 offset  20 constrained 0
  (   9) dim  1 offset  21 constrained 0
  (  10) dim  1 offset  22 constrained 0
  (  11) dim  1 offset  23 constrained 0
  (  12) dim  1 offset  24 constrained 0
  (  13) dim  1 offset  25 constrained 0
  (  14) dim  0 offset  26
  (  15) dim  0 offset  28
  (  16) dim  0 offset  30
  (  17) dim  0 offset  32
  (  18) dim  0 offset  34
  (  19) dim  0 offset  36
  (  20) dim  0 offset  38
  (  21) dim  0 offset  40
  (  22) dim  0 offset  42
  (  23) dim  0 offset  44
  (  24) dim  0 offset  46
  (  25) dim  0 offset  48
  (  26) dim  0 offset  50
  (  27) dim  0 offset  52
  (  28) dim  0 offset  54
  (  29) dim  0 offset  56
  (  30) dim  0 offset  58
  (  31) dim  0 offset  60
  (  32) dim  0 offset  62
  (  33) dim  0 offset  62
  (  34) dim  0 offset  62
  (  35) dim  0 offset  62
  (  36) dim  0 offset  62
  (  37) dim  0 offset  62
  (  38) dim  0 offset  62
  (  39) dim  0 offset  62
  (  40) dim  0 offset  62
  (  41) dim  0 offset  62
  (  42) dim  0 offset  62
  (  43) dim  0 offset  62
  (  44) dim  0 offset  62
  (  45) dim  0 offset  62
  (  46) dim  0 offset  62
  (  47) dim  0 offset  62
  (  48) dim  0 offset  62
  (  49) dim  0 offset  62
  (  50) dim  0 offset  62
  field 1 with 3 components
Process 0:
  (   0) dim  3 offset   0
  (   1) dim  3 offset   3
  (   2) dim  3 offset   6
  (   3) dim  3 offset   9
  (   4) dim  3 offset  12
  (   5) dim  3 offset  15
  (   6) dim  0 offset  19
  (   7) dim  0 offset  20
  (   8) dim  0 offset  21
  (   9) dim  0 offset  22
  (  10) dim  0 offset  23
  (  11) dim  0 offset  24
  (  12) dim  0 offset  25
  (  13) dim  0 offset  26
  (  14) dim  0 offset  26
  (  15) dim  0 offset  28
  (  16) dim  0 offset  30
  (  17) dim  0 offset  32
  (  18) dim  0 offset  34
  (  19) dim  0 offset  36
  (  20) dim  0 offset  38
  (  21) dim  0 offset  40
  (  22) dim  0 offset  42
  (  23) dim  0 offset  44
  (  24) dim  0 offset  46
  (  25) dim  0 offset  48
  (  26) dim  0 offset  50
  (  27) dim  0 offset  52
  (  28) dim  0 offset  54
  (  29) dim  0 offset  56
  (  30) dim  0 offset  58
  (  31) dim  0 offset  60
  (  32) dim  0 offset  62
  (  33) dim  0 offset  62
  (  34) dim  0 offset  62
  (  35) dim  0 offset  62
  (  36) dim  0 offset  62
  (  37) dim  0 offset  62
  (  38) dim  0 offset  62
  (  39) dim  0 offset  62
  (  40) dim  0 offset  62
  (  41) dim  0 offset  62
  (  42) dim  0 offset  62
  (  43) dim  0 offset  62
  (  44) dim  0 offset  62
  (  45) dim  0 offset  62
  (  46) dim  0 offset  62
  (  47) dim  0 offset  62
  (  48) dim  0 offset  62
  (  49) dim  0 offset  62
  (  50) dim  0 offset  62
  field 2 with 2 components
Process 0:
  (   0) dim  0 offset   3
  (   1) dim  0 offset   6
  (   2) dim  0 offset   9
  (   3) dim  0 offset  12
  (   4) dim  0 offset  15
  (   5) dim  0 offset  18
  (   6) dim  0 offset  19
  (   7) dim  0 offset  20
  (   8) dim  0 offset  21
  (   9) dim  0 offset  22
  (  10) dim  0 offset  23
  (  11) dim  0 offset  24
  (  12) dim  0 offset  25
  (  13) dim  0 offset  26
  (  14) dim  2 offset  26
  (  15) dim  2 offset  28
  (  16) dim  2 offset  30
  (  17) dim  2 offset  32
  (  18) dim  2 offset  34
  (  19) dim  2 offset  36
  (  20) dim  2 offset  38
  (  21) dim  2 offset  40
  (  22) dim  2 offset  42
  (  23) dim  2 offset  44
  (  24) dim  2 offset  46
  (  25) dim  2 offset  48
  (  26) dim  2 offset  50
  (  27) dim  2 offset  52
  (  28) dim  2 offset  54
  (  29) dim  2 offset  56
  (  30) dim  2 offset  58
  (  31) dim  2 offset  60
  (  32) dim  0 offset  62
  (  33) dim  0 offset  62
  (  34) dim  0 offset  62
  (  35) dim  0 offset  62
  (  36) dim  0 offset  62
  (  37) dim  0 offset  62
  (  38) dim  0 offset  62
  (  39) dim  0 offset  62
  (  40) dim  0 offset  62
  (  41) dim  0 offset  62
  (  42) dim  0 offset  62
  (  43) dim  0 offset  62
  (  44) dim  0 offset  62
  (  45) dim  0 offset  62
  (  46) dim  0 offset  62
  (  47) dim  0 offset  62
  (  48) dim  0 offset  62
  (  49) dim  0 offset  62
  (  50) dim  0 offset  62
[0] DMPlexGetFieldType_Internal(): Could not classify VTK output type of section field 2 "w"
