[Bbar#M] A 3D test

ProcessLib/SmallDeformation/Tests.cmake

@@ -303,6 +303,21 @@ AddTest(
     axisymmetric_bbar_ts_1_t_1.000000.vtu axisymmetric_bbar_ts_1_t_1.000000.vtu analytic_eps epsilon 1.e-9 5.e-7
 )
 
+AddTest(
+    NAME Mechanics_Simple3DBbar
+    PATH Mechanics/Simple3DBbar
+    EXECUTABLE ogs
+    EXECUTABLE_ARGS simple_3d_bbar.prj
+    WRAPPER time
+    TESTER vtkdiff
+    REQUIREMENTS NOT OGS_USE_MPI
+    RUNTIME 1
+    DIFF_DATA
+    simple_3d_bbar_ts_1_t_1.000000.vtu simple_3d_bbar_ts_1_t_1.000000.vtu analytic_displacement displacement 1.e-10 1.e-10
+    simple_3d_bbar_ts_1_t_1.000000.vtu simple_3d_bbar_ts_1_t_1.000000.vtu analytic_sigma sigma 6.e-5 5.e-7
+    simple_3d_bbar_ts_1_t_1.000000.vtu simple_3d_bbar_ts_1_t_1.000000.vtu analytic_eps epsilon 1.e-9 5.e-7
+)
+
 if(NOT OGS_USE_PETSC)
     NotebookTest(NOTEBOOKFILE Mechanics/CooksMembrane/CooksMembraneBbar.ipynb RUNTIME 1)
     NotebookTest(NOTEBOOKFILE Mechanics/Linear/SimpleMechanics.ipynb RUNTIME 5)

Tests/Data/Mechanics/Simple3DBbar/analytical_solution.py

+import pyvista as pv
+
+mesh = pv.read("simple_3d_bbar_ts_1_t_1.000000.vtu")
+E = 1.0e10
+nv = 0.499
+p = 1e7
+
+sigma_analytical = mesh.point_data["sigma"]
+sigma_analytical[:, 0] = 0.0
+sigma_analytical[:, 1] = 0.0
+sigma_analytical[:, 2] = -p
+sigma_analytical[:, 3] = 0.0
+sigma_analytical[:, 4] = 0.0
+sigma_analytical[:, 5] = 0.0
+mesh.point_data["analytic_sigma"] = sigma_analytical
+
+eps_analytical = mesh.point_data["epsilon"]
+eps_1 = -p / E
+eps_0 = -eps_1 * nv
+eps_analytical[:, 0] = eps_0
+eps_analytical[:, 1] = eps_0
+eps_analytical[:, 2] = eps_1
+eps_analytical[:, 3] = 0.0
+eps_analytical[:, 4] = 0.0
+eps_analytical[:, 5] = 0.0
+mesh.point_data["analytic_eps"] = eps_analytical
+
+
+u_analytical = mesh.point_data["displacement"]
+u_analytical[:, 0] = eps_0 * mesh.points[:, 0]
+u_analytical[:, 1] = eps_0 * mesh.points[:, 1]
+u_analytical[:, 2] = eps_1 * mesh.points[:, 2]
+mesh.point_data["analytic_displacement"] = u_analytical
+
+mesh.save("simple_3d_bbar_ts_1_t_1.000000.vtu")
+# mesh.plot(scalars="analytic_displacement")

Tests/Data/Mechanics/Simple3DBbar/bottom.vtu

+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="63" NumberOfCells="40">
+      <PointData>
+        <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="0" RangeMax="62">
+          +AEAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAKQAAAAAAAAAqAAAAAAAAACsAAAAAAAAALAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAzAAAAAAAAADQAAAAAAAAANQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA=
+        </DataArray>
+      </PointData>
+      <CellData>
+        <DataArray type="Int32" Name="MaterialIDs" format="binary" RangeMin="0" RangeMax="0">
+          oAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
+        </DataArray>
+        <DataArray type="UInt64" Name="bulk_element_ids" format="binary" RangeMin="0" RangeMax="39">
+          QAEAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAAA==
+        </DataArray>
+        <DataArray type="UInt64" Name="bulk_face_ids" format="binary" RangeMin="0" RangeMax="0">
+          QAEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA==
+        </DataArray>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="1.004987562112089">
+          6AUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACamZmZmZmpPwAAAAAAAAAAAAAAAAAAAACamZmZmZm5PwAAAAAAAAAAAAAAAAAAAAAzMzMzMzPDPwAAAAAAAAAAAAAAAAAAAACamZmZmZnJPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADQPwAAAAAAAAAAAAAAAAAAAAAzMzMzMzPTPwAAAAAAAAAAAAAAAAAAAABmZmZmZmbWPwAAAAAAAAAAAAAAAAAAAACamZmZmZnZPwAAAAAAAAAAAAAAAAAAAADNzMzMzMzcPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADgPwAAAAAAAAAAAAAAAAAAAACamZmZmZnhPwAAAAAAAAAAAAAAAAAAAAAzMzMzMzPjPwAAAAAAAAAAAAAAAAAAAADNzMzMzMzkPwAAAAAAAAAAAAAAAAAAAABmZmZmZmbmPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADoPwAAAAAAAAAAAAAAAAAAAACamZmZmZnpPwAAAAAAAAAAAAAAAAAAAAAzMzMzMzPrPwAAAAAAAAAAAAAAAAAAAADNzMzMzMzsPwAAAAAAAAAAAAAAAAAAAABmZmZmZmbuPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADwPwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAJqZmZmZmak/AAAAAAAAAACamZmZmZmpP5qZmZmZmak/AAAAAAAAAACamZmZmZm5P5qZmZmZmak/AAAAAAAAAAAzMzMzMzPDP5qZmZmZmak/AAAAAAAAAACamZmZmZnJP5qZmZmZmak/AAAAAAAAAAAAAAAAAADQP5qZmZmZmak/AAAAAAAAAAAzMzMzMzPTP5qZmZmZmak/AAAAAAAAAABmZmZmZmbWP5qZmZmZmak/AAAAAAAAAACamZmZmZnZP5qZmZmZmak/AAAAAAAAAADNzMzMzMzcP5qZmZmZmak/AAAAAAAAAAAAAAAAAADgP5qZmZmZmak/AAAAAAAAAACamZmZmZnhP5qZmZmZmak/AAAAAAAAAAAzMzMzMzPjP5qZmZmZmak/AAAAAAAAAADNzMzMzMzkP5qZmZmZmak/AAAAAAAAAABmZmZmZmbmP5qZmZmZmak/AAAAAAAAAAAAAAAAAADoP5qZmZmZmak/AAAAAAAAAACamZmZmZnpP5qZmZmZmak/AAAAAAAAAAAzMzMzMzPrP5qZmZmZmak/AAAAAAAAAADNzMzMzMzsP5qZmZmZmak/AAAAAAAAAABmZmZmZmbuP5qZmZmZmak/AAAAAAAAAAAAAAAAAADwP5qZmZmZmak/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmbk/AAAAAAAAAACamZmZmZmpP5qZmZmZmbk/AAAAAAAAAACamZmZmZm5P5qZmZmZmbk/AAAAAAAAAAAzMzMzMzPDP5qZmZmZmbk/AAAAAAAAAACamZmZmZnJP5qZmZmZmbk/AAAAAAAAAAAAAAAAAADQP5qZmZmZmbk/AAAAAAAAAAAzMzMzMzPTP5qZmZmZmbk/AAAAAAAAAABmZmZmZmbWP5qZmZmZmbk/AAAAAAAAAACamZmZmZnZP5qZmZmZmbk/AAAAAAAAAADNzMzMzMzcP5qZmZmZmbk/AAAAAAAAAAAAAAAAAADgP5qZmZmZmbk/AAAAAAAAAACamZmZmZnhP5qZmZmZmbk/AAAAAAAAAAAzMzMzMzPjP5qZmZmZmbk/AAAAAAAAAADNzMzMzMzkP5qZmZmZmbk/AAAAAAAAAABmZmZmZmbmP5qZmZmZmbk/AAAAAAAAAAAAAAAAAADoP5qZmZmZmbk/AAAAAAAAAACamZmZmZnpP5qZmZmZmbk/AAAAAAAAAAAzMzMzMzPrP5qZmZmZmbk/AAAAAAAAAADNzMzMzMzsP5qZmZmZmbk/AAAAAAAAAABmZmZmZmbuP5qZmZmZmbk/AAAAAAAAAAAAAAAAAADwP5qZmZmZmbk/AAAAAAAAAAA=
+          <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+            <Value index="0">
+              0
+            </Value>
+            <Value index="1">
+              1.0049875621
+            </Value>
+          </InformationKey>
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="62">
+          AAUAAAAAAAAAAAAAAAAAABUAAAAAAAAAFgAAAAAAAAABAAAAAAAAAAEAAAAAAAAAFgAAAAAAAAAXAAAAAAAAAAIAAAAAAAAAAgAAAAAAAAAXAAAAAAAAABgAAAAAAAAAAwAAAAAAAAADAAAAAAAAABgAAAAAAAAAGQAAAAAAAAAEAAAAAAAAAAQAAAAAAAAAGQAAAAAAAAAaAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAaAAAAAAAAABsAAAAAAAAABgAAAAAAAAAGAAAAAAAAABsAAAAAAAAAHAAAAAAAAAAHAAAAAAAAAAcAAAAAAAAAHAAAAAAAAAAdAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAdAAAAAAAAAB4AAAAAAAAACQAAAAAAAAAJAAAAAAAAAB4AAAAAAAAAHwAAAAAAAAAKAAAAAAAAAAoAAAAAAAAAHwAAAAAAAAAgAAAAAAAAAAsAAAAAAAAACwAAAAAAAAAgAAAAAAAAACEAAAAAAAAADAAAAAAAAAAMAAAAAAAAACEAAAAAAAAAIgAAAAAAAAANAAAAAAAAAA0AAAAAAAAAIgAAAAAAAAAjAAAAAAAAAA4AAAAAAAAADgAAAAAAAAAjAAAAAAAAACQAAAAAAAAADwAAAAAAAAAPAAAAAAAAACQAAAAAAAAAJQAAAAAAAAAQAAAAAAAAABAAAAAAAAAAJQAAAAAAAAAmAAAAAAAAABEAAAAAAAAAEQAAAAAAAAAmAAAAAAAAACcAAAAAAAAAEgAAAAAAAAASAAAAAAAAACcAAAAAAAAAKAAAAAAAAAATAAAAAAAAABMAAAAAAAAAKAAAAAAAAAApAAAAAAAAABQAAAAAAAAAFQAAAAAAAAAqAAAAAAAAACsAAAAAAAAAFgAAAAAAAAAWAAAAAAAAACsAAAAAAAAALAAAAAAAAAAXAAAAAAAAABcAAAAAAAAALAAAAAAAAAAtAAAAAAAAABgAAAAAAAAAGAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAAGQAAAAAAAAAZAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAaAAAAAAAAABoAAAAAAAAALwAAAAAAAAAwAAAAAAAAABsAAAAAAAAAGwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAHAAAAAAAAAAcAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAdAAAAAAAAAB0AAAAAAAAAMgAAAAAAAAAzAAAAAAAAAB4AAAAAAAAAHgAAAAAAAAAzAAAAAAAAADQAAAAAAAAAHwAAAAAAAAAfAAAAAAAAADQAAAAAAAAANQAAAAAAAAAgAAAAAAAAACAAAAAAAAAANQAAAAAAAAA2AAAAAAAAACEAAAAAAAAAIQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAIgAAAAAAAAAiAAAAAAAAADcAAAAAAAAAOAAAAAAAAAAjAAAAAAAAACMAAAAAAAAAOAAAAAAAAAA5AAAAAAAAACQAAAAAAAAAJAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAJQAAAAAAAAAlAAAAAAAAADoAAAAAAAAAOwAAAAAAAAAmAAAAAAAAACYAAAAAAAAAOwAAAAAAAAA8AAAAAAAAACcAAAAAAAAAJwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAKAAAAAAAAAAoAAAAAAAAAD0AAAAAAAAAPgAAAAAAAAApAAAAAAAAAA==
+        </DataArray>
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="4" RangeMax="160">
+          QAEAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAA==
+        </DataArray>
+        <DataArray type="UInt8" Name="types" format="binary" RangeMin="9" RangeMax="9">
+          KAAAAAAAAAAJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJ
+        </DataArray>
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>

Tests/Data/Mechanics/Simple3DBbar/geometry_origin.vtu

+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="1"                    NumberOfCells="1"                   >
+      <PointData>
+        <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0"                    RangeMax="0"                    offset="0"                   />
+      </PointData>
+      <CellData>
+        <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0"                    RangeMax="0"                    offset="24"                  />
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="0"                    offset="48"                  />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="92"                  />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="116"                 />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="140"                 />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _CAAAAAAAAAAAAAAAAAAAAA==CAAAAAAAAAAAAAAAAAAAAA==GAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA=CAAAAAAAAAAAAAAAAAAAAA==CAAAAAAAAAABAAAAAAAAAA==AQAAAAAAAAAB
+  </AppendedData>
+</VTKFile>

Tests/Data/Mechanics/Simple3DBbar/simple_3d_bbar.prj

+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+    <meshes>
+        <mesh>hex_brick.vtu</mesh>
+        <mesh>x0.vtu</mesh>
+        <mesh>y0.vtu</mesh>
+        <mesh>top.vtu</mesh>
+        <mesh>bottom.vtu</mesh>
+    </meshes>    
+    <processes>
+        <process>
+            <name>SD</name>
+            <type>SMALL_DEFORMATION</type>
+            <integration_order>3</integration_order>
+            <use_b_bar>true</use_b_bar>
+            <constitutive_relation>
+                <type>LinearElasticIsotropic</type>
+                <youngs_modulus>E</youngs_modulus>
+                <poissons_ratio>nu</poissons_ratio>
+            </constitutive_relation>
+            <specific_body_force>0 0 0</specific_body_force>
+            <process_variables>
+                <process_variable>displacement</process_variable>
+            </process_variables>
+            <secondary_variables>
+                <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+                <secondary_variable internal_name="sigma" output_name="sigma"/>
+            </secondary_variables>
+        </process>
+    </processes>
+    <time_loop>
+        <processes>
+            <process ref="SD">
+                <nonlinear_solver>basic_newton</nonlinear_solver>
+                <convergence_criterion>
+                    <type>PerComponentDeltaX</type>
+                    <norm_type>NORM2</norm_type>
+                    <reltols> 1.0e-13 1.0e-13 1e-13</reltols>
+                </convergence_criterion>
+                <time_discretization>
+                    <type>BackwardEuler</type>
+                </time_discretization>
+                <time_stepping>
+                    <type>FixedTimeStepping</type>
+                    <t_initial> 0.0 </t_initial>
+                    <t_end>1</t_end>
+                    <timesteps>
+                        <pair>
+                            <repeat>1</repeat>
+                            <delta_t>1</delta_t>
+                        </pair>
+                    </timesteps>
+                </time_stepping>
+            </process>
+        </processes>
+        <output>
+            <type>VTK</type>
+            <prefix>simple_3d_bbar</prefix>
+            <timesteps>
+                <pair>
+                    <repeat> 1 </repeat>
+                    <each_steps> 1 </each_steps>
+                </pair>
+            </timesteps>
+            <variables>
+                <variable>displacement</variable>
+                <variable>epsilon</variable>
+                <variable>sigma</variable>
+            </variables>
+            <suffix>_ts_{:timestep}_t_{:time}</suffix>
+        </output>
+    </time_loop>
+    <media>
+        <medium>
+            <phases>
+                <phase>
+                    <type>Solid</type>
+                    <properties>
+                        <property>
+                            <name>density</name>
+                            <type>Parameter</type>
+                            <parameter_name>rho</parameter_name>
+                        </property>
+                    </properties>
+                </phase>
+            </phases>
+        </medium>
+    </media>
+    <parameters>
+        <parameter>
+            <name>E</name>
+            <type>Constant</type>
+            <value>1.0E+10</value>
+        </parameter>
+        <parameter>
+            <name>nu</name>
+            <type>Constant</type>
+            <value>0.499</value>
+        </parameter>
+        <parameter>
+            <name>rho</name>
+            <type>Constant</type>
+            <value>0.0</value>
+        </parameter>
+        <parameter>
+            <name>zero</name>
+            <type>Constant</type>
+            <value>0.0</value>
+        </parameter>
+        <parameter>
+            <name>u0</name>
+            <type>Constant</type>
+            <values>0 0 0</values>
+        </parameter>
+        <parameter>
+            <name>load</name>
+            <type>Constant</type>
+            <values>-10e+6</values>
+        </parameter>
+    </parameters>
+    <process_variables>
+        <process_variable>
+            <name>displacement</name>
+            <components>3</components>
+            <order>1</order>
+            <initial_condition>u0</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <mesh>bottom</mesh>
+                    <type>Dirichlet</type>
+                    <component>2</component>
+                    <parameter>zero</parameter>
+                </boundary_condition>
+                <boundary_condition>
+                    <mesh>x0</mesh>
+                    <type>Dirichlet</type>
+                    <component>0</component>
+                    <parameter>zero</parameter>
+                </boundary_condition>
+                <boundary_condition>
+                    <mesh>y0</mesh>
+                    <type>Dirichlet</type>
+                    <component>1</component>
+                    <parameter>zero</parameter>
+                </boundary_condition>
+                <boundary_condition>
+                    <mesh>top</mesh>
+                    <type>Neumann</type>
+                    <component>2</component>
+                    <parameter>load</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+        </process_variable>
+    </process_variables>
+    <nonlinear_solvers>
+        <nonlinear_solver>
+            <name>basic_newton</name>
+            <type>Newton</type>
+            <max_iter>10</max_iter>
+            <linear_solver>general_linear_solver</linear_solver>
+        </nonlinear_solver>
+    </nonlinear_solvers>
+    <linear_solvers>
+        <linear_solver>
+            <name>general_linear_solver</name>
+            <eigen>
+                <solver_type>SparseLU</solver_type>
+                <scaling>true</scaling>
+            </eigen>
+        </linear_solver>
+    </linear_solvers>
+</OpenGeoSysProject>

Tests/Data/Mechanics/Simple3DBbar/top.vtu

+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="63" NumberOfCells="40">
+      <PointData>
+        <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="1260" RangeMax="1322">
+          +AEAAAAAAADsBAAAAAAAAO0EAAAAAAAA7gQAAAAAAADvBAAAAAAAAPAEAAAAAAAA8QQAAAAAAADyBAAAAAAAAPMEAAAAAAAA9AQAAAAAAAD1BAAAAAAAAPYEAAAAAAAA9wQAAAAAAAD4BAAAAAAAAPkEAAAAAAAA+gQAAAAAAAD7BAAAAAAAAPwEAAAAAAAA/QQAAAAAAAD+BAAAAAAAAP8EAAAAAAAAAAUAAAAAAAABBQAAAAAAAAIFAAAAAAAAAwUAAAAAAAAEBQAAAAAAAAUFAAAAAAAABgUAAAAAAAAHBQAAAAAAAAgFAAAAAAAACQUAAAAAAAAKBQAAAAAAAAsFAAAAAAAADAUAAAAAAAANBQAAAAAAAA4FAAAAAAAADwUAAAAAAAAQBQAAAAAAABEFAAAAAAAAEgUAAAAAAAATBQAAAAAAABQFAAAAAAAAFQUAAAAAAAAWBQAAAAAAABcFAAAAAAAAGAUAAAAAAAAZBQAAAAAAABoFAAAAAAAAGwUAAAAAAAAcBQAAAAAAAB0FAAAAAAAAHgUAAAAAAAAfBQAAAAAAACAFAAAAAAAAIQUAAAAAAAAiBQAAAAAAACMFAAAAAAAAJAUAAAAAAAAlBQAAAAAAACYFAAAAAAAAJwUAAAAAAAAoBQAAAAAAACkFAAAAAAAAKgUAAAAAAAA=
+        </DataArray>
+      </PointData>
+      <CellData>
+        <DataArray type="Int32" Name="MaterialIDs" format="binary" RangeMin="0" RangeMax="0">
+          oAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
+        </DataArray>
+        <DataArray type="UInt64" Name="bulk_element_ids" format="binary" RangeMin="760" RangeMax="799">
+          QAEAAAAAAAD4AgAAAAAAAPkCAAAAAAAA+gIAAAAAAAD7AgAAAAAAAPwCAAAAAAAA/QIAAAAAAAD+AgAAAAAAAP8CAAAAAAAAAAMAAAAAAAABAwAAAAAAAAIDAAAAAAAAAwMAAAAAAAAEAwAAAAAAAAUDAAAAAAAABgMAAAAAAAAHAwAAAAAAAAgDAAAAAAAACQMAAAAAAAAKAwAAAAAAAAsDAAAAAAAADAMAAAAAAAANAwAAAAAAAA4DAAAAAAAADwMAAAAAAAAQAwAAAAAAABEDAAAAAAAAEgMAAAAAAAATAwAAAAAAABQDAAAAAAAAFQMAAAAAAAAWAwAAAAAAABcDAAAAAAAAGAMAAAAAAAAZAwAAAAAAABoDAAAAAAAAGwMAAAAAAAAcAwAAAAAAAB0DAAAAAAAAHgMAAAAAAAAfAwAAAAAAAA==
+        </DataArray>
+        <DataArray type="UInt64" Name="bulk_face_ids" format="binary" RangeMin="5" RangeMax="5">
+          QAEAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFAAAAAAAAAA==
+        </DataArray>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="1" RangeMax="1.4177446878757824">
+          6AUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA8D+amZmZmZmpPwAAAAAAAAAAAAAAAAAA8D+amZmZmZm5PwAAAAAAAAAAAAAAAAAA8D8zMzMzMzPDPwAAAAAAAAAAAAAAAAAA8D+amZmZmZnJPwAAAAAAAAAAAAAAAAAA8D8AAAAAAADQPwAAAAAAAAAAAAAAAAAA8D8zMzMzMzPTPwAAAAAAAAAAAAAAAAAA8D9mZmZmZmbWPwAAAAAAAAAAAAAAAAAA8D+amZmZmZnZPwAAAAAAAAAAAAAAAAAA8D/NzMzMzMzcPwAAAAAAAAAAAAAAAAAA8D8AAAAAAADgPwAAAAAAAAAAAAAAAAAA8D+amZmZmZnhPwAAAAAAAAAAAAAAAAAA8D8zMzMzMzPjPwAAAAAAAAAAAAAAAAAA8D/NzMzMzMzkPwAAAAAAAAAAAAAAAAAA8D9mZmZmZmbmPwAAAAAAAAAAAAAAAAAA8D8AAAAAAADoPwAAAAAAAAAAAAAAAAAA8D+amZmZmZnpPwAAAAAAAAAAAAAAAAAA8D8zMzMzMzPrPwAAAAAAAAAAAAAAAAAA8D/NzMzMzMzsPwAAAAAAAAAAAAAAAAAA8D9mZmZmZmbuPwAAAAAAAAAAAAAAAAAA8D8AAAAAAADwPwAAAAAAAAAAAAAAAAAA8D8AAAAAAAAAAJqZmZmZmak/AAAAAAAA8D+amZmZmZmpP5qZmZmZmak/AAAAAAAA8D+amZmZmZm5P5qZmZmZmak/AAAAAAAA8D8zMzMzMzPDP5qZmZmZmak/AAAAAAAA8D+amZmZmZnJP5qZmZmZmak/AAAAAAAA8D8AAAAAAADQP5qZmZmZmak/AAAAAAAA8D8zMzMzMzPTP5qZmZmZmak/AAAAAAAA8D9mZmZmZmbWP5qZmZmZmak/AAAAAAAA8D+amZmZmZnZP5qZmZmZmak/AAAAAAAA8D/NzMzMzMzcP5qZmZmZmak/AAAAAAAA8D8AAAAAAADgP5qZmZmZmak/AAAAAAAA8D+amZmZmZnhP5qZmZmZmak/AAAAAAAA8D8zMzMzMzPjP5qZmZmZmak/AAAAAAAA8D/NzMzMzMzkP5qZmZmZmak/AAAAAAAA8D9mZmZmZmbmP5qZmZmZmak/AAAAAAAA8D8AAAAAAADoP5qZmZmZmak/AAAAAAAA8D+amZmZmZnpP5qZmZmZmak/AAAAAAAA8D8zMzMzMzPrP5qZmZmZmak/AAAAAAAA8D/NzMzMzMzsP5qZmZmZmak/AAAAAAAA8D9mZmZmZmbuP5qZmZmZmak/AAAAAAAA8D8AAAAAAADwP5qZmZmZmak/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmbk/AAAAAAAA8D+amZmZmZmpP5qZmZmZmbk/AAAAAAAA8D+amZmZmZm5P5qZmZmZmbk/AAAAAAAA8D8zMzMzMzPDP5qZmZmZmbk/AAAAAAAA8D+amZmZmZnJP5qZmZmZmbk/AAAAAAAA8D8AAAAAAADQP5qZmZmZmbk/AAAAAAAA8D8zMzMzMzPTP5qZmZmZmbk/AAAAAAAA8D9mZmZmZmbWP5qZmZmZmbk/AAAAAAAA8D+amZmZmZnZP5qZmZmZmbk/AAAAAAAA8D/NzMzMzMzcP5qZmZmZmbk/AAAAAAAA8D8AAAAAAADgP5qZmZmZmbk/AAAAAAAA8D+amZmZmZnhP5qZmZmZmbk/AAAAAAAA8D8zMzMzMzPjP5qZmZmZmbk/AAAAAAAA8D/NzMzMzMzkP5qZmZmZmbk/AAAAAAAA8D9mZmZmZmbmP5qZmZmZmbk/AAAAAAAA8D8AAAAAAADoP5qZmZmZmbk/AAAAAAAA8D+amZmZmZnpP5qZmZmZmbk/AAAAAAAA8D8zMzMzMzPrP5qZmZmZmbk/AAAAAAAA8D/NzMzMzMzsP5qZmZmZmbk/AAAAAAAA8D9mZmZmZmbuP5qZmZmZmbk/AAAAAAAA8D8AAAAAAADwP5qZmZmZmbk/AAAAAAAA8D8=
+          <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+            <Value index="0">
+              1
+            </Value>
+            <Value index="1">
+              1.4177446879
+            </Value>
+          </InformationKey>
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="62">
+          AAUAAAAAAAAAAAAAAAAAAAEAAAAAAAAAFgAAAAAAAAAVAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAAXAAAAAAAAABYAAAAAAAAAAgAAAAAAAAADAAAAAAAAABgAAAAAAAAAFwAAAAAAAAADAAAAAAAAAAQAAAAAAAAAGQAAAAAAAAAYAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAaAAAAAAAAABkAAAAAAAAABQAAAAAAAAAGAAAAAAAAABsAAAAAAAAAGgAAAAAAAAAGAAAAAAAAAAcAAAAAAAAAHAAAAAAAAAAbAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAdAAAAAAAAABwAAAAAAAAACAAAAAAAAAAJAAAAAAAAAB4AAAAAAAAAHQAAAAAAAAAJAAAAAAAAAAoAAAAAAAAAHwAAAAAAAAAeAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAgAAAAAAAAAB8AAAAAAAAACwAAAAAAAAAMAAAAAAAAACEAAAAAAAAAIAAAAAAAAAAMAAAAAAAAAA0AAAAAAAAAIgAAAAAAAAAhAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAjAAAAAAAAACIAAAAAAAAADgAAAAAAAAAPAAAAAAAAACQAAAAAAAAAIwAAAAAAAAAPAAAAAAAAABAAAAAAAAAAJQAAAAAAAAAkAAAAAAAAABAAAAAAAAAAEQAAAAAAAAAmAAAAAAAAACUAAAAAAAAAEQAAAAAAAAASAAAAAAAAACcAAAAAAAAAJgAAAAAAAAASAAAAAAAAABMAAAAAAAAAKAAAAAAAAAAnAAAAAAAAABMAAAAAAAAAFAAAAAAAAAApAAAAAAAAACgAAAAAAAAAFQAAAAAAAAAWAAAAAAAAACsAAAAAAAAAKgAAAAAAAAAWAAAAAAAAABcAAAAAAAAALAAAAAAAAAArAAAAAAAAABcAAAAAAAAAGAAAAAAAAAAtAAAAAAAAACwAAAAAAAAAGAAAAAAAAAAZAAAAAAAAAC4AAAAAAAAALQAAAAAAAAAZAAAAAAAAABoAAAAAAAAALwAAAAAAAAAuAAAAAAAAABoAAAAAAAAAGwAAAAAAAAAwAAAAAAAAAC8AAAAAAAAAGwAAAAAAAAAcAAAAAAAAADEAAAAAAAAAMAAAAAAAAAAcAAAAAAAAAB0AAAAAAAAAMgAAAAAAAAAxAAAAAAAAAB0AAAAAAAAAHgAAAAAAAAAzAAAAAAAAADIAAAAAAAAAHgAAAAAAAAAfAAAAAAAAADQAAAAAAAAAMwAAAAAAAAAfAAAAAAAAACAAAAAAAAAANQAAAAAAAAA0AAAAAAAAACAAAAAAAAAAIQAAAAAAAAA2AAAAAAAAADUAAAAAAAAAIQAAAAAAAAAiAAAAAAAAADcAAAAAAAAANgAAAAAAAAAiAAAAAAAAACMAAAAAAAAAOAAAAAAAAAA3AAAAAAAAACMAAAAAAAAAJAAAAAAAAAA5AAAAAAAAADgAAAAAAAAAJAAAAAAAAAAlAAAAAAAAADoAAAAAAAAAOQAAAAAAAAAlAAAAAAAAACYAAAAAAAAAOwAAAAAAAAA6AAAAAAAAACYAAAAAAAAAJwAAAAAAAAA8AAAAAAAAADsAAAAAAAAAJwAAAAAAAAAoAAAAAAAAAD0AAAAAAAAAPAAAAAAAAAAoAAAAAAAAACkAAAAAAAAAPgAAAAAAAAA9AAAAAAAAAA==
+        </DataArray>
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="4" RangeMax="160">
+          QAEAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAA==
+        </DataArray>
+        <DataArray type="UInt8" Name="types" format="binary" RangeMin="9" RangeMax="9">
+          KAAAAAAAAAAJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJ
+        </DataArray>
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>

Tests/Data/Mechanics/Simple3DBbar/x0.vtu

+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="63" NumberOfCells="40">
+      <PointData>
+        <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="0" RangeMax="1302">
+          +AEAAAAAAAAAAAAAAAAAABUAAAAAAAAAKgAAAAAAAAA/AAAAAAAAAFQAAAAAAAAAaQAAAAAAAAB+AAAAAAAAAJMAAAAAAAAAqAAAAAAAAAC9AAAAAAAAANIAAAAAAAAA5wAAAAAAAAD8AAAAAAAAABEBAAAAAAAAJgEAAAAAAAA7AQAAAAAAAFABAAAAAAAAZQEAAAAAAAB6AQAAAAAAAI8BAAAAAAAApAEAAAAAAAC5AQAAAAAAAM4BAAAAAAAA4wEAAAAAAAD4AQAAAAAAAA0CAAAAAAAAIgIAAAAAAAA3AgAAAAAAAEwCAAAAAAAAYQIAAAAAAAB2AgAAAAAAAIsCAAAAAAAAoAIAAAAAAAC1AgAAAAAAAMoCAAAAAAAA3wIAAAAAAAD0AgAAAAAAAAkDAAAAAAAAHgMAAAAAAAAzAwAAAAAAAEgDAAAAAAAAXQMAAAAAAAByAwAAAAAAAIcDAAAAAAAAnAMAAAAAAACxAwAAAAAAAMYDAAAAAAAA2wMAAAAAAADwAwAAAAAAAAUEAAAAAAAAGgQAAAAAAAAvBAAAAAAAAEQEAAAAAAAAWQQAAAAAAABuBAAAAAAAAIMEAAAAAAAAmAQAAAAAAACtBAAAAAAAAMIEAAAAAAAA1wQAAAAAAADsBAAAAAAAAAEFAAAAAAAAFgUAAAAAAAA=
+        </DataArray>
+      </PointData>
+      <CellData>
+        <DataArray type="Int32" Name="MaterialIDs" format="binary" RangeMin="0" RangeMax="0">
+          oAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
+        </DataArray>
+        <DataArray type="UInt64" Name="bulk_element_ids" format="binary" RangeMin="0" RangeMax="780">
+          QAEAAAAAAAAAAAAAAAAAABQAAAAAAAAAKAAAAAAAAAA8AAAAAAAAAFAAAAAAAAAAZAAAAAAAAAB4AAAAAAAAAIwAAAAAAAAAoAAAAAAAAAC0AAAAAAAAAMgAAAAAAAAA3AAAAAAAAADwAAAAAAAAAAQBAAAAAAAAGAEAAAAAAAAsAQAAAAAAAEABAAAAAAAAVAEAAAAAAABoAQAAAAAAAHwBAAAAAAAAkAEAAAAAAACkAQAAAAAAALgBAAAAAAAAzAEAAAAAAADgAQAAAAAAAPQBAAAAAAAACAIAAAAAAAAcAgAAAAAAADACAAAAAAAARAIAAAAAAABYAgAAAAAAAGwCAAAAAAAAgAIAAAAAAACUAgAAAAAAAKgCAAAAAAAAvAIAAAAAAADQAgAAAAAAAOQCAAAAAAAA+AIAAAAAAAAMAwAAAAAAAA==
+        </DataArray>
+        <DataArray type="UInt64" Name="bulk_face_ids" format="binary" RangeMin="4" RangeMax="4">
+          QAEAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAEAAAAAAAAAA==
+        </DataArray>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="1.004987562112089">
+          6AUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAJqZmZmZmak/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmbk/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAmpmZmZmZqT8AAAAAAAAAAJqZmZmZmak/mpmZmZmZqT8AAAAAAAAAAJqZmZmZmbk/mpmZmZmZqT8AAAAAAAAAAAAAAAAAAAAAmpmZmZmZuT8AAAAAAAAAAJqZmZmZmak/mpmZmZmZuT8AAAAAAAAAAJqZmZmZmbk/mpmZmZmZuT8AAAAAAAAAAAAAAAAAAAAAMzMzMzMzwz8AAAAAAAAAAJqZmZmZmak/MzMzMzMzwz8AAAAAAAAAAJqZmZmZmbk/MzMzMzMzwz8AAAAAAAAAAAAAAAAAAAAAmpmZmZmZyT8AAAAAAAAAAJqZmZmZmak/mpmZmZmZyT8AAAAAAAAAAJqZmZmZmbk/mpmZmZmZyT8AAAAAAAAAAAAAAAAAAAAAAAAAAAAA0D8AAAAAAAAAAJqZmZmZmak/AAAAAAAA0D8AAAAAAAAAAJqZmZmZmbk/AAAAAAAA0D8AAAAAAAAAAAAAAAAAAAAAMzMzMzMz0z8AAAAAAAAAAJqZmZmZmak/MzMzMzMz0z8AAAAAAAAAAJqZmZmZmbk/MzMzMzMz0z8AAAAAAAAAAAAAAAAAAAAAZmZmZmZm1j8AAAAAAAAAAJqZmZmZmak/ZmZmZmZm1j8AAAAAAAAAAJqZmZmZmbk/ZmZmZmZm1j8AAAAAAAAAAAAAAAAAAAAAmpmZmZmZ2T8AAAAAAAAAAJqZmZmZmak/mpmZmZmZ2T8AAAAAAAAAAJqZmZmZmbk/mpmZmZmZ2T8AAAAAAAAAAAAAAAAAAAAAzczMzMzM3D8AAAAAAAAAAJqZmZmZmak/zczMzMzM3D8AAAAAAAAAAJqZmZmZmbk/zczMzMzM3D8AAAAAAAAAAAAAAAAAAAAAAAAAAAAA4D8AAAAAAAAAAJqZmZmZmak/AAAAAAAA4D8AAAAAAAAAAJqZmZmZmbk/AAAAAAAA4D8AAAAAAAAAAAAAAAAAAAAAmpmZmZmZ4T8AAAAAAAAAAJqZmZmZmak/mpmZmZmZ4T8AAAAAAAAAAJqZmZmZmbk/mpmZmZmZ4T8AAAAAAAAAAAAAAAAAAAAAMzMzMzMz4z8AAAAAAAAAAJqZmZmZmak/MzMzMzMz4z8AAAAAAAAAAJqZmZmZmbk/MzMzMzMz4z8AAAAAAAAAAAAAAAAAAAAAzczMzMzM5D8AAAAAAAAAAJqZmZmZmak/zczMzMzM5D8AAAAAAAAAAJqZmZmZmbk/zczMzMzM5D8AAAAAAAAAAAAAAAAAAAAAZmZmZmZm5j8AAAAAAAAAAJqZmZmZmak/ZmZmZmZm5j8AAAAAAAAAAJqZmZmZmbk/ZmZmZmZm5j8AAAAAAAAAAAAAAAAAAAAAAAAAAAAA6D8AAAAAAAAAAJqZmZmZmak/AAAAAAAA6D8AAAAAAAAAAJqZmZmZmbk/AAAAAAAA6D8AAAAAAAAAAAAAAAAAAAAAmpmZmZmZ6T8AAAAAAAAAAJqZmZmZmak/mpmZmZmZ6T8AAAAAAAAAAJqZmZmZmbk/mpmZmZmZ6T8AAAAAAAAAAAAAAAAAAAAAMzMzMzMz6z8AAAAAAAAAAJqZmZmZmak/MzMzMzMz6z8AAAAAAAAAAJqZmZmZmbk/MzMzMzMz6z8AAAAAAAAAAAAAAAAAAAAAzczMzMzM7D8AAAAAAAAAAJqZmZmZmak/zczMzMzM7D8AAAAAAAAAAJqZmZmZmbk/zczMzMzM7D8AAAAAAAAAAAAAAAAAAAAAZmZmZmZm7j8AAAAAAAAAAJqZmZmZmak/ZmZmZmZm7j8AAAAAAAAAAJqZmZmZmbk/ZmZmZmZm7j8AAAAAAAAAAAAAAAAAAAAAAAAAAAAA8D8AAAAAAAAAAJqZmZmZmak/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmbk/AAAAAAAA8D8=
+          <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+            <Value index="0">
+              0
+            </Value>
+            <Value index="1">
+              1.0049875621
+            </Value>
+          </InformationKey>
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="62">
+          AAUAAAAAAAABAAAAAAAAAAAAAAAAAAAAAwAAAAAAAAAEAAAAAAAAAAIAAAAAAAAAAQAAAAAAAAAEAAAAAAAAAAUAAAAAAAAABAAAAAAAAAADAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAFAAAAAAAAAAQAAAAAAAAABwAAAAAAAAAIAAAAAAAAAAcAAAAAAAAABgAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACAAAAAAAAAAHAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAKAAAAAAAAAAkAAAAAAAAADAAAAAAAAAANAAAAAAAAAAsAAAAAAAAACgAAAAAAAAANAAAAAAAAAA4AAAAAAAAADQAAAAAAAAAMAAAAAAAAAA8AAAAAAAAAEAAAAAAAAAAOAAAAAAAAAA0AAAAAAAAAEAAAAAAAAAARAAAAAAAAABAAAAAAAAAADwAAAAAAAAASAAAAAAAAABMAAAAAAAAAEQAAAAAAAAAQAAAAAAAAABMAAAAAAAAAFAAAAAAAAAATAAAAAAAAABIAAAAAAAAAFQAAAAAAAAAWAAAAAAAAABQAAAAAAAAAEwAAAAAAAAAWAAAAAAAAABcAAAAAAAAAFgAAAAAAAAAVAAAAAAAAABgAAAAAAAAAGQAAAAAAAAAXAAAAAAAAABYAAAAAAAAAGQAAAAAAAAAaAAAAAAAAABkAAAAAAAAAGAAAAAAAAAAbAAAAAAAAABwAAAAAAAAAGgAAAAAAAAAZAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAcAAAAAAAAABsAAAAAAAAAHgAAAAAAAAAfAAAAAAAAAB0AAAAAAAAAHAAAAAAAAAAfAAAAAAAAACAAAAAAAAAAHwAAAAAAAAAeAAAAAAAAACEAAAAAAAAAIgAAAAAAAAAgAAAAAAAAAB8AAAAAAAAAIgAAAAAAAAAjAAAAAAAAACIAAAAAAAAAIQAAAAAAAAAkAAAAAAAAACUAAAAAAAAAIwAAAAAAAAAiAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAlAAAAAAAAACQAAAAAAAAAJwAAAAAAAAAoAAAAAAAAACYAAAAAAAAAJQAAAAAAAAAoAAAAAAAAACkAAAAAAAAAKAAAAAAAAAAnAAAAAAAAACoAAAAAAAAAKwAAAAAAAAApAAAAAAAAACgAAAAAAAAAKwAAAAAAAAAsAAAAAAAAACsAAAAAAAAAKgAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALAAAAAAAAAArAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAuAAAAAAAAAC0AAAAAAAAAMAAAAAAAAAAxAAAAAAAAAC8AAAAAAAAALgAAAAAAAAAxAAAAAAAAADIAAAAAAAAAMQAAAAAAAAAwAAAAAAAAADMAAAAAAAAANAAAAAAAAAAyAAAAAAAAADEAAAAAAAAANAAAAAAAAAA1AAAAAAAAADQAAAAAAAAAMwAAAAAAAAA2AAAAAAAAADcAAAAAAAAANQAAAAAAAAA0AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA3AAAAAAAAADYAAAAAAAAAOQAAAAAAAAA6AAAAAAAAADgAAAAAAAAANwAAAAAAAAA6AAAAAAAAADsAAAAAAAAAOgAAAAAAAAA5AAAAAAAAADwAAAAAAAAAPQAAAAAAAAA7AAAAAAAAADoAAAAAAAAAPQAAAAAAAAA+AAAAAAAAAA==
+        </DataArray>
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="4" RangeMax="160">
+          QAEAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAA==
+        </DataArray>
+        <DataArray type="UInt8" Name="types" format="binary" RangeMin="9" RangeMax="9">
+          KAAAAAAAAAAJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJ
+        </DataArray>
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>