/*========================================================================= Program: Visualization Toolkit Module: $RCSfile$ Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen All rights reserved. See Copyright.txt or http://www.kitware.com/Copyright.htm for details. This software is distributed WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the above copyright notice for more information. =========================================================================*/ // Testings of the vtkReebGraph features. // Julien Tierny , 2010 #include "vtkActor.h" #include "vtkAreaContourSpectrumFilter.h" #include "vtkCamera.h" #include "vtkDataSetAttributes.h" #include "vtkDoubleArray.h" #include "vtkEdgeListIterator.h" #include "vtkIdList.h" #include "vtkLight.h" #include "vtkObjectFactory.h" #include "vtkPointData.h" #include "vtkPolyData.h" #include "vtkPolyDataMapper.h" #include "vtkPolyDataToReebGraphFilter.h" #include "vtkProperty.h" #include "vtkPNGWriter.h" #include "vtkReebGraph.h" #include "vtkReebGraphSurfaceSkeletonFilter.h" #include "vtkReebGraphSimplificationFilter.h" #include "vtkReebGraphSimplificationMetric.h" #include "vtkReebGraphToJoinSplitTreeFilter.h" #include "vtkReebGraphVolumeSkeletonFilter.h" #include "vtkRenderer.h" #include "vtkRenderWindow.h" #include "vtkRenderWindowInteractor.h" #include "vtkSphereSource.h" #include "vtkTable.h" #include "vtkTriangle.h" #include "vtkUnstructuredGrid.h" #include "vtkUnstructuredGridToReebGraphFilter.h" #include "vtkVariantArray.h" #include "vtkVolumeContourSpectrumFilter.h" #include "vtkWindowToImageFilter.h" #include // loading code generated automatically. // put whatever mesh you like here. int LoadVolumeMesh(vtkUnstructuredGrid *vMesh) { vMesh->Allocate(); vtkPoints *points = vtkPoints::New(); points->InsertNextPoint( 0, 0, 0); points->InsertNextPoint( 0.632981, 0.774167, 0); points->InsertNextPoint( 1.26596, 1.54833, 0); points->InsertNextPoint( 1.89894, 2.3225, 0); points->InsertNextPoint( -0.774167, 0.632981, 0); points->InsertNextPoint( -0.141186, 1.40715, 0); points->InsertNextPoint( 0.491796, 2.18132, 0); points->InsertNextPoint( 1.12478, 2.95548, 0); points->InsertNextPoint( -1.54833, 1.26596, 0); points->InsertNextPoint( -0.915353, 2.04013, 0); points->InsertNextPoint( -0.282372, 2.8143, 0); points->InsertNextPoint( 0.35061, 3.58846, 0); points->InsertNextPoint( -2.3225, 1.89894, 0); points->InsertNextPoint( -1.68952, 2.67311, 0); points->InsertNextPoint( -1.05654, 3.44728, 0); points->InsertNextPoint( -0.423557, 4.22145, 0); points->InsertNextPoint( -3.09667, 2.53193, 0); points->InsertNextPoint( -2.46369, 3.30609, 0); points->InsertNextPoint( -1.83071, 4.08026, 0); points->InsertNextPoint( -1.19772, 4.85443, 0); points->InsertNextPoint( -3.87084, 3.16491, 0); points->InsertNextPoint( -3.23785, 3.93907, 0); points->InsertNextPoint( -2.60487, 4.71324, 0); points->InsertNextPoint( -1.97189, 5.48741, 0); points->InsertNextPoint( -4.645, 3.79789, 0); points->InsertNextPoint( -4.01202, 4.57205, 0); points->InsertNextPoint( -3.37904, 5.34622, 0); points->InsertNextPoint( -2.74606, 6.12039, 0); points->InsertNextPoint( -5.41917, 4.43087, 0); points->InsertNextPoint( -4.78619, 5.20504, 0); points->InsertNextPoint( -4.15321, 5.9792, 0); points->InsertNextPoint( -3.52023, 6.75337, 0); points->InsertNextPoint( 0, 0, 1); points->InsertNextPoint( 0.632981, 0.774167, 1); points->InsertNextPoint( 1.26596, 1.54833, 1); points->InsertNextPoint( 1.89894, 2.3225, 1); points->InsertNextPoint( -0.774167, 0.632981, 1); points->InsertNextPoint( -0.141186, 1.40715, 1); points->InsertNextPoint( 0.491796, 2.18132, 1); points->InsertNextPoint( 1.12478, 2.95548, 1); points->InsertNextPoint( -1.54833, 1.26596, 1); points->InsertNextPoint( -0.915353, 2.04013, 1); points->InsertNextPoint( -0.282372, 2.8143, 1); points->InsertNextPoint( 0.35061, 3.58846, 1); points->InsertNextPoint( -2.3225, 1.89894, 1); points->InsertNextPoint( -1.68952, 2.67311, 1); points->InsertNextPoint( -1.05654, 3.44728, 1); points->InsertNextPoint( -0.423557, 4.22145, 1); points->InsertNextPoint( -3.09667, 2.53193, 1); points->InsertNextPoint( -2.46369, 3.30609, 1); points->InsertNextPoint( -1.83071, 4.08026, 1); points->InsertNextPoint( -1.19772, 4.85443, 1); points->InsertNextPoint( -3.87084, 3.16491, 1); points->InsertNextPoint( -3.23785, 3.93907, 1); points->InsertNextPoint( -2.60487, 4.71324, 1); points->InsertNextPoint( -1.97189, 5.48741, 1); points->InsertNextPoint( -4.645, 3.79789, 1); points->InsertNextPoint( -4.01202, 4.57205, 1); points->InsertNextPoint( -3.37904, 5.34622, 1); points->InsertNextPoint( -2.74606, 6.12039, 1); points->InsertNextPoint( -5.41917, 4.43087, 1); points->InsertNextPoint( -4.78619, 5.20504, 1); points->InsertNextPoint( -4.15321, 5.9792, 1); points->InsertNextPoint( -3.52023, 6.75337, 1); points->InsertNextPoint( -0.0352964, 0.351787, 0.25); points->InsertNextPoint( -0.264135, 0.86182, 0.25); points->InsertNextPoint( -0.228838, 0.510032, 0.5); points->InsertNextPoint( 0.0876524, 0.897116, 0.5); points->InsertNextPoint( -0.105889, 1.05536, 0.75); points->InsertNextPoint( 0.122949, 0.545329, 0.75); points->InsertNextPoint( 0.597685, 1.12595, 0.25); points->InsertNextPoint( 0.368847, 1.63599, 0.25); points->InsertNextPoint( 0.404143, 1.2842, 0.5); points->InsertNextPoint( 0.720634, 1.67128, 0.5); points->InsertNextPoint( 0.527092, 1.82953, 0.75); points->InsertNextPoint( 0.75593, 1.3195, 0.75); points->InsertNextPoint( 1.23067, 1.90012, 0.25); points->InsertNextPoint( 1.00183, 2.41015, 0.25); points->InsertNextPoint( 1.03712, 2.05837, 0.5); points->InsertNextPoint( 1.35362, 2.44545, 0.5); points->InsertNextPoint( 1.16007, 2.6037, 0.75); points->InsertNextPoint( 1.38891, 2.09366, 0.75); points->InsertNextPoint( -0.809464, 0.984768, 0.25); points->InsertNextPoint( -1.0383, 1.4948, 0.25); points->InsertNextPoint( -1.00301, 1.14301, 0.5); points->InsertNextPoint( -0.686515, 1.5301, 0.5); points->InsertNextPoint( -0.880056, 1.68834, 0.75); points->InsertNextPoint( -0.651218, 1.17831, 0.75); points->InsertNextPoint( 0.456499, 2.5331, 0.25); points->InsertNextPoint( 0.227661, 3.04313, 0.25); points->InsertNextPoint( 0.262957, 2.69135, 0.5); points->InsertNextPoint( 0.579448, 3.07843, 0.5); points->InsertNextPoint( 0.385906, 3.23668, 0.75); points->InsertNextPoint( 0.614744, 2.72664, 0.75); points->InsertNextPoint( -1.58363, 1.61775, 0.25); points->InsertNextPoint( -1.81247, 2.12778, 0.25); points->InsertNextPoint( -1.77717, 1.776, 0.5); points->InsertNextPoint( -1.46068, 2.16308, 0.5); points->InsertNextPoint( -1.65422, 2.32132, 0.75); points->InsertNextPoint( -1.42539, 1.81129, 0.75); points->InsertNextPoint( -0.950649, 2.39192, 0.25); points->InsertNextPoint( -1.17949, 2.90195, 0.25); points->InsertNextPoint( -1.14419, 2.55016, 0.5); points->InsertNextPoint( -0.8277, 2.93725, 0.5); points->InsertNextPoint( -1.02124, 3.09549, 0.75); points->InsertNextPoint( -0.792404, 2.58546, 0.75); points->InsertNextPoint( -0.317668, 3.16608, 0.25); points->InsertNextPoint( -0.546506, 3.67612, 0.25); points->InsertNextPoint( -0.51121, 3.32433, 0.5); points->InsertNextPoint( -0.194719, 3.71141, 0.5); points->InsertNextPoint( -0.388261, 3.86966, 0.75); points->InsertNextPoint( -0.159423, 3.35963, 0.75); points->InsertNextPoint( -2.3578, 2.25073, 0.25); points->InsertNextPoint( -2.58664, 2.76076, 0.25); points->InsertNextPoint( -2.55134, 2.40898, 0.5); points->InsertNextPoint( -2.23485, 2.79606, 0.5); points->InsertNextPoint( -2.42839, 2.95431, 0.75); points->InsertNextPoint( -2.19955, 2.44427, 0.75); points->InsertNextPoint( -1.09184, 3.79907, 0.25); points->InsertNextPoint( -1.32067, 4.3091, 0.25); points->InsertNextPoint( -1.28538, 3.95731, 0.5); points->InsertNextPoint( -0.968886, 4.34439, 0.5); points->InsertNextPoint( -1.16243, 4.50264, 0.75); points->InsertNextPoint( -0.93359, 3.99261, 0.75); points->InsertNextPoint( -3.13196, 2.88371, 0.25); points->InsertNextPoint( -3.3608, 3.39374, 0.25); points->InsertNextPoint( -3.32551, 3.04196, 0.5); points->InsertNextPoint( -3.00902, 3.42904, 0.5); points->InsertNextPoint( -3.20256, 3.58729, 0.75); points->InsertNextPoint( -2.97372, 3.07725, 0.75); points->InsertNextPoint( -2.49898, 3.65788, 0.25); points->InsertNextPoint( -2.72782, 4.16791, 0.25); points->InsertNextPoint( -2.69253, 3.81612, 0.5); points->InsertNextPoint( -2.37603, 4.20321, 0.5); points->InsertNextPoint( -2.56958, 4.36145, 0.75); points->InsertNextPoint( -2.34074, 3.85142, 0.75); points->InsertNextPoint( -1.866, 4.43205, 0.25); points->InsertNextPoint( -2.09484, 4.94208, 0.25); points->InsertNextPoint( -2.05954, 4.59029, 0.5); points->InsertNextPoint( -1.74305, 4.97738, 0.5); points->InsertNextPoint( -1.9366, 5.13562, 0.75); points->InsertNextPoint( -1.70776, 4.62559, 0.75); points->InsertNextPoint( -3.90613, 3.51669, 0.25); points->InsertNextPoint( -4.13497, 4.02673, 0.25); points->InsertNextPoint( -4.09967, 3.67494, 0.5); points->InsertNextPoint( -3.78318, 4.06202, 0.5); points->InsertNextPoint( -3.97672, 4.22027, 0.75); points->InsertNextPoint( -3.74789, 3.71024, 0.75); points->InsertNextPoint( -2.64017, 5.06503, 0.25); points->InsertNextPoint( -2.86901, 5.57506, 0.25); points->InsertNextPoint( -2.83371, 5.22327, 0.5); points->InsertNextPoint( -2.51722, 5.61036, 0.5); points->InsertNextPoint( -2.71076, 5.7686, 0.75); points->InsertNextPoint( -2.48192, 5.25857, 0.75); points->InsertNextPoint( -4.6803, 4.14967, 0.25); points->InsertNextPoint( -4.90914, 4.65971, 0.25); points->InsertNextPoint( -4.87384, 4.30792, 0.5); points->InsertNextPoint( -4.55735, 4.695, 0.5); points->InsertNextPoint( -4.75089, 4.85325, 0.75); points->InsertNextPoint( -4.52205, 4.34322, 0.75); points->InsertNextPoint( -4.04732, 4.92384, 0.25); points->InsertNextPoint( -4.27616, 5.43387, 0.25); points->InsertNextPoint( -4.24086, 5.08209, 0.5); points->InsertNextPoint( -3.92437, 5.46917, 0.5); points->InsertNextPoint( -4.11791, 5.62742, 0.75); points->InsertNextPoint( -3.88907, 5.11738, 0.75); points->InsertNextPoint( -3.41434, 5.69801, 0.25); points->InsertNextPoint( -3.64317, 6.20804, 0.25); points->InsertNextPoint( -3.60788, 5.85625, 0.5); points->InsertNextPoint( -3.29139, 6.24334, 0.5); points->InsertNextPoint( -3.48493, 6.40158, 0.75); points->InsertNextPoint( -3.25609, 5.89155, 0.75); vMesh->SetPoints(points); points->Delete(); vtkIdType *vertexIds = (vtkIdType *) malloc(sizeof(vtkIdType)*4); vertexIds[0] = 0; vertexIds[1] = 1; vertexIds[2] = 32; vertexIds[3] = 64; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 4; vertexIds[1] = 1; vertexIds[2] = 36; vertexIds[3] = 65; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 4; vertexIds[1] = 1; vertexIds[2] = 32; vertexIds[3] = 66; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 1; vertexIds[1] = 33; vertexIds[2] = 5; vertexIds[3] = 67; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 33; vertexIds[1] = 5; vertexIds[2] = 36; vertexIds[3] = 68; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 1; vertexIds[1] = 33; vertexIds[2] = 32; vertexIds[3] = 69; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 1; vertexIds[1] = 2; vertexIds[2] = 33; vertexIds[3] = 70; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 5; vertexIds[1] = 2; vertexIds[2] = 37; vertexIds[3] = 71; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 5; vertexIds[1] = 2; vertexIds[2] = 33; vertexIds[3] = 72; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 2; vertexIds[1] = 34; vertexIds[2] = 6; vertexIds[3] = 73; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 34; vertexIds[1] = 6; vertexIds[2] = 37; vertexIds[3] = 74; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 2; vertexIds[1] = 34; vertexIds[2] = 33; vertexIds[3] = 75; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 2; vertexIds[1] = 3; vertexIds[2] = 34; vertexIds[3] = 76; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 6; vertexIds[1] = 3; vertexIds[2] = 38; vertexIds[3] = 77; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 6; vertexIds[1] = 3; vertexIds[2] = 34; vertexIds[3] = 78; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 3; vertexIds[1] = 35; vertexIds[2] = 7; vertexIds[3] = 79; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 35; vertexIds[1] = 7; vertexIds[2] = 38; vertexIds[3] = 80; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 3; vertexIds[1] = 35; vertexIds[2] = 34; vertexIds[3] = 81; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 4; vertexIds[1] = 5; vertexIds[2] = 36; vertexIds[3] = 82; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 8; vertexIds[1] = 5; vertexIds[2] = 40; vertexIds[3] = 83; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 8; vertexIds[1] = 5; vertexIds[2] = 36; vertexIds[3] = 84; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 5; vertexIds[1] = 37; vertexIds[2] = 9; vertexIds[3] = 85; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 37; vertexIds[1] = 9; vertexIds[2] = 40; vertexIds[3] = 86; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 5; vertexIds[1] = 37; vertexIds[2] = 36; vertexIds[3] = 87; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 6; vertexIds[1] = 7; vertexIds[2] = 38; vertexIds[3] = 88; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 10; vertexIds[1] = 7; vertexIds[2] = 42; vertexIds[3] = 89; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 10; vertexIds[1] = 7; vertexIds[2] = 38; vertexIds[3] = 90; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 7; vertexIds[1] = 39; vertexIds[2] = 11; vertexIds[3] = 91; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 39; vertexIds[1] = 11; vertexIds[2] = 42; vertexIds[3] = 92; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 7; vertexIds[1] = 39; vertexIds[2] = 38; vertexIds[3] = 93; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 8; vertexIds[1] = 9; vertexIds[2] = 40; vertexIds[3] = 94; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 12; vertexIds[1] = 9; vertexIds[2] = 44; vertexIds[3] = 95; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 12; vertexIds[1] = 9; vertexIds[2] = 40; vertexIds[3] = 96; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 9; vertexIds[1] = 41; vertexIds[2] = 13; vertexIds[3] = 97; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 41; vertexIds[1] = 13; vertexIds[2] = 44; vertexIds[3] = 98; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 9; vertexIds[1] = 41; vertexIds[2] = 40; vertexIds[3] = 99; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 9; vertexIds[1] = 10; vertexIds[2] = 41; vertexIds[3] = 100; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 13; vertexIds[1] = 10; vertexIds[2] = 45; vertexIds[3] = 101; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 13; vertexIds[1] = 10; vertexIds[2] = 41; vertexIds[3] = 102; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 10; vertexIds[1] = 42; vertexIds[2] = 14; vertexIds[3] = 103; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 42; vertexIds[1] = 14; vertexIds[2] = 45; vertexIds[3] = 104; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 10; vertexIds[1] = 42; vertexIds[2] = 41; vertexIds[3] = 105; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 10; vertexIds[1] = 11; vertexIds[2] = 42; vertexIds[3] = 106; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 14; vertexIds[1] = 11; vertexIds[2] = 46; vertexIds[3] = 107; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 14; vertexIds[1] = 11; vertexIds[2] = 42; vertexIds[3] = 108; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 11; vertexIds[1] = 43; vertexIds[2] = 15; vertexIds[3] = 109; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 43; vertexIds[1] = 15; vertexIds[2] = 46; vertexIds[3] = 110; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 11; vertexIds[1] = 43; vertexIds[2] = 42; vertexIds[3] = 111; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 12; vertexIds[1] = 13; vertexIds[2] = 44; vertexIds[3] = 112; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 16; vertexIds[1] = 13; vertexIds[2] = 48; vertexIds[3] = 113; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 16; vertexIds[1] = 13; vertexIds[2] = 44; vertexIds[3] = 114; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 13; vertexIds[1] = 45; vertexIds[2] = 17; vertexIds[3] = 115; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 45; vertexIds[1] = 17; vertexIds[2] = 48; vertexIds[3] = 116; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 13; vertexIds[1] = 45; vertexIds[2] = 44; vertexIds[3] = 117; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 14; vertexIds[1] = 15; vertexIds[2] = 46; vertexIds[3] = 118; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 18; vertexIds[1] = 15; vertexIds[2] = 50; vertexIds[3] = 119; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 18; vertexIds[1] = 15; vertexIds[2] = 46; vertexIds[3] = 120; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 15; vertexIds[1] = 47; vertexIds[2] = 19; vertexIds[3] = 121; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 47; vertexIds[1] = 19; vertexIds[2] = 50; vertexIds[3] = 122; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 15; vertexIds[1] = 47; vertexIds[2] = 46; vertexIds[3] = 123; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 16; vertexIds[1] = 17; vertexIds[2] = 48; vertexIds[3] = 124; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 20; vertexIds[1] = 17; vertexIds[2] = 52; vertexIds[3] = 125; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 20; vertexIds[1] = 17; vertexIds[2] = 48; vertexIds[3] = 126; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 17; vertexIds[1] = 49; vertexIds[2] = 21; vertexIds[3] = 127; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 49; vertexIds[1] = 21; vertexIds[2] = 52; vertexIds[3] = 128; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 17; vertexIds[1] = 49; vertexIds[2] = 48; vertexIds[3] = 129; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 17; vertexIds[1] = 18; vertexIds[2] = 49; vertexIds[3] = 130; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 21; vertexIds[1] = 18; vertexIds[2] = 53; vertexIds[3] = 131; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 21; vertexIds[1] = 18; vertexIds[2] = 49; vertexIds[3] = 132; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 18; vertexIds[1] = 50; vertexIds[2] = 22; vertexIds[3] = 133; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 50; vertexIds[1] = 22; vertexIds[2] = 53; vertexIds[3] = 134; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 18; vertexIds[1] = 50; vertexIds[2] = 49; vertexIds[3] = 135; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 18; vertexIds[1] = 19; vertexIds[2] = 50; vertexIds[3] = 136; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 22; vertexIds[1] = 19; vertexIds[2] = 54; vertexIds[3] = 137; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 22; vertexIds[1] = 19; vertexIds[2] = 50; vertexIds[3] = 138; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 19; vertexIds[1] = 51; vertexIds[2] = 23; vertexIds[3] = 139; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 51; vertexIds[1] = 23; vertexIds[2] = 54; vertexIds[3] = 140; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 19; vertexIds[1] = 51; vertexIds[2] = 50; vertexIds[3] = 141; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 20; vertexIds[1] = 21; vertexIds[2] = 52; vertexIds[3] = 142; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 24; vertexIds[1] = 21; vertexIds[2] = 56; vertexIds[3] = 143; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 24; vertexIds[1] = 21; vertexIds[2] = 52; vertexIds[3] = 144; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 21; vertexIds[1] = 53; vertexIds[2] = 25; vertexIds[3] = 145; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 53; vertexIds[1] = 25; vertexIds[2] = 56; vertexIds[3] = 146; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 21; vertexIds[1] = 53; vertexIds[2] = 52; vertexIds[3] = 147; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 22; vertexIds[1] = 23; vertexIds[2] = 54; vertexIds[3] = 148; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 26; vertexIds[1] = 23; vertexIds[2] = 58; vertexIds[3] = 149; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 26; vertexIds[1] = 23; vertexIds[2] = 54; vertexIds[3] = 150; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 23; vertexIds[1] = 55; vertexIds[2] = 27; vertexIds[3] = 151; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 55; vertexIds[1] = 27; vertexIds[2] = 58; vertexIds[3] = 152; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 23; vertexIds[1] = 55; vertexIds[2] = 54; vertexIds[3] = 153; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 24; vertexIds[1] = 25; vertexIds[2] = 56; vertexIds[3] = 154; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 28; vertexIds[1] = 25; vertexIds[2] = 60; vertexIds[3] = 155; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 28; vertexIds[1] = 25; vertexIds[2] = 56; vertexIds[3] = 156; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 25; vertexIds[1] = 57; vertexIds[2] = 29; vertexIds[3] = 157; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 57; vertexIds[1] = 29; vertexIds[2] = 60; vertexIds[3] = 158; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 25; vertexIds[1] = 57; vertexIds[2] = 56; vertexIds[3] = 159; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 25; vertexIds[1] = 26; vertexIds[2] = 57; vertexIds[3] = 160; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 29; vertexIds[1] = 26; vertexIds[2] = 61; vertexIds[3] = 161; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 29; vertexIds[1] = 26; vertexIds[2] = 57; vertexIds[3] = 162; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 26; vertexIds[1] = 58; vertexIds[2] = 30; vertexIds[3] = 163; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 58; vertexIds[1] = 30; vertexIds[2] = 61; vertexIds[3] = 164; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 26; vertexIds[1] = 58; vertexIds[2] = 57; vertexIds[3] = 165; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 26; vertexIds[1] = 27; vertexIds[2] = 58; vertexIds[3] = 166; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 30; vertexIds[1] = 27; vertexIds[2] = 62; vertexIds[3] = 167; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 30; vertexIds[1] = 27; vertexIds[2] = 58; vertexIds[3] = 168; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 27; vertexIds[1] = 59; vertexIds[2] = 31; vertexIds[3] = 169; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 59; vertexIds[1] = 31; vertexIds[2] = 62; vertexIds[3] = 170; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 27; vertexIds[1] = 59; vertexIds[2] = 58; vertexIds[3] = 171; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 4; vertexIds[1] = 1; vertexIds[2] = 0; vertexIds[3] = 64; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 4; vertexIds[1] = 32; vertexIds[2] = 0; vertexIds[3] = 64; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 32; vertexIds[1] = 4; vertexIds[2] = 1; vertexIds[3] = 64; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 5; vertexIds[1] = 1; vertexIds[2] = 4; vertexIds[3] = 65; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 5; vertexIds[1] = 36; vertexIds[2] = 4; vertexIds[3] = 65; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 36; vertexIds[1] = 5; vertexIds[2] = 1; vertexIds[3] = 65; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 36; vertexIds[1] = 1; vertexIds[2] = 4; vertexIds[3] = 66; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 36; vertexIds[1] = 32; vertexIds[2] = 4; vertexIds[3] = 66; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 32; vertexIds[1] = 36; vertexIds[2] = 1; vertexIds[3] = 66; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 36; vertexIds[1] = 33; vertexIds[2] = 1; vertexIds[3] = 67; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 36; vertexIds[1] = 5; vertexIds[2] = 1; vertexIds[3] = 67; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 5; vertexIds[1] = 36; vertexIds[2] = 33; vertexIds[3] = 67; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 37; vertexIds[1] = 5; vertexIds[2] = 33; vertexIds[3] = 68; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 37; vertexIds[1] = 36; vertexIds[2] = 33; vertexIds[3] = 68; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 36; vertexIds[1] = 37; vertexIds[2] = 5; vertexIds[3] = 68; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 36; vertexIds[1] = 33; vertexIds[2] = 1; vertexIds[3] = 69; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 36; vertexIds[1] = 32; vertexIds[2] = 1; vertexIds[3] = 69; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 32; vertexIds[1] = 36; vertexIds[2] = 33; vertexIds[3] = 69; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 5; vertexIds[1] = 2; vertexIds[2] = 1; vertexIds[3] = 70; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 5; vertexIds[1] = 33; vertexIds[2] = 1; vertexIds[3] = 70; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 33; vertexIds[1] = 5; vertexIds[2] = 2; vertexIds[3] = 70; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 6; vertexIds[1] = 2; vertexIds[2] = 5; vertexIds[3] = 71; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 6; vertexIds[1] = 37; vertexIds[2] = 5; vertexIds[3] = 71; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 37; vertexIds[1] = 6; vertexIds[2] = 2; vertexIds[3] = 71; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 37; vertexIds[1] = 2; vertexIds[2] = 5; vertexIds[3] = 72; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 37; vertexIds[1] = 33; vertexIds[2] = 5; vertexIds[3] = 72; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 33; vertexIds[1] = 37; vertexIds[2] = 2; vertexIds[3] = 72; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 37; vertexIds[1] = 34; vertexIds[2] = 2; vertexIds[3] = 73; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 37; vertexIds[1] = 6; vertexIds[2] = 2; vertexIds[3] = 73; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 6; vertexIds[1] = 37; vertexIds[2] = 34; vertexIds[3] = 73; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 6; vertexIds[2] = 34; vertexIds[3] = 74; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 37; vertexIds[2] = 34; vertexIds[3] = 74; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 37; vertexIds[1] = 38; vertexIds[2] = 6; vertexIds[3] = 74; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 37; vertexIds[1] = 34; vertexIds[2] = 2; vertexIds[3] = 75; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 37; vertexIds[1] = 33; vertexIds[2] = 2; vertexIds[3] = 75; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 33; vertexIds[1] = 37; vertexIds[2] = 34; vertexIds[3] = 75; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 6; vertexIds[1] = 3; vertexIds[2] = 2; vertexIds[3] = 76; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 6; vertexIds[1] = 34; vertexIds[2] = 2; vertexIds[3] = 76; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 34; vertexIds[1] = 6; vertexIds[2] = 3; vertexIds[3] = 76; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 7; vertexIds[1] = 3; vertexIds[2] = 6; vertexIds[3] = 77; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 7; vertexIds[1] = 38; vertexIds[2] = 6; vertexIds[3] = 77; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 7; vertexIds[2] = 3; vertexIds[3] = 77; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 3; vertexIds[2] = 6; vertexIds[3] = 78; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 34; vertexIds[2] = 6; vertexIds[3] = 78; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 34; vertexIds[1] = 38; vertexIds[2] = 3; vertexIds[3] = 78; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 35; vertexIds[2] = 3; vertexIds[3] = 79; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 7; vertexIds[2] = 3; vertexIds[3] = 79; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 7; vertexIds[1] = 38; vertexIds[2] = 35; vertexIds[3] = 79; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 39; vertexIds[1] = 7; vertexIds[2] = 35; vertexIds[3] = 80; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 39; vertexIds[1] = 38; vertexIds[2] = 35; vertexIds[3] = 80; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 39; vertexIds[2] = 7; vertexIds[3] = 80; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 35; vertexIds[2] = 3; vertexIds[3] = 81; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 34; vertexIds[2] = 3; vertexIds[3] = 81; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 34; vertexIds[1] = 38; vertexIds[2] = 35; vertexIds[3] = 81; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 8; vertexIds[1] = 5; vertexIds[2] = 4; vertexIds[3] = 82; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 8; vertexIds[1] = 36; vertexIds[2] = 4; vertexIds[3] = 82; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 36; vertexIds[1] = 8; vertexIds[2] = 5; vertexIds[3] = 82; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 9; vertexIds[1] = 5; vertexIds[2] = 8; vertexIds[3] = 83; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 9; vertexIds[1] = 40; vertexIds[2] = 8; vertexIds[3] = 83; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 40; vertexIds[1] = 9; vertexIds[2] = 5; vertexIds[3] = 83; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 40; vertexIds[1] = 5; vertexIds[2] = 8; vertexIds[3] = 84; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 40; vertexIds[1] = 36; vertexIds[2] = 8; vertexIds[3] = 84; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 36; vertexIds[1] = 40; vertexIds[2] = 5; vertexIds[3] = 84; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 40; vertexIds[1] = 37; vertexIds[2] = 5; vertexIds[3] = 85; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 40; vertexIds[1] = 9; vertexIds[2] = 5; vertexIds[3] = 85; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 9; vertexIds[1] = 40; vertexIds[2] = 37; vertexIds[3] = 85; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 41; vertexIds[1] = 9; vertexIds[2] = 37; vertexIds[3] = 86; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 41; vertexIds[1] = 40; vertexIds[2] = 37; vertexIds[3] = 86; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 40; vertexIds[1] = 41; vertexIds[2] = 9; vertexIds[3] = 86; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 40; vertexIds[1] = 37; vertexIds[2] = 5; vertexIds[3] = 87; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 40; vertexIds[1] = 36; vertexIds[2] = 5; vertexIds[3] = 87; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 36; vertexIds[1] = 40; vertexIds[2] = 37; vertexIds[3] = 87; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 10; vertexIds[1] = 7; vertexIds[2] = 6; vertexIds[3] = 88; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 10; vertexIds[1] = 38; vertexIds[2] = 6; vertexIds[3] = 88; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 10; vertexIds[2] = 7; vertexIds[3] = 88; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 11; vertexIds[1] = 7; vertexIds[2] = 10; vertexIds[3] = 89; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 11; vertexIds[1] = 42; vertexIds[2] = 10; vertexIds[3] = 89; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 42; vertexIds[1] = 11; vertexIds[2] = 7; vertexIds[3] = 89; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 42; vertexIds[1] = 7; vertexIds[2] = 10; vertexIds[3] = 90; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 42; vertexIds[1] = 38; vertexIds[2] = 10; vertexIds[3] = 90; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 42; vertexIds[2] = 7; vertexIds[3] = 90; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 42; vertexIds[1] = 39; vertexIds[2] = 7; vertexIds[3] = 91; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 42; vertexIds[1] = 11; vertexIds[2] = 7; vertexIds[3] = 91; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 11; vertexIds[1] = 42; vertexIds[2] = 39; vertexIds[3] = 91; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 43; vertexIds[1] = 11; vertexIds[2] = 39; vertexIds[3] = 92; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 43; vertexIds[1] = 42; vertexIds[2] = 39; vertexIds[3] = 92; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 42; vertexIds[1] = 43; vertexIds[2] = 11; vertexIds[3] = 92; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 42; vertexIds[1] = 39; vertexIds[2] = 7; vertexIds[3] = 93; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 42; vertexIds[1] = 38; vertexIds[2] = 7; vertexIds[3] = 93; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 38; vertexIds[1] = 42; vertexIds[2] = 39; vertexIds[3] = 93; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 12; vertexIds[1] = 9; vertexIds[2] = 8; vertexIds[3] = 94; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 12; vertexIds[1] = 40; vertexIds[2] = 8; vertexIds[3] = 94; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 40; vertexIds[1] = 12; vertexIds[2] = 9; vertexIds[3] = 94; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 13; vertexIds[1] = 9; vertexIds[2] = 12; vertexIds[3] = 95; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 13; vertexIds[1] = 44; vertexIds[2] = 12; vertexIds[3] = 95; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 44; vertexIds[1] = 13; vertexIds[2] = 9; vertexIds[3] = 95; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 44; vertexIds[1] = 9; vertexIds[2] = 12; vertexIds[3] = 96; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 44; vertexIds[1] = 40; vertexIds[2] = 12; vertexIds[3] = 96; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 40; vertexIds[1] = 44; vertexIds[2] = 9; vertexIds[3] = 96; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 44; vertexIds[1] = 41; vertexIds[2] = 9; vertexIds[3] = 97; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 44; vertexIds[1] = 13; vertexIds[2] = 9; vertexIds[3] = 97; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 13; vertexIds[1] = 44; vertexIds[2] = 41; vertexIds[3] = 97; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 45; vertexIds[1] = 13; vertexIds[2] = 41; vertexIds[3] = 98; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 45; vertexIds[1] = 44; vertexIds[2] = 41; vertexIds[3] = 98; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 44; vertexIds[1] = 45; vertexIds[2] = 13; vertexIds[3] = 98; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 44; vertexIds[1] = 41; vertexIds[2] = 9; vertexIds[3] = 99; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 44; vertexIds[1] = 40; vertexIds[2] = 9; vertexIds[3] = 99; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 40; vertexIds[1] = 44; vertexIds[2] = 41; vertexIds[3] = 99; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 13; vertexIds[1] = 10; vertexIds[2] = 9; vertexIds[3] = 100; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 13; vertexIds[1] = 41; vertexIds[2] = 9; vertexIds[3] = 100; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 41; vertexIds[1] = 13; vertexIds[2] = 10; vertexIds[3] = 100; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 14; vertexIds[1] = 10; vertexIds[2] = 13; vertexIds[3] = 101; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 14; vertexIds[1] = 45; vertexIds[2] = 13; vertexIds[3] = 101; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 45; vertexIds[1] = 14; vertexIds[2] = 10; vertexIds[3] = 101; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 45; vertexIds[1] = 10; vertexIds[2] = 13; vertexIds[3] = 102; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 45; vertexIds[1] = 41; vertexIds[2] = 13; vertexIds[3] = 102; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 41; vertexIds[1] = 45; vertexIds[2] = 10; vertexIds[3] = 102; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 45; vertexIds[1] = 42; vertexIds[2] = 10; vertexIds[3] = 103; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 45; vertexIds[1] = 14; vertexIds[2] = 10; vertexIds[3] = 103; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 14; vertexIds[1] = 45; vertexIds[2] = 42; vertexIds[3] = 103; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 14; vertexIds[2] = 42; vertexIds[3] = 104; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 45; vertexIds[2] = 42; vertexIds[3] = 104; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 45; vertexIds[1] = 46; vertexIds[2] = 14; vertexIds[3] = 104; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 45; vertexIds[1] = 42; vertexIds[2] = 10; vertexIds[3] = 105; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 45; vertexIds[1] = 41; vertexIds[2] = 10; vertexIds[3] = 105; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 41; vertexIds[1] = 45; vertexIds[2] = 42; vertexIds[3] = 105; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 14; vertexIds[1] = 11; vertexIds[2] = 10; vertexIds[3] = 106; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 14; vertexIds[1] = 42; vertexIds[2] = 10; vertexIds[3] = 106; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 42; vertexIds[1] = 14; vertexIds[2] = 11; vertexIds[3] = 106; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 15; vertexIds[1] = 11; vertexIds[2] = 14; vertexIds[3] = 107; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 15; vertexIds[1] = 46; vertexIds[2] = 14; vertexIds[3] = 107; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 15; vertexIds[2] = 11; vertexIds[3] = 107; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 11; vertexIds[2] = 14; vertexIds[3] = 108; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 42; vertexIds[2] = 14; vertexIds[3] = 108; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 42; vertexIds[1] = 46; vertexIds[2] = 11; vertexIds[3] = 108; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 43; vertexIds[2] = 11; vertexIds[3] = 109; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 15; vertexIds[2] = 11; vertexIds[3] = 109; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 15; vertexIds[1] = 46; vertexIds[2] = 43; vertexIds[3] = 109; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 47; vertexIds[1] = 15; vertexIds[2] = 43; vertexIds[3] = 110; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 47; vertexIds[1] = 46; vertexIds[2] = 43; vertexIds[3] = 110; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 47; vertexIds[2] = 15; vertexIds[3] = 110; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 43; vertexIds[2] = 11; vertexIds[3] = 111; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 42; vertexIds[2] = 11; vertexIds[3] = 111; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 42; vertexIds[1] = 46; vertexIds[2] = 43; vertexIds[3] = 111; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 16; vertexIds[1] = 13; vertexIds[2] = 12; vertexIds[3] = 112; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 16; vertexIds[1] = 44; vertexIds[2] = 12; vertexIds[3] = 112; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 44; vertexIds[1] = 16; vertexIds[2] = 13; vertexIds[3] = 112; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 17; vertexIds[1] = 13; vertexIds[2] = 16; vertexIds[3] = 113; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 17; vertexIds[1] = 48; vertexIds[2] = 16; vertexIds[3] = 113; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 48; vertexIds[1] = 17; vertexIds[2] = 13; vertexIds[3] = 113; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 48; vertexIds[1] = 13; vertexIds[2] = 16; vertexIds[3] = 114; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 48; vertexIds[1] = 44; vertexIds[2] = 16; vertexIds[3] = 114; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 44; vertexIds[1] = 48; vertexIds[2] = 13; vertexIds[3] = 114; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 48; vertexIds[1] = 45; vertexIds[2] = 13; vertexIds[3] = 115; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 48; vertexIds[1] = 17; vertexIds[2] = 13; vertexIds[3] = 115; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 17; vertexIds[1] = 48; vertexIds[2] = 45; vertexIds[3] = 115; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 49; vertexIds[1] = 17; vertexIds[2] = 45; vertexIds[3] = 116; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 49; vertexIds[1] = 48; vertexIds[2] = 45; vertexIds[3] = 116; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 48; vertexIds[1] = 49; vertexIds[2] = 17; vertexIds[3] = 116; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 48; vertexIds[1] = 45; vertexIds[2] = 13; vertexIds[3] = 117; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 48; vertexIds[1] = 44; vertexIds[2] = 13; vertexIds[3] = 117; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 44; vertexIds[1] = 48; vertexIds[2] = 45; vertexIds[3] = 117; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 18; vertexIds[1] = 15; vertexIds[2] = 14; vertexIds[3] = 118; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 18; vertexIds[1] = 46; vertexIds[2] = 14; vertexIds[3] = 118; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 18; vertexIds[2] = 15; vertexIds[3] = 118; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 19; vertexIds[1] = 15; vertexIds[2] = 18; vertexIds[3] = 119; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 19; vertexIds[1] = 50; vertexIds[2] = 18; vertexIds[3] = 119; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 50; vertexIds[1] = 19; vertexIds[2] = 15; vertexIds[3] = 119; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 50; vertexIds[1] = 15; vertexIds[2] = 18; vertexIds[3] = 120; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 50; vertexIds[1] = 46; vertexIds[2] = 18; vertexIds[3] = 120; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 50; vertexIds[2] = 15; vertexIds[3] = 120; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 50; vertexIds[1] = 47; vertexIds[2] = 15; vertexIds[3] = 121; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 50; vertexIds[1] = 19; vertexIds[2] = 15; vertexIds[3] = 121; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 19; vertexIds[1] = 50; vertexIds[2] = 47; vertexIds[3] = 121; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 51; vertexIds[1] = 19; vertexIds[2] = 47; vertexIds[3] = 122; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 51; vertexIds[1] = 50; vertexIds[2] = 47; vertexIds[3] = 122; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 50; vertexIds[1] = 51; vertexIds[2] = 19; vertexIds[3] = 122; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 50; vertexIds[1] = 47; vertexIds[2] = 15; vertexIds[3] = 123; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 50; vertexIds[1] = 46; vertexIds[2] = 15; vertexIds[3] = 123; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 46; vertexIds[1] = 50; vertexIds[2] = 47; vertexIds[3] = 123; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 20; vertexIds[1] = 17; vertexIds[2] = 16; vertexIds[3] = 124; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 20; vertexIds[1] = 48; vertexIds[2] = 16; vertexIds[3] = 124; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 48; vertexIds[1] = 20; vertexIds[2] = 17; vertexIds[3] = 124; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 21; vertexIds[1] = 17; vertexIds[2] = 20; vertexIds[3] = 125; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 21; vertexIds[1] = 52; vertexIds[2] = 20; vertexIds[3] = 125; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 52; vertexIds[1] = 21; vertexIds[2] = 17; vertexIds[3] = 125; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 52; vertexIds[1] = 17; vertexIds[2] = 20; vertexIds[3] = 126; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 52; vertexIds[1] = 48; vertexIds[2] = 20; vertexIds[3] = 126; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 48; vertexIds[1] = 52; vertexIds[2] = 17; vertexIds[3] = 126; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 52; vertexIds[1] = 49; vertexIds[2] = 17; vertexIds[3] = 127; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 52; vertexIds[1] = 21; vertexIds[2] = 17; vertexIds[3] = 127; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 21; vertexIds[1] = 52; vertexIds[2] = 49; vertexIds[3] = 127; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 53; vertexIds[1] = 21; vertexIds[2] = 49; vertexIds[3] = 128; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 53; vertexIds[1] = 52; vertexIds[2] = 49; vertexIds[3] = 128; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 52; vertexIds[1] = 53; vertexIds[2] = 21; vertexIds[3] = 128; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 52; vertexIds[1] = 49; vertexIds[2] = 17; vertexIds[3] = 129; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 52; vertexIds[1] = 48; vertexIds[2] = 17; vertexIds[3] = 129; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 48; vertexIds[1] = 52; vertexIds[2] = 49; vertexIds[3] = 129; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 21; vertexIds[1] = 18; vertexIds[2] = 17; vertexIds[3] = 130; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 21; vertexIds[1] = 49; vertexIds[2] = 17; vertexIds[3] = 130; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 49; vertexIds[1] = 21; vertexIds[2] = 18; vertexIds[3] = 130; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 22; vertexIds[1] = 18; vertexIds[2] = 21; vertexIds[3] = 131; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 22; vertexIds[1] = 53; vertexIds[2] = 21; vertexIds[3] = 131; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 53; vertexIds[1] = 22; vertexIds[2] = 18; vertexIds[3] = 131; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 53; vertexIds[1] = 18; vertexIds[2] = 21; vertexIds[3] = 132; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 53; vertexIds[1] = 49; vertexIds[2] = 21; vertexIds[3] = 132; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 49; vertexIds[1] = 53; vertexIds[2] = 18; vertexIds[3] = 132; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 53; vertexIds[1] = 50; vertexIds[2] = 18; vertexIds[3] = 133; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 53; vertexIds[1] = 22; vertexIds[2] = 18; vertexIds[3] = 133; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 22; vertexIds[1] = 53; vertexIds[2] = 50; vertexIds[3] = 133; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 22; vertexIds[2] = 50; vertexIds[3] = 134; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 53; vertexIds[2] = 50; vertexIds[3] = 134; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 53; vertexIds[1] = 54; vertexIds[2] = 22; vertexIds[3] = 134; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 53; vertexIds[1] = 50; vertexIds[2] = 18; vertexIds[3] = 135; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 53; vertexIds[1] = 49; vertexIds[2] = 18; vertexIds[3] = 135; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 49; vertexIds[1] = 53; vertexIds[2] = 50; vertexIds[3] = 135; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 22; vertexIds[1] = 19; vertexIds[2] = 18; vertexIds[3] = 136; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 22; vertexIds[1] = 50; vertexIds[2] = 18; vertexIds[3] = 136; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 50; vertexIds[1] = 22; vertexIds[2] = 19; vertexIds[3] = 136; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 23; vertexIds[1] = 19; vertexIds[2] = 22; vertexIds[3] = 137; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 23; vertexIds[1] = 54; vertexIds[2] = 22; vertexIds[3] = 137; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 23; vertexIds[2] = 19; vertexIds[3] = 137; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 19; vertexIds[2] = 22; vertexIds[3] = 138; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 50; vertexIds[2] = 22; vertexIds[3] = 138; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 50; vertexIds[1] = 54; vertexIds[2] = 19; vertexIds[3] = 138; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 51; vertexIds[2] = 19; vertexIds[3] = 139; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 23; vertexIds[2] = 19; vertexIds[3] = 139; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 23; vertexIds[1] = 54; vertexIds[2] = 51; vertexIds[3] = 139; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 55; vertexIds[1] = 23; vertexIds[2] = 51; vertexIds[3] = 140; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 55; vertexIds[1] = 54; vertexIds[2] = 51; vertexIds[3] = 140; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 55; vertexIds[2] = 23; vertexIds[3] = 140; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 51; vertexIds[2] = 19; vertexIds[3] = 141; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 50; vertexIds[2] = 19; vertexIds[3] = 141; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 50; vertexIds[1] = 54; vertexIds[2] = 51; vertexIds[3] = 141; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 24; vertexIds[1] = 21; vertexIds[2] = 20; vertexIds[3] = 142; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 24; vertexIds[1] = 52; vertexIds[2] = 20; vertexIds[3] = 142; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 52; vertexIds[1] = 24; vertexIds[2] = 21; vertexIds[3] = 142; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 25; vertexIds[1] = 21; vertexIds[2] = 24; vertexIds[3] = 143; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 25; vertexIds[1] = 56; vertexIds[2] = 24; vertexIds[3] = 143; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 56; vertexIds[1] = 25; vertexIds[2] = 21; vertexIds[3] = 143; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 56; vertexIds[1] = 21; vertexIds[2] = 24; vertexIds[3] = 144; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 56; vertexIds[1] = 52; vertexIds[2] = 24; vertexIds[3] = 144; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 52; vertexIds[1] = 56; vertexIds[2] = 21; vertexIds[3] = 144; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 56; vertexIds[1] = 53; vertexIds[2] = 21; vertexIds[3] = 145; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 56; vertexIds[1] = 25; vertexIds[2] = 21; vertexIds[3] = 145; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 25; vertexIds[1] = 56; vertexIds[2] = 53; vertexIds[3] = 145; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 57; vertexIds[1] = 25; vertexIds[2] = 53; vertexIds[3] = 146; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 57; vertexIds[1] = 56; vertexIds[2] = 53; vertexIds[3] = 146; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 56; vertexIds[1] = 57; vertexIds[2] = 25; vertexIds[3] = 146; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 56; vertexIds[1] = 53; vertexIds[2] = 21; vertexIds[3] = 147; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 56; vertexIds[1] = 52; vertexIds[2] = 21; vertexIds[3] = 147; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 52; vertexIds[1] = 56; vertexIds[2] = 53; vertexIds[3] = 147; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 26; vertexIds[1] = 23; vertexIds[2] = 22; vertexIds[3] = 148; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 26; vertexIds[1] = 54; vertexIds[2] = 22; vertexIds[3] = 148; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 26; vertexIds[2] = 23; vertexIds[3] = 148; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 27; vertexIds[1] = 23; vertexIds[2] = 26; vertexIds[3] = 149; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 27; vertexIds[1] = 58; vertexIds[2] = 26; vertexIds[3] = 149; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 58; vertexIds[1] = 27; vertexIds[2] = 23; vertexIds[3] = 149; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 58; vertexIds[1] = 23; vertexIds[2] = 26; vertexIds[3] = 150; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 58; vertexIds[1] = 54; vertexIds[2] = 26; vertexIds[3] = 150; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 58; vertexIds[2] = 23; vertexIds[3] = 150; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 58; vertexIds[1] = 55; vertexIds[2] = 23; vertexIds[3] = 151; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 58; vertexIds[1] = 27; vertexIds[2] = 23; vertexIds[3] = 151; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 27; vertexIds[1] = 58; vertexIds[2] = 55; vertexIds[3] = 151; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 59; vertexIds[1] = 27; vertexIds[2] = 55; vertexIds[3] = 152; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 59; vertexIds[1] = 58; vertexIds[2] = 55; vertexIds[3] = 152; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 58; vertexIds[1] = 59; vertexIds[2] = 27; vertexIds[3] = 152; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 58; vertexIds[1] = 55; vertexIds[2] = 23; vertexIds[3] = 153; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 58; vertexIds[1] = 54; vertexIds[2] = 23; vertexIds[3] = 153; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 54; vertexIds[1] = 58; vertexIds[2] = 55; vertexIds[3] = 153; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 28; vertexIds[1] = 25; vertexIds[2] = 24; vertexIds[3] = 154; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 28; vertexIds[1] = 56; vertexIds[2] = 24; vertexIds[3] = 154; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 56; vertexIds[1] = 28; vertexIds[2] = 25; vertexIds[3] = 154; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 29; vertexIds[1] = 25; vertexIds[2] = 28; vertexIds[3] = 155; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 29; vertexIds[1] = 60; vertexIds[2] = 28; vertexIds[3] = 155; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 60; vertexIds[1] = 29; vertexIds[2] = 25; vertexIds[3] = 155; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 60; vertexIds[1] = 25; vertexIds[2] = 28; vertexIds[3] = 156; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 60; vertexIds[1] = 56; vertexIds[2] = 28; vertexIds[3] = 156; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 56; vertexIds[1] = 60; vertexIds[2] = 25; vertexIds[3] = 156; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 60; vertexIds[1] = 57; vertexIds[2] = 25; vertexIds[3] = 157; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 60; vertexIds[1] = 29; vertexIds[2] = 25; vertexIds[3] = 157; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 29; vertexIds[1] = 60; vertexIds[2] = 57; vertexIds[3] = 157; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 61; vertexIds[1] = 29; vertexIds[2] = 57; vertexIds[3] = 158; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 61; vertexIds[1] = 60; vertexIds[2] = 57; vertexIds[3] = 158; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 60; vertexIds[1] = 61; vertexIds[2] = 29; vertexIds[3] = 158; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 60; vertexIds[1] = 57; vertexIds[2] = 25; vertexIds[3] = 159; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 60; vertexIds[1] = 56; vertexIds[2] = 25; vertexIds[3] = 159; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 56; vertexIds[1] = 60; vertexIds[2] = 57; vertexIds[3] = 159; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 29; vertexIds[1] = 26; vertexIds[2] = 25; vertexIds[3] = 160; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 29; vertexIds[1] = 57; vertexIds[2] = 25; vertexIds[3] = 160; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 57; vertexIds[1] = 29; vertexIds[2] = 26; vertexIds[3] = 160; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 30; vertexIds[1] = 26; vertexIds[2] = 29; vertexIds[3] = 161; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 30; vertexIds[1] = 61; vertexIds[2] = 29; vertexIds[3] = 161; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 61; vertexIds[1] = 30; vertexIds[2] = 26; vertexIds[3] = 161; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 61; vertexIds[1] = 26; vertexIds[2] = 29; vertexIds[3] = 162; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 61; vertexIds[1] = 57; vertexIds[2] = 29; vertexIds[3] = 162; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 57; vertexIds[1] = 61; vertexIds[2] = 26; vertexIds[3] = 162; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 61; vertexIds[1] = 58; vertexIds[2] = 26; vertexIds[3] = 163; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 61; vertexIds[1] = 30; vertexIds[2] = 26; vertexIds[3] = 163; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 30; vertexIds[1] = 61; vertexIds[2] = 58; vertexIds[3] = 163; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 62; vertexIds[1] = 30; vertexIds[2] = 58; vertexIds[3] = 164; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 62; vertexIds[1] = 61; vertexIds[2] = 58; vertexIds[3] = 164; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 61; vertexIds[1] = 62; vertexIds[2] = 30; vertexIds[3] = 164; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 61; vertexIds[1] = 58; vertexIds[2] = 26; vertexIds[3] = 165; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 61; vertexIds[1] = 57; vertexIds[2] = 26; vertexIds[3] = 165; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 57; vertexIds[1] = 61; vertexIds[2] = 58; vertexIds[3] = 165; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 30; vertexIds[1] = 27; vertexIds[2] = 26; vertexIds[3] = 166; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 30; vertexIds[1] = 58; vertexIds[2] = 26; vertexIds[3] = 166; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 58; vertexIds[1] = 30; vertexIds[2] = 27; vertexIds[3] = 166; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 31; vertexIds[1] = 27; vertexIds[2] = 30; vertexIds[3] = 167; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 31; vertexIds[1] = 62; vertexIds[2] = 30; vertexIds[3] = 167; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 62; vertexIds[1] = 31; vertexIds[2] = 27; vertexIds[3] = 167; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 62; vertexIds[1] = 27; vertexIds[2] = 30; vertexIds[3] = 168; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 62; vertexIds[1] = 58; vertexIds[2] = 30; vertexIds[3] = 168; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 58; vertexIds[1] = 62; vertexIds[2] = 27; vertexIds[3] = 168; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 62; vertexIds[1] = 59; vertexIds[2] = 27; vertexIds[3] = 169; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 62; vertexIds[1] = 31; vertexIds[2] = 27; vertexIds[3] = 169; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 31; vertexIds[1] = 62; vertexIds[2] = 59; vertexIds[3] = 169; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 63; vertexIds[1] = 31; vertexIds[2] = 59; vertexIds[3] = 170; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 63; vertexIds[1] = 62; vertexIds[2] = 59; vertexIds[3] = 170; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 62; vertexIds[1] = 63; vertexIds[2] = 31; vertexIds[3] = 170; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 62; vertexIds[1] = 59; vertexIds[2] = 27; vertexIds[3] = 171; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 62; vertexIds[1] = 58; vertexIds[2] = 27; vertexIds[3] = 171; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); vertexIds[0] = 58; vertexIds[1] = 62; vertexIds[2] = 59; vertexIds[3] = 171; vMesh->InsertNextCell(VTK_TETRA, 4, vertexIds); free(vertexIds); return 1; } // loading code generated automatically. // put whatever mesh you like here. int LoadSurfaceMesh(vtkPolyData *sMesh) { sMesh->Allocate(); vtkPoints *points = vtkPoints::New(); points->InsertNextPoint( -0.774167, 0.632981, 0); points->InsertNextPoint( -0.141186, 1.40715, 0); points->InsertNextPoint( 0.632981, 0.774167, 0); points->InsertNextPoint( 0, 0, 0); points->InsertNextPoint( -0.774167, 0.632981, 1); points->InsertNextPoint( -0.141186, 1.40715, 1); points->InsertNextPoint( 0.632981, 0.774167, 1); points->InsertNextPoint( 0, 0, 1); points->InsertNextPoint( 0.491796, 2.18132, 0); points->InsertNextPoint( 1.26596, 1.54833, 0); points->InsertNextPoint( 0.491796, 2.18132, 1); points->InsertNextPoint( 1.26596, 1.54833, 1); points->InsertNextPoint( 1.12478, 2.95548, 0); points->InsertNextPoint( 1.89894, 2.3225, 0); points->InsertNextPoint( 1.12478, 2.95548, 1); points->InsertNextPoint( 1.89894, 2.3225, 1); points->InsertNextPoint( -1.54833, 1.26596, 0); points->InsertNextPoint( -0.915353, 2.04013, 0); points->InsertNextPoint( -1.54833, 1.26596, 1); points->InsertNextPoint( -0.915353, 2.04013, 1); points->InsertNextPoint( -0.282372, 2.8143, 0); points->InsertNextPoint( 0.35061, 3.58846, 0); points->InsertNextPoint( -0.282372, 2.8143, 1); points->InsertNextPoint( 0.35061, 3.58846, 1); points->InsertNextPoint( -2.3225, 1.89894, 0); points->InsertNextPoint( -1.68952, 2.67311, 0); points->InsertNextPoint( -2.3225, 1.89894, 1); points->InsertNextPoint( -1.68952, 2.67311, 1); points->InsertNextPoint( -1.05654, 3.44728, 0); points->InsertNextPoint( -1.05654, 3.44728, 1); points->InsertNextPoint( -0.423557, 4.22145, 0); points->InsertNextPoint( -0.423557, 4.22145, 1); points->InsertNextPoint( -3.09667, 2.53193, 0); points->InsertNextPoint( -2.46369, 3.30609, 0); points->InsertNextPoint( -3.09667, 2.53193, 1); points->InsertNextPoint( -2.46369, 3.30609, 1); points->InsertNextPoint( -1.83071, 4.08026, 0); points->InsertNextPoint( -1.19772, 4.85443, 0); points->InsertNextPoint( -1.83071, 4.08026, 1); points->InsertNextPoint( -1.19772, 4.85443, 1); points->InsertNextPoint( -3.87084, 3.16491, 0); points->InsertNextPoint( -3.23785, 3.93907, 0); points->InsertNextPoint( -3.87084, 3.16491, 1); points->InsertNextPoint( -3.23785, 3.93907, 1); points->InsertNextPoint( -2.60487, 4.71324, 0); points->InsertNextPoint( -2.60487, 4.71324, 1); points->InsertNextPoint( -1.97189, 5.48741, 0); points->InsertNextPoint( -1.97189, 5.48741, 1); points->InsertNextPoint( -4.645, 3.79789, 0); points->InsertNextPoint( -4.01202, 4.57205, 0); points->InsertNextPoint( -4.645, 3.79789, 1); points->InsertNextPoint( -4.01202, 4.57205, 1); points->InsertNextPoint( -3.37904, 5.34622, 0); points->InsertNextPoint( -2.74606, 6.12039, 0); points->InsertNextPoint( -3.37904, 5.34622, 1); points->InsertNextPoint( -2.74606, 6.12039, 1); points->InsertNextPoint( -5.41917, 4.43087, 0); points->InsertNextPoint( -4.78619, 5.20504, 0); points->InsertNextPoint( -5.41917, 4.43087, 1); points->InsertNextPoint( -4.78619, 5.20504, 1); points->InsertNextPoint( -4.15321, 5.9792, 0); points->InsertNextPoint( -4.15321, 5.9792, 1); points->InsertNextPoint( -3.52023, 6.75337, 0); points->InsertNextPoint( -3.52023, 6.75337, 1); points->InsertNextPoint( -0.0941238, 0.938099, 0); points->InsertNextPoint( -0.0470619, 0.469049, 0); points->InsertNextPoint( -0.0941238, 0.938099, 1); points->InsertNextPoint( -0.0470619, 0.469049, 1); points->InsertNextPoint( -0.516111, 0.421988, 0.333333); points->InsertNextPoint( -0.258056, 0.210994, 0.666667); points->InsertNextPoint( 0.210994, 0.258056, 0.333333); points->InsertNextPoint( 0.421988, 0.516111, 0.666667); points->InsertNextPoint( 0.538857, 1.71227, 0); points->InsertNextPoint( 0.585919, 1.24322, 0); points->InsertNextPoint( 0.538857, 1.71227, 1); points->InsertNextPoint( 0.585919, 1.24322, 1); points->InsertNextPoint( 0.069808, 1.6652, 0.666667); points->InsertNextPoint( 0.280802, 1.92326, 0.333333); points->InsertNextPoint( 0.843975, 1.03222, 0.333333); points->InsertNextPoint( 1.05497, 1.29028, 0.666667); points->InsertNextPoint( 1.17184, 2.48643, 0); points->InsertNextPoint( 1.2189, 2.01738, 0); points->InsertNextPoint( 1.17184, 2.48643, 1); points->InsertNextPoint( 1.2189, 2.01738, 1); points->InsertNextPoint( 1.47696, 1.80639, 0.333333); points->InsertNextPoint( 1.68795, 2.06445, 0.666667); points->InsertNextPoint( 1.38283, 2.74449, 0.666667); points->InsertNextPoint( 1.64089, 2.5335, 0.333333); points->InsertNextPoint( -0.868291, 1.57108, 0); points->InsertNextPoint( -0.821229, 1.10203, 0); points->InsertNextPoint( -0.868291, 1.57108, 1); points->InsertNextPoint( -0.821229, 1.10203, 1); points->InsertNextPoint( -1.29028, 1.05497, 0.333333); points->InsertNextPoint( -1.03222, 0.843975, 0.666667); points->InsertNextPoint( -0.657297, 1.82914, 0.666667); points->InsertNextPoint( -0.399241, 1.61814, 0.333333); points->InsertNextPoint( 0.397672, 3.11941, 0); points->InsertNextPoint( 0.444734, 2.65036, 0); points->InsertNextPoint( 0.397672, 3.11941, 1); points->InsertNextPoint( 0.444734, 2.65036, 1); points->InsertNextPoint( -0.0243159, 2.6033, 0.333333); points->InsertNextPoint( 0.23374, 2.39231, 0.666667); points->InsertNextPoint( 0.608665, 3.37747, 0.666667); points->InsertNextPoint( 0.866721, 3.16648, 0.333333); points->InsertNextPoint( -1.64246, 2.20406, 0); points->InsertNextPoint( -1.5954, 1.73501, 0); points->InsertNextPoint( -1.64246, 2.20406, 1); points->InsertNextPoint( -1.5954, 1.73501, 1); points->InsertNextPoint( -2.06445, 1.68795, 0.333333); points->InsertNextPoint( -1.80639, 1.47696, 0.666667); points->InsertNextPoint( -1.00948, 2.97823, 0); points->InsertNextPoint( -0.962415, 2.50918, 0); points->InsertNextPoint( -1.00948, 2.97823, 1); points->InsertNextPoint( -0.962415, 2.50918, 1); points->InsertNextPoint( -1.47853, 2.93117, 0.666667); points->InsertNextPoint( -1.26753, 3.18922, 0.333333); points->InsertNextPoint( -0.704359, 2.29819, 0.333333); points->InsertNextPoint( -0.493365, 2.55624, 0.666667); points->InsertNextPoint( -0.376495, 3.7524, 0); points->InsertNextPoint( -0.329433, 3.28335, 0); points->InsertNextPoint( -0.376495, 3.7524, 1); points->InsertNextPoint( -0.329433, 3.28335, 1); points->InsertNextPoint( -0.165502, 4.01045, 0.666667); points->InsertNextPoint( 0.0925541, 3.79946, 0.333333); points->InsertNextPoint( -2.41663, 2.83704, 0); points->InsertNextPoint( -2.36956, 2.36799, 0); points->InsertNextPoint( -2.41663, 2.83704, 1); points->InsertNextPoint( -2.36956, 2.36799, 1); points->InsertNextPoint( -2.83861, 2.32093, 0.333333); points->InsertNextPoint( -2.58056, 2.10994, 0.666667); points->InsertNextPoint( -2.20563, 3.0951, 0.666667); points->InsertNextPoint( -1.94758, 2.8841, 0.333333); points->InsertNextPoint( -1.15066, 4.38538, 0); points->InsertNextPoint( -1.1036, 3.91633, 0); points->InsertNextPoint( -1.15066, 4.38538, 1); points->InsertNextPoint( -1.1036, 3.91633, 1); points->InsertNextPoint( -1.57265, 3.86927, 0.333333); points->InsertNextPoint( -1.31459, 3.65827, 0.666667); points->InsertNextPoint( -0.939669, 4.64343, 0.666667); points->InsertNextPoint( -0.681613, 4.43244, 0.333333); points->InsertNextPoint( -3.19079, 3.47002, 0); points->InsertNextPoint( -3.14373, 3.00097, 0); points->InsertNextPoint( -3.19079, 3.47002, 1); points->InsertNextPoint( -3.14373, 3.00097, 1); points->InsertNextPoint( -3.61278, 2.95391, 0.333333); points->InsertNextPoint( -3.35472, 2.74292, 0.666667); points->InsertNextPoint( -2.55781, 4.24419, 0); points->InsertNextPoint( -2.51075, 3.77514, 0); points->InsertNextPoint( -2.55781, 4.24419, 1); points->InsertNextPoint( -2.51075, 3.77514, 1); points->InsertNextPoint( -3.02686, 4.19713, 0.666667); points->InsertNextPoint( -2.81587, 4.45518, 0.333333); points->InsertNextPoint( -2.25269, 3.56415, 0.333333); points->InsertNextPoint( -2.0417, 3.8222, 0.666667); points->InsertNextPoint( -1.92483, 5.01836, 0); points->InsertNextPoint( -1.87777, 4.54931, 0); points->InsertNextPoint( -1.92483, 5.01836, 1); points->InsertNextPoint( -1.87777, 4.54931, 1); points->InsertNextPoint( -1.71384, 5.27641, 0.666667); points->InsertNextPoint( -1.45578, 5.06542, 0.333333); points->InsertNextPoint( -3.96496, 4.10301, 0); points->InsertNextPoint( -3.9179, 3.63396, 0); points->InsertNextPoint( -3.96496, 4.10301, 1); points->InsertNextPoint( -3.9179, 3.63396, 1); points->InsertNextPoint( -4.38695, 3.58689, 0.333333); points->InsertNextPoint( -4.12889, 3.3759, 0.666667); points->InsertNextPoint( -3.75397, 4.36106, 0.666667); points->InsertNextPoint( -3.49591, 4.15007, 0.333333); points->InsertNextPoint( -2.699, 5.65134, 0); points->InsertNextPoint( -2.65193, 5.18229, 0); points->InsertNextPoint( -2.699, 5.65134, 1); points->InsertNextPoint( -2.65193, 5.18229, 1); points->InsertNextPoint( -3.12098, 5.13523, 0.333333); points->InsertNextPoint( -2.86293, 4.92423, 0.666667); points->InsertNextPoint( -2.488, 5.9094, 0.666667); points->InsertNextPoint( -2.22995, 5.6984, 0.333333); points->InsertNextPoint( -4.73913, 4.73599, 0); points->InsertNextPoint( -4.69206, 4.26694, 0); points->InsertNextPoint( -4.73913, 4.73599, 1); points->InsertNextPoint( -4.69206, 4.26694, 1); points->InsertNextPoint( -5.20818, 4.68892, 0.666667); points->InsertNextPoint( -4.99718, 4.94698, 0.333333); points->InsertNextPoint( -5.16111, 4.21988, 0.333333); points->InsertNextPoint( -4.90306, 4.00888, 0.666667); points->InsertNextPoint( -4.10615, 5.51015, 0); points->InsertNextPoint( -4.05908, 5.0411, 0); points->InsertNextPoint( -4.10615, 5.51015, 1); points->InsertNextPoint( -4.05908, 5.0411, 1); points->InsertNextPoint( -4.57519, 5.46309, 0.666667); points->InsertNextPoint( -4.3642, 5.72115, 0.333333); points->InsertNextPoint( -3.80103, 4.83011, 0.333333); points->InsertNextPoint( -3.59003, 5.08817, 0.666667); points->InsertNextPoint( -3.47316, 6.28432, 0); points->InsertNextPoint( -3.4261, 5.81527, 0); points->InsertNextPoint( -3.47316, 6.28432, 1); points->InsertNextPoint( -3.4261, 5.81527, 1); points->InsertNextPoint( -3.94221, 6.23726, 0.666667); points->InsertNextPoint( -3.73122, 6.49531, 0.333333); points->InsertNextPoint( -3.26217, 6.54238, 0.666667); points->InsertNextPoint( -3.00411, 6.33138, 0.333333); sMesh->SetPoints(points); points->Delete(); vtkIdType *vertexIds = (vtkIdType *) malloc(sizeof(vtkIdType)*3); vertexIds[0] = 0; vertexIds[1] = 1; vertexIds[2] = 64; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 0; vertexIds[1] = 2; vertexIds[2] = 65; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 4; vertexIds[1] = 5; vertexIds[2] = 66; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 4; vertexIds[1] = 6; vertexIds[2] = 67; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 4; vertexIds[1] = 0; vertexIds[2] = 68; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 4; vertexIds[1] = 3; vertexIds[2] = 69; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 2; vertexIds[1] = 3; vertexIds[2] = 70; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 2; vertexIds[1] = 7; vertexIds[2] = 71; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 1; vertexIds[1] = 8; vertexIds[2] = 72; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 1; vertexIds[1] = 9; vertexIds[2] = 73; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 5; vertexIds[1] = 10; vertexIds[2] = 74; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 5; vertexIds[1] = 11; vertexIds[2] = 75; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 1; vertexIds[1] = 5; vertexIds[2] = 76; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 1; vertexIds[1] = 10; vertexIds[2] = 77; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 9; vertexIds[1] = 2; vertexIds[2] = 78; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 9; vertexIds[1] = 6; vertexIds[2] = 79; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 8; vertexIds[1] = 12; vertexIds[2] = 80; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 8; vertexIds[1] = 13; vertexIds[2] = 81; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 10; vertexIds[1] = 14; vertexIds[2] = 82; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 10; vertexIds[1] = 15; vertexIds[2] = 83; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 13; vertexIds[1] = 9; vertexIds[2] = 84; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 13; vertexIds[1] = 11; vertexIds[2] = 85; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 12; vertexIds[1] = 14; vertexIds[2] = 86; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 12; vertexIds[1] = 15; vertexIds[2] = 87; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 16; vertexIds[1] = 17; vertexIds[2] = 88; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 16; vertexIds[1] = 1; vertexIds[2] = 89; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 18; vertexIds[1] = 19; vertexIds[2] = 90; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 18; vertexIds[1] = 5; vertexIds[2] = 91; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 18; vertexIds[1] = 16; vertexIds[2] = 92; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 18; vertexIds[1] = 0; vertexIds[2] = 93; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 17; vertexIds[1] = 19; vertexIds[2] = 94; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 17; vertexIds[1] = 5; vertexIds[2] = 95; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 20; vertexIds[1] = 21; vertexIds[2] = 96; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 20; vertexIds[1] = 12; vertexIds[2] = 97; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 22; vertexIds[1] = 23; vertexIds[2] = 98; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 22; vertexIds[1] = 14; vertexIds[2] = 99; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 22; vertexIds[1] = 20; vertexIds[2] = 100; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 22; vertexIds[1] = 8; vertexIds[2] = 101; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 21; vertexIds[1] = 23; vertexIds[2] = 102; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 21; vertexIds[1] = 14; vertexIds[2] = 103; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 24; vertexIds[1] = 25; vertexIds[2] = 104; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 24; vertexIds[1] = 17; vertexIds[2] = 105; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 26; vertexIds[1] = 27; vertexIds[2] = 106; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 26; vertexIds[1] = 19; vertexIds[2] = 107; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 26; vertexIds[1] = 24; vertexIds[2] = 108; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 26; vertexIds[1] = 16; vertexIds[2] = 109; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 25; vertexIds[1] = 28; vertexIds[2] = 110; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 25; vertexIds[1] = 20; vertexIds[2] = 111; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 27; vertexIds[1] = 29; vertexIds[2] = 112; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 27; vertexIds[1] = 22; vertexIds[2] = 113; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 25; vertexIds[1] = 27; vertexIds[2] = 114; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 25; vertexIds[1] = 29; vertexIds[2] = 115; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 20; vertexIds[1] = 17; vertexIds[2] = 116; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 20; vertexIds[1] = 19; vertexIds[2] = 117; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 28; vertexIds[1] = 30; vertexIds[2] = 118; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 28; vertexIds[1] = 21; vertexIds[2] = 119; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 29; vertexIds[1] = 31; vertexIds[2] = 120; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 29; vertexIds[1] = 23; vertexIds[2] = 121; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 30; vertexIds[1] = 31; vertexIds[2] = 122; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 30; vertexIds[1] = 23; vertexIds[2] = 123; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 32; vertexIds[1] = 33; vertexIds[2] = 124; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 32; vertexIds[1] = 25; vertexIds[2] = 125; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 34; vertexIds[1] = 35; vertexIds[2] = 126; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 34; vertexIds[1] = 27; vertexIds[2] = 127; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 34; vertexIds[1] = 32; vertexIds[2] = 128; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 34; vertexIds[1] = 24; vertexIds[2] = 129; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 33; vertexIds[1] = 35; vertexIds[2] = 130; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 33; vertexIds[1] = 27; vertexIds[2] = 131; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 36; vertexIds[1] = 37; vertexIds[2] = 132; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 36; vertexIds[1] = 30; vertexIds[2] = 133; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 38; vertexIds[1] = 39; vertexIds[2] = 134; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 38; vertexIds[1] = 31; vertexIds[2] = 135; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 38; vertexIds[1] = 36; vertexIds[2] = 136; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 38; vertexIds[1] = 28; vertexIds[2] = 137; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 37; vertexIds[1] = 39; vertexIds[2] = 138; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 37; vertexIds[1] = 31; vertexIds[2] = 139; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 40; vertexIds[1] = 41; vertexIds[2] = 140; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 40; vertexIds[1] = 33; vertexIds[2] = 141; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 42; vertexIds[1] = 43; vertexIds[2] = 142; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 42; vertexIds[1] = 35; vertexIds[2] = 143; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 42; vertexIds[1] = 40; vertexIds[2] = 144; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 42; vertexIds[1] = 32; vertexIds[2] = 145; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 41; vertexIds[1] = 44; vertexIds[2] = 146; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 41; vertexIds[1] = 36; vertexIds[2] = 147; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 43; vertexIds[1] = 45; vertexIds[2] = 148; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 43; vertexIds[1] = 38; vertexIds[2] = 149; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 41; vertexIds[1] = 43; vertexIds[2] = 150; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 41; vertexIds[1] = 45; vertexIds[2] = 151; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 36; vertexIds[1] = 33; vertexIds[2] = 152; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 36; vertexIds[1] = 35; vertexIds[2] = 153; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 44; vertexIds[1] = 46; vertexIds[2] = 154; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 44; vertexIds[1] = 37; vertexIds[2] = 155; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 45; vertexIds[1] = 47; vertexIds[2] = 156; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 45; vertexIds[1] = 39; vertexIds[2] = 157; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 46; vertexIds[1] = 47; vertexIds[2] = 158; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 46; vertexIds[1] = 39; vertexIds[2] = 159; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 48; vertexIds[1] = 49; vertexIds[2] = 160; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 48; vertexIds[1] = 41; vertexIds[2] = 161; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 50; vertexIds[1] = 51; vertexIds[2] = 162; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 50; vertexIds[1] = 43; vertexIds[2] = 163; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 50; vertexIds[1] = 48; vertexIds[2] = 164; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 50; vertexIds[1] = 40; vertexIds[2] = 165; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 49; vertexIds[1] = 51; vertexIds[2] = 166; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 49; vertexIds[1] = 43; vertexIds[2] = 167; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 52; vertexIds[1] = 53; vertexIds[2] = 168; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 52; vertexIds[1] = 46; vertexIds[2] = 169; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 54; vertexIds[1] = 55; vertexIds[2] = 170; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 54; vertexIds[1] = 47; vertexIds[2] = 171; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 54; vertexIds[1] = 52; vertexIds[2] = 172; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 54; vertexIds[1] = 44; vertexIds[2] = 173; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 53; vertexIds[1] = 55; vertexIds[2] = 174; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 53; vertexIds[1] = 47; vertexIds[2] = 175; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 56; vertexIds[1] = 57; vertexIds[2] = 176; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 56; vertexIds[1] = 49; vertexIds[2] = 177; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 58; vertexIds[1] = 59; vertexIds[2] = 178; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 58; vertexIds[1] = 51; vertexIds[2] = 179; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 56; vertexIds[1] = 58; vertexIds[2] = 180; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 56; vertexIds[1] = 59; vertexIds[2] = 181; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 58; vertexIds[1] = 56; vertexIds[2] = 182; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 58; vertexIds[1] = 48; vertexIds[2] = 183; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 57; vertexIds[1] = 60; vertexIds[2] = 184; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 57; vertexIds[1] = 52; vertexIds[2] = 185; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 59; vertexIds[1] = 61; vertexIds[2] = 186; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 59; vertexIds[1] = 54; vertexIds[2] = 187; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 57; vertexIds[1] = 59; vertexIds[2] = 188; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 57; vertexIds[1] = 61; vertexIds[2] = 189; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 52; vertexIds[1] = 49; vertexIds[2] = 190; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 52; vertexIds[1] = 51; vertexIds[2] = 191; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 60; vertexIds[1] = 62; vertexIds[2] = 192; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 60; vertexIds[1] = 53; vertexIds[2] = 193; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 61; vertexIds[1] = 63; vertexIds[2] = 194; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 61; vertexIds[1] = 55; vertexIds[2] = 195; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 60; vertexIds[1] = 61; vertexIds[2] = 196; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 60; vertexIds[1] = 63; vertexIds[2] = 197; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 62; vertexIds[1] = 63; vertexIds[2] = 198; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 62; vertexIds[1] = 55; vertexIds[2] = 199; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 2; vertexIds[1] = 64; vertexIds[2] = 1; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 0; vertexIds[1] = 64; vertexIds[2] = 2; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 3; vertexIds[1] = 65; vertexIds[2] = 2; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 0; vertexIds[1] = 65; vertexIds[2] = 3; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 6; vertexIds[1] = 66; vertexIds[2] = 5; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 4; vertexIds[1] = 66; vertexIds[2] = 6; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 7; vertexIds[1] = 67; vertexIds[2] = 6; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 4; vertexIds[1] = 67; vertexIds[2] = 7; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 3; vertexIds[1] = 68; vertexIds[2] = 0; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 4; vertexIds[1] = 68; vertexIds[2] = 3; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 7; vertexIds[1] = 69; vertexIds[2] = 3; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 4; vertexIds[1] = 69; vertexIds[2] = 7; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 7; vertexIds[1] = 70; vertexIds[2] = 3; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 2; vertexIds[1] = 70; vertexIds[2] = 7; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 6; vertexIds[1] = 71; vertexIds[2] = 7; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 2; vertexIds[1] = 71; vertexIds[2] = 6; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 9; vertexIds[1] = 72; vertexIds[2] = 8; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 1; vertexIds[1] = 72; vertexIds[2] = 9; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 2; vertexIds[1] = 73; vertexIds[2] = 9; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 1; vertexIds[1] = 73; vertexIds[2] = 2; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 11; vertexIds[1] = 74; vertexIds[2] = 10; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 5; vertexIds[1] = 74; vertexIds[2] = 11; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 6; vertexIds[1] = 75; vertexIds[2] = 11; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 5; vertexIds[1] = 75; vertexIds[2] = 6; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 10; vertexIds[1] = 76; vertexIds[2] = 5; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 1; vertexIds[1] = 76; vertexIds[2] = 10; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 8; vertexIds[1] = 77; vertexIds[2] = 10; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 1; vertexIds[1] = 77; vertexIds[2] = 8; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 6; vertexIds[1] = 78; vertexIds[2] = 2; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 9; vertexIds[1] = 78; vertexIds[2] = 6; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 11; vertexIds[1] = 79; vertexIds[2] = 6; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 9; vertexIds[1] = 79; vertexIds[2] = 11; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 13; vertexIds[1] = 80; vertexIds[2] = 12; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 8; vertexIds[1] = 80; vertexIds[2] = 13; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 9; vertexIds[1] = 81; vertexIds[2] = 13; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 8; vertexIds[1] = 81; vertexIds[2] = 9; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 15; vertexIds[1] = 82; vertexIds[2] = 14; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 10; vertexIds[1] = 82; vertexIds[2] = 15; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 11; vertexIds[1] = 83; vertexIds[2] = 15; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 10; vertexIds[1] = 83; vertexIds[2] = 11; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 11; vertexIds[1] = 84; vertexIds[2] = 9; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 13; vertexIds[1] = 84; vertexIds[2] = 11; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 15; vertexIds[1] = 85; vertexIds[2] = 11; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 13; vertexIds[1] = 85; vertexIds[2] = 15; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 15; vertexIds[1] = 86; vertexIds[2] = 14; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 12; vertexIds[1] = 86; vertexIds[2] = 15; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 13; vertexIds[1] = 87; vertexIds[2] = 15; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 12; vertexIds[1] = 87; vertexIds[2] = 13; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 1; vertexIds[1] = 88; vertexIds[2] = 17; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 16; vertexIds[1] = 88; vertexIds[2] = 1; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 0; vertexIds[1] = 89; vertexIds[2] = 1; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 16; vertexIds[1] = 89; vertexIds[2] = 0; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 5; vertexIds[1] = 90; vertexIds[2] = 19; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 18; vertexIds[1] = 90; vertexIds[2] = 5; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 4; vertexIds[1] = 91; vertexIds[2] = 5; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 18; vertexIds[1] = 91; vertexIds[2] = 4; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 0; vertexIds[1] = 92; vertexIds[2] = 16; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 18; vertexIds[1] = 92; vertexIds[2] = 0; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 4; vertexIds[1] = 93; vertexIds[2] = 0; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 18; vertexIds[1] = 93; vertexIds[2] = 4; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 5; vertexIds[1] = 94; vertexIds[2] = 19; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 17; vertexIds[1] = 94; vertexIds[2] = 5; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 1; vertexIds[1] = 95; vertexIds[2] = 5; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 17; vertexIds[1] = 95; vertexIds[2] = 1; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 12; vertexIds[1] = 96; vertexIds[2] = 21; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 20; vertexIds[1] = 96; vertexIds[2] = 12; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 8; vertexIds[1] = 97; vertexIds[2] = 12; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 20; vertexIds[1] = 97; vertexIds[2] = 8; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 14; vertexIds[1] = 98; vertexIds[2] = 23; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 22; vertexIds[1] = 98; vertexIds[2] = 14; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 10; vertexIds[1] = 99; vertexIds[2] = 14; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 22; vertexIds[1] = 99; vertexIds[2] = 10; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 8; vertexIds[1] = 100; vertexIds[2] = 20; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 22; vertexIds[1] = 100; vertexIds[2] = 8; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 10; vertexIds[1] = 101; vertexIds[2] = 8; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 22; vertexIds[1] = 101; vertexIds[2] = 10; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 14; vertexIds[1] = 102; vertexIds[2] = 23; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 21; vertexIds[1] = 102; vertexIds[2] = 14; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 12; vertexIds[1] = 103; vertexIds[2] = 14; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 21; vertexIds[1] = 103; vertexIds[2] = 12; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 17; vertexIds[1] = 104; vertexIds[2] = 25; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 24; vertexIds[1] = 104; vertexIds[2] = 17; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 16; vertexIds[1] = 105; vertexIds[2] = 17; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 24; vertexIds[1] = 105; vertexIds[2] = 16; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 19; vertexIds[1] = 106; vertexIds[2] = 27; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 26; vertexIds[1] = 106; vertexIds[2] = 19; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 18; vertexIds[1] = 107; vertexIds[2] = 19; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 26; vertexIds[1] = 107; vertexIds[2] = 18; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 16; vertexIds[1] = 108; vertexIds[2] = 24; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 26; vertexIds[1] = 108; vertexIds[2] = 16; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 18; vertexIds[1] = 109; vertexIds[2] = 16; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 26; vertexIds[1] = 109; vertexIds[2] = 18; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 20; vertexIds[1] = 110; vertexIds[2] = 28; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 25; vertexIds[1] = 110; vertexIds[2] = 20; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 17; vertexIds[1] = 111; vertexIds[2] = 20; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 25; vertexIds[1] = 111; vertexIds[2] = 17; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 22; vertexIds[1] = 112; vertexIds[2] = 29; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 27; vertexIds[1] = 112; vertexIds[2] = 22; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 19; vertexIds[1] = 113; vertexIds[2] = 22; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 27; vertexIds[1] = 113; vertexIds[2] = 19; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 29; vertexIds[1] = 114; vertexIds[2] = 27; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 25; vertexIds[1] = 114; vertexIds[2] = 29; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 28; vertexIds[1] = 115; vertexIds[2] = 29; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 25; vertexIds[1] = 115; vertexIds[2] = 28; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 19; vertexIds[1] = 116; vertexIds[2] = 17; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 20; vertexIds[1] = 116; vertexIds[2] = 19; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 22; vertexIds[1] = 117; vertexIds[2] = 19; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 20; vertexIds[1] = 117; vertexIds[2] = 22; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 21; vertexIds[1] = 118; vertexIds[2] = 30; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 28; vertexIds[1] = 118; vertexIds[2] = 21; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 20; vertexIds[1] = 119; vertexIds[2] = 21; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 28; vertexIds[1] = 119; vertexIds[2] = 20; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 23; vertexIds[1] = 120; vertexIds[2] = 31; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 29; vertexIds[1] = 120; vertexIds[2] = 23; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 22; vertexIds[1] = 121; vertexIds[2] = 23; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 29; vertexIds[1] = 121; vertexIds[2] = 22; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 23; vertexIds[1] = 122; vertexIds[2] = 31; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 30; vertexIds[1] = 122; vertexIds[2] = 23; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 21; vertexIds[1] = 123; vertexIds[2] = 23; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 30; vertexIds[1] = 123; vertexIds[2] = 21; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 25; vertexIds[1] = 124; vertexIds[2] = 33; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 32; vertexIds[1] = 124; vertexIds[2] = 25; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 24; vertexIds[1] = 125; vertexIds[2] = 25; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 32; vertexIds[1] = 125; vertexIds[2] = 24; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 27; vertexIds[1] = 126; vertexIds[2] = 35; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 34; vertexIds[1] = 126; vertexIds[2] = 27; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 26; vertexIds[1] = 127; vertexIds[2] = 27; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 34; vertexIds[1] = 127; vertexIds[2] = 26; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 24; vertexIds[1] = 128; vertexIds[2] = 32; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 34; vertexIds[1] = 128; vertexIds[2] = 24; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 26; vertexIds[1] = 129; vertexIds[2] = 24; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 34; vertexIds[1] = 129; vertexIds[2] = 26; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 27; vertexIds[1] = 130; vertexIds[2] = 35; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 33; vertexIds[1] = 130; vertexIds[2] = 27; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 25; vertexIds[1] = 131; vertexIds[2] = 27; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 33; vertexIds[1] = 131; vertexIds[2] = 25; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 30; vertexIds[1] = 132; vertexIds[2] = 37; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 36; vertexIds[1] = 132; vertexIds[2] = 30; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 28; vertexIds[1] = 133; vertexIds[2] = 30; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 36; vertexIds[1] = 133; vertexIds[2] = 28; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 31; vertexIds[1] = 134; vertexIds[2] = 39; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 38; vertexIds[1] = 134; vertexIds[2] = 31; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 29; vertexIds[1] = 135; vertexIds[2] = 31; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 38; vertexIds[1] = 135; vertexIds[2] = 29; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 28; vertexIds[1] = 136; vertexIds[2] = 36; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 38; vertexIds[1] = 136; vertexIds[2] = 28; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 29; vertexIds[1] = 137; vertexIds[2] = 28; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 38; vertexIds[1] = 137; vertexIds[2] = 29; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 31; vertexIds[1] = 138; vertexIds[2] = 39; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 37; vertexIds[1] = 138; vertexIds[2] = 31; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 30; vertexIds[1] = 139; vertexIds[2] = 31; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 37; vertexIds[1] = 139; vertexIds[2] = 30; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 33; vertexIds[1] = 140; vertexIds[2] = 41; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 40; vertexIds[1] = 140; vertexIds[2] = 33; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 32; vertexIds[1] = 141; vertexIds[2] = 33; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 40; vertexIds[1] = 141; vertexIds[2] = 32; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 35; vertexIds[1] = 142; vertexIds[2] = 43; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 42; vertexIds[1] = 142; vertexIds[2] = 35; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 34; vertexIds[1] = 143; vertexIds[2] = 35; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 42; vertexIds[1] = 143; vertexIds[2] = 34; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 32; vertexIds[1] = 144; vertexIds[2] = 40; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 42; vertexIds[1] = 144; vertexIds[2] = 32; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 34; vertexIds[1] = 145; vertexIds[2] = 32; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 42; vertexIds[1] = 145; vertexIds[2] = 34; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 36; vertexIds[1] = 146; vertexIds[2] = 44; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 41; vertexIds[1] = 146; vertexIds[2] = 36; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 33; vertexIds[1] = 147; vertexIds[2] = 36; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 41; vertexIds[1] = 147; vertexIds[2] = 33; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 38; vertexIds[1] = 148; vertexIds[2] = 45; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 43; vertexIds[1] = 148; vertexIds[2] = 38; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 35; vertexIds[1] = 149; vertexIds[2] = 38; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 43; vertexIds[1] = 149; vertexIds[2] = 35; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 45; vertexIds[1] = 150; vertexIds[2] = 43; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 41; vertexIds[1] = 150; vertexIds[2] = 45; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 44; vertexIds[1] = 151; vertexIds[2] = 45; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 41; vertexIds[1] = 151; vertexIds[2] = 44; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 35; vertexIds[1] = 152; vertexIds[2] = 33; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 36; vertexIds[1] = 152; vertexIds[2] = 35; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 38; vertexIds[1] = 153; vertexIds[2] = 35; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 36; vertexIds[1] = 153; vertexIds[2] = 38; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 37; vertexIds[1] = 154; vertexIds[2] = 46; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 44; vertexIds[1] = 154; vertexIds[2] = 37; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 36; vertexIds[1] = 155; vertexIds[2] = 37; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 44; vertexIds[1] = 155; vertexIds[2] = 36; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 39; vertexIds[1] = 156; vertexIds[2] = 47; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 45; vertexIds[1] = 156; vertexIds[2] = 39; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 38; vertexIds[1] = 157; vertexIds[2] = 39; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 45; vertexIds[1] = 157; vertexIds[2] = 38; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 39; vertexIds[1] = 158; vertexIds[2] = 47; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 46; vertexIds[1] = 158; vertexIds[2] = 39; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 37; vertexIds[1] = 159; vertexIds[2] = 39; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 46; vertexIds[1] = 159; vertexIds[2] = 37; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 41; vertexIds[1] = 160; vertexIds[2] = 49; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 48; vertexIds[1] = 160; vertexIds[2] = 41; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 40; vertexIds[1] = 161; vertexIds[2] = 41; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 48; vertexIds[1] = 161; vertexIds[2] = 40; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 43; vertexIds[1] = 162; vertexIds[2] = 51; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 50; vertexIds[1] = 162; vertexIds[2] = 43; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 42; vertexIds[1] = 163; vertexIds[2] = 43; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 50; vertexIds[1] = 163; vertexIds[2] = 42; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 40; vertexIds[1] = 164; vertexIds[2] = 48; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 50; vertexIds[1] = 164; vertexIds[2] = 40; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 42; vertexIds[1] = 165; vertexIds[2] = 40; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 50; vertexIds[1] = 165; vertexIds[2] = 42; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 43; vertexIds[1] = 166; vertexIds[2] = 51; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 49; vertexIds[1] = 166; vertexIds[2] = 43; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 41; vertexIds[1] = 167; vertexIds[2] = 43; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 49; vertexIds[1] = 167; vertexIds[2] = 41; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 46; vertexIds[1] = 168; vertexIds[2] = 53; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 52; vertexIds[1] = 168; vertexIds[2] = 46; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 44; vertexIds[1] = 169; vertexIds[2] = 46; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 52; vertexIds[1] = 169; vertexIds[2] = 44; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 47; vertexIds[1] = 170; vertexIds[2] = 55; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 54; vertexIds[1] = 170; vertexIds[2] = 47; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 45; vertexIds[1] = 171; vertexIds[2] = 47; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 54; vertexIds[1] = 171; vertexIds[2] = 45; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 44; vertexIds[1] = 172; vertexIds[2] = 52; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 54; vertexIds[1] = 172; vertexIds[2] = 44; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 45; vertexIds[1] = 173; vertexIds[2] = 44; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 54; vertexIds[1] = 173; vertexIds[2] = 45; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 47; vertexIds[1] = 174; vertexIds[2] = 55; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 53; vertexIds[1] = 174; vertexIds[2] = 47; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 46; vertexIds[1] = 175; vertexIds[2] = 47; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 53; vertexIds[1] = 175; vertexIds[2] = 46; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 49; vertexIds[1] = 176; vertexIds[2] = 57; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 56; vertexIds[1] = 176; vertexIds[2] = 49; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 48; vertexIds[1] = 177; vertexIds[2] = 49; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 56; vertexIds[1] = 177; vertexIds[2] = 48; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 51; vertexIds[1] = 178; vertexIds[2] = 59; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 58; vertexIds[1] = 178; vertexIds[2] = 51; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 50; vertexIds[1] = 179; vertexIds[2] = 51; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 58; vertexIds[1] = 179; vertexIds[2] = 50; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 59; vertexIds[1] = 180; vertexIds[2] = 58; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 56; vertexIds[1] = 180; vertexIds[2] = 59; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 57; vertexIds[1] = 181; vertexIds[2] = 59; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 56; vertexIds[1] = 181; vertexIds[2] = 57; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 48; vertexIds[1] = 182; vertexIds[2] = 56; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 58; vertexIds[1] = 182; vertexIds[2] = 48; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 50; vertexIds[1] = 183; vertexIds[2] = 48; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 58; vertexIds[1] = 183; vertexIds[2] = 50; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 52; vertexIds[1] = 184; vertexIds[2] = 60; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 57; vertexIds[1] = 184; vertexIds[2] = 52; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 49; vertexIds[1] = 185; vertexIds[2] = 52; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 57; vertexIds[1] = 185; vertexIds[2] = 49; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 54; vertexIds[1] = 186; vertexIds[2] = 61; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 59; vertexIds[1] = 186; vertexIds[2] = 54; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 51; vertexIds[1] = 187; vertexIds[2] = 54; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 59; vertexIds[1] = 187; vertexIds[2] = 51; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 61; vertexIds[1] = 188; vertexIds[2] = 59; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 57; vertexIds[1] = 188; vertexIds[2] = 61; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 60; vertexIds[1] = 189; vertexIds[2] = 61; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 57; vertexIds[1] = 189; vertexIds[2] = 60; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 51; vertexIds[1] = 190; vertexIds[2] = 49; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 52; vertexIds[1] = 190; vertexIds[2] = 51; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 54; vertexIds[1] = 191; vertexIds[2] = 51; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 52; vertexIds[1] = 191; vertexIds[2] = 54; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 53; vertexIds[1] = 192; vertexIds[2] = 62; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 60; vertexIds[1] = 192; vertexIds[2] = 53; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 52; vertexIds[1] = 193; vertexIds[2] = 53; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 60; vertexIds[1] = 193; vertexIds[2] = 52; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 55; vertexIds[1] = 194; vertexIds[2] = 63; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 61; vertexIds[1] = 194; vertexIds[2] = 55; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 54; vertexIds[1] = 195; vertexIds[2] = 55; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 61; vertexIds[1] = 195; vertexIds[2] = 54; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 63; vertexIds[1] = 196; vertexIds[2] = 61; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 60; vertexIds[1] = 196; vertexIds[2] = 63; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 62; vertexIds[1] = 197; vertexIds[2] = 63; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 60; vertexIds[1] = 197; vertexIds[2] = 62; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 55; vertexIds[1] = 198; vertexIds[2] = 63; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 62; vertexIds[1] = 198; vertexIds[2] = 55; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 53; vertexIds[1] = 199; vertexIds[2] = 55; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); vertexIds[0] = 62; vertexIds[1] = 199; vertexIds[2] = 53; sMesh->InsertNextCell(VTK_TRIANGLE, 3, vertexIds); free(vertexIds); return 1; } class AreaSimplificationMetric : public vtkReebGraphSimplificationMetric{ public: vtkTypeMacro(AreaSimplificationMetric, vtkReebGraphSimplificationMetric); static AreaSimplificationMetric* New(); double ComputeMetric(vtkDataSet *mesh, vtkDataArray *scalarField, vtkIdType startCriticalPoint, vtkAbstractArray *vertexList, vtkIdType endCriticalPoint); }; vtkStandardNewMacro(AreaSimplificationMetric); double AreaSimplificationMetric::ComputeMetric(vtkDataSet *mesh, vtkDataArray *scalarField, vtkIdType startCriticalPoint, vtkAbstractArray* vertexList, vtkIdType endCriticalPoint) { // In this example, the metric algorithm just evaluates the area of the // surface region corresponding to the arc of the Reeb graph passed as an // argument. // As a result, the arcs corresponding to small surface regions (below the // threshold specified to the simplificatin filter) will be // simplified in priority in the surface simplification algorithm. double fieldLowerBound = scalarField->GetComponent(startCriticalPoint,0), fieldUpperBound = scalarField->GetComponent(endCriticalPoint,0); double cumulativeArea = 0; std::map visitedTriangles; for(int i = 0; i < vertexList->GetNumberOfTuples(); i++) { int vId = vertexList->GetVariantValue(i).ToInt(); vtkIdList *starTriangleList = vtkIdList::New(); mesh->GetPointCells(vId, starTriangleList); for(int j = 0; j < starTriangleList->GetNumberOfIds(); j++) { vtkIdType tId = starTriangleList->GetId(j); vtkTriangle *t = vtkTriangle::SafeDownCast(mesh->GetCell(tId)); std::map::iterator tIt = visitedTriangles.find(tId); if(tIt == visitedTriangles.end()) { if((scalarField->GetComponent(t->GetPointIds()->GetId(0), 0) <= fieldUpperBound) &&(scalarField->GetComponent(t->GetPointIds()->GetId(1), 0) <= fieldUpperBound) &&(scalarField->GetComponent(t->GetPointIds()->GetId(2), 0) <= fieldUpperBound) &&(scalarField->GetComponent(t->GetPointIds()->GetId(0), 0) >= fieldLowerBound) &&(scalarField->GetComponent(t->GetPointIds()->GetId(1), 0) >= fieldLowerBound) &&(scalarField->GetComponent(t->GetPointIds()->GetId(2), 0) >= fieldLowerBound)) { // the triangle fully maps inside the arc function interval cumulativeArea += t->ComputeArea(); } visitedTriangles[tId] = true; } } starTriangleList->Delete(); } return cumulativeArea/(this->UpperBound - this->LowerBound); } int DisplayReebGraph(vtkReebGraph *g) { vtkDataArray *vertexInfo = vtkDataArray::SafeDownCast( g->GetVertexData()->GetAbstractArray("Vertex Ids")); if(!vertexInfo) return 1; vtkVariantArray *edgeInfo = vtkVariantArray::SafeDownCast( g->GetEdgeData()->GetAbstractArray("Vertex Ids")); if(!edgeInfo) return 2; cout << " Reeb graph nodes:" << endl; for(int i = 0; i < vertexInfo->GetNumberOfTuples(); i++) cout << " Node #" << i << ") VertexMeshId: " << ((int) *(vertexInfo->GetTuple(i))) << endl; cout << " Reeb graph arcs:" << endl; vtkEdgeListIterator *eIt = vtkEdgeListIterator::New(); g->GetEdges(eIt); do{ vtkEdgeType e = eIt->Next(); vtkAbstractArray *deg2NodeList = edgeInfo->GetPointer(e.Id)->ToArray(); cout << " Arc #" << e.Id << ": " << *(vertexInfo->GetTuple(e.Source)) << " -> " << *(vertexInfo->GetTuple(e.Target)) << " (" << deg2NodeList->GetNumberOfTuples() << " degree-2 nodes)" << endl; }while(eIt->HasNext()); eIt->Delete(); return 0; } int DisplaySurfaceSkeleton(vtkPolyData *surfaceMesh, vtkTable *skeleton) { // Rendering setting vtkRenderer *renderer = vtkRenderer::New(); vtkRenderWindow *renderWindow = vtkRenderWindow::New(); renderWindow->AddRenderer(renderer); renderWindow->SetSize(400, 400); vtkRenderWindowInteractor *windowInteractor = vtkRenderWindowInteractor::New(); windowInteractor->SetRenderWindow(renderWindow); vtkPolyDataMapper *surfaceMapper = vtkPolyDataMapper::New(); surfaceMapper->SetInput(surfaceMesh); vtkActor *surfaceActor = vtkActor::New(); surfaceActor->SetMapper(surfaceMapper); surfaceActor->GetProperty()->SetOpacity(0.3); renderer->AddActor(surfaceActor); // Actual display of the skeleton vtkSphereSource *nodeSphere = vtkSphereSource::New(); nodeSphere->SetThetaResolution(50); nodeSphere->SetPhiResolution(20); nodeSphere->SetRadius(0.03); vtkPolyDataMapper *sphereMapper = vtkPolyDataMapper::New(); sphereMapper->SetInputConnection(nodeSphere->GetOutputPort()); // 2 nodes per arc of the skeleton vtkActor **nodeActors = (vtkActor **) malloc(sizeof(vtkActor *)* 2*skeleton->GetNumberOfColumns()); int sampleId = 0; double *point = (double *) malloc(sizeof(double)*3); vtkPolyData *embeddedSkeleton = vtkPolyData::New(); embeddedSkeleton->Allocate(); vtkPoints *skeletonSamples = vtkPoints::New(); skeletonSamples->SetNumberOfPoints( skeleton->GetNumberOfColumns()*skeleton->GetNumberOfRows()); for(int i = 0; i < skeleton->GetNumberOfColumns(); i++) { vtkDoubleArray *arc = vtkDoubleArray::SafeDownCast(skeleton->GetColumn(i)); // critical point at the origin of the arc arc->GetTupleValue(0, point); nodeActors[2*i] = vtkActor::New(); nodeActors[2*i]->SetMapper(sphereMapper); nodeActors[2*i]->GetProperty()->SetColor(0, 0, 1); nodeActors[2*i]->SetPosition(point); renderer->AddActor(nodeActors[2*i]); arc->GetTupleValue(arc->GetNumberOfTuples() - 1, point); nodeActors[2*i + 1] = vtkActor::New(); nodeActors[2*i + 1]->SetMapper(sphereMapper); nodeActors[2*i + 1]->GetProperty()->SetColor(0, 0, 1); nodeActors[2*i + 1]->SetPosition(point); renderer->AddActor(nodeActors[2*i + 1]); // now add the samples to the skeleton polyData int initialSampleId = sampleId; for(int j = 0; j < arc->GetNumberOfTuples(); j++) { arc->GetTupleValue(j, point); skeletonSamples->SetPoint(sampleId, point); sampleId++; } for(int j = 1; j < arc->GetNumberOfTuples(); j++) { vtkIdType samplePair[2]; samplePair[0] = j - 1 + initialSampleId; samplePair[1] = j + initialSampleId; embeddedSkeleton->InsertNextCell(VTK_LINE, 2, samplePair); } } embeddedSkeleton->SetPoints(skeletonSamples); free(point); skeletonSamples->Delete(); vtkPolyDataMapper *lineMapper = vtkPolyDataMapper::New(); lineMapper->SetInput(embeddedSkeleton); vtkActor *skeletonActor = vtkActor::New(); skeletonActor->SetMapper(lineMapper); skeletonActor->GetProperty()->SetColor(0, 1, 0); renderer->AddActor(skeletonActor); windowInteractor->Initialize(); // Interactive mode //windowInteractor->Start(); // Testing mode windowInteractor->Render(); vtkWindowToImageFilter *windowToImage = vtkWindowToImageFilter::New(); windowToImage->SetInput(renderWindow); vtkPNGWriter *pngWriter = vtkPNGWriter::New(); pngWriter->SetInputConnection(windowToImage->GetOutputPort()); pngWriter->SetFileName("TestReebGraph2D.png"); pngWriter->Write(); pngWriter->Delete(); windowToImage->Delete(); // End of testing mode skeletonActor->Delete(); lineMapper->Delete(); for(int i = 0; i < 2*skeleton->GetNumberOfColumns(); i++) nodeActors[i]->Delete(); embeddedSkeleton->Delete(); free(nodeActors); sphereMapper->Delete(); nodeSphere->Delete(); surfaceActor->Delete(); surfaceMapper->Delete(); windowInteractor->Delete(); renderWindow->Delete(); renderer->Delete(); return 0; } int DisplayVolumeSkeleton(vtkUnstructuredGrid* vtkNotUsed(volumeMesh), vtkTable *skeleton) { // no rendering of the mesh this time, just the skeleton. // one could just display the bounding surface of the volume (if the volume is // embeddable in IR^3). // Rendering setting vtkRenderer *renderer = vtkRenderer::New(); vtkRenderWindow *renderWindow = vtkRenderWindow::New(); renderWindow->AddRenderer(renderer); renderWindow->SetSize(400, 400); vtkRenderWindowInteractor *windowInteractor = vtkRenderWindowInteractor::New(); windowInteractor->SetRenderWindow(renderWindow); // Actual display of the skeleton vtkSphereSource *nodeSphere = vtkSphereSource::New(); nodeSphere->SetThetaResolution(50); nodeSphere->SetPhiResolution(20); nodeSphere->SetRadius(0.03); vtkPolyDataMapper *sphereMapper = vtkPolyDataMapper::New(); sphereMapper->SetInputConnection(nodeSphere->GetOutputPort()); // 2 nodes per arc of the skeleton vtkActor **nodeActors = (vtkActor **) malloc(sizeof(vtkActor *)* 2*skeleton->GetNumberOfColumns()); int sampleId = 0; double *point = (double *) malloc(sizeof(double)*3); vtkPolyData *embeddedSkeleton = vtkPolyData::New(); embeddedSkeleton->Allocate(); vtkPoints *skeletonSamples = vtkPoints::New(); skeletonSamples->SetNumberOfPoints( skeleton->GetNumberOfColumns()*skeleton->GetNumberOfRows()); for(int i = 0; i < skeleton->GetNumberOfColumns(); i++) { vtkDoubleArray *arc = vtkDoubleArray::SafeDownCast(skeleton->GetColumn(i)); // critical point at the origin of the arc arc->GetTupleValue(0, point); nodeActors[2*i] = vtkActor::New(); nodeActors[2*i]->SetMapper(sphereMapper); nodeActors[2*i]->GetProperty()->SetColor(0, 0, 1); nodeActors[2*i]->SetPosition(point); renderer->AddActor(nodeActors[2*i]); arc->GetTupleValue(arc->GetNumberOfTuples() - 1, point); nodeActors[2*i + 1] = vtkActor::New(); nodeActors[2*i + 1]->SetMapper(sphereMapper); nodeActors[2*i + 1]->GetProperty()->SetColor(0, 0, 1); nodeActors[2*i + 1]->SetPosition(point); renderer->AddActor(nodeActors[2*i + 1]); // now add the samples to the skeleton polyData int initialSampleId = sampleId; for(int j = 0; j < arc->GetNumberOfTuples(); j++) { arc->GetTupleValue(j, point); skeletonSamples->SetPoint(sampleId, point); sampleId++; } for(int j = 1; j < arc->GetNumberOfTuples(); j++) { vtkIdType samplePair[2]; samplePair[0] = j - 1 + initialSampleId; samplePair[1] = j + initialSampleId; embeddedSkeleton->InsertNextCell(VTK_LINE, 2, samplePair); } } embeddedSkeleton->SetPoints(skeletonSamples); free(point); skeletonSamples->Delete(); vtkPolyDataMapper *lineMapper = vtkPolyDataMapper::New(); lineMapper->SetInput(embeddedSkeleton); vtkActor *skeletonActor = vtkActor::New(); skeletonActor->SetMapper(lineMapper); skeletonActor->GetProperty()->SetColor(0, 1, 0); renderer->AddActor(skeletonActor); windowInteractor->Initialize(); // Interactive mode //windowInteractor->Start(); // Testing mode windowInteractor->Render(); vtkWindowToImageFilter *windowToImage = vtkWindowToImageFilter::New(); windowToImage->SetInput(renderWindow); vtkPNGWriter *pngWriter = vtkPNGWriter::New(); pngWriter->SetInputConnection(windowToImage->GetOutputPort()); pngWriter->SetFileName("TestReebGraph3D.png"); pngWriter->Write(); pngWriter->Delete(); windowToImage->Delete(); // End of testing mode skeletonActor->Delete(); lineMapper->Delete(); for(int i = 0; i < 2*skeleton->GetNumberOfColumns(); i++) nodeActors[i]->Delete(); embeddedSkeleton->Delete(); free(nodeActors); sphereMapper->Delete(); nodeSphere->Delete(); windowInteractor->Delete(); renderWindow->Delete(); renderer->Delete(); return 0; } int TestReebGraph(int vtkNotUsed(argc), char* vtkNotUsed(argv)[] ) { int errorCode; cout << endl << "Reeb Graph Tests ========================== Surface Mesh Tests" << endl; // Loading the mesh vtkPolyData *surfaceMesh = vtkPolyData::New(); LoadSurfaceMesh(surfaceMesh); // Attaching a height scalar field to it vtkDoubleArray *surfaceScalarField = vtkDoubleArray::New(); surfaceScalarField->SetNumberOfTuples(surfaceMesh->GetNumberOfPoints()); for(vtkIdType vId = 0; vId < surfaceMesh->GetNumberOfPoints(); vId++) { double *p = (double *) malloc(sizeof(double)*3); surfaceMesh->GetPoint(vId, p); double scalarValue = p[1]; // add a bit of noise for the split tree test if(vId == 2) scalarValue -= 10*scalarValue; surfaceScalarField->SetTuple1(vId, scalarValue); free(p); } surfaceMesh->GetPointData()->SetScalars(surfaceScalarField); cout << " Test 2D.1 Reeb graph computation... " << endl; vtkPolyDataToReebGraphFilter *surfaceReebGraphFilter = vtkPolyDataToReebGraphFilter::New(); surfaceReebGraphFilter->SetInput(surfaceMesh); surfaceReebGraphFilter->Update(); vtkReebGraph *surfaceReebGraph = surfaceReebGraphFilter->GetOutput(); cout << " Test 2D.1 "; if(surfaceReebGraph->GetNumberOfEdges() == 12) cout << "OK!" << endl; else { cout << "Failed!" << endl; return EXIT_FAILURE; } cout << " Test 2D.2 Customized Reeb graph simplification... " << endl; vtkReebGraphSimplificationFilter *surfaceSimplification = vtkReebGraphSimplificationFilter::New(); AreaSimplificationMetric *metric = AreaSimplificationMetric::New(); metric->SetLowerBound(0); // determining the maximum area double globalArea = 0; for(int i = 0; i < surfaceMesh->GetNumberOfCells(); i++) { vtkTriangle *t = vtkTriangle::SafeDownCast(surfaceMesh->GetCell(i)); globalArea += t->ComputeArea(); } metric->SetUpperBound(globalArea); surfaceSimplification->SetSimplificationMetric(metric); surfaceSimplification->SetInput(surfaceReebGraph); surfaceSimplification->SetSimplificationThreshold(0.01); surfaceSimplification->Update(); vtkReebGraph *simplifiedSurfaceReebGraph = surfaceSimplification->GetOutput(); metric->Delete(); cout << " Test 2D.2 "; if(simplifiedSurfaceReebGraph->GetNumberOfEdges() == 12) cout << "OK!" << endl; else { cout << "Failed!" << endl; return EXIT_FAILURE; } cout << " Test 2D.3 Reeb graph traversal..." << endl; errorCode = DisplayReebGraph(simplifiedSurfaceReebGraph); cout << " Test 2D.3 "; if(!errorCode) cout << "OK!" << endl; else { cout << "Failed! (code " << errorCode << ")" << endl; return EXIT_FAILURE; } cout << " Test 2D.4 Reeb graph based surface skeleton... " << endl; vtkReebGraphSurfaceSkeletonFilter *surfaceSkeletonFilter = vtkReebGraphSurfaceSkeletonFilter::New(); surfaceSkeletonFilter->SetInput(0, surfaceMesh); surfaceSkeletonFilter->SetInput(1, simplifiedSurfaceReebGraph); surfaceSkeletonFilter->SetNumberOfSamples(5); surfaceSkeletonFilter->Update(); vtkTable *surfaceSkeleton = surfaceSkeletonFilter->GetOutput(); errorCode = DisplaySurfaceSkeleton(surfaceMesh, surfaceSkeleton); cout << " Test 2D.4 "; if(surfaceSkeleton->GetNumberOfColumns() == 12) cout << "OK!" << endl; else { cout << "Failed!" << endl; return EXIT_FAILURE; } cout << " Test 2D.5 Area contour spectrum..." << endl; vtkAreaContourSpectrumFilter *areaSpectrumFilter = vtkAreaContourSpectrumFilter::New(); areaSpectrumFilter->SetInput(0, surfaceMesh); areaSpectrumFilter->SetInput(1, simplifiedSurfaceReebGraph); areaSpectrumFilter->SetArcId(0); areaSpectrumFilter->Update(); vtkTable *areaSpectrum = areaSpectrumFilter->GetOutput(); cout << " Test 2D.5 "; if(areaSpectrum->GetNumberOfRows() == 100) cout << "OK!" << endl; else { cout << "Failed!" << endl; return EXIT_FAILURE; } cout << " Test 2D.6 Reeb graph to split tree filter..." << endl; cout << " Not currently tested..." << endl; // vtkReebGraphToJoinSplitTreeFilter *splitTreeFilter = // vtkReebGraphToJoinSplitTreeFilter::New(); // splitTreeFilter->SetInput(0, surfaceMesh); // splitTreeFilter->SetInput(1, simplifiedSurfaceReebGraph); // splitTreeFilter->SetIsSplitTree(true); // splitTreeFilter->Update(); // vtkReebGraph *splitTree = splitTreeFilter->GetOutput(); // DisplayReebGraph(splitTree); // cout << " Test 2D.6 "; // if(splitTree->GetNumberOfEdges() == 3) // cout << "OK!" << endl; // else // { // cout << "Failed!" << endl; // return EXIT_FAILURE; // } // splitTreeFilter->Delete(); areaSpectrumFilter->Delete(); surfaceSkeletonFilter->Delete(); surfaceSimplification->Delete(); surfaceReebGraphFilter->Delete(); surfaceScalarField->Delete(); surfaceMesh->Delete(); cout << endl << endl << endl << "Reeb Graph Tests ========================== Volume Mesh Tests" << endl; // Loading the mesh vtkUnstructuredGrid *volumeMesh = vtkUnstructuredGrid::New(); LoadVolumeMesh(volumeMesh); // Attaching a height scalar field to it vtkDoubleArray *volumeScalarField = vtkDoubleArray::New(); volumeScalarField->SetNumberOfTuples(volumeMesh->GetNumberOfPoints()); for(vtkIdType vId = 0; vId < volumeMesh->GetNumberOfPoints(); vId++) { double *p = (double *) malloc(sizeof(double)*3); volumeMesh->GetPoint(vId, p); volumeScalarField->SetTuple1(vId, p[1]); free(p); } volumeMesh->GetPointData()->SetScalars(volumeScalarField); cout << " Test 3D.1 Reeb graph computation... " << endl; vtkUnstructuredGridToReebGraphFilter *volumeReebGraphFilter = vtkUnstructuredGridToReebGraphFilter::New(); volumeReebGraphFilter->SetInput(volumeMesh); volumeReebGraphFilter->Update(); vtkReebGraph *volumeReebGraph = volumeReebGraphFilter->GetOutput(); cout << " Test 3D.1 "; if(volumeReebGraph->GetNumberOfEdges() == 10) cout << "OK!" << endl; else { cout << "Failed!" << endl; return EXIT_FAILURE; } cout << " Test 3D.2 Customized Reeb graph simplification... " << endl; // in this example, we don't define any custom simplification metric and use // the default one (persistence). vtkReebGraphSimplificationFilter *volumeSimplification = vtkReebGraphSimplificationFilter::New(); volumeSimplification->SetInput(volumeReebGraph); volumeSimplification->SetSimplificationThreshold(0.05); volumeSimplification->Update(); vtkReebGraph *simplifiedVolumeReebGraph = volumeSimplification->GetOutput(); cout << " Test 3D.2 "; if(simplifiedVolumeReebGraph->GetNumberOfEdges() == 10) cout << "OK!" << endl; else { cout << "Failed!" << endl; return EXIT_FAILURE; } cout << " Test 3D.3 Reeb graph traversal..." << endl; errorCode = DisplayReebGraph(simplifiedVolumeReebGraph); cout << " Test 3D.3 "; if(!errorCode) cout << "OK!" << endl; else { cout << "Failed! (code " << errorCode << ")" << endl; return EXIT_FAILURE; } cout << " Test 3D.4 Reeb graph based volume skeleton... " << endl; vtkReebGraphVolumeSkeletonFilter *volumeSkeletonFilter = vtkReebGraphVolumeSkeletonFilter::New(); volumeSkeletonFilter->SetInput(0, volumeMesh); volumeSkeletonFilter->SetInput(1, simplifiedVolumeReebGraph); volumeSkeletonFilter->Update(); vtkTable *volumeSkeleton = volumeSkeletonFilter->GetOutput(); errorCode = DisplayVolumeSkeleton(volumeMesh, volumeSkeleton); cout << " Test 3D.4 "; if(volumeSkeleton->GetNumberOfColumns() == 10) cout << "OK!" << endl; else { cout << "Failed!" << endl; return EXIT_FAILURE; } cout << " Test 3D.5 Volume contour spectrum..." << endl; vtkVolumeContourSpectrumFilter *volumeSpectrumFilter = vtkVolumeContourSpectrumFilter::New(); volumeSpectrumFilter->SetInput(0, volumeMesh); volumeSpectrumFilter->SetInput(1, simplifiedVolumeReebGraph); volumeSpectrumFilter->SetArcId(0); volumeSpectrumFilter->Update(); vtkTable *volumeSpectrum = volumeSpectrumFilter->GetOutput(); cout << " Test 3D.5 "; if(volumeSpectrum->GetNumberOfRows() == 100) cout << "OK!" << endl; else { cout << "Failed!" << endl; return EXIT_FAILURE; } cout << " Test 3D.6 Reeb graph to join tree filter..." << endl; cout << " Not currently tested..." << endl; // vtkReebGraphToJoinSplitTreeFilter *joinTreeFilter = // vtkReebGraphToJoinSplitTreeFilter::New(); // joinTreeFilter->SetInput(0, volumeMesh); // joinTreeFilter->SetInput(1, simplifiedVolumeReebGraph); // joinTreeFilter->SetIsSplitTree(true); // joinTreeFilter->Update(); // vtkReebGraph *joinTree = joinTreeFilter->GetOutput(); // DisplayReebGraph(joinTree); // cout << " Test 3D.6 "; // if(joinTree->GetNumberOfEdges() == 1) // cout << "OK!" << endl; // else // { // cout << "Failed!" << endl; // return EXIT_FAILURE; // } // joinTreeFilter->Delete(); volumeSpectrumFilter->Delete(); volumeSkeletonFilter->Delete(); volumeSimplification->Delete(); volumeReebGraphFilter->Delete(); volumeScalarField->Delete(); volumeMesh->Delete(); return 0; }