diff --git a/ProcessLib/LIE/SmallDeformation/Tests.cmake b/ProcessLib/LIE/SmallDeformation/Tests.cmake
index 0dffe32ccd639db97c669f88238eeb3d1440f0f0..50924d6ef4482ed101d106382c942903afde1005 100644
--- a/ProcessLib/LIE/SmallDeformation/Tests.cmake
+++ b/ProcessLib/LIE/SmallDeformation/Tests.cmake
@@ -152,3 +152,7 @@ AddTest(
     expected_sfrac_pcs_0_ts_1_t_1.000000.vtu sfrac_pcs_0_ts_1_t_1.000000.vtu displacement_jump4 displacement_jump4 1e-16 1e-16
     expected_sfrac_pcs_0_ts_1_t_1.000000.vtu sfrac_pcs_0_ts_1_t_1.000000.vtu f_stress_n f_stress_n 1e-6 5e-6
 )
+
+if (NOT OGS_USE_MPI)
+    OgsTest(PROJECTFILE LIE/Mechanics/elastic_push_pull_two_fractures.prj RUNTIME 1)
+endif()
diff --git a/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures.prj b/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures.prj
new file mode 100644
index 0000000000000000000000000000000000000000..a181c4e800bde3647097844412d642beb48cd6a4
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures.prj
@@ -0,0 +1,305 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+    <mesh>five_element_two_fractures.vtu</mesh>
+    <geometry>five_element_two_fractures.gml</geometry>
+    <processes>
+        <process>
+            <name>SD</name>
+            <type>SMALL_DEFORMATION_WITH_LIE</type>
+            <integration_order>2</integration_order>
+            <dimension>2</dimension>
+            <process_variables>
+                <process_variable>displacement</process_variable>
+                <process_variable>displacement_jump1</process_variable>
+                <process_variable>displacement_jump2</process_variable>
+            </process_variables>
+            <constitutive_relation>
+                <type>LinearElasticIsotropic</type>
+                <youngs_modulus>E</youngs_modulus>
+                <poissons_ratio>nu</poissons_ratio>
+            </constitutive_relation>
+            <fracture_model>
+                <type>LinearElasticIsotropic</type>
+                <normal_stiffness>Kn</normal_stiffness>
+                <shear_stiffness>Ks</shear_stiffness>
+                <penalty_aperture_cutoff>1e-5</penalty_aperture_cutoff>
+                <tension_cutoff>0</tension_cutoff>
+            </fracture_model>
+            <fracture_properties>
+                <material_id>0</material_id>
+                <initial_aperture>aperture0</initial_aperture>
+            </fracture_properties>
+            <fracture_properties>
+                <material_id>1</material_id>
+                <initial_aperture>aperture0</initial_aperture>
+            </fracture_properties>
+            <secondary_variables>
+                <secondary_variable internal_name="sigma_xx" output_name="sigma_xx"/>
+                <secondary_variable internal_name="sigma_yy" output_name="sigma_yy"/>
+                <secondary_variable internal_name="sigma_xy" output_name="sigma_xy"/>
+            </secondary_variables>
+        </process>
+    </processes>
+    <time_loop>
+        <processes>
+            <process ref="SD">
+                <nonlinear_solver>basic_newton</nonlinear_solver>
+                <convergence_criterion>
+                    <type>DeltaX</type>
+                    <norm_type>NORM2</norm_type>
+                    <abstol>1e-14</abstol>
+                </convergence_criterion>
+                <time_discretization>
+                    <type>BackwardEuler</type>
+                </time_discretization>
+                <time_stepping>
+                    <type>FixedTimeStepping</type>
+                    <t_initial>0</t_initial>
+                    <t_end>10</t_end>
+                    <timesteps>
+                        <pair>
+                            <repeat>6</repeat>
+                            <delta_t>0.5</delta_t>
+                        </pair>
+                    </timesteps>
+                </time_stepping>
+            </process>
+        </processes>
+        <output>
+            <type>VTK</type>
+            <prefix>elastic_push_pull_two_fractures</prefix>
+            <suffix>_t_{:gtime}</suffix>
+            <timesteps>
+                <pair>
+                    <repeat>100</repeat>
+                    <each_steps>1</each_steps>
+                </pair>
+            </timesteps>
+            <variables>
+                <variable>displacement</variable>
+                <variable>displacement_jump1</variable>
+                <variable>displacement_jump2</variable>
+                <variable>sigma_xx</variable>
+                <variable>sigma_yy</variable>
+                <variable>sigma_xy</variable>
+            </variables>
+        </output>
+    </time_loop>
+    <parameters>
+        <parameter>
+            <name>E</name>
+            <type>Constant</type>
+            <value>200e9</value>
+        </parameter>
+        <parameter>
+            <name>nu</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+        <!-- Fracture properties -->
+        <parameter>
+            <name>Kn</name>
+            <type>Constant</type>
+            <value>20e9</value>
+        </parameter>
+        <parameter>
+            <name>Ks</name>
+            <type>Constant</type>
+            <value>20e9</value>
+        </parameter>
+        <parameter>
+            <name>friction_angle</name>
+            <type>Constant</type>
+            <value>22</value>
+        </parameter>
+        <parameter>
+            <name>dilatancy_angle</name>
+            <type>Constant</type>
+            <value>10</value>
+        </parameter>
+        <parameter>
+            <name>cohesion</name>
+            <type>Constant</type>
+            <value>1e3</value>
+        </parameter>
+        <parameter>
+            <name>displacement0</name>
+            <type>Constant</type>
+            <values>0 0</values>
+        </parameter>
+        <parameter>
+            <name>aperture0</name>
+            <type>Constant</type>
+            <value>1e-5</value>
+        </parameter>
+        <parameter>
+            <name>force</name>
+            <type>Constant</type>
+            <value>1e5</value>
+        </parameter>
+        <parameter>
+            <name>push_pull_right</name>
+            <type>CurveScaled</type>
+            <curve>time_ramp_right</curve>
+            <parameter>force</parameter>
+        </parameter>
+        <parameter>
+            <name>push_pull_top</name>
+            <type>CurveScaled</type>
+            <curve>time_ramp_top</curve>
+            <parameter>force</parameter>
+        </parameter>
+        <parameter>
+            <name>zero_u</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+    </parameters>
+    <curves>
+        <curve>
+            <name>time_ramp_right</name>
+            <coords>0 1  3 4 5  7 8 9 10</coords>
+            <values>0 1 -1 0 0  0 0 1 0</values>
+        </curve>
+        <curve>
+            <name>time_ramp_top</name>
+            <coords>0 1  3 4 5  7 8 9 10</coords>
+            <values>0 0  0 0 1 -1 0 1 0</values>
+        </curve>
+    </curves>
+    <process_variables>
+        <process_variable>
+            <name>displacement</name>
+            <components>2</components>
+            <order>1</order>
+            <initial_condition>displacement0</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <geometrical_set>five_element_two_fractures</geometrical_set>
+                    <geometry>PLY_LEFT</geometry>
+                    <type>Dirichlet</type>
+                    <component>0</component>
+                    <parameter>zero_u</parameter>
+                </boundary_condition>
+                <boundary_condition>
+                    <geometrical_set>five_element_two_fractures</geometrical_set>
+                    <geometry>PLY_BOTTOM</geometry>
+                    <type>Dirichlet</type>
+                    <component>1</component>
+                    <parameter>zero_u</parameter>
+                </boundary_condition>
+                <boundary_condition>
+                    <geometrical_set>five_element_two_fractures</geometrical_set>
+                    <geometry>PLY_RIGHT</geometry>
+                    <type>Neumann</type>
+                    <component>0</component>
+                    <parameter>push_pull_right</parameter>
+                </boundary_condition>
+                <boundary_condition>
+                    <geometrical_set>five_element_two_fractures</geometrical_set>
+                    <geometry>PLY_TOP</geometry>
+                    <type>Neumann</type>
+                    <component>1</component>
+                    <parameter>push_pull_top</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+        </process_variable>
+        <process_variable>
+            <name>displacement_jump1</name>
+            <components>2</components>
+            <order>1</order>
+            <initial_condition>displacement0</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <geometrical_set>five_element_two_fractures</geometrical_set>
+                    <geometry>FRACTURE1_LEFT</geometry>
+                    <type>Dirichlet</type>
+                    <component>0</component>
+                    <parameter>zero_u</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+        </process_variable>
+        <process_variable>
+            <name>displacement_jump2</name>
+            <components>2</components>
+            <order>1</order>
+            <initial_condition>displacement0</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <geometrical_set>five_element_two_fractures</geometrical_set>
+                    <geometry>FRACTURE2_LEFT</geometry>
+                    <type>Dirichlet</type>
+                    <component>0</component>
+                    <parameter>zero_u</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+        </process_variable>
+    </process_variables>
+    <nonlinear_solvers>
+        <nonlinear_solver>
+            <name>basic_newton</name>
+            <type>Newton</type>
+            <max_iter>25</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>1</scaling>
+            </eigen>
+        </linear_solver>
+    </linear_solvers>
+    <test_definition>
+        <vtkdiff>
+            <regex>elastic_push_pull_two_fractures_t_.*.vtu</regex>
+            <field>displacement</field>
+            <absolute_tolerance>1e-16</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+        <vtkdiff>
+            <regex>elastic_push_pull_two_fractures_t_.*.vtu</regex>
+            <field>displacement_jump1</field>
+            <absolute_tolerance>1e-16</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+        <vtkdiff>
+            <regex>elastic_push_pull_two_fractures_t_.*.vtu</regex>
+            <field>displacement_jump2</field>
+            <absolute_tolerance>1e-16</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+        <vtkdiff>
+            <regex>elastic_push_pull_two_fractures_t_.*.vtu</regex>
+            <field>sigma_xx</field>
+            <absolute_tolerance>1e-10</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+        <vtkdiff>
+            <regex>elastic_push_pull_two_fractures_t_.*.vtu</regex>
+            <field>sigma_yy</field>
+            <absolute_tolerance>1e-11</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+        <vtkdiff>
+            <regex>elastic_push_pull_two_fractures_t_.*.vtu</regex>
+            <field>sigma_xy</field>
+            <absolute_tolerance>1e-10</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+        <vtkdiff>
+            <regex>elastic_push_pull_two_fractures_t_.*.vtu</regex>
+            <field>f_stress_n</field>
+            <absolute_tolerance>1e-16</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+        <vtkdiff>
+            <regex>elastic_push_pull_two_fractures_t_.*.vtu</regex>
+            <field>f_stress_s</field>
+            <absolute_tolerance>1e-11</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+    </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures_t_0.vtu b/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures_t_0.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..86f473aa04c18976edf71dc661e6ff5e85e3612e
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures_t_0.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <FieldData>
+      <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45"                   RangeMax="103"                  offset="0"                   />
+    </FieldData>
+    <Piece NumberOfPoints="8"                    NumberOfCells="5"                   >
+      <PointData>
+        <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="0"                    offset="84"                  />
+        <DataArray type="Float64" Name="displacement_jump1" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="0"                    offset="144"                 />
+        <DataArray type="Float64" Name="displacement_jump2" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="0"                    offset="204"                 />
+        <DataArray type="Float64" Name="sigma_xx" format="appended" RangeMin="0"                    RangeMax="0"                    offset="264"                 />
+        <DataArray type="Float64" Name="sigma_xy" format="appended" RangeMin="0"                    RangeMax="0"                    offset="324"                 />
+        <DataArray type="Float64" Name="sigma_yy" format="appended" RangeMin="0"                    RangeMax="0"                    offset="384"                 />
+      </PointData>
+      <CellData>
+        <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0"                    RangeMax="2"                    offset="444"                 />
+        <DataArray type="Float64" Name="aperture" format="appended" RangeMin="0"                    RangeMax="1e-05"                offset="512"                 />
+        <DataArray type="Float64" Name="f_stress_n" format="appended" RangeMin="0"                    RangeMax="0"                    offset="588"                 />
+        <DataArray type="Float64" Name="f_stress_s" format="appended" RangeMin="0"                    RangeMax="0"                    offset="648"                 />
+        <DataArray type="Float64" Name="levelset1" format="appended" RangeMin="-0.5"                 RangeMax="0.5"                  offset="708"                 />
+        <DataArray type="Float64" Name="levelset2" format="appended" RangeMin="-0.5"                 RangeMax="0.5"                  offset="776"                 />
+        <DataArray type="Float64" Name="strain_xx" format="appended" RangeMin="0"                    RangeMax="0"                    offset="844"                 />
+        <DataArray type="Float64" Name="strain_xy" format="appended" RangeMin="0"                    RangeMax="0"                    offset="904"                 />
+        <DataArray type="Float64" Name="strain_yy" format="appended" RangeMin="0"                    RangeMax="0"                    offset="964"                 />
+        <DataArray type="Float64" Name="strain_zz" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1024"                />
+        <DataArray type="Float64" Name="stress_xx" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1084"                />
+        <DataArray type="Float64" Name="stress_xy" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1144"                />
+        <DataArray type="Float64" Name="stress_yy" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1204"                />
+        <DataArray type="Float64" Name="stress_zz" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1264"                />
+        <DataArray type="Float64" Name="w_n" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1324"                />
+        <DataArray type="Float64" Name="w_s" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1384"                />
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="3.1622776602"         offset="1444"                />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="1536"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="1640"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="1716"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPWM9I1NjHVTTdLMrBMTrE0NgAALlgEuQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAADAAAAAAAAAA=eF5jYKAMAAAAQAABAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAADAAAAAAAAAA=eF5jYKAMAAAAQAABAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAADAAAAAAAAAA=eF5jYKAMAAAAQAABAQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAEQAAAAAAAAA=eF5jYGBgYARiJiQMAABUAAg=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFgAAAAAAAAA=eF77mPG4Y+uPJ3YfoTQDDgAAbJcLJw==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAEgAAAAAAAAA=eF5jYMAGHthD6f0wEQAjjgK/AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAEAAAAAAAAAA=eF5jYMAHHthD6f0ADZ4Cvw==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAIwAAAAAAAAA=eF5jYMAHPtgTEHfALo9LHN08GB+XPRwOqHyYOoQ4ALH4CCs=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAKwAAAAAAAAA=eF5jZYAAFijNBKWZGVABI5RmhdIw9TB1MD5MHmYOG5o6mDg7lAYADkAAOQ==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFgAAAAAAAAA=eF5jYoAAFijNAaV5oLQAlAYAAvgAKw==AQAAAAAAAAAAgAAAAAAAAAUAAAAAAAAADQAAAAAAAAA=eF5jZubk5AQAAFYAIg==
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures_t_1.vtu b/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures_t_1.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..0c41b33f86024d766a9fa0ef9cd8fb0878caa7d6
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures_t_1.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <FieldData>
+      <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45"                   RangeMax="103"                  offset="0"                   />
+    </FieldData>
+    <Piece NumberOfPoints="8"                    NumberOfCells="5"                   >
+      <PointData>
+        <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="5e-07"                offset="84"                  />
+        <DataArray type="Float64" Name="displacement_jump1" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="3.0440279142e-23"     offset="220"                 />
+        <DataArray type="Float64" Name="displacement_jump2" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="2.2054896457e-22"     offset="312"                 />
+        <DataArray type="Float64" Name="sigma_xx" format="appended" RangeMin="100000"               RangeMax="100000"               offset="408"                 />
+        <DataArray type="Float64" Name="sigma_xy" format="appended" RangeMin="-2.2249097855e-11"    RangeMax="1.28272497e-11"       offset="500"                 />
+        <DataArray type="Float64" Name="sigma_yy" format="appended" RangeMin="-1.6604511801e-12"    RangeMax="4.9901190925e-12"     offset="644"                 />
+      </PointData>
+      <CellData>
+        <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0"                    RangeMax="2"                    offset="784"                 />
+        <DataArray type="Float64" Name="aperture" format="appended" RangeMin="0"                    RangeMax="1e-05"                offset="852"                 />
+        <DataArray type="Float64" Name="f_stress_n" format="appended" RangeMin="0"                    RangeMax="1.0602290566e-27"     offset="928"                 />
+        <DataArray type="Float64" Name="f_stress_s" format="appended" RangeMin="-1.3864254603e-13"    RangeMax="2.1907966211e-12"     offset="1000"                />
+        <DataArray type="Float64" Name="levelset1" format="appended" RangeMin="-0.5"                 RangeMax="0.5"                  offset="1080"                />
+        <DataArray type="Float64" Name="levelset2" format="appended" RangeMin="-0.5"                 RangeMax="0.5"                  offset="1148"                />
+        <DataArray type="Float64" Name="strain_xx" format="appended" RangeMin="0"                    RangeMax="5e-07"                offset="1216"                />
+        <DataArray type="Float64" Name="strain_xy" format="appended" RangeMin="-4.8577953294e-23"    RangeMax="0"                    offset="1292"                />
+        <DataArray type="Float64" Name="strain_yy" format="appended" RangeMin="0"                    RangeMax="1.3737014523e-23"     offset="1384"                />
+        <DataArray type="Float64" Name="strain_zz" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1472"                />
+        <DataArray type="Float64" Name="stress_xx" format="appended" RangeMin="0"                    RangeMax="100000"               offset="1532"                />
+        <DataArray type="Float64" Name="stress_xy" format="appended" RangeMin="-9.7155906588e-12"    RangeMax="0"                    offset="1604"                />
+        <DataArray type="Float64" Name="stress_yy" format="appended" RangeMin="0"                    RangeMax="2.7474029045e-12"     offset="1700"                />
+        <DataArray type="Float64" Name="stress_zz" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1788"                />
+        <DataArray type="Float64" Name="w_n" format="appended" RangeMin="0"                    RangeMax="5.2897245787e-38"     offset="1848"                />
+        <DataArray type="Float64" Name="w_s" format="appended" RangeMin="-6.9321273015e-24"    RangeMax="1.0953983106e-22"     offset="1920"                />
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="3.1622776602"         offset="2000"                />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="2092"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="2196"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="2272"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPWM9I1NjHVTTdLMrBMTrE0NgAALlgEuQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAARAAAAAAAAAA=eF5jYEAFvW+3Lvh+bIEdjN8H5WdobVjQfsfXGiZuA+HvhvGtEmo0ft6x3A3THwHhw9XLHdR6xMYeDJc3gfCtAXcRKuY=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAIgAAAAAAAAA=eF5jYKAOENhQb1KZ6GAd4LA16/ghA2sFCH83IX0ArjQJxA==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAJwAAAAAAAAA=eF5jYMAP2rl1P87tLdjNwrWGdftOu90w8QecYL41Pr3EAABrhgqxAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAIwAAAAAAAAA=eF779x8IMn84/ILSDCCQ9cPhD5o4E1T8L5T/DUoDAH5dKtY=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAASwAAAAAAAAA=eF4BQAC//3yfQeHD1Uk9AliFNrL7pL1cRFm0goq1vfJcxYAZNaw9V1MGl5a3oj2+zrqA/Z6rvZ6ddv39XKY9DaQAdpB2uL1FRCI8AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAARwAAAAAAAAA=eF773hQUduFgou3Me/rP4z9Ntf0xZzOv2f4eWx/Z7R23egr27jl36snlL862uiea9nFtmmp7aU35fUaz2r2vIbQtAIddIrM=AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAEQAAAAAAAAA=eF5jYGBgYARiJiQMAABUAAg=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFgAAAAAAAAA=eF77mPG4Y+uPJ3YfoTQDDgAAbJcLJw==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFQAAAAAAAAA=eF5jYACChiArBjAIhdKYAAAyUgGcAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAGwAAAAAAAAA=eF5zn2pzQr7Zea9X8WYXG9dmWwYcAADUQAZ4AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAEgAAAAAAAAA=eF5jYMAGHthD6f0wEQAjjgK/AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAEAAAAAAAAAA=eF5jYMAHHthD6f0ADZ4Cvw==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFgAAAAAAAAA=eF5jYEAFfW+3Lvh+bIEdOg0A0A8QQg==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAIwAAAAAAAAA=eF5jYEAD9zrlObLsdkuUFQeZKbHvblEw/SsU67sbAGdICKw=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAIAAAAAAAAAA=eF5jYEAFqxvcs9NmG1izvM7ie6ikYw0R5bAAAGnxBls=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAEwAAAAAAAAA=eF5jYMACsn44wGhGKA0AMaoE6A==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAJQAAAAAAAAA=eF5jYEAFX/YU5IgqTtvbcyvm988rCXs9tl3QSYlduhcArbYN2g==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAIAAAAAAAAAA=eF5jYEAFf/R+fZyp1WHrzxdx8kRFiy1ENNEKAI1yCLg=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFAAAAAAAAAA=eF5jYAABAwswxWAEpTEBABhMANM=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAGwAAAAAAAAA=eF4LcNiadfyQwu52bt2Pc3sTrBlwAAD1yweJAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAIwAAAAAAAAA=eF5jYMAHPtgTEHfALo9LHN08GB+XPRwOqHyYOoQ4ALH4CCs=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAKwAAAAAAAAA=eF5jZYAAFijNBKWZGVABI5RmhdIw9TB1MD5MHmYOG5o6mDg7lAYADkAAOQ==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFgAAAAAAAAA=eF5jYoAAFijNAaV5oLQAlAYAAvgAKw==AQAAAAAAAAAAgAAAAAAAAAUAAAAAAAAADQAAAAAAAAA=eF5jZubk5AQAAFYAIg==
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures_t_2.vtu b/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures_t_2.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..5382171334293d76ac48e4481f77a224f66a37af
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures_t_2.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <FieldData>
+      <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45"                   RangeMax="103"                  offset="0"                   />
+    </FieldData>
+    <Piece NumberOfPoints="8"                    NumberOfCells="5"                   >
+      <PointData>
+        <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="3.7384045184e-23"     offset="84"                  />
+        <DataArray type="Float64" Name="displacement_jump1" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="3.893521669e-23"      offset="248"                 />
+        <DataArray type="Float64" Name="displacement_jump2" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="5.5026260157e-23"     offset="340"                 />
+        <DataArray type="Float64" Name="sigma_xx" format="appended" RangeMin="-3.7093386428e-11"    RangeMax="2.5692053265e-11"     offset="436"                 />
+        <DataArray type="Float64" Name="sigma_xy" format="appended" RangeMin="-6.5123768933e-12"    RangeMax="4.8096420343e-12"     offset="572"                 />
+        <DataArray type="Float64" Name="sigma_yy" format="appended" RangeMin="-3.9131442357e-12"    RangeMax="3.7109697464e-12"     offset="716"                 />
+      </PointData>
+      <CellData>
+        <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0"                    RangeMax="2"                    offset="860"                 />
+        <DataArray type="Float64" Name="aperture" format="appended" RangeMin="0"                    RangeMax="1e-05"                offset="928"                 />
+        <DataArray type="Float64" Name="f_stress_n" format="appended" RangeMin="-2.3728936029e-27"    RangeMax="3.3053271929e-28"     offset="1004"                />
+        <DataArray type="Float64" Name="f_stress_s" format="appended" RangeMin="-4.6918958783e-13"    RangeMax="3.8918866227e-13"     offset="1080"                />
+        <DataArray type="Float64" Name="levelset1" format="appended" RangeMin="-0.5"                 RangeMax="0.5"                  offset="1160"                />
+        <DataArray type="Float64" Name="levelset2" format="appended" RangeMin="-0.5"                 RangeMax="0.5"                  offset="1228"                />
+        <DataArray type="Float64" Name="strain_xx" format="appended" RangeMin="-4.9269023769e-23"    RangeMax="6.3125781892e-24"     offset="1296"                />
+        <DataArray type="Float64" Name="strain_xy" format="appended" RangeMin="-5.5413769302e-23"    RangeMax="7.3587798427e-24"     offset="1388"                />
+        <DataArray type="Float64" Name="strain_yy" format="appended" RangeMin="-9.8580549082e-24"    RangeMax="0"                    offset="1484"                />
+        <DataArray type="Float64" Name="strain_zz" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1580"                />
+        <DataArray type="Float64" Name="stress_xx" format="appended" RangeMin="-7.2759576142e-12"    RangeMax="0"                    offset="1640"                />
+        <DataArray type="Float64" Name="stress_xy" format="appended" RangeMin="-1.108275386e-11"     RangeMax="1.4717559685e-12"     offset="1712"                />
+        <DataArray type="Float64" Name="stress_yy" format="appended" RangeMin="-1.9716109816e-12"    RangeMax="0"                    offset="1808"                />
+        <DataArray type="Float64" Name="stress_zz" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1904"                />
+        <DataArray type="Float64" Name="w_n" format="appended" RangeMin="-1.1754943508e-37"    RangeMax="1.6530389308e-38"     offset="1964"                />
+        <DataArray type="Float64" Name="w_s" format="appended" RangeMin="-2.3459479391e-23"    RangeMax="1.9459433114e-23"     offset="2040"                />
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="3.1622776602"         offset="2120"                />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="2212"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="2316"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="2392"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPWM9I1NjHVTTdLMrBMTrE0NgAALlgEuQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAWAAAAAAAAAA=eF5jYEAFb6T3T+zfYW8N478ImiV4Sstpd4LPN96rd6R3w8QfeIP5cHVF++0WLl2vultZXolFc6KudQWED5cXMOSaU95gaR3xIif3+3Inay8IfzcAx1EjhA==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAJAAAAAAAAAA=eF5jYKAOOCMixnrl0lcrpaMCiult7rsbRMH8XYT0AQDaQwrdAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAJgAAAAAAAAA=eF5jYMAPtP2kiqJjfKwDrjAbr0twtIaJa0D4u/HpJQYAAKKtCGQ=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAARQAAAAAAAAA=eF4LUtMpsb29dG8QlE7avFy81XPOXmFn3+7V9ntsa1IPOV3fnGu72ifiRdW27r17ZBTF81KO7D0fn7dIcsm+vQC6ux7RAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAASwAAAAAAAAA=eF4BQAC//wAAAAAAICo6d99So/6Pgb3q6kN57EqavSj2JDdIcHI9CrgFbZfOg70ShEgDj0yUvc2lxO4sJ5U94cpz10iknL1udh4BAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAASwAAAAAAAAA=eF4BQAC//yoV9QjPNZG9ADE1SajuIL1UwnYZ0ASIvW52pGmCdHM9k/bXCLzNh71Qn1MGLr9svaIctTQuUpA96kR06OH2kL22GB87AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAEQAAAAAAAAA=eF5jYGBgYARiJiQMAABUAAg=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFgAAAAAAAAA=eF77mPG4Y+uPJ3YfoTQDDgAAbJcLJw==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFwAAAAAAAAA=eF5jYAACAysrEMXQkL6LAQcAAEAGAkY=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAGwAAAAAAAAA=eF7jX9H38n9StO0G1a/tok0JexlwAAAK0AgtAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAEgAAAAAAAAA=eF5jYMAGHthD6f0wEQAjjgK/AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAEAAAAAAAAAA=eF5jYMAHHthD6f0ADZ4Cvw==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAJAAAAAAAAAA=eF5jYEAFPhvc6u2VNXZ7Vn15OL1Nztr8cY4M8wnf3QCBkQqKAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAJQAAAAAAAAA=eF5jYEAFFxa9rlipfG5XhMHq43n7A3Z/MX5heveUojUAs4UNow==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAJQAAAAAAAAA=eF5jYEAFURxi6tevqu/+9/Fj3OuDSrsb7gpYr/a+swsAmuUNKQ==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFAAAAAAAAAA=eF5jYMAGJuyF0A1QesFeACR3A+g=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAJQAAAAAAAAA=eF5jYEAFO17acD16q7p3W5yhBmf8ir07QniljZ9U2gIAk40K1A==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAJQAAAAAAAAA=eF5jYEAFE+rbI9eEN+7duvdU4jbv6r0JHE/Ly6aZ7AUAorQMyA==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFgAAAAAAAAA=eF5jYACCBjELBjBw2cGAAwAANKwByw==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAGwAAAAAAAAA=eF5TOiqgmN5mbq3tJ1UUHWOzmwEHAAClIgUrAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAIwAAAAAAAAA=eF5jYMAHPtgTEHfALo9LHN08GB+XPRwOqHyYOoQ4ALH4CCs=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAKwAAAAAAAAA=eF5jZYAAFijNBKWZGVABI5RmhdIw9TB1MD5MHmYOG5o6mDg7lAYADkAAOQ==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFgAAAAAAAAA=eF5jYoAAFijNAaV5oLQAlAYAAvgAKw==AQAAAAAAAAAAgAAAAAAAAAUAAAAAAAAADQAAAAAAAAA=eF5jZubk5AQAAFYAIg==
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures_t_3.vtu b/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures_t_3.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..06acaa60b80f6a0e7ee983b13179a83ebe604514
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures_t_3.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <FieldData>
+      <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45"                   RangeMax="103"                  offset="0"                   />
+    </FieldData>
+    <Piece NumberOfPoints="8"                    NumberOfCells="5"                   >
+      <PointData>
+        <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="5e-07"                offset="84"                  />
+        <DataArray type="Float64" Name="displacement_jump1" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="9.7209865587e-23"     offset="224"                 />
+        <DataArray type="Float64" Name="displacement_jump2" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="3.3374502509e-22"     offset="316"                 />
+        <DataArray type="Float64" Name="sigma_xx" format="appended" RangeMin="-100000"              RangeMax="-100000"              offset="412"                 />
+        <DataArray type="Float64" Name="sigma_xy" format="appended" RangeMin="-1.198190553e-11"     RangeMax="2.0436821284e-11"     offset="508"                 />
+        <DataArray type="Float64" Name="sigma_yy" format="appended" RangeMin="-4.5616845589e-12"    RangeMax="2.3561010176e-12"     offset="644"                 />
+      </PointData>
+      <CellData>
+        <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0"                    RangeMax="2"                    offset="784"                 />
+        <DataArray type="Float64" Name="aperture" format="appended" RangeMin="0"                    RangeMax="1e-05"                offset="852"                 />
+        <DataArray type="Float64" Name="f_stress_n" format="appended" RangeMin="-2.1204581132e-27"    RangeMax="0"                    offset="928"                 />
+        <DataArray type="Float64" Name="f_stress_s" format="appended" RangeMin="-3.1373041282e-12"    RangeMax="0"                    offset="1000"                />
+        <DataArray type="Float64" Name="levelset1" format="appended" RangeMin="-0.5"                 RangeMax="0.5"                  offset="1080"                />
+        <DataArray type="Float64" Name="levelset2" format="appended" RangeMin="-0.5"                 RangeMax="0.5"                  offset="1148"                />
+        <DataArray type="Float64" Name="strain_xx" format="appended" RangeMin="-5e-07"               RangeMax="0"                    offset="1216"                />
+        <DataArray type="Float64" Name="strain_xy" format="appended" RangeMin="-6.4144543038e-24"    RangeMax="3.8222716351e-23"     offset="1292"                />
+        <DataArray type="Float64" Name="strain_yy" format="appended" RangeMin="-1.3737014523e-23"    RangeMax="0"                    offset="1388"                />
+        <DataArray type="Float64" Name="strain_zz" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1476"                />
+        <DataArray type="Float64" Name="stress_xx" format="appended" RangeMin="-100000"              RangeMax="0"                    offset="1536"                />
+        <DataArray type="Float64" Name="stress_xy" format="appended" RangeMin="-1.2828908608e-12"    RangeMax="7.6445432701e-12"     offset="1612"                />
+        <DataArray type="Float64" Name="stress_yy" format="appended" RangeMin="-2.7474029045e-12"    RangeMax="0"                    offset="1708"                />
+        <DataArray type="Float64" Name="stress_zz" format="appended" RangeMin="0"                    RangeMax="0"                    offset="1796"                />
+        <DataArray type="Float64" Name="w_n" format="appended" RangeMin="-1.0579449157e-37"    RangeMax="0"                    offset="1856"                />
+        <DataArray type="Float64" Name="w_s" format="appended" RangeMin="-1.5686520641e-22"    RangeMax="0"                    offset="1928"                />
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="3.1622776602"         offset="2008"                />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="2100"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="2204"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="2280"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPWM9I1NjHVTTdLMrBMTrE0NgAALlgEuQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAARwAAAAAAAAA=eF5jYEAFvW+3Lvh+bME+GL8Pyi/ofBVz09N7N0xcA8K3hvE/fNBwND2vbN0DVR/wEcyHq+d4a3ct+V88XN4Ewt8NAA6lL5k=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAIgAAAAAAAAA=eF5jYKAOyFjebqvxzsu64OQD50syUdY+EP5uQvoA7soLIA==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAJgAAAAAAAAA=eF5jYMAPXOrlD+VuKbf+zOu79b5xojVM/BWEvxufXmIAAJZbC08=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAJQAAAAAAAAA=eF77/x8IMn8c+AelGRmAIOvHgV9QPoxmgor/h/K/QWkA830u3A==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAARAAAAAAAAAA=eF5jYACBcCswxRC2S2O6KsNLi7W2H7pmi6+WXbhXJa9UQMGoa+8Dv/nVnyY32HrnrZ5oFrlq7zrrZ9olFdtsAWBZF+E=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAARgAAAAAAAAA=eF47EeTWl55VsNd5VgYfg8CUvUkXp2asa+rZG9BWbXdauMBWRJhj5oGtmXuNZTq8fkpN3luZdudw+c4W204IvRcAqs4euQ==AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAEQAAAAAAAAA=eF5jYGBgYARiJiQMAABUAAg=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFgAAAAAAAAA=eF77mPG4Y+uPJ3YfoTQDDgAAbJcLJw==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFAAAAAAAAAA=eF5jYACBhF1giiEVSmMCAEFOAjo=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAGwAAAAAAAAA=eF7jOhTLpJGWt9fVbO2mKTO69zLgAADgTgczAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAEgAAAAAAAAA=eF5jYMAGHthD6f0wEQAjjgK/AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAEAAAAAAAAAA=eF5jYMAHHthD6f0ADZ4Cvw==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFgAAAAAAAAA=eF5jYEAFfW+3Lvh+bME+dBoA3Y8Rwg==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAJQAAAAAAAAA=eF5jYEAFDxi9t21lkd/t8JC7QKjgx67JhjNOrZVytwYAhj8KtQ==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAIAAAAAAAAAA=eF5jYEAFqxrcs9NmG+x+8SqL76GSzm6IqMAOAIWaCMU=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFgAAAAAAAAA=eF5jYEAFTCAi68cBBijNCKUBP1oGag==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAJQAAAAAAAAA=eF5jYEAFe9wVyhZMLNvboV5q8eR44F6FxvZTdecX2AIAk0sMLw==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAIAAAAAAAAAA=eF5jYEAFv/R+fZyp1bHXgi/i5ImKlr0Q0bRdAJlcCiQ=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAADAAAAAAAAAA=eF5jYCAOAAAAKAABAQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFAAAAAAAAAA=eF5jYAABmx1gisEJSmMCADiEAe8=AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAGwAAAAAAAAA=eF4rOPnA+ZKM126XevlDuVvSdzPgAAALrgg3AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAIwAAAAAAAAA=eF5jYMAHPtgTEHfALo9LHN08GB+XPRwOqHyYOoQ4ALH4CCs=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAKwAAAAAAAAA=eF5jZYAAFijNBKWZGVABI5RmhdIw9TB1MD5MHmYOG5o6mDg7lAYADkAAOQ==AQAAAAAAAAAAgAAAAAAAACgAAAAAAAAAFgAAAAAAAAA=eF5jYoAAFijNAaV5oLQAlAYAAvgAKw==AQAAAAAAAAAAgAAAAAAAAAUAAAAAAAAADQAAAAAAAAA=eF5jZubk5AQAAFYAIg==
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/LIE/Mechanics/five_element_two_fractures.gml b/Tests/Data/LIE/Mechanics/five_element_two_fractures.gml
new file mode 100644
index 0000000000000000000000000000000000000000..7d96b7cfbd3a4690f469b40cc8d546a296ef41b2
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/five_element_two_fractures.gml
@@ -0,0 +1,32 @@
+<?xml version="1.0" encoding="utf-8"?>
+<OpenGeoSysGLI>
+    <name>five_element_two_fractures</name>
+    <points>
+        <point id="0" x="0" y="0" z="0"/>
+        <point id="1" x="0" y="1" z="0" name="FRACTURE1_LEFT"/>
+        <point id="2" x="0" y="2" z="0" name="FRACTURE2_LEFT"/>
+        <point id="3" x="0" y="3" z="0"/>
+        <point id="4" x="1" y="0" z="0"/>
+        <point id="5" x="1" y="1" z="0" name="FRACTURE1_RIGHT"/>
+        <point id="6" x="1" y="2" z="0" name="FRACTURE2_RIGHT"/>
+        <point id="7" x="1" y="3" z="0"/>
+    </points>
+    <polylines>
+        <polyline id="0" name="PLY_BOTTOM">
+            <pnt>0</pnt>
+            <pnt>4</pnt>
+        </polyline>
+        <polyline id="1" name="PLY_TOP">
+            <pnt>3</pnt>
+            <pnt>7</pnt>
+        </polyline>
+        <polyline id="2" name="PLY_LEFT">
+            <pnt>0</pnt>
+            <pnt>3</pnt>
+        </polyline>
+        <polyline id="3" name="PLY_RIGHT">
+            <pnt>4</pnt>
+            <pnt>7</pnt>
+        </polyline>
+    </polylines>
+</OpenGeoSysGLI>
diff --git a/Tests/Data/LIE/Mechanics/five_element_two_fractures.vtu b/Tests/Data/LIE/Mechanics/five_element_two_fractures.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..1d234b1a3e89f18a8e1c641a72681d9390f2e8ff
Binary files /dev/null and b/Tests/Data/LIE/Mechanics/five_element_two_fractures.vtu differ