Refactor/oduran/speed up integration tests

src/coreComponents/integrationTests/fluidFlowTests/testMixedDimHydrostaticEquilibrium.cpp

@@ -45,9 +45,7 @@ INSTANTIATE_TEST_SUITE_P(
       "DFN_5_fractures_hex_binarized.vtu",
       "DFN_5_fractures_tet_binarized.vtu"
       ),
-    ::testing::Values(
-      std::make_tuple( 1, 1, 1 )
-      )
+    ::testing::Values( 1 )
     )
   );
 

src/coreComponents/integrationTests/fluidFlowTests/testMixedDimHydrostaticEquilibriumFixture.hpp

@@ -60,10 +60,10 @@ extern CommandLineOptions g_commandLineOptions;
  *
  * Parametrized with std::tuple<std::string, std::tuple<int,int,int>>:
  *  - std::string:                Mesh file name (e.g. "fractured_mesh_hex_DFN_1.vtu")
- *  - tuple<int, int, int>:       Number of MPI partitions in x, y, z directions
+ *  - int:       Number of MPI partitions
  */
 class MixedDimHydrostaticEquilibriumTest
-  : public ::testing::TestWithParam< std::tuple< std::string, std::tuple< int, int, int > > >
+  : public ::testing::TestWithParam< std::tuple< std::string, int > >
 {
 protected:
   void SetUp() override
@@ -240,10 +240,7 @@ class MixedDimHydrostaticEquilibriumTest
 TEST_P( MixedDimHydrostaticEquilibriumTest, Run )
 {
   std::string const & meshFileName = std::get< 0 >( GetParam() );
-  std::tuple< int, int, int > const & partitions = std::get< 1 >( GetParam() );
-  int const xPartitions = std::get< 0 >( partitions );
-  int const yPartitions = std::get< 1 >( partitions );
-  int const zPartitions = std::get< 2 >( partitions );
+  int const n_partitions = std::get< 1 >( GetParam() );
 
   std::string nodeSetNames = "{ f1_node_set }";
   if( meshFileName.find( "_DFN_123.vtu" ) != std::string::npos )
@@ -292,10 +289,7 @@ TEST_P( MixedDimHydrostaticEquilibriumTest, Run )
   std::string::size_type const extPos = baseName.rfind( ".vtu" );
   if( extPos != std::string::npos )
     baseName.erase( extPos );
-  std::string const xmlName = "test_hydrostatic_" + baseName
-                              + "_" + std::to_string( xPartitions )
-                              + "x" + std::to_string( yPartitions )
-                              + "x" + std::to_string( zPartitions );
+  std::string const xmlName = "test_hydrostatic_" + baseName + "_" + std::to_string( n_partitions );
   std::string const xmlPath = testBinaryDir + "/" + xmlName + ".xml";
   // Only rank 0 writes the XML file; all ranks then barrier before reading it.
   if( MpiWrapper::commRank( MPI_COMM_GEOS ) == 0 )
@@ -309,9 +303,9 @@ TEST_P( MixedDimHydrostaticEquilibriumTest, Run )
     auto options = std::make_unique< CommandLineOptions >( g_commandLineOptions );
     options->inputFileNames.push_back( xmlPath );
     options->problemName = xmlName;
-    options->xPartitionsOverride = xPartitions;
-    options->yPartitionsOverride = yPartitions;
-    options->zPartitionsOverride = zPartitions;
+    options->xPartitionsOverride = n_partitions;
+    options->yPartitionsOverride = 0;
+    options->zPartitionsOverride = 0;
     options->overridePartitionNumbers = true;
 
     GeosxState state( std::move( options ) );

src/coreComponents/integrationTests/fluidFlowTests/testMixedDimHydrostaticEquilibrium_mpi.cpp

@@ -33,38 +33,15 @@ INSTANTIATE_TEST_SUITE_P(
   MixedDimHydrostaticEquilibriumTest,
   ::testing::Combine(
     ::testing::Values(
-      // Flat tet meshes - single and triple fracture only
-      "fractured_mesh_tet_DFN_1.vtu",
-      "fractured_mesh_tet_DFN_123.vtu",
-
-      // Wavy tet meshes - single and triple fracture only
-      "fractured_wavy_mesh_tet_DFN_1.vtu",
-      "fractured_wavy_mesh_tet_DFN_123.vtu",
-
       // Full span hex meshes - single and triple fracture only
       "fractured_full_span_mesh_hex_DFN_1.vtu",
       "fractured_full_span_mesh_hex_DFN_123.vtu",
 
       // Full span tet meshes - single and triple fracture only
       "fractured_full_span_mesh_tet_DFN_1.vtu",
-      "fractured_full_span_mesh_tet_DFN_123.vtu",
-
-      // T-shaped wavy meshes
-      "t_shaped_wavy_mesh_hex_DFN_t1t2.vtu",
-      "t_shaped_wavy_mesh_tet_DFN_t1t2.vtu",
-
-      // Y-shaped wavy meshes
-      "y_shaped_wavy_mesh_hex_DFN_y1y2y3.vtu",
-      "y_shaped_wavy_mesh_tet_DFN_y1y2y3.vtu",
-
-      // 5-fracture DFN market meshes
-      "DFN_5_fractures_hex_binarized.vtu",
-      "DFN_5_fractures_tet_binarized.vtu"
+      "fractured_full_span_mesh_tet_DFN_123.vtu"
       ),
-    ::testing::Values(
-      std::make_tuple( 1, 4, 1 ),
-      std::make_tuple( 2, 1, 2 )
-      )
+    ::testing::Values( 4 )
     )
   );
 

src/coreComponents/integrationTests/fluidFlowTests/testMixedDimSinglePhaseFlow.cpp

@@ -33,9 +33,7 @@ INSTANTIATE_TEST_SUITE_P(
       "fractured_mesh_hex_DFN_123.vtu"
       ),
     ::testing::Bool(),
-    ::testing::Values(
-      std::make_tuple( 1, 1, 1 )
-      )
+    ::testing::Values( 1 )
     )
   );
 

src/coreComponents/integrationTests/fluidFlowTests/testMixedDimSinglePhaseFlowFixture.hpp

@@ -56,7 +56,7 @@ extern CommandLineOptions g_commandLineOptions;
  * - bool: Whether to run the flow solver (true) or just check initialization (false).
  * - tuple<int, int, int>: Number of partitions in the x, y, z directions.
  */
-class MixedDimSinglePhaseFlowTest : public ::testing::TestWithParam< std::tuple< std::string, bool, std::tuple< int, int, int > > >
+class MixedDimSinglePhaseFlowTest : public ::testing::TestWithParam< std::tuple< std::string, bool, int > >
 {
 protected:
   void SetUp() override
@@ -181,10 +181,7 @@ TEST_P( MixedDimSinglePhaseFlowTest, Run )
 {
   std::string const & meshFileName = std::get< 0 >( GetParam() );
   bool const runSolver = std::get< 1 >( GetParam() );
-  std::tuple< int, int, int > const & partitions = std::get< 2 >( GetParam() );
-  int const xPartitions = std::get< 0 >( partitions );
-  int const yPartitions = std::get< 1 >( partitions );
-  int const zPartitions = std::get< 2 >( partitions );
+  int const n_partitions = std::get< 2 >( GetParam() );
 
   std::string xmlPath = testBinaryDir + "/test_mixed_dim_single_phase_flow.xml";
 
@@ -214,9 +211,9 @@ TEST_P( MixedDimSinglePhaseFlowTest, Run )
     options->inputFileNames.push_back( xmlPath );
     options->problemName = "test_mixed_dim_single_phase_flow";
 
-    options->xPartitionsOverride = xPartitions;
-    options->yPartitionsOverride = yPartitions;
-    options->zPartitionsOverride = zPartitions;
+    options->xPartitionsOverride = n_partitions;
+    options->yPartitionsOverride = 0;
+    options->zPartitionsOverride = 0;
     options->overridePartitionNumbers = true;
 
     GeosxState state( std::move( options ) );

src/coreComponents/integrationTests/fluidFlowTests/testMixedDimSinglePhaseFlow_mpi.cpp

@@ -28,7 +28,7 @@ CommandLineOptions g_commandLineOptions;
  * Combinations:
  * - Meshes: 1 fracture, 2 fractures, 3 fractures
  * - Solver: Run solver (true) vs Initialization only (false)
- * - Partitioning: 1x1x4, 1x2x2 (4 ranks)
+ * - Partitioning: 4 ranks
  */
 INSTANTIATE_TEST_SUITE_P(
   MixedDimPartitionedFlowCases,
@@ -38,10 +38,7 @@ INSTANTIATE_TEST_SUITE_P(
       "fractured_mesh_hex_DFN_123.vtu"
       ),
     ::testing::Bool(),
-    ::testing::Values(
-      std::make_tuple( 1, 1, 4 ),
-      std::make_tuple( 1, 2, 2 )
-      )
+    ::testing::Values( 4 )
     )
   );
 

src/coreComponents/integrationTests/mechanicTest/testFEMConsistency.cpp

@@ -23,16 +23,16 @@
 CommandLineOptions g_commandLineOptions;
 
 /**
- * @brief Serial execution test cases (single rank, partition 1x1x1).
+ * @brief Serial execution test cases (single rank, partition number 1).
  */
 INSTANTIATE_TEST_SUITE_P(
   FractureStressCasesSerial,
   ConsistencyTest,
   ::testing::Combine(
     ::testing::Values(
-      // Flat hex meshes - single and triple fracture only
-      "fractured_mesh_hex_DFN_1.vtu",
-      "fractured_mesh_hex_DFN_123.vtu",
+      // Wavy tet meshes - single and triple fracture only
+      "fractured_wavy_mesh_tet_DFN_1.vtu",
+      "fractured_wavy_mesh_tet_DFN_123.vtu",
 
       // Wavy hex meshes - single and triple fracture only
       "fractured_wavy_mesh_hex_DFN_1.vtu",
@@ -61,9 +61,7 @@ INSTANTIATE_TEST_SUITE_P(
     ::testing::Values( -1.0e6 ),     // s_xx
     ::testing::Values( -0.5e6 ),     // s_yy
     ::testing::Values( -2.0e6 ),     // s_zz
-    ::testing::Values(
-      std::make_tuple( 1, 1, 1 )
-      )
+    ::testing::Values( 1 )
     )
   );
 

src/coreComponents/integrationTests/mechanicTest/testFEMConsistencyFixture.hpp

@@ -73,7 +73,7 @@ extern CommandLineOptions g_commandLineOptions;
  * - real64: Applied stress component ZZ (s_zz).
  * - tuple<int, int, int>: Number of partitions in the x, y, z directions.
  */
-class ConsistencyTest : public ::testing::TestWithParam< std::tuple< std::string, real64, real64, real64, std::tuple< int, int, int > > >
+class ConsistencyTest : public ::testing::TestWithParam< std::tuple< std::string, real64, real64, real64, int > >
 {
 protected:
   void SetUp() override
@@ -111,8 +111,8 @@ class ConsistencyTest : public ::testing::TestWithParam< std::tuple< std::string
       tolJumpT="1.e-4" 
       discretization="FE1" 
       targetRegions="{ Region, Fracture }" 
-      logLevel="1">
-      <NonlinearSolverParameters newtonTol="1.0e-7" newtonMaxIter="20" logLevel="1"/>
+      logLevel="0">
+      <NonlinearSolverParameters newtonTol="1.0e-7" newtonMaxIter="20" logLevel="0"/>
       <LinearSolverParameters directParallel="0"/>
     </SolidMechanicsAugmentedLagrangianContact>
     <SurfaceGenerator name="SurfaceGen" targetRegions="{ Region, Fracture }" fractureRegion="Fracture" initialRockToughness="10.0e9"/>
@@ -172,10 +172,7 @@ TEST_P( ConsistencyTest, Run )
   real64 const s_xx = std::get< 1 >( params );
   real64 const s_yy = std::get< 2 >( params );
   real64 const s_zz = std::get< 3 >( params );
-  std::tuple< int, int, int > const & partitions = std::get< 4 >( params );
-  int const xPartitions = std::get< 0 >( partitions );
-  int const yPartitions = std::get< 1 >( partitions );
-  int const zPartitions = std::get< 2 >( partitions );
+  int const nPartitions = std::get< 4 >( params );
   std::string xmlPath = testBinaryDir + "/test_fem_consistency_" + meshFileName + ".xml";
 
   std::string nodeSetNames = "{ f1_node_set }";
@@ -229,9 +226,9 @@ TEST_P( ConsistencyTest, Run )
   options->inputFileNames.push_back( xmlPath );
   options->problemName = "test_fem_consistency";
 
-  options->xPartitionsOverride = xPartitions;
-  options->yPartitionsOverride = yPartitions;
-  options->zPartitionsOverride = zPartitions;
+  options->xPartitionsOverride = nPartitions;
+  options->yPartitionsOverride = 0;
+  options->zPartitionsOverride = 0;
   options->overridePartitionNumbers = true;
 
   // Scoped state to ensure cleanup

src/coreComponents/integrationTests/mechanicTest/testFEMConsistency_mpi.cpp

@@ -30,7 +30,7 @@ CommandLineOptions g_commandLineOptions;
  * 2. s_xx (real64): Applied stress component XX.
  * 3. s_yy (real64): Applied stress component YY.
  * 4. s_zz (real64): Applied stress component ZZ.
- * 5. Partitioning (tuple<int, int, int>): Number of partitions in x, y, z directions.
+ * 5. Partitioning int: Number of partitions 4.
  *
  * The test cases cover a subset of meshes with various partitioning schemes for 4 ranks.
  */
@@ -43,37 +43,14 @@ INSTANTIATE_TEST_SUITE_P(
       "fractured_mesh_hex_DFN_1.vtu",
       "fractured_mesh_hex_DFN_123.vtu",
 
-      // Wavy hex meshes - single and triple fracture only
-      "fractured_wavy_mesh_hex_DFN_1.vtu",
-      "fractured_wavy_mesh_hex_DFN_123.vtu",
-
-      // Full span hex meshes - single and triple fracture only
-      "fractured_full_span_mesh_hex_DFN_1.vtu",
-      "fractured_full_span_mesh_hex_DFN_123.vtu",
-
       // Full span tet meshes - single and triple fracture only
       "fractured_full_span_mesh_tet_DFN_1.vtu",
-      "fractured_full_span_mesh_tet_DFN_123.vtu",
-
-      // T-shaped wavy meshes
-      "t_shaped_wavy_mesh_hex_DFN_t1t2.vtu",
-      "t_shaped_wavy_mesh_tet_DFN_t1t2.vtu",
-
-      // Y-shaped wavy meshes
-      "y_shaped_wavy_mesh_hex_DFN_y1y2y3.vtu",
-      "y_shaped_wavy_mesh_tet_DFN_y1y2y3.vtu",
-
-      // 5-fracture DFN market meshes
-      "DFN_5_fractures_hex_binarized.vtu",
-      "DFN_5_fractures_tet_binarized.vtu"
+      "fractured_full_span_mesh_tet_DFN_123.vtu"
       ),
     ::testing::Values( -1.0e6 ),     // s_xx
     ::testing::Values( -0.5e6 ),     // s_yy
     ::testing::Values( -2.0e6 ),     // s_zz
-    ::testing::Values(
-      std::make_tuple( 2, 2, 1 ),
-      std::make_tuple( 4, 1, 1 )
-      )
+    ::testing::Values( 4 )
     )
   );
 

src/coreComponents/integrationTests/meshTests/dfn_market/fractured_full_span_mesh_hex_DFN_1.vtu

@@ -1,15 +1,17 @@
-<?xml version='1.0' encoding='utf-8'?>
+<?xml version="1.0"?>
 <VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
   <UnstructuredGrid>
-    <Piece NumberOfPoints="27" NumberOfCells="8">
-      <PointData>
-        <DataArray type="UInt8" Name="f1_node_set" format="binary" RangeMin="0" RangeMax="1">AQAAAACAAAAbAAAADQAAAA==eJxjYGRgwIEAAJkACg==</DataArray>
+    <Piece NumberOfPoints="125" NumberOfCells="64">
+      <PointData Scalars="f1_node_set">
+        <DataArray type="UInt8" Name="f1_node_set" format="binary" RangeMin="0" RangeMax="1">
+          AQAAAACAAAB9AAAADwAAAA==eJxjYGBkAIKBIQAGpAAa
+        </DataArray>
       </PointData>
       <CellData>
       </CellData>
       <Points>
         <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="1.7320508075688772">
-          AQAAAACAAACIAgAATwAAAA==eJyF0cENACAIQ1E2c/8tHIERPHFBn3Ah/tjQQsSv9nrzBJeu3p0nuObUv84TXL7kT/OkV++5ptzS9VzTnjRH+6qafMvHdAfp734AQnk+Ow==
+          AQAAAACAAAC4CwAADgEAAA==eJyN1bkVwzAMA1Bvlv27jKBSpUZIxQbJD8RGz/ABiwSg5/lX79dvfAHfwA9w8cx14gv4Bn6A5/eSN/EFfAM/wPO95E18Ad/AD/C8n7yJL+Ab+AGepbmq7+qL/rvpSGvqKnkTT101nYondZW8iaeumg+0j9RVPpd46qr5TH1KXSVv4qmr5mPNQX6W3+SHqaYbzVV9V19aDolPa+oqeW9zTjoVT+oqeW9zVD7QPlJX+dxtTstn6pPyWnk61XKnzUF+lt+mmi7bnDVvzWOq/XfTkdbUVTsnlXPSqXhSV+0cVo7KB9qHzuOpdm41n6lPyuuplmvNx5qD/DzVdN9yQnPWvKdaX1oOie97/QAgkVK8
           <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
             <Value index="0">
               0
@@ -21,16 +23,16 @@
         </DataArray>
       </Points>
       <Cells>
-        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="26">
-          AQAAAACAAAAAAgAAZwAAAA==eJyFjzcOwDAMA504vTen/f+hGXJaCBjWcouOBJ37L4MF9LCFHRzhIF4OS8kxr4eT5Hj5r2ElPfa/wFk8620kxzzrXSUntW+DOzxhEC+2z7wDXpKT2hfk/4G3eLF95lnvKzkfjh4DQQ==
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="124">
+          AQAAAACAAAAAEAAAUgIAAA==eJyFk0VXEGAUBbE7EMEWAxUUsTsAuxsD7O5u/fEuvLOZc77z3uauZs59i9vR8f8mJacnpyWXJZcne5OrxU1OzpAHbkVyjTxwU5Iz5YFbmVwrD9zU5Cx54FYl18nDn/Sdk5yd5E/69iXXi6PvXHng6LtBHjj6zpMHjr4b5YHj7/nywPH3Jnn4k74LkwuS/EnfgWS/OPp2ygNH383ywNF3kTxw9N0iDxx/d8kDx9+D8vAnfbuTi5P8Sd+h5FZx9O2RB46+2+SBo+8SeeDou10eOP5eKg8cf++Qp9r3zuSu5P7kPnGtfcPtTh6Qp9o33J7kQXmqfcPtTR6Sp9o3f9L3SPKwOPp633D0PSoPXGvfcPQdlqfaNxx/j8hT7Zs/6XssOSqutW84+h6Xp9o3HH1PyFPtG24keVKeat/8Sd/TyVPiWvuGo+8Zeap9w9H3rDzVvuH4+5w81b7PJy8kryaviGvtG+5i8po81b7hLiWvy1PtG+5y8oY81b75k743k2PiWvuGo+8teap9w9H3tjzVvuH4+4481b75k74TyXFxrX3D0feuPNW+4eh7T55q33D8fV+eat/8Sd+HyQfiWvuGo+8jeap9w9H3sTzVvuH4+4k81b6fJp8lXydfiWvtG+558o081b7hXiTfylPtG+5l8p081b75k74fku/FtfYNR9+P8lT7hqPvJ3mqfcPx92d5qn3zJ32/Jr+Ia+0bjr7f5Kn2DUff7/JU+4bj7x/yVPvmT/r+Sv4U19o3HH1/y1PtG46+f+Sp9g3H33/l+QdSt3wB
         </DataArray>
-        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="8" RangeMax="64">
-          AQAAAACAAABAAAAAHgAAAA==eJzjYIAAASgtAaUVoLQGlDaA0hZQ2gFKAwAeQAEh
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="8" RangeMax="512">
+          AQAAAACAAAAAAgAAegAAAA==eJwtxREAgwAAALB+CsPwGIbhMQzDY3gMwzAMj2EYhmEYhmEYHtpkYXCLHPvlxKkzv527cOmPK39du3Hrzr1/Hjx68uzFqzfvPnz6cvC4Cx059suJU2d+O3fh0h9X/rp249ade/88ePTk2YtXb959+PTl4Hn3BzUDHyM=
         </DataArray>
         <DataArray type="UInt8" Name="types" format="binary" RangeMin="12" RangeMax="12">
-          AQAAAACAAAAIAAAACwAAAA==eJzj4YEAAAG4AGE=
+          AQAAAACAAABAAAAADAAAAA==eJzj4aEMAABhwAMB
         </DataArray>
       </Cells>
     </Piece>
   </UnstructuredGrid>
-</VTKFile>
+</VTKFile>

src/coreComponents/integrationTests/meshTests/dfn_market/fractured_full_span_mesh_hex_DFN_12.vtu

@@ -1,16 +1,20 @@
-<?xml version='1.0' encoding='utf-8'?>
+<?xml version="1.0"?>
 <VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
   <UnstructuredGrid>
-    <Piece NumberOfPoints="27" NumberOfCells="8">
-      <PointData>
-        <DataArray type="UInt8" Name="f1_node_set" format="binary" RangeMin="0" RangeMax="1">AQAAAACAAAAbAAAADQAAAA==eJxjYGRgwIEAAJkACg==</DataArray>
-        <DataArray type="UInt8" Name="f2_node_set" format="binary" RangeMin="0" RangeMax="1">AQAAAACAAAAbAAAAEgAAAA==eJxjYGBgZGRkAAM0BgAAmQAK</DataArray>
+    <Piece NumberOfPoints="125" NumberOfCells="64">
+      <PointData Scalars="f1_node_set">
+        <DataArray type="UInt8" Name="f1_node_set" format="binary" RangeMin="0" RangeMax="1">
+          AQAAAACAAAB9AAAADwAAAA==eJxjYGBkAIKBIQAGpAAa
+        </DataArray>
+        <DataArray type="UInt8" Name="f2_node_set" format="binary" RangeMin="0" RangeMax="1">
+          AQAAAACAAAB9AAAAEwAAAA==eJxjYIABRhBgwAJoJwEABqQAGg==
+        </DataArray>
       </PointData>
       <CellData>
       </CellData>
       <Points>
         <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="1.7320508075688772">
-          AQAAAACAAACIAgAATwAAAA==eJyF0cENACAIQ1E2c/8tHIERPHFBn3Ah/tjQQsSv9nrzBJeu3p0nuObUv84TXL7kT/OkV++5ptzS9VzTnjRH+6qafMvHdAfp734AQnk+Ow==
+          AQAAAACAAAC4CwAADgEAAA==eJyN1bkVwzAMA1Bvlv27jKBSpUZIxQbJD8RGz/ABiwSg5/lX79dvfAHfwA9w8cx14gv4Bn6A5/eSN/EFfAM/wPO95E18Ad/AD/C8n7yJL+Ab+AGepbmq7+qL/rvpSGvqKnkTT101nYondZW8iaeumg+0j9RVPpd46qr5TH1KXSVv4qmr5mPNQX6W3+SHqaYbzVV9V19aDolPa+oqeW9zTjoVT+oqeW9zVD7QPlJX+dxtTstn6pPyWnk61XKnzUF+lt+mmi7bnDVvzWOq/XfTkdbUVTsnlXPSqXhSV+0cVo7KB9qHzuOpdm41n6lPyuuplmvNx5qD/DzVdN9yQnPWvKdaX1oOie97/QAgkVK8
           <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
             <Value index="0">
               0
@@ -22,16 +26,16 @@
         </DataArray>
       </Points>
       <Cells>
-        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="26">
-          AQAAAACAAAAAAgAAZwAAAA==eJyFjzcOwDAMA504vTen/f+hGXJaCBjWcouOBJ37L4MF9LCFHRzhIF4OS8kxr4eT5Hj5r2ElPfa/wFk8620kxzzrXSUntW+DOzxhEC+2z7wDXpKT2hfk/4G3eLF95lnvKzkfjh4DQQ==
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="124">
+          AQAAAACAAAAAEAAAUgIAAA==eJyFk0VXEGAUBbE7EMEWAxUUsTsAuxsD7O5u/fEuvLOZc77z3uauZs59i9vR8f8mJacnpyWXJZcne5OrxU1OzpAHbkVyjTxwU5Iz5YFbmVwrD9zU5Cx54FYl18nDn/Sdk5yd5E/69iXXi6PvXHng6LtBHjj6zpMHjr4b5YHj7/nywPH3Jnn4k74LkwuS/EnfgWS/OPp2ygNH383ywNF3kTxw9N0iDxx/d8kDx9+D8vAnfbuTi5P8Sd+h5FZx9O2RB46+2+SBo+8SeeDou10eOP5eKg8cf++Qp9r3zuSu5P7kPnGtfcPtTh6Qp9o33J7kQXmqfcPtTR6Sp9o3f9L3SPKwOPp633D0PSoPXGvfcPQdlqfaNxx/j8hT7Zs/6XssOSqutW84+h6Xp9o3HH1PyFPtG24keVKeat/8Sd/TyVPiWvuGo+8Zeap9w9H3rDzVvuH4+5w81b7PJy8kryaviGvtG+5i8po81b7hLiWvy1PtG+5y8oY81b75k743k2PiWvuGo+8teap9w9H3tjzVvuH4+4481b75k74TyXFxrX3D0feuPNW+4eh7T55q33D8fV+eat/8Sd+HyQfiWvuGo+8jeap9w9H3sTzVvuH4+4k81b6fJp8lXydfiWvtG+558o081b7hXiTfylPtG+5l8p081b75k74fku/FtfYNR9+P8lT7hqPvJ3mqfcPx92d5qn3zJ32/Jr+Ia+0bjr7f5Kn2DUff7/JU+4bj7x/yVPvmT/r+Sv4U19o3HH1/y1PtG46+f+Sp9g3H33/l+QdSt3wB
         </DataArray>
-        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="8" RangeMax="64">
-          AQAAAACAAABAAAAAHgAAAA==eJzjYIAAASgtAaUVoLQGlDaA0hZQ2gFKAwAeQAEh
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="8" RangeMax="512">
+          AQAAAACAAAAAAgAAegAAAA==eJwtxREAgwAAALB+CsPwGIbhMQzDY3gMwzAMj2EYhmEYhmEYHtpkYXCLHPvlxKkzv527cOmPK39du3Hrzr1/Hjx68uzFqzfvPnz6cvC4Cx059suJU2d+O3fh0h9X/rp249ade/88ePTk2YtXb959+PTl4Hn3BzUDHyM=
         </DataArray>
         <DataArray type="UInt8" Name="types" format="binary" RangeMin="12" RangeMax="12">
-          AQAAAACAAAAIAAAACwAAAA==eJzj4YEAAAG4AGE=
+          AQAAAACAAABAAAAADAAAAA==eJzj4aEMAABhwAMB
         </DataArray>
       </Cells>
     </Piece>
   </UnstructuredGrid>
-</VTKFile>
+</VTKFile>