diff --git a/Applications/Utils/ModelPreparation/PartitionMesh/NodeWiseMeshPartitioner.cpp b/Applications/Utils/ModelPreparation/PartitionMesh/NodeWiseMeshPartitioner.cpp
index 7c86dd8bb00169e1c053a3981b5a17e759ae3009..6ddd8939294d305ad48c209660a3f4e899dcace8 100644
--- a/Applications/Utils/ModelPreparation/PartitionMesh/NodeWiseMeshPartitioner.cpp
+++ b/Applications/Utils/ModelPreparation/PartitionMesh/NodeWiseMeshPartitioner.cpp
@@ -337,12 +337,14 @@ std::size_t copyNodePropertyVectorValues(
{
auto const& nodes = p.nodes;
auto const nnodes = nodes.size();
+ auto const n_components = pv.getNumberOfComponents();
for (std::size_t i = 0; i < nnodes; ++i)
{
const auto global_id = nodes[i]->getID();
- partitioned_pv[offset + i] = pv[global_id];
+ std::copy_n(&pv[n_components * global_id], n_components,
+ &partitioned_pv[offset + n_components * i]);
}
- return nnodes;
+ return n_components * nnodes;
}
/// Copies the properties from global property vector \c pv to the
@@ -356,19 +358,22 @@ std::size_t copyCellPropertyVectorValues(
MeshLib::PropertyVector<T>& partitioned_pv)
{
std::size_t const n_regular(p.regular_elements.size());
+ auto const n_components = pv.getNumberOfComponents();
for (std::size_t i = 0; i < n_regular; ++i)
{
const auto id = p.regular_elements[i]->getID();
- partitioned_pv[offset + i] = pv[id];
+ std::copy_n(&pv[n_components * id], n_components,
+ &partitioned_pv[offset + n_components * i]);
}
std::size_t const n_ghost(p.ghost_elements.size());
for (std::size_t i = 0; i < n_ghost; ++i)
{
const auto id = p.ghost_elements[i]->getID();
- partitioned_pv[offset + n_regular + i] = pv[id];
+ std::copy_n(&pv[n_components * id], n_components,
+ &partitioned_pv[offset + n_components * (n_regular + i)]);
}
- return n_regular + n_ghost;
+ return n_components * (n_regular + n_ghost);
}
template <typename T>
diff --git a/MeshLib/IO/MPI_IO/NodePartitionedMeshReader.cpp b/MeshLib/IO/MPI_IO/NodePartitionedMeshReader.cpp
index ac0d28fc87465376b704d958023d022a631f38f8..b8db24bd7fce0028f996d9d76e1a858ce71997e9 100644
--- a/MeshLib/IO/MPI_IO/NodePartitionedMeshReader.cpp
+++ b/MeshLib/IO/MPI_IO/NodePartitionedMeshReader.cpp
@@ -303,7 +303,7 @@ void NodePartitionedMeshReader::readDomainSpecificPartOfPropertyVectors(
"{:d}.",
_mpi_rank, global_offset,
global_offset +
- pvpmd.offset * vec_pvmd[i]->data_type_size_in_bytes);
+ pvpmd.offset * vec_pvmd[i]->number_of_components * vec_pvmd[i]->data_type_size_in_bytes);
if (vec_pvmd[i]->is_int_type)
{
if (vec_pvmd[i]->is_data_type_signed)
diff --git a/MeshLib/IO/MPI_IO/NodePartitionedMeshReader.h b/MeshLib/IO/MPI_IO/NodePartitionedMeshReader.h
index b954ef8033df6c6eb932f65b5ce3b17e3a90c3fd..e7fae45575d882cc350cad80d6560f8f4fa654ad 100644
--- a/MeshLib/IO/MPI_IO/NodePartitionedMeshReader.h
+++ b/MeshLib/IO/MPI_IO/NodePartitionedMeshReader.h
@@ -192,7 +192,8 @@ private:
pv->resize(pvpmd.number_of_tuples * pvmd.number_of_components);
// jump to the place for reading the specific part of the
// PropertyVector
- is.seekg(global_offset + pvpmd.offset * sizeof(T));
+ is.seekg(global_offset +
+ pvpmd.offset * pvmd.number_of_components * sizeof(T));
// read the values
unsigned long const number_of_bytes = pvmd.data_type_size_in_bytes *
pvpmd.number_of_tuples *
diff --git a/ProcessLib/SmallDeformation/Tests.cmake b/ProcessLib/SmallDeformation/Tests.cmake
index 044b006646583613498af2a0b7905fe927e0db47..57fcdd54a009f1ad990bb2917b52e0e36bec8efe 100644
--- a/ProcessLib/SmallDeformation/Tests.cmake
+++ b/ProcessLib/SmallDeformation/Tests.cmake
@@ -39,6 +39,7 @@ if (NOT OGS_USE_MPI)
OgsTest(PROJECTFILE Mechanics/InitialStates/into_initial_state.prj)
OgsTest(PROJECTFILE Mechanics/InitialStates/equilibrium_restart.prj)
OgsTest(PROJECTFILE Mechanics/InitialStates/non_equilibrium_initial_state.prj)
+ OgsTest(PROJECTFILE Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise.prj)
OgsTest(PROJECTFILE Mechanics/MohrCoulombAbboSloan/slope.prj RUNTIME 50)
set_tests_properties(ogs-Mechanics/MohrCoulombAbboSloan/slope PROPERTIES WILL_FAIL true)
endif()
@@ -50,6 +51,7 @@ endif()
if (OGS_USE_MPI)
# OgsTest(WRAPPER mpirun -np 4 PROJECTFILE Mechanics/Linear/disc_with_hole.prj)
+ OgsTest(WRAPPER mpirun -np 2 PROJECTFILE Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise.prj)
endif()
if (OGS_USE_MFRONT)
diff --git a/Tests/Data/Mechanics/InitialStates/equilibrium_restart.prj b/Tests/Data/Mechanics/InitialStates/equilibrium_restart.prj
index 5edd7f8e76bb99dd4e8b76dafea0a551b7e371e4..f14ed6a93175129b530216e360fb0d8b9f5ec719 100644
--- a/Tests/Data/Mechanics/InitialStates/equilibrium_restart.prj
+++ b/Tests/Data/Mechanics/InitialStates/equilibrium_restart.prj
@@ -2,10 +2,10 @@
<OpenGeoSysProject>
<meshes>
<mesh>soil_column_q8_with_ic.vtu</mesh>
- <mesh>soil_column_bottom.vtu</mesh>
- <mesh>soil_column_top.vtu</mesh>
- <mesh>soil_column_left.vtu</mesh>
- <mesh>soil_column_right.vtu</mesh>
+ <mesh>soil_column_q8_bottom.vtu</mesh>
+ <mesh>soil_column_q8_top.vtu</mesh>
+ <mesh>soil_column_q8_left.vtu</mesh>
+ <mesh>soil_column_q8_right.vtu</mesh>
</meshes>
<processes>
<process>
@@ -113,7 +113,7 @@
<parameter>maximum_top_pressure</parameter>
</parameter>
<parameter>
- <mesh>soil_column_right</mesh>
+ <mesh>soil_column_q8_right</mesh>
<name>side_load</name>
<type>Function</type>
<expression>-0.8*2200.0*9.81*(100.0-y)</expression>
@@ -135,25 +135,25 @@
<boundary_conditions>
<!-- fixed boundaries -->
<boundary_condition>
- <mesh>soil_column_bottom</mesh>
+ <mesh>soil_column_q8_bottom</mesh>
<type>Dirichlet</type>
<component>1</component>
<parameter>Zero</parameter>
</boundary_condition>
<boundary_condition>
- <mesh>soil_column_left</mesh>
+ <mesh>soil_column_q8_left</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>Zero</parameter>
</boundary_condition>
<boundary_condition>
- <mesh>soil_column_top</mesh>
+ <mesh>soil_column_q8_top</mesh>
<type>Neumann</type>
<component>1</component>
<parameter>top_load</parameter>
</boundary_condition>
<boundary_condition>
- <mesh>soil_column_right</mesh>
+ <mesh>soil_column_q8_right</mesh>
<type>Neumann</type>
<component>0</component>
<parameter>side_load</parameter>
diff --git a/Tests/Data/Mechanics/InitialStates/into_initial_state.prj b/Tests/Data/Mechanics/InitialStates/into_initial_state.prj
index 02d9ba23669f3a93c7a728454c4b9741c88d119b..f94af1ff6b86f9977bb36caba876d4431dd5312b 100644
--- a/Tests/Data/Mechanics/InitialStates/into_initial_state.prj
+++ b/Tests/Data/Mechanics/InitialStates/into_initial_state.prj
@@ -2,10 +2,10 @@
<OpenGeoSysProject>
<meshes>
<mesh>soil_column_q8.vtu</mesh>
- <mesh>soil_column_bottom.vtu</mesh>
- <mesh>soil_column_top.vtu</mesh>
- <mesh>soil_column_left.vtu</mesh>
- <mesh>soil_column_right.vtu</mesh>
+ <mesh>soil_column_q8_bottom.vtu</mesh>
+ <mesh>soil_column_q8_top.vtu</mesh>
+ <mesh>soil_column_q8_left.vtu</mesh>
+ <mesh>soil_column_q8_right.vtu</mesh>
</meshes>
<processes>
<process>
@@ -113,7 +113,7 @@
<parameter>maximum_top_pressure</parameter>
</parameter>
<parameter>
- <mesh>soil_column_right</mesh>
+ <mesh>soil_column_q8_right</mesh>
<name>side_load</name>
<type>Function</type>
<expression>-0.8*2200.0*9.81*(100.0-y)</expression>
@@ -135,25 +135,25 @@
<boundary_conditions>
<!-- fixed boundaries -->
<boundary_condition>
- <mesh>soil_column_bottom</mesh>
+ <mesh>soil_column_q8_bottom</mesh>
<type>Dirichlet</type>
<component>1</component>
<parameter>Zero</parameter>
</boundary_condition>
<boundary_condition>
- <mesh>soil_column_left</mesh>
+ <mesh>soil_column_q8_left</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>Zero</parameter>
</boundary_condition>
<boundary_condition>
- <mesh>soil_column_top</mesh>
+ <mesh>soil_column_q8_top</mesh>
<type>Neumann</type>
<component>1</component>
<parameter>top_load</parameter>
</boundary_condition>
<boundary_condition>
- <mesh>soil_column_right</mesh>
+ <mesh>soil_column_q8_right</mesh>
<type>Neumann</type>
<component>0</component>
<parameter>side_load</parameter>
diff --git a/Tests/Data/Mechanics/InitialStates/non_equilibrium_initial_state.prj b/Tests/Data/Mechanics/InitialStates/non_equilibrium_initial_state.prj
index aeac9d4bf56b37520f34ab6204f2ecf2d3ebf70b..566ad1347c9735dc2a0987b30cc0a50057774764 100644
--- a/Tests/Data/Mechanics/InitialStates/non_equilibrium_initial_state.prj
+++ b/Tests/Data/Mechanics/InitialStates/non_equilibrium_initial_state.prj
@@ -2,10 +2,10 @@
<OpenGeoSysProject>
<meshes>
<mesh>soil_column_q8.vtu</mesh>
- <mesh>soil_column_bottom.vtu</mesh>
- <mesh>soil_column_top.vtu</mesh>
- <mesh>soil_column_left.vtu</mesh>
- <mesh>soil_column_right.vtu</mesh>
+ <mesh>soil_column_q8_bottom.vtu</mesh>
+ <mesh>soil_column_q8_top.vtu</mesh>
+ <mesh>soil_column_q8_left.vtu</mesh>
+ <mesh>soil_column_q8_right.vtu</mesh>
</meshes>
<processes>
<process>
@@ -124,7 +124,7 @@
<parameter>maximum_top_pressure</parameter>
</parameter>
<parameter>
- <mesh>soil_column_right</mesh>
+ <mesh>soil_column_q8_right</mesh>
<name>side_load</name>
<type>Function</type>
<expression>-0.8*2200.0*9.81*(100.0-y)</expression>
@@ -146,25 +146,25 @@
<boundary_conditions>
<!-- fixed boundaries -->
<boundary_condition>
- <mesh>soil_column_bottom</mesh>
+ <mesh>soil_column_q8_bottom</mesh>
<type>Dirichlet</type>
<component>1</component>
<parameter>Zero</parameter>
</boundary_condition>
<boundary_condition>
- <mesh>soil_column_left</mesh>
+ <mesh>soil_column_q8_left</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>Zero</parameter>
</boundary_condition>
<boundary_condition>
- <mesh>soil_column_top</mesh>
+ <mesh>soil_column_q8_top</mesh>
<type>Neumann</type>
<component>1</component>
<parameter>top_load</parameter>
</boundary_condition>
<boundary_condition>
- <mesh>soil_column_right</mesh>
+ <mesh>soil_column_q8_right</mesh>
<type>Neumann</type>
<component>0</component>
<parameter>Zero</parameter>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_bottom.vtu b/Tests/Data/Mechanics/InitialStates/soil_column_bottom.vtu
index c68af3d102c880bb315b922f13990afaf485978d..43ed5d0dee03a13ee77d4452fecccbfacba517cb 100644
--- a/Tests/Data/Mechanics/InitialStates/soil_column_bottom.vtu
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_bottom.vtu
@@ -1,33 +1,24 @@
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
<UnstructuredGrid>
- <Piece NumberOfPoints="5" NumberOfCells="2">
+ <Piece NumberOfPoints="3" NumberOfCells="2" >
<PointData>
- <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="0" RangeMax="124">
- KAAAAHwAAAAAAAAAAAAAAAAAAAACAAAAAAAAAAEAAAAAAAAAUwAAAAAAAAA=
- </DataArray>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="2" offset="0" />
</PointData>
<CellData>
- <DataArray type="UInt64" Name="bulk_element_ids" format="binary" RangeMin="0" RangeMax="1">
- EAAAAAAAAAAAAAAAAQAAAAAAAAA=
- </DataArray>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="1" offset="40" />
</CellData>
<Points>
- <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="10">
- eAAAAAAAAAAAAB5AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACRAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABRAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAARAAAAAAAAAAAAAAAAAAAAAAA==
- </DataArray>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="10" offset="68" />
</Points>
<Cells>
- <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="4">
- MAAAAAEAAAAAAAAAAwAAAAAAAAAEAAAAAAAAAAMAAAAAAAAAAgAAAAAAAAAAAAAAAAAAAA==
- </DataArray>
- <DataArray type="Int64" Name="offsets" format="binary" RangeMin="3" RangeMax="6">
- EAAAAAMAAAAAAAAABgAAAAAAAAA=
- </DataArray>
- <DataArray type="UInt8" Name="types" format="binary" RangeMin="21" RangeMax="21">
- AgAAABUV
- </DataArray>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="172" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="220" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="248" />
</Cells>
</Piece>
</UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _GAAAAAIAAAAAAAAAAAAAAAAAAAABAAAAAAAAAA==EAAAAAAAAAAAAAAAAQAAAAAAAAA=SAAAAAAAAAAAACRAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABRAAAAAAAAAAAAAAAAAAAAAAA==IAAAAAEAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAAAAAAAAAAAEAAAAAIAAAAAAAAABAAAAAAAAAA=AgAAAAMD
+ </AppendedData>
</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_cell_properties_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_cell_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..25a56b8b2e51801712a3ef30c75f5458926d326a
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_cell_properties_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_cell_properties_val2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_cell_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..05a13f7ad50b5b90d969e11282aed55d2d053567
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_cell_properties_val2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_msh_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_msh_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..9903ca1abe50deb4864cee34e170631c96618d2e
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_msh_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_msh_ele2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_msh_ele2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..7d0392c3f5170652bdda7bdf7255b1f8af19a804
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_msh_ele2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_msh_ele_g2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_msh_ele_g2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_msh_nod2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_msh_nod2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..397dfc238d27388fca6619bd85df6f3316326b09
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_msh_nod2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_node_properties_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_node_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..02dba888ca799ac0ebd79efd5517f11823592525
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_node_properties_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_node_properties_val2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_node_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..e8a95df72230c122826eacfbcccf0f364eb4855e
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_bottom_partitioned_node_properties_val2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_left.vtu b/Tests/Data/Mechanics/InitialStates/soil_column_left.vtu
index fe15564751da543a537179d5d0e8d8b7921c2d7c..750cac0f03b6e1cb3a7db1199aef3ba128aface5 100644
--- a/Tests/Data/Mechanics/InitialStates/soil_column_left.vtu
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_left.vtu
@@ -1,33 +1,24 @@
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
<UnstructuredGrid>
- <Piece NumberOfPoints="41" NumberOfCells="20">
+ <Piece NumberOfPoints="20" NumberOfCells="18" >
<PointData>
- <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="0" RangeMax="82">
- SAEAAAAAAAAAAAAAQAAAAAAAAAADAAAAAAAAAEEAAAAAAAAAQwAAAAAAAABEAAAAAAAAAD8AAAAAAAAADwAAAAAAAABCAAAAAAAAABIAAAAAAAAARgAAAAAAAABFAAAAAAAAABUAAAAAAAAARwAAAAAAAAAYAAAAAAAAAEgAAAAAAAAADAAAAAAAAAA8AAAAAAAAAAkAAAAAAAAAOQAAAAAAAABQAAAAAAAAAFIAAAAAAAAAJAAAAAAAAAAGAAAAAAAAADYAAAAAAAAATQAAAAAAAAAzAAAAAAAAAFEAAAAAAAAAIQAAAAAAAABKAAAAAAAAADAAAAAAAAAATwAAAAAAAAAtAAAAAAAAAB4AAAAAAAAATgAAAAAAAAAqAAAAAAAAACcAAAAAAAAATAAAAAAAAAAbAAAAAAAAAEsAAAAAAAAASQAAAAAAAAA=
- </DataArray>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="60" offset="0" />
</PointData>
<CellData>
- <DataArray type="UInt64" Name="bulk_element_ids" format="binary" RangeMin="0" RangeMax="38">
- oAAAACYAAAAAAAAAJAAAAAAAAAAiAAAAAAAAACAAAAAAAAAAHgAAAAAAAAAcAAAAAAAAABoAAAAAAAAAGAAAAAAAAAAWAAAAAAAAABQAAAAAAAAAEgAAAAAAAAAQAAAAAAAAAA4AAAAAAAAADAAAAAAAAAAKAAAAAAAAAAgAAAAAAAAABgAAAAAAAAAEAAAAAAAAAAIAAAAAAAAAAAAAAAAAAAA=
- </DataArray>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="38" offset="220" />
</CellData>
<Points>
- <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="100">
- 2AMAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAHkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAFEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAKUAAAAAAAAAAAAAAAAAAAAAAAAAAAACANkAAAAAAAAAAAAAAAAAAAAAAAAAAAACAO0AAAAAAAAAAAAAAAAAAAAAAAAAAAAAABEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAOUAAAAAAAAAAAAAAAAAAAAAAAAAAAACAMUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPkAAAAAAAAAAAAAAAAAAAAAAAAAAAADAQkAAAAAAAAAAAAAAAAAAAAAAAAAAAABAQEAAAAAAAAAAAAAAAAAAAAAAAAAAAACAQUAAAAAAAAAAAAAAAAAAAAAAAAAAAABARUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAREAAAAAAAAAAAAAAAAAAAAAAAAAAAADAR0AAAAAAAAAAAAAAAAAAAAAAAAAAAAAANEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAWUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAALkAAAAAAAAAAAAAAAAAAAAAAAAAAAADAV0AAAAAAAAAAAAAAAAAAAAAAAAAAAADgVUAAAAAAAAAAAAAAAAAAAAAAAAAAAABgWEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAATkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAJEAAAAAAAAAAAAAAAAAAAAAAAAAAAACAVkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAgUkAAAAAAAAAAAAAAAAAAAAAAAAAAAABAVUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAgV0AAAAAAAAAAAAAAAAAAAAAAAAAAAACAS0AAAAAAAAAAAAAAAAAAAAAAAAAAAADATEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAVEAAAAAAAAAAAAAAAAAAAAAAAAAAAACgVEAAAAAAAAAAAAAAAAAAAAAAAAAAAADAUkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAASUAAAAAAAAAAAAAAAAAAAAAAAAAAAABgU0AAAAAAAAAAAAAAAAAAAAAAAAAAAACAUUAAAAAAAAAAAAAAAAAAAAAAAAAAAABAUEAAAAAAAAAAAAAAAAAAAAAAAAAAAADgUEAAAAAAAAAAAAAAAAAAAAAAAAAAAACARkAAAAAAAAAAAAAAAAAAAAAAAAAAAABAT0AAAAAAAAAAAAAAAAAAAAAAAAAAAABASkAAAAAAAAAAAA==
- </DataArray>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="100" offset="420" />
</Points>
<Cells>
- <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="40">
- 4AEAABEAAAAAAAAAEwAAAAAAAAAVAAAAAAAAABMAAAAAAAAAGAAAAAAAAAAbAAAAAAAAABgAAAAAAAAAGgAAAAAAAAAUAAAAAAAAABoAAAAAAAAAHgAAAAAAAAAfAAAAAAAAAB4AAAAAAAAAIAAAAAAAAAAiAAAAAAAAACAAAAAAAAAAIwAAAAAAAAAZAAAAAAAAACMAAAAAAAAAJAAAAAAAAAAlAAAAAAAAACQAAAAAAAAAFgAAAAAAAAAnAAAAAAAAABYAAAAAAAAAHAAAAAAAAAAdAAAAAAAAABwAAAAAAAAAIQAAAAAAAAAoAAAAAAAAACEAAAAAAAAAJgAAAAAAAAAPAAAAAAAAACYAAAAAAAAADgAAAAAAAAANAAAAAAAAAA4AAAAAAAAADAAAAAAAAAAKAAAAAAAAAAwAAAAAAAAACQAAAAAAAAALAAAAAAAAAAkAAAAAAAAABwAAAAAAAAAFAAAAAAAAAAcAAAAAAAAAEAAAAAAAAAAEAAAAAAAAABAAAAAAAAAAEgAAAAAAAAAIAAAAAAAAABIAAAAAAAAAFwAAAAAAAAADAAAAAAAAABcAAAAAAAAAAgAAAAAAAAABAAAAAAAAAAIAAAAAAAAAAAAAAAAAAAAGAAAAAAAAAA==
- </DataArray>
- <DataArray type="Int64" Name="offsets" format="binary" RangeMin="3" RangeMax="60">
- oAAAAAMAAAAAAAAABgAAAAAAAAAJAAAAAAAAAAwAAAAAAAAADwAAAAAAAAASAAAAAAAAABUAAAAAAAAAGAAAAAAAAAAbAAAAAAAAAB4AAAAAAAAAIQAAAAAAAAAkAAAAAAAAACcAAAAAAAAAKgAAAAAAAAAtAAAAAAAAADAAAAAAAAAAMwAAAAAAAAA2AAAAAAAAADkAAAAAAAAAPAAAAAAAAAA=
- </DataArray>
- <DataArray type="UInt8" Name="types" format="binary" RangeMin="21" RangeMax="21">
- FAAAABUVFRUVFRUVFRUVFRUVFRUVFRUV
- </DataArray>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1068" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="1460" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1660" />
</Cells>
</Piece>
</UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _oAAAADMAAAAAAAAAAAAAAAAAAAAGAAAAAAAAADAAAAAAAAAANgAAAAAAAAAtAAAAAAAAACoAAAAAAAAAJwAAAAAAAAAkAAAAAAAAADwAAAAAAAAAIQAAAAAAAAA5AAAAAAAAAB4AAAAAAAAAGAAAAAAAAAAVAAAAAAAAABIAAAAAAAAADwAAAAAAAAAMAAAAAAAAAAkAAAAAAAAAAwAAAAAAAAA=kAAAACYAAAAAAAAAJAAAAAAAAAAiAAAAAAAAACAAAAAAAAAAHgAAAAAAAAAcAAAAAAAAABoAAAAAAAAAGAAAAAAAAAAWAAAAAAAAABQAAAAAAAAADgAAAAAAAAAMAAAAAAAAAAoAAAAAAAAACAAAAAAAAAAGAAAAAAAAAAQAAAAAAAAAAgAAAAAAAAAAAAAAAAAAAA==4AEAAAAAAAAAAAAAAAAAAABAVUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAJEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAVEAAAAAAAAAAAAAAAAAAAAAAAAAAAACAVkAAAAAAAAAAAAAAAAAAAAAAAAAAAADAUkAAAAAAAAAAAAAAAAAAAAAAAAAAAACAUUAAAAAAAAAAAAAAAAAAAAAAAAAAAABAUEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAATkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAWUAAAAAAAAAAAAAAAAAAAAAAAAAAAACAS0AAAAAAAAAAAAAAAAAAAAAAAAAAAADAV0AAAAAAAAAAAAAAAAAAAAAAAAAAAAAASUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAREAAAAAAAAAAAAAAAAAAAAAAAAAAAACAQUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAOUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAANEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAALkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAFEAAAAAAAAAAAA==IAEAAAkAAAAAAAAACwAAAAAAAAALAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAwAAAAAAAAADAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAHAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAKAAAAAAAAAAoAAAAAAAAADAAAAAAAAAANAAAAAAAAAA4AAAAAAAAADgAAAAAAAAAPAAAAAAAAAA8AAAAAAAAAEAAAAAAAAAAQAAAAAAAAABEAAAAAAAAAEQAAAAAAAAASAAAAAAAAABIAAAAAAAAAAgAAAAAAAAACAAAAAAAAABMAAAAAAAAAEwAAAAAAAAABAAAAAAAAAA==kAAAAAIAAAAAAAAABAAAAAAAAAAGAAAAAAAAAAgAAAAAAAAACgAAAAAAAAAMAAAAAAAAAA4AAAAAAAAAEAAAAAAAAAASAAAAAAAAABQAAAAAAAAAFgAAAAAAAAAYAAAAAAAAABoAAAAAAAAAHAAAAAAAAAAeAAAAAAAAACAAAAAAAAAAIgAAAAAAAAAkAAAAAAAAAA==EgAAAAMDAwMDAwMDAwMDAwMDAwMDAw==
+ </AppendedData>
</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_cell_properties_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_cell_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..fbe6218eccd824e69ecbc4dfd08f92dd5de0d139
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_cell_properties_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_cell_properties_val2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_cell_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..b1ab73110ec05ba7bb9374150b9e14aa5533b81d
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_cell_properties_val2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_msh_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_msh_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..56d689bd170b1c603f3b7265a99a480827f4e74a
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_msh_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_msh_ele2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_msh_ele2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..f680466d6335dfc9face022190804916181d461b
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_msh_ele2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_msh_ele_g2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_msh_ele_g2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_msh_nod2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_msh_nod2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..713af52b8468caab27878bd1dbc7f12c819d7e2d
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_msh_nod2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_node_properties_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_node_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..fa2afacbecea341917b930cb069307d0ee170821
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_node_properties_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_node_properties_val2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_node_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..56eda696176859e3739f45a76ee6da57069886c8
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_left_partitioned_node_properties_val2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma.vtu b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..378e50a7823ad1af6581bea161e76386dc0e1151
--- /dev/null
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma.vtu
@@ -0,0 +1,23 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="63" NumberOfCells="40" >
+ <PointData>
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="nonequilibrium_stress" NumberOfComponents="4" format="appended" RangeMin="555634.11785" RangeMax="3100953.7425" offset="0" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="100.49875621" offset="324" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="584" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="888" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1048" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAACAAAAABQAA4QAAAA==eJxjZGBgOLLe8iCQasi65wyiUQAjEfI9cRZg+eXc5MmH8JqD5Z9YOJElL7PHFCyvmOZIlvyTLBOwfNxkB7Lk10gag8XNG+3Jki85YQgW33jblix5m3IDsLi2mQ1Z8hw+emDxpROsyJKfsVIbLC7/2oIseQ12TbD4TFdzsuR3pKiBxYUXmJIl73FIGSze+8uYLPkb8opgcfZQI7LkL/ySAYs3rjcgS/7DVnGw+DdhPbLkBQqFwOKvTmiRJZ8QwQUWv1erTpa8w9z/B0DsS0YqZMk3CD8Fyx97rkCyPABN16JHAQAAAACAAADoBQAAsAAAAA==eJx9zjsKwmAQReEswiXc2srGRn5BBBHEiA/cTZaWpdnMbWbmeLp8JDczTf/ajN4FTt/5ObvA6T9+L7vA3RZ2sgvc7WAnu8DdHnayC9wdYCe7wKPl2O8UF7g7wT3ZBe57znBPdoG7C9yTXeC+5wr3ZBe4u8E92QUejXu/U1zg0TL3O8UFHq2Pfqe4wN2z3yku8Gi8+p3iAo+Wd79TXODR+ul3igvcffud4ir+AxBKO7k=AQAAAACAAAAABQAA0QAAAA==eJx1i0dSAgAAA+kq0iyIgqKggA3BhqLC/5/FwewlM8llJztJofCfolgRy+ZLYjXs6HtizTy//bCj18UD8/wOw47eFBvm+bXCjt4R2+b5HYUd/UQ8Ns/vNOzoZ2LXPL9e2NEvxHPz/PphR78UB+b5XYUd/Vocmud3E3b0sTgyz+827OgT8c48v2nY0e/FmXl+D2FHfxIfzfN7Djv6izg3z28RdvRXcWme31vY0T/Ed/P8PsOO/iWuzPP7Djv6j7g2z+837Ogb8c88v63tdiHvE2E=AQAAAACAAABAAQAAZgAAAA==eJwtxRECwgAAAMAgCAbBIBgMgiAYBMFgMBgMBkEQBIPBIAgGwWAwGAyCINhTg+7k1qu/jQNvHXrnyLH3PvjoxCefnTpz7sKlK1989c13127c+uGnO7/ce/DoybPf/vjrxT+ciwzRAQAAAACAAAAoAAAADAAAAA==eJzj5CQOAAAc/AFp
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise.prj b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise.prj
new file mode 100644
index 0000000000000000000000000000000000000000..e9b5672abde1723e6c6b84bfb7e13b280a7e598c
--- /dev/null
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise.prj
@@ -0,0 +1,218 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <meshes>
+ <mesh>soil_column_nonequilibrium_sigma.vtu</mesh>
+ <mesh>soil_column_bottom.vtu</mesh>
+ <mesh>soil_column_top.vtu</mesh>
+ <mesh>soil_column_left.vtu</mesh>
+ <mesh>soil_column_right.vtu</mesh>
+ </meshes>
+ <processes>
+ <process>
+ <name>SD</name>
+ <type>SMALL_DEFORMATION</type>
+ <integration_order>3</integration_order>
+ <constitutive_relation>
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <solid_density>rho_sr</solid_density>
+ <specific_body_force>0 0</specific_body_force>
+ <process_variables>
+ <process_variable>displacement</process_variable>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="sigma" output_name="sigma"/>
+ <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <initial_stress>nonequilibrium_stress</initial_stress>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="SD">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <compensate_non_equilibrium_initial_residuum>true</compensate_non_equilibrium_initial_residuum>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1e-12</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>1</t_initial>
+ <t_end>2</t_end>
+ <timesteps>
+ <pair>
+ <repeat>2</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>soil_column_nonequilibrium_sigma_elementwise</prefix>
+ <timesteps>
+ <pair>
+ <repeat>2</repeat>
+ <each_steps>1</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>displacement</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ </variables>
+ <suffix>_ts_{:timestep}</suffix>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <mesh>soil_column_nonequilibrium_sigma</mesh>
+ <name>nonequilibrium_stress</name>
+ <type>MeshElement</type>
+ <field_name>nonequilibrium_stress</field_name>
+ <!--
+ <expression>2200*9.81*0.8*(y-100)</expression>
+ <expression>-0.5*1e6+2200*9.81*(y-100)</expression>
+ <expression>0</expression>
+ <expression>0</expression>
+ -->
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1e9</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>.0</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>2200.0</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0</values>
+ </parameter>
+ <parameter>
+ <name>Dirichlet_left</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>Zero</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>maximum_top_pressure</name>
+ <type>Constant</type>
+ <value>-1e6</value>
+ </parameter>
+ <parameter>
+ <name>top_load</name>
+ <type>CurveScaled</type>
+ <curve>top_load_ramp</curve>
+ <parameter>maximum_top_pressure</parameter>
+ </parameter>
+ <parameter>
+ <mesh>soil_column_right</mesh>
+ <name>side_load</name>
+ <type>Function</type>
+ <expression>-0.8*2200.0*9.81*(100.0-y)</expression>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>top_load_ramp</name>
+ <coords>0.0 1.0 2.0</coords>
+ <values>0.0 0.0 0.5</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>2</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <!-- fixed boundaries -->
+ <boundary_condition>
+ <mesh>soil_column_bottom</mesh>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>Zero</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>soil_column_left</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>Zero</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>soil_column_top</mesh>
+ <type>Neumann</type>
+ <component>1</component>
+ <parameter>top_load</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>soil_column_right</mesh>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Zero</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>
+ <lis>-i cg -p jacobi -tol 1e-16 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ <precon_type>DIAGONAL</precon_type>
+ <!--solver_type>CG</solver_type>
+ <precon_type>DIAGONAL</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-16</error_tolerance-->
+ </eigen>
+ <petsc>
+ <prefix>sd</prefix>
+ <parameters>-sd_ksp_type cg -sd_pc_type bjacobi -sd_ksp_rtol 1e-16 -sd_ksp_max_it 10000</parameters>
+ </petsc>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <file>soil_column_nonequilibrium_sigma_elementwise_ts_1.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-15</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>soil_column_nonequilibrium_sigma_elementwise_ts_1.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-7</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise_ts_1.pvtu b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise_ts_1.pvtu
new file mode 100644
index 0000000000000000000000000000000000000000..e2f1dfe2db6595d2b5980803a896931108dd9ba7
--- /dev/null
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise_ts_1.pvtu
@@ -0,0 +1,43 @@
+<?xml version="1.0"?>
+<VTKFile type="PUnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
+ <PUnstructuredGrid GhostLevel="1">
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="97" format="binary" RangeMin="34" RangeMax="125">
+ AQAAAACAAABhAAAAVgAAAA==eF5VjDEOgDAIAP/C7KZTv2IMQUXDADS0Dqbp3+3qeHfJNRCrfAdVccPsg5Ai6C2Q1vaLHicHpHkCI2VIUORWQskwzKM7B/qFh2t2Y6tjsPStfwt7JC0=
+ </DataArray>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="40" format="binary" RangeMin="45" RangeMax="121">
+ AQAAAACAAAAoAAAAMAAAAA==eF4z0zPWM9Q1MtdNNzNOsTS3MEs110vJLCqp1DMyMDIwMDMyMDQxNTI0BwDHbAmw
+ </DataArray>
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="198" format="binary" RangeMin="2184903.3215410206" RangeMax="3531311.24553486">
+ AQAAAACAAADAGAAAzgIAAA==eF7t1s1r03Acx/FOKFNBfICCilAZCMpgOHyYghWfNj14FLx5UWQMD0MG3vQf2MWD6MGjPQjDucPAgx0UVhA9iEJg4EBxm8KE6dqkTxaiSZrkx2ffXz6tzSaCub74vrN8k3TpSiQSs5Mn86Zt36s/u5hPNI/Xnx6uvNmVzXVFeoG4ofW+dwtL/Y8n/9m+5w8ypE/cIJ5dt35zP7H0x6+dyFu/fWB+SNyf3gta9w/fTdKXvaD1v91/vjg7env0feZP59l+sc/22+79i7N/ZduAe31jW+T5/77Wm/sP9meBs/2q728rfc088w3p73t5zPWp44Pg07loN4gvEY+vb5G+7AZx1g/3h+4fur8fvdN+u479cpvzqi+OHHV95fqFwBfmR24dmP6a07m631bcJH10Ns98vfv+EZfj/uJ2dv5yB/MTe4645++9fz5w7/vqY07nzf1mot3Setx9k/TRsS972C+K/bekbwVeEudt7fVjX3P9ZL71fqWD/tirfvf6hmfOBd6Yyyzfvbqc0Tmb93z7jO8m6aNjX/awXyV9dOzLru+njB3Dc7VGMF8HV3+/pHn1+XS8TPro2Jc97Fsd9nFe7Z+6c9idz347K+5Pcm8+qXX/8N0kfXTso2N/lfTRsY+O/aI4bwT7K5Hz/xA9HXiF9NGxL3vYr5N+A9x7/7YG8z8j+psv97nn/7x7TT8X7RbxZNCvkr7sFvGwXyd92S3iYb8GvvfSRM+Nnv3k/CniSe3+sS/vN0U87JdI3yL9YkT/0dNe19NDZwJXf98cL4Orz6c0r74fjldJHx376NivkT56u30T3Ps+79a6//9xk8bV59fxCumjYx8d+3h/sC/fv7CP7vcdP9h9KP/dtu0XX07Dfuvu+SVXvx8k95/PjeqXSR8d++jq++V4CfzJzQ9Tg+mdWlefP8crpL9K+ujYR8c+7g/7JukXNf1ftHn7uw==
+ <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+ <Value index="0">
+ 2184903.3215
+ </Value>
+ <Value index="1">
+ 3531311.2455
+ </Value>
+ </InformationKey>
+ </DataArray>
+ </FieldData>
+ <PPointData>
+ <PDataArray type="Float64" Name="MaterialForces" NumberOfComponents="2"/>
+ <PDataArray type="Float64" Name="NodalForces" NumberOfComponents="2"/>
+ <PDataArray type="Float64" Name="displacement" NumberOfComponents="2"/>
+ <PDataArray type="Float64" Name="epsilon" NumberOfComponents="4"/>
+ <PDataArray type="Float64" Name="sigma" NumberOfComponents="4"/>
+ </PPointData>
+ <PCellData>
+ <PDataArray type="Float64" Name="nonequilibrium_stress" NumberOfComponents="4"/>
+ <PDataArray type="Float64" Name="principal_stress_values" NumberOfComponents="3"/>
+ <PDataArray type="Float64" Name="principal_stress_vector_1" NumberOfComponents="3"/>
+ <PDataArray type="Float64" Name="principal_stress_vector_2" NumberOfComponents="3"/>
+ <PDataArray type="Float64" Name="principal_stress_vector_3" NumberOfComponents="3"/>
+ </PCellData>
+ <PPoints>
+ <PDataArray type="Float64" Name="Points" NumberOfComponents="3"/>
+ </PPoints>
+ <Piece Source="soil_column_nonequilibrium_sigma_elementwise_ts_1_0.vtu"/>
+ <Piece Source="soil_column_nonequilibrium_sigma_elementwise_ts_1_1.vtu"/>
+ </PUnstructuredGrid>
+</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise_ts_1.vtu b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise_ts_1.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..a167d7837da46b534b9c0862568665f2b3f50421
--- /dev/null
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise_ts_1.vtu
@@ -0,0 +1,37 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="97" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45" RangeMax="103" offset="140" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="360" format="appended" RangeMin="1054838.5056" RangeMax="3531311.2455" offset="204" />
+ </FieldData>
+ <Piece NumberOfPoints="63" NumberOfCells="40" >
+ <PointData>
+ <DataArray type="Float64" Name="MaterialForces" NumberOfComponents="2" format="appended" RangeMin="134.8875" RangeMax="3880.30625" offset="2068" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="509010.45198" RangeMax="15521225" offset="2740" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0.05" offset="3320" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.0005" RangeMax="0.0005" offset="4160" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="1054838.5056" RangeMax="3531311.2455" offset="5772" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="nonequilibrium_stress" NumberOfComponents="4" format="appended" RangeMin="555634.11785" RangeMax="3100953.7425" offset="7148" />
+ <DataArray type="Float64" Name="principal_stress_values" NumberOfComponents="3" format="appended" RangeMin="1054838.5056" RangeMax="3531311.2455" offset="7472" />
+ <DataArray type="Float64" Name="principal_stress_vector_1" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1" offset="7956" />
+ <DataArray type="Float64" Name="principal_stress_vector_2" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1" offset="8072" />
+ <DataArray type="Float64" Name="principal_stress_vector_3" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1" offset="8188" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="100.49875621" offset="8244" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="8540" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="8844" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="9004" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAACAAABhAAAAVgAAAA==eF5VjDEOgDAIAP/C7KZTv2IMQUXDADS0Dqbp3+3qeHfJNRCrfAdVccPsg5Ai6C2Q1vaLHicHpHkCI2VIUORWQskwzKM7B/qFh2t2Y6tjsPStfwt7JC0=AQAAAACAAAAUAAAAHAAAAA==eF4z0zPWM9C1NLLQTU82Tk1KMTRPNQUAMLgFHQ==AQAAAACAAAAALQAAYwUAAA==eF7VmE+IG3UUxyfJTtwIoaUS1FosItpKoXRd+w/NaNW2Hjx4EPSkBS2U4kUoeLN7VPDioVRRDx56EKvFw0IVE1xREEWqh8VFChXrKGjpbpFpdc1OzO83eTPDm/f9vXSTLPRdJuyH7+f3e+83M7vZkud5X53ZO9fpdo8vf3xwzuvXt7+cvPzdhhPNkpPPKvxryLf/eDGcevfMkP75If14fxqncq9/ruXmFyDX/P3+YX4wfmrIfMLfeG6P5bvPHxDnh/ks5FTD+THX/P3zG7P/ApxfY379kYV/L0Ge5Bch7yw0/3z1mZ8hv3j+6Ev3zP6hnF+Y8hXFzzn3c57f/9P13ZYfq8nzk3h+fxJfa3+s+Dnnfs61/Wvrk98DnGq1/XP/qOdL7xfyd2Run79Nn++0+U927Wc8ao2Kxzc45/1RufsPtfkqfF6Zr3p+Ke9cp5/KnU/6++3oA3Z+l194nPOWwmE+eT67a+b/T/Fz3n9/wPxg/mjI9btK/xHsX5sf9/cuM67+uopf5n7bfD59+7Rdf9ubjxX8iFN/bh6mvKv4Oed+zvNleO8yU8z7AeLJ7//l1C/nk/vXsf+WBzi9v5CfCvl5Xl5/EXLuX1H2zzm9n4i7zvfYN1O2vyPtR8XzkTjdn24eQZ6vUfhLSn51PHm+JJ4vnE/eH3j/yf09xHwUf5Ser7m/rn++mT8W84n/oVd22P5P/bVPPB+J58uRT/2e4uc8X4586jfzKfKwNRxPzt/wWOR+gDi9X8t9vsI4FeU9sb/Mj3jZybP9y/1lfs75/tF8zOfJJ7db/utt4+O9y8w4ubb+anlF4WWFlzx5/1Tk55z7ER/UX1L8iJO/AvyGv/XBNvv33eYDj4h5iVOVFU75ruLnnIr8iFPe9O/yc05FfpBvUd4Dec9z9Rc2ez/vnhygf83POffLXJ+v1+fxKvdv5rf1pvvmzPXT34PC/Nw8hJyKuCfmM7/MQ8ipiPf67551+GUeQk5l5relxzuy356PY337/trqyJMfre/2h6kfrT8qv/kZmv/ZF++15z+9rzk2bvZ3f47nz1fiWp5zvj73a/tDvAJ4/vmT9pf//oz799vkryp+wFO/zDN/TeZB2ckzf53xpL91bZN/4su77Xw+eudBMe/i5v4+CLmf+uuMU5Gfc+5HnPw1xc859yPuOfN+W8vHTp58/zL+ScXPOVUMOPdr52c8px1+cL6Byf+0+S7Lt1zdW9i/+bwgcLr/SoAn+Uab/BOKn3Pu5zz5/83G1F9X/JxzP+fkj/v+mpwPyM85zbfk5FGT/JOKn3PKlxRO/qri55z70f7N5x+WN9n8+0/tEc//nMBp/jHglF8rf1nxc879KF9x+v3AzZO/fyRORXxC8cs883PO/fL8Mn9N8fM8lZnf4uytlt/x4S7ut8/nksCpfMBpf+SfUPyccz/Ip35z3ejwc879IG/n68qXFU7rLwM+Kj/i4/ab+a1/eYN9Pk74OwvzmwCcquLOt8lfVfyccz/Ip/4ayJMfcfJznr//DK+DPOL0/Jbd/qAEOPfL+dH55fPL/K7zO/TszTa/7tB0YX5lhXtO3mjf2P7k/jG8KvPAfC9/Huaj5rU+R3nyg/UH9qM8+WuKH/FrCjf5h9/rfmHWf/2zqcL8rip8SeHkryp+xJcUTv6a4kd8SeGX3P0FV/rc/H56rcCjZuTkfkB+3h/VFcC5H8wn9cdsfe6X+8v8qP/ve/z4LaHtr9zYUZjf2wo/rHDym/UrDr/Mo9ZhwKnQ/qnQ/qnQ/qn+cfbXCO7sfZ7p8b8Zp+83+/sczCcgvzn/Yv+Zn3Pul/OZX55v5o/F/UVN5Kf6H6rwkSw=AQAAAACAAADwAwAA5QEAAA==eF5dkzFoE2EUxzM6WizB7g4OQoMBcfnwQUUyWDhcOpTKgfFyisNRCjmCQ5DSphIwhIYmKjEajTHGkEHFIh95g4hCwBscoji00KFDh4AgEYX0u3svIV9ueXz3+/P//v/HXaPb7T42yviWJoSCJy3eB+cWvOP3n3gmS6XS7tc3aNy0LOuXgaRfk25wNjFFHByLeDgajf7+2MRQYlLvyIFF+lnisMd+u+qec63XWNH8TRljfZv42P+aPysNdOxJ/V1xNkH6JeJgsP5oOBw+yb/CtuafFRn27xMHj/VF1efHeh3bluYvgf0fEIcRv6D6nEq+xL7mb4s072eeOHjMv6k+V2/Vpvq60uU8XeLj/fhjY/kFgqYvSJP9bxOHHPP/qk9n8TlGtP3fE8j6f8ThJOvzqs/fS1UEra8rB9z3IXFA3veC6hM9/wxzU/rTnP8KcfDixA9UnztnnmJay+OKkE36feLgMV9ToxauoHdjUp+SvTjpXeJgsv+sojhTxg9ani2R4fwniEOaeczP9+cR1rV9OiLO+S8ThyLrU/59P0t4qOnvyybvc5M4mMwbal9znR20tb62iHDfOvHx99ZT/bPVAl6c+n6+sP934jD6//ycnwvruMKT9K5cDc7b6NB7uM78GHpiy0U=AQAAAACAAADwAwAAnwEAAA==eF5dk7FLw0AUh1Nbi4MtGUQkFIybiENBBKkl5EihRRCKrViUUgfBpc51cMjo3+BU6OJYnHQw9mUW6SidOolj0aVU5Xy5u4a763K8fB/vfu8lTRiGP8pUyS+l6093TWLwn8/rK5LgHH4EX0D4VWkRE89ixwThO7y2QXAYU0r3kaewbveaZCJq4Qe85j5ymPOo/vg7I0PNt6X+yBX/vHGq+/6m5COHoeSPHhpsTil/8CjlQa70r2dOdN+fSj5ykO9/uzzW5jWdPSkP8rg/HrQc1khe6W/G+4yqSlhT8oe5I2Ir/pqyf+TsecST7KySvpLHDvKSj1yZt9071N5vzjGkPMjBF37aiPZ7oPVfCrpiPynOwRX9Ir+eKWv7mTzP8yQ5h76UpxJ65FbxV5wb4Sc4j/PP2OmSsZanKuUvdlyY73uR5Sto867G+fl8BeiKfll23w75VPs7F6J/ivP4+8f/EW01tsi94i8H28L/5hymgr+z2iK2nMdPB9fCz7I8VryfDebPXLznxbNK3M/vOriXgWfV4ZXzgcHqEvwDHJFgeg==AQAAAACAAADwAwAAYgIAAA==eF5jYEAFs2zcxeMNDPbA+CfiqgU42azhfBjIFlnn/rAqZX8Pe/CW246qezKhfIsFXrUNV9T3FEL5MPW5YH7Jfjb1Wwue5bzYDePH/BPfPnvG6915UD5M/RK51teBO+r2+1fuiV52/MluGH9DcFFi3Ntf1sugfJj6fLD+lv02G5i2/f3GvScPym9I2zNvfcsjaxgfpr5zJgh07v/NwS3cLqJlA+OLzD7qxsKoZdMN5aO6p2//LiZ9X/PJxjZLoXyJqVZG7XpucD5M/YPAHUCRifsbbvetZ1Tg3QPjpzXJ/HjromMD48PUQ9w3Zf8FZfkt2t+eW8P4ElPXb8p9pLQHxuf4f6luC5OBzd9riwtsuabv969Z2/dM4cfuP1A+//2PYTwHv+2GycPM7wX7Z+b+x0Y3105mVYT6b+b+EwnPTgSraNr0QPkw9dIx/Ye+aszZzyt38uRsOa09MH7M9kVHzwarw/kw9avA/p+3X2+q2NQtotJ7YHzHB8o+hZJqcD5M/Vyg664vXrD/5OfZrGFLdWxgfA4/d9dlcQZwPkz9E3B4Ldy/2CZmqyaj0h4Y/8fZLt6qXwpwPky9jjEILN7vFR61cepkLRsYf0bZ0nXZB0zgfJj6YnD4Ltmfw/gnXUVVew+Mf9vFtmZPiCqcD1O/+6sGMASW7p8hGuJ85ZXIHhh/Sr6D95oYJTgfpp75Oig+lu9nKOAyDQyPsobxX4QmPFWpEtkDl4cCL5B3Alfsl5XOaP07Q9cGxt/UvCrsZbIZnA9TPwkcfSv3T7D7Ed1UrLQHxj9xJCNYMFIKzgcAHYSeXA==AQAAAACAAADgBwAApwQAAA==eF5l1WlQk0cYB3COEFIgHAkRKOHKcIkhohyh3ZVBUYqMUipIvFoEFEeUwlAGsXgAqWIprVJAW8EDVIo0UDniwa5Cua0HFAQUuWODIAJJUG7oTMcv3ff9+pv/s8+zu7Mvi7HOT7g8h8ZuKtq2hntXq334dJPx1wdTqpEqnF8TUqBCo4RDC5ZeeHw7uJUkjB1Lf43eEH6NJs7mn2pHqQ5XuTeMVRQ/Z7K5trh/FBzt/8ojPUtBqd/skLZsxBpCsRcsAmNYk+gt4V2T1wakMTPI1nHMzSi3g9p/4e9/lbvWg8p/0hrsQAkYJ/xmWXWW2E4Dp9sEFHaP9FD8ddh5i+HrTLxX/baA9vkhzwnCD1se3+abwsHcfkmutVQFFIRHWPAtM+Xa+OP/XBOS8x8XjNTeaJYh+UvNx62JCqQk/Eb9O/PLPxnjqgG/TJlkljLf2HrzNuG+CVCr9F2VNKmGpwjPt47zD+BzYfH3CzUHE9Uh6ZWGq0ODsngwx+3wKoeuYTBJuEa0yOdAzySIu29WssNVjsi8aahQRRM44CrFFb+oTg6cJvz6d5suPmhUgjz3lLmlhBVwmXAZcNm1mPYEHTsaLskyMaH0V+q/66GX+TDQ6YthXojQg0uEZ/mVFDY8NsPGddpvunXV4Bzhho75+UU7adBk/eLG91d7wXvCn57gnbrkaQMv57uKz2gx4CLhReN3K4CBOXTStK/LlLxAZB7a03NebGPBfbxFlSRWHy8Qrs2ImfvkVy7sM1j5NjJpmrK/Sj2BHWszB3bNbpcktMoAmc+WHktaG8uAp2/9Qk8VPUHk/D7MgFtHUjpRy/Na5pkCBuX8NWYcv5Aa6eCGTzlxCiAH5P54neXUM+8PIg6KjPAeMKKcX7LlvR+9vtGC/rSw/QuvtLBa8f89iV+/Ep7vBP5T8iObAwwwWR8V/nA9vk4FeHEemx7tNsZ0Ip+ZgHXvD1lBz+D9NZp0BqW+Y+qehVKXJVTu3NrKEykRmefOivbk87lY42LwAzGSAz3Ck+d5q3c7C7B1u5VFVRGbcn9sFqMkGQ9s8R0nzxIsNIMriPw0u33ib38WDEwNXFfdaAnJ/h6F9/k25TpD7XLrxIrEQaBOeIK6l9/0OQ8c3cIYbAoZQh8RXiP9LK9DYoP3OZpdsmHPI7J/k15PLcslLiwbE915/rM+JH1Db1h6qmQNvmpwWrciZAjMk++PUe5ej5M83Pcqe83JA9TzaxHbN30bMooy7Be6ojXZmEZ4ZGjP4bztLrAob8I0KFgDaxF+Qtwd1H3PEDaDzEGRrIEy/8sP9b8UMNqkbnOU+xtZaeQeEGGPy5u9zdhlSsr+eD8Ly1jsqEU6jcqagV1yxCDc5m6lzlPgCHNstxbJDo0Dsr7TxYHsLdlrcRS/hSc+oQvJ+5MsgxOXAlWoI62hhF5mAsn5Ry6nTD0tZUCroo3l6b+x4Qz5/t+efuYudcJ5Ow0LtjzUh7OEq95tKHtSzYH0iAv9Yh86ZX3W+R0Zc3qu0F02VODFEEIVkT82LAL6TDvsX0Abr3huCsn/y5U/ks52h/Ggre2foodCO2xA1I/XbVdsiWfjfwGMUQxrAQAAAACAAADgBwAA9AMAAA==eF511X1ME2ccB/DrOy9GkK6JghZ1sDJfAiIO2vU6qnWdCc7gdJgNRjcbmcLG1EB8jQeBuLFpZFvEbSFTDDPRRQfRqdvu1utEYQPFTl4MKuoG4ru8lLlrr909j9eHy5O1/zyX59Pv7+nzu+euaoIgzp0wsnwwuIs7bmcJ8TNLRkZQefW0XHQO84aMtd15+mYyEAwGfxP8Keb23soU3/xmmE/RGFm/8L0DSyfc0z1cmqZpI4kw7rQNJRRccdEh92H+8dncr3TubuhPTmXCfNemV5GfM5r8VZ0NZMh5zB27M3u2DN4Mm6eW9Ref/uY+LROuzzhfguvHHVyKvLzN3pKqGYN54DzmW3N+VuYvHqSF6eDp/8knGXpiPCVuWJ/SLmJBH1d02JDfu5nfuSB5FNYH7sO85tbKH1aYO8iQ4/ljm/tO+swc9NfcC6F/xi1Brq19m6dzZJaQ85i369dXPSy6RIbzvUOFZTv2qxhwHbsxnQX7bDNM+La/kzIO7bhMhhyMv0t8PJHn26hjsD8xggewfMPdpneI6haY701cAOc1qxcjH2ncVzhP85QOOaijlnjfT+wH/aO3UT6AecKqsvaeSMIC+p9dnwrnbZVW5FZH8bSdC/8l5aKD/BKJryNKzzfFDiPH8zJvcH1Jw7P+d3Lz4PmiTmQjn5pW+0uPZxQ5GCskbt5Y2viJezCsU2Mv3mn1noT7d6yZA+tvP/sK8nfzZNEO+xDsb6HgAcw9L9TtT6120QHR8Xz/9Lvf67m/SCUBng8DC/o0rrUgv23ebJ0ZeQruHzio45X41C6b0V/kQa7A8t1H3qpIauahU9pkFowffWhGPu3P/F8HZlwTz38yq8S844AtbhnJIddgfnGLalN7/CQmAM/XbDh/r9WEPFB5/WhXroLhRY/G3Gpf93pm5DAZcjnmNyqqiKwvouHzk12fCH+f83kj8n3b5iqylDpGJroG8/yCL+ueKyfC+mw64ch7hI/0w/sTz0aANXdmIn85a67pQlksfP4cgoP7fF3icVHF9XKTjJGLLsfybhNXVLK7lVYRoH86FoxrehchP3j4YXVfyWSLQnQl5gXfUjWfGx6L7x8drC/1DrG+mgD7i4HznvQM5O9fq/mx6aqSkYuuxnym5UH8lZpDdMhVmDuPb3iTqVXD/lNaFVw/Z0868nKiTvfd8imM8L8V3CW4CvOVdvWnEfxl2is6np8RNeuNq/5B8hF0rwuco5Y7acg3nLnobeyPYEbCeN52+o9i6xg5LjpY/7zEt84xrY5iplhWQR9wgfe7XDfhWa4UbuRrBeMRncPcOT3yn6N6mSWXAPsfcEUJo0Lih5cbUm36Efo/ZT2wQw==AQAAAACAAAAABQAA4QAAAA==eF5jZGBgOLLe8iCQasi65wyiUQAjEfI9cRZg+eXc5MmH8JqD5Z9YOJElL7PHFCyvmOZIlvyTLBOwfNxkB7Lk10gag8XNG+3Jki85YQgW33jblix5m3IDsLi2mQ1Z8hw+emDxpROsyJKfsVIbLC7/2oIseQ12TbD4TFdzsuR3pKiBxYUXmJIl73FIGSze+8uYLPkb8opgcfZQI7LkL/ySAYs3rjcgS/7DVnGw+DdhPbLkBQqFwOKvTmiRJZ8QwQUWv1erTpa8w9z/B0DsS0YqZMk3CD8Fyx97rkCyPABN16JHAQAAAACAAADAAwAAWAEAAA==eF778/9//a917gcZGRgYjqy3PMgABX/wiJvfcQOL98RZwMX/4hD/BxQv4XQ7yARkh/CaY4gzoIn/B4pvNHMFq5fZY4piL0icEYv4u2QXsPiTLBMU87GJg8zXnugM5q+RNCYoDgQNGfucwPySE4ZEiS997Qjm25QboJgPEmfEIv5IAqKew0cPxRyQOBOaOMhf8m4OYHNmrNRGUQ8SZ8AUZ9j5zA6sXoNdk6A4yD5jR1uw+I4UNbg4Kw5xDiBeO9saLO5xSBkuzgYVZ0ATZwdi9W+Q9HRDXhEuzolDnAWIFwZA0tOFXzIY4kxYxKVXm4Hd82GrOHI4/5cCijOgiYP8NZXV9CBIn0ChEIp7QOKMaOIgPn+C8UFmIJ0QwYViL0gcmB7/xyOJg8zv3GV48CtQ3H7u/wPI5oDEP6CJg/QziRocnAlKGMJP4eIgwAwUlwfSDUjiAMNC6MU=AQAAAACAAADAAwAARAAAAA==eF5jYEAGH+wZsIJR8VHxgRGX7/392lb6kg2x6mkt3tBcr/B3/kmi3fPE8+Sq8vSjRKvHJV6vfkab8/3pPejiABlPQC8=AQAAAACAAADAAwAARAAAAA==eF5jYACBD/YMWMGo+Kj44BCX7/392lb60h5i1dNavKG5XuHv/JNEu+eJ58lV5elHiVaPS7xe/Yw25/vTNjARAPpuQS8=AQAAAACAAADAAwAAGAAAAA==eF5jYMAGPthjFR4VHxUfFR9U4gDqgi9ZAQAAAACAAADoBQAAzAAAAA==eF6N0DsKAkEURNFZhEvo2MjERFoQQQRR8YO76aXN0kzqIty2xcrm4JQ1b5p+ZVG/exk48Xs82+n593/4nZ3n0a6lnB4779vJSk6PnR47WcvpsdNjJxs5PXZ67EnbytPTeRk42cnZY2ePPWl7OXvs7LGTg5w9dvbYk3aUs8fOHjs5ydljZ489qWd5ejovA0/aRc532fkuezJf5enpvAyc3OTcx8597Em9y7mPnfvYk/aQcx8797En81POfezcx05ecu5j5z4ffwMQSju5AQAAAACAAAAABQAA0QAAAA==eF51i0dSAgAQBAEREQQUBDERjGTMGf7/LA5OX6Zq9tI1Xb2Fwv8VxbK4Z74k7oeOfSBWzPNXDR27Jh6a568eOnZDPDLPXzN07GOxZZ6/k9CxO2LbPH+noWP3xK55/s5Cxz4X++b5uwgd+0q8NM/fdejYQ3Fgnr9R6Ng34tj8SLwNHftevDPP30Po2BPx0Tx/09Cx5+LMPH+L0LFX4tI8f+vQsZ/FJ/P8vYSO/Sa+mufvPXTsT/HDPH9foWP/iN/m+fsNHXsj/pnnb2vdDiHvE2E=AQAAAACAAABAAQAAZgAAAA==eF4txRECwgAAAMAgCAbBIBgMgiAYBMFgMBgMBkEQBIPBIAgGwWAwGAyCINhTg+7k1qu/jQNvHXrnyLH3PvjoxCefnTpz7sKlK1989c13127c+uGnO7/ce/DoybPf/vjrxT+ciwzRAQAAAACAAAAoAAAADAAAAA==eF7j5CQOAAAc/AFp
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise_ts_1_0.vtu b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise_ts_1_0.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..d647911fc780ec5c073d014c31fb9f5cd0d0058f
--- /dev/null
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise_ts_1_0.vtu
@@ -0,0 +1,37 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="97" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="40" format="appended" RangeMin="45" RangeMax="121" offset="140" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="198" format="appended" RangeMin="2184903.3215" RangeMax="3531311.2455" offset="228" />
+ </FieldData>
+ <Piece NumberOfPoints="36" NumberOfCells="22" >
+ <PointData>
+ <DataArray type="Float64" Name="MaterialForces" NumberOfComponents="2" format="appended" RangeMin="95.379865947" RangeMax="3880.30625" offset="1212" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="539550" RangeMax="15521225" offset="1628" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0.0275" offset="1976" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.0005" RangeMax="0.0005" offset="2480" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="2118783.7139" RangeMax="3531311.2455" offset="3452" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="nonequilibrium_stress" NumberOfComponents="4" format="appended" RangeMin="1731663.2249" RangeMax="3100953.7425" offset="4268" />
+ <DataArray type="Float64" Name="principal_stress_values" NumberOfComponents="3" format="appended" RangeMin="2184903.3215" RangeMax="3531311.2455" offset="4468" />
+ <DataArray type="Float64" Name="principal_stress_vector_1" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1" offset="4724" />
+ <DataArray type="Float64" Name="principal_stress_vector_2" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1" offset="4776" />
+ <DataArray type="Float64" Name="principal_stress_vector_3" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1" offset="4828" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="55.901699437" offset="4880" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5072" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5272" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5384" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAACAAABhAAAAVgAAAA==eF5VjDEOgDAIAP/C7KZTv2IMQUXDADS0Dqbp3+3qeHfJNRCrfAdVccPsg5Ai6C2Q1vaLHicHpHkCI2VIUORWQskwzKM7B/qFh2t2Y6tjsPStfwt7JC0=AQAAAACAAAAoAAAAMAAAAA==eF4z0zPWM9Q1MtdNNzNOsTS3MEs110vJLCqp1DMyMDIwMDMyMDQxNTI0BwDHbAmwAQAAAACAAADAGAAAzgIAAA==eF7t1s1r03Acx/FOKFNBfICCilAZCMpgOHyYghWfNj14FLx5UWQMD0MG3vQf2MWD6MGjPQjDucPAgx0UVhA9iEJg4EBxm8KE6dqkTxaiSZrkx2ffXz6tzSaCub74vrN8k3TpSiQSs5Mn86Zt36s/u5hPNI/Xnx6uvNmVzXVFeoG4ofW+dwtL/Y8n/9m+5w8ypE/cIJ5dt35zP7H0x6+dyFu/fWB+SNyf3gta9w/fTdKXvaD1v91/vjg7env0feZP59l+sc/22+79i7N/ZduAe31jW+T5/77Wm/sP9meBs/2q728rfc088w3p73t5zPWp44Pg07loN4gvEY+vb5G+7AZx1g/3h+4fur8fvdN+u479cpvzqi+OHHV95fqFwBfmR24dmP6a07m631bcJH10Ns98vfv+EZfj/uJ2dv5yB/MTe4645++9fz5w7/vqY07nzf1mot3Setx9k/TRsS972C+K/bekbwVeEudt7fVjX3P9ZL71fqWD/tirfvf6hmfOBd6Yyyzfvbqc0Tmb93z7jO8m6aNjX/awXyV9dOzLru+njB3Dc7VGMF8HV3+/pHn1+XS8TPro2Jc97Fsd9nFe7Z+6c9idz347K+5Pcm8+qXX/8N0kfXTso2N/lfTRsY+O/aI4bwT7K5Hz/xA9HXiF9NGxL3vYr5N+A9x7/7YG8z8j+psv97nn/7x7TT8X7RbxZNCvkr7sFvGwXyd92S3iYb8GvvfSRM+Nnv3k/CniSe3+sS/vN0U87JdI3yL9YkT/0dNe19NDZwJXf98cL4Orz6c0r74fjldJHx376NivkT56u30T3Ps+79a6//9xk8bV59fxCumjYx8d+3h/sC/fv7CP7vcdP9h9KP/dtu0XX07Dfuvu+SVXvx8k95/PjeqXSR8d++jq++V4CfzJzQ9Tg+mdWlefP8crpL9K+ujYR8c+7g/7JukXNf1ftHn7uw==AQAAAACAAABAAgAAJQEAAA==eF6LO3PmzJyAeQfUzoJpBwYwyNirAeavc0iAyDtoQuV1Zs6cufPk2gMCaUBwJ+AARH2C7YtUED/hgDFE3iEhFSJ/3MjY+PPuNQcWpCOrb7L9kA5RfwEi7/ABKp8NtE9n3eoDBanI6rv2GqRB1GdB5B0MoPZzgugFq9DM79r7AaqeCyLvUABVv/r////zJ608wIDi/oq9GlD3rITIOwRAzfMB+udmy4oDB1DUT4Kr94PIOzRA5W8A/SNcvvxAQAay+vlw97+DyDskQOWrgP7xyVyG5t+CvRug5rdB5B0UoOZrAam26KUHHqQgq5+8d0UKRL0GRB4enjD7E6DmnQL6b7/vEnj4QflQdT22aUD//LRffCAFQsPNeQHVbwDlb4DSAJf6DJk=AQAAAACAAABAAgAA8QAAAA==eF5jZGBouMUb4Pj5/3/5nbNiHRkgoOETmJ/ryAiRP/gJKs8ElPzkEe8oAKRtygUOwtRD+AoHofJgcZA8iJ+zGGIukvp9P/7//28NVQ+UP/gAzIeof/Y3Gq/5QHm4ehYgPyEiyvEDlA9TL4GkHiiP4p5bmyIw1L9Acg9QHsU9Ibzh6O63Q1YPlEcx/1x6qOMGFPMV7BKg7mGAyB+8gCTvcSgYTT2DXQbUfKj8wQ1Q85mB9CGZIEcFFPcI7POAms8CkYe7ByL/Ye8MJPfalAcchOmH8h0noNgvsE+DAeHenMW+jjD3skD4BydA9QMARHiy7A==AQAAAACAAABAAgAAaAEAAA==eF5jYEAF4VqK6XOOSuyB8a2aZKreiCrC+TBQJ7LO/WFVyn6Z2VsbA/hV9jRA+TP1c1kem6vsqYfyYeoh8iX7nbKO+Apc0NrTBOU3FZ1563/aYE8zlA9Tv1eu9XXgjrr9dvNnf758nnvPASj/1p7wrcFzf+0+BOWjmt+yX1tHQWaVsqINjC/jtPLshX1aNo1QPkz9gpkg0LnfJl38ZZKnos1CKD813lHyVaMunA9TfxBsX9/+R3qvsw91idjsh/InONy3iXBVtIHJw9R/CNwBFJm43+XhMafQaQo276H8oNpjXD7zLGxg8jD1kPCcsr9J9uozEyYxm3oo/8QbvvriH5Zw/jz9kNP3v6na8F5fXGDLNX3/H53Aq4mLdGz4oHzz98eXBXDbwPla/F8TVYL1beaB/TNz/y0+du2noVY2c6F8mP0weSsho6rYk4J79GL6D33VmAOXh/HT5I8lJlpo7NGF8gEUguYKAQAAAACAAACABAAAxQIAAA==eF6rsbjFcLT6we5TKz9e8k1y2M8ABSu6Es8eWzzfOiJu27wf+dd2n0STLz+9ds7Hy3t2N00vKD934fju42jyLj8uaX8x22IdN+O/S7/FSwzz387snmhe+XN3w+89UWXHLmDIB+wXj1M4Mst6j2XrFcMnF6zR7VdaZyNbV/LTmqtHXPps9p3dp9Hk7zWrHa8K22U93e/1ny/TTmG4n/mYYNjDpre7n1w7/KjRdy6G+W56LUyeGv93Gxj+UH9wYxOGfoNoTr7PBxdYa29duHBj7SFrdPcXLY2p3Pn3qXXRwfmrXxpfxzB/RtSbByvLv1rreXCuksh7juH/e7MfMiaZc+7piwrUKQi/iiF/7jjX8cLfv3aL5O5/67x1g/VZNHmt/ZOC/jic2b3F6JeXm+uT3efQ5B22mN+JrGDbc+Vf8qEo3ie7j6LJLzbOD/018fvun6fOWuZffYyhf8HX8snssw7udv15SvJN6u/dJ9Dk7+RyS91UkbLZuvfs3or1H3ZfQJN3mpTWxVX4Z3cbh1O76JPru8+gyTdHz9dfr6eyp3by36Pifz5ixO+tR/NNWstM9hTYvOua+ucrhvwqxXMWQg+k96Q8vecv/Pfn7lto8h1RL5UUOmVshKYdULkyf9HuS2jy6fNt/pfOeWP9RiuHLevbFWt086fuDG/S/GC0x8tkh0e6wjNrdPNPvBRO2irivMf9wUx23+RbGP7/fV7qhiiX1h6l3vsJDy2WWl9Bkw+dq7Q5Pl9lT9fhi/GtG+9bo+sPuBoiE7xSdk+i/antustf7EbXf8hI0zj/gsaerre+4gLr1u2+jSbPMserwrfJeM+X8/zmezd+x/D/r6iq120Wont4O5c+6Sz6YH0PTX6pGQfP7FbNPezKi3/MrvhqfQ1N3qNoC0NfueYenXv+C9/feGaNHr9uZpv0eHzU9wAAwHvj0g==AQAAAACAAACABAAAUAIAAA==eF5jYWBgOLLe8uD3///rf61zP8gABSemzbdoLK+1ZYLKf0GTd7pweJv0zVl7/wPBYaD8JzT5K7OSFHsvT7JlBLI12C0PfgOqm+GKkDf0rfKL3nl9LwNU/guavEPZMvuutxW2/4Di6uxg81HkL4VJvPqy5TrY/A9bzcHmXy1yQ8hb/UrPaZ9hywCV/4Im/23LsnPsNaf3wuQ/ocmviFCbWsdydy/I/ztSzMD6hRa4wuUPiD72Xm5WDzYfm3zDtX+vxKqP2ILCZztQ/jOavIbXeV+bf5fB+huETcH6/c+6IPTfiDXhrnm+Fyb/GU0+7t9jMX7Fm+Dwr8ei/0dskoZdwKq9oPDxOGR88CtQvueXM1ze9b5EtPfbJ2DzQfKf0eSvmdzLnfL5Klz+C5r8hIhQnYkpC8DyAoVGYPmT6gj5X3nbgowjf9jiki+6taf2f80tnPoV8k513DzNtA8U/jfkDcHxyxbqBJf/VCw3+QCf8D4GqPwXNPkwq4LPv7/+3AuTR9ev0R+wLHfyL3D6dpirD5Z3bnKEy2dV8U3MX3oKnL5A8l/R5Ce08r4zXiy4DxT+9kD572jyf29YudgeVQS7/8IvHbD+hvUOcPmQVXOnuzlwgN0Pkv+GJr+L8eM2CQ5msPnnseg32dSx5ur832D/JURoHXwDlK/eaQ+XT+Ofk6Alw74PJv8eTZ5t/6QOP20hsPtA8m/R5BsdVRPyD3+G5g91cP74KmwHl5d7cvfnkyMQ80Hy79Dkj9lu+H5uO0L+BZq8nVTygeJNrPsA8ljNCA==AQAAAACAAADAAgAAggAAAA==eF5jZGBgOLLe8iCQasi65wyiUQAjEfI9cRZg+eXc5MmH8JqD5Z9YOJElL7PHFCyvmOZIlvyTLBOwfNxkB7Lk10gag8XNG+3Jki85YQgW33jblix5m3IDsLi2mQ1Z8hw+emDxpROsyJKfsVIbLC7/2oIseQ12TbD4TFdzkuUBN99XPQ==AQAAAACAAAAQAgAArQAAAA==eF778v9//a917geZGBgYjqy3PMgABV/wiJvfcQOL98RZwMW/QsUZ0cRB6ks4IeIhvOYExUHmbDRzBYvL7DElKA4y512yC1j8SZYJinps4iD12hOdweJrJI0JioPMydjnBBYvOWEIF/+BQxxkztLXjmBxm3IDuPg3HOIgcx5JOILDk8NHDy7+HYc4yBx5NwewOTNWaqPYi0380////3c8swOLa7BrEhQHAGqfopQ=AQAAAACAAAAQAgAAFQAAAA==eF5jYEAGH+wZsIJR8ZEiDgCkaBoLAQAAAACAAAAQAgAAFQAAAA==eF5jYACBD/YMWMGo+EgTBwB0xxoLAQAAAACAAAAQAgAAFAAAAA==eF5jYMAGPthjFR4VH/biANQJGgs=AQAAAACAAABgAwAAfAAAAA==eF6NzjEKgDAQRNEcwiOktrKxEQURxEvl6Db7EX5ccbo8ksmU8pVhffeaOPE7znZ6/v7DPTvnbNcop8fOezuZ5PTY6bGTWU6PnR47WeT02OmxR9omj57Oa+Jkl7PHzh57pB1y9tjZk9wvp5x/7cTOv5c80jl7Hr8Bt4EcsQ==AQAAAACAAADAAgAAggAAAA==eF51jEkSggAQA0FAWVQEBRdwQf//Rw+kL6mauXSlK5kkWS8VczEzvxGLoEfeiVvz7MqgR67Fyjy7JuiRD+LePLtj0COfxNY8uy7okc9ib57dJeiRR3Ewz+4a9Mh38Wae3SPokWdxMs/uaT2If4svkT/08B8jfhG/5vn7s94fxaMGAg==AQAAAACAAACwAAAAQQAAAA==eF4txRECgDAAAMAgCAaDIBgMBkEQBMEgGARBEPT/lwTdyfXdb3Bw9OjJydnFsxev3ry7+nDz6cu3H7/+AP2wA/U=AQAAAACAAAAWAAAACwAAAA==eF7j5MQGAAj7AMc=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise_ts_1_1.vtu b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise_ts_1_1.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..663c1f3c0f73afd15439113271e11a999e254a7a
--- /dev/null
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_elementwise_ts_1_1.vtu
@@ -0,0 +1,37 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="97" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="40" format="appended" RangeMin="45" RangeMax="121" offset="140" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="189" format="appended" RangeMin="1054838.5056" RangeMax="2317663.5559" offset="228" />
+ </FieldData>
+ <Piece NumberOfPoints="35" NumberOfCells="21" >
+ <PointData>
+ <DataArray type="Float64" Name="MaterialForces" NumberOfComponents="2" format="appended" RangeMin="134.8875" RangeMax="2733.7625" offset="1468" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="509010.45198" RangeMax="5269775" offset="1896" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0.0225" RangeMax="0.05" offset="2300" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.0005" RangeMax="0.0005" offset="2804" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="1054838.5056" RangeMax="2384275.1716" offset="3764" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="nonequilibrium_stress" NumberOfComponents="4" format="appended" RangeMin="555634.11785" RangeMax="1867755.1655" offset="4596" />
+ <DataArray type="Float64" Name="principal_stress_values" NumberOfComponents="3" format="appended" RangeMin="1054838.5056" RangeMax="2317663.5559" offset="4800" />
+ <DataArray type="Float64" Name="principal_stress_vector_1" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1" offset="5112" />
+ <DataArray type="Float64" Name="principal_stress_vector_2" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1" offset="5272" />
+ <DataArray type="Float64" Name="principal_stress_vector_3" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1" offset="5432" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="45" RangeMax="100.49875621" offset="5484" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5688" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5884" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5992" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAACAAABhAAAAVgAAAA==eF5VjDEOgDAIAP/C7KZTv2IMQUXDADS0Dqbp3+3qeHfJNRCrfAdVccPsg5Ai6C2Q1vaLHicHpHkCI2VIUORWQskwzKM7B/qFh2t2Y6tjsPStfwt7JC0=AQAAAACAAAAoAAAAMAAAAA==eF4z0zPWM9Q1MtdNNzNOsTS3MEs110vJLCqp1DMyMDIwMDMyMDQxNTI0BwDHbAmwAQAAAACAAACgFwAAjgMAAA==eF692M1rE0EYx/Ft2rS2tfZFo1DEUMWiCGojiiKp9aAIXgQPepKiIiJe/AeM3vwDFAUVFbyIL+ChUIRdGo8qqIdIQQ9CrWIpTU2api6xMTPJkwm/zLNPe2j3Ysi389nM7Fui4zjO6MX+ZLpYLMaOxpNOZfty7ubLhQvNnsP0d9/vzrzvCXsNgT3n0vgZwceOvr0bf6U+P/Vp6O75V3seJRur47PM+BAzvtLjNH5K8LGjz/Sq/0Pwfwv+L+i0qfEn3m7T839x/3C17/48MTnwIO1ynbYQ0yvHZ5DGpwUfO/rY0Z8WfOzoYyc/uE/Ggz9fef358REveP2Mb+/Gt3fj5wW/IPh+gD8e7dPz758/ZF0/vufiXKeNelbw7d342NGfF3x7Nz529PPQy/e3jfr6tXVaX8cJnL/HzQ99+/il+3h80J8T/BnGbyy9/uRv1vN/fOog9JS+fvmec0tvFz+y3fi+4Nu78e3d+P8E396Nb+9hj5vfxLcrV7ePtFT7IvTy/a+THY9+QfCxo48dfVw/9NV7TwJ8bv1VT49s0vfX3ucHqr0wHp+6fqaN7bW+rdP5u1r+rOBjRx87+mnokVTX5fG/a/X1PVvqeev+e72mSsfx5f1Hqz6OX66PHX1f8IvMePJxPHX1uutaT7Kj9O+d8P669QsF9rD2bb12U73VOt749m587LWb6iHreOPbu/Glru4vt2v61tN+f/jkFn3/s3U6P2n/PnTayMeOPnbJp+c3+WnBzwj+NOOr4zJ8tk33dcP7oKfctsCec9sDe9gjPy/49m58e195P1TpPvTy9e/r3zeqq7/rtIx3GJ+uX/Jx/HJ97OivEfxuwW9n/J7S66GHxTG1PrfeDFQ7PR+bmE6++t55JKCTnxV87OhznfyC4GOn84d8ew97i8z+n176+vpYtE8/H2zzp+dLgxM0v8iSfezL9XOCvyD4Gcb/UHo/sX5yTJ2njZG9deufYDqtb0zo5Ddbe8pNBPacGwvsxm8QfOy0xYSuthuW+dH3B5+ZP63/XKWrvwvVfb7y7xvldwg+dvSxo98q+BuY8eTb1z+iz697z3ap32+J6PGhuvOH6+X7V1H//4Ct0/Gn8XnBx44+dvQzgo8dfezkq76jZaf+/jD6cxCOn6/3b+vkc732+bga/rzgY0cfO/nUs9DL97duttP3K+p5wf8j+NjRx44+rh/6c4KfYfz/+2E1ig==AQAAAACAAAAwAgAALgEAAA==eF5VkrFLw1AQh4s4dHBxkbo4OTkVdXB5eJObJYOD4JIheVroX+AUQcShSEGLBFFTG6Xg4uT26IHSoWRwliIZHZ0Ep/iS+w3JW47L93Hv945Msywb78b8c6i1njmswzD82x5yE32tOF31XfQue8LpBXxuw571Bw78su+qni/+vHBiLTxOkqSzOmBHl/0LQ1r8Z+EUge/Y8rQUsXNU9k9UA3lawilFnrf8PYt3XPdKfjAwe5j/JZxczFvJ8/3ecKMy31Mz+AvCiZDHze/7DLlWyd9Tm8jTFk4fyNOx+1oeX3O9vB+KTOyJfy6ctuC/2/d3h31OK/NvTYD5r8JpH3mPbZn0T9lHFZ9Mq+iveCTfyQFfs/Xs4BHepZF6b6byH1CKvY2Qz8W9TeT5B4vf/Ac=AQAAAACAAAAwAgAAHAEAAA==eF5jYmBgsCkPcJzw//9/63KBg0AuQ85iX8cLSHwGBoF9GmB1CgdZIPIHJ4D5EPmECC/HDyjqOfZxIKkHyh9sgMqD+CG87o4KSPoZGGTsIPYrHGSCyB8sgMqzAWmPQ86OG1DMZ2hYgKQeKH/wAVSeGazPwTEAxXw9Ox+oe/6A1TkcTIDK/wPy4yOsHCv+IZsvsS8Haj4P2P1WB1dA1bOA7TN2vIFivtI+C6j5rBD5gwFQ9wiA9Ws5LkBxv4wdLHzege3XOqjBCDHvDpgv5RiB5l8fqHuEwOZJHayAyr/7B1L/ywHo7wPOUi5Q9bz7WMD8kINnwOp/HRCFyoPC55BMkKMCivsF9nlA3QPhf9g7Ayl8genjIEw9AACTtwg=AQAAAACAAAAwAgAAaAEAAA==eF5jYICAeTPBYD+Uy6AX03/oq8ac/VZCRlWxJwX3wPhp8scSEy009uhC+TD1e+RaXwfumLf/wtLoPXuLL1jD+A17ivx29z2B82HqlxXYcl1fvGD/pJvGr1g2yu+B8RN4f7BPUfpkDePD1L8L3AE0YeF+DqWlc3YwKu6B8Qt6V/vdTJCH82HqTY1BYPH+mv0z1v3L5LUxg/IPKTFN3yuuAufD1NeKrHN/WLVkv8fVqfwzTB5Zw/gWkxh6pW9/2g3jw9Qf/aoBDIGl+wNWr74ucFDIBsYPYZjOu1pbDs6Hqee5vhjoo+X7H/DskDln9Nsaxj9g0nq72O+GNS+UD1MfAvJO4Ir9p69fZVeVUbWB8d2j+ic8iDCE82Hq54Cjb+V+B5eWvYuWsu2B8RmK5K33V0rumQ3lz9MPOX3/m6oNxL7p+//oBF5NXKRjwwfla/F/TVQJ1reBpYdbfOzaT0OtbOZC+QCHIusBAQAAAACAAABgBAAAvAIAAA==eF5l01tIk2EYB/C1iYqHtKGG81Tm2pZbmprpvq/ZYZSZUyGEmLiTCUGRhwsTywzTykQv0uwwFUXH0GShWUHfY+VAqDnaUNGh4jzU1OmcmlObuXXTTa/v7Y//+xzg2dD5nOrr3iKGOtaGBLIzn0n/nl1YsvQwwR/cjrRty4tt2CjiyYW9pNrbLPCuUvyoKlzFphBXxLt7yStZwJ5Kb7UaTJgW8QvxPce9UhmweLjWysixYJOI22i53wLv0GB4n9lwvWp6T38RH5VJi+wl7Cu1W/+H6yTGES+Nn0vKtlNxO1lL6ttxEuj/Dbuyd/1ZLDxY26Ihe5BgAvGK6n5jIC0IoGHCjxFAgp+I1/CfaixcLtDnBydPYxRA56dy2ikOuomQjBVFrOUrsGXEyyybfJlfEMw0WXMHTBTchHhXJ1wj97DAqcqJG0h3xdH5Cos6Tgo3mJAwm+ljrnHZk195TpaWP2dBvn+yeb+WjNsQ1w5K08Y72DCrbAwRuf7GNhGfkww5h+MSYXEsQ2wdXSfQ/p2X6gv4TB5MRYbdChX6AFp/Tjaacq6AD+03ogyiuz6whbhaHPvqsqcXnln9mv6lngJofY+Qvma11BU/1NWsrexdIdYRz1KNeDtk0bho+KaKlxewZ/8VnfyVA99ZwLz/xlsnIsMC4lfvtZmzsUQwptrKP51dxcid/3uGoU6syDkPWFpNrrLJF99G8qEveOxZEQO4uw11HDfSnv2WtDxTbwiOAi9Mw+Pq1IQDcXn143Xq8gkQlJpEK+95OFqfJp+QDMxzQC985E69Eo3bkXyMb5nnlDISgisfhGf9cmChSL61OYpuNibAxemXboKccUKP5Hd0NIO/xzEIrzFKZhIU2AjimU3hb8V5ESBN0nzgKBcI1NUxrNg8PROeWAQHfVUqAr0fl8aUYkF5LPwFNNe+Lw==AQAAAACAAABgBAAAXgIAAA==eF5jYmBgSIjQOvj2////1TvtDzJAQaOjakL+4c97QewPW9UPvgPKfxW2g8sfs93w/dx29n0w+U9o8nJP7v58cgQh/wJN3k4q+UDxJlawfIOwKlh/fp4NXP60sf2E9pU/9sLk36DJy/xZf1fB/6QtTP4dmvzcplrfLWs+2jIB2QKFSgc/AOVfnrCCy1fKvE97PZfdjhEq/xFNvqfhldEc259g+0HyX9DkW8+9OXM4V2IfC5DtMFce7P5kZUu4vNSN2rrpC46A9YPkv6PJn9wXL+Np8XPvH6C4PVD+C5q8y8PyJfbvIOGXECEFdv/dWnO4fLDQsadzbdn3/QWKxwPlv6LJ14REd068yr7vD1T+D5r8jFPLc19s5IKGv+hBkDnhN0zh8n8WbZJZVSkOl/+JJt/GXqfDLyYD9j9I/iuavM3nX5q3OJX3cTGA/M9/ECj9/6KRCcJ+tp3JtnKvbNmh8n/Q5G9MMvJ2vPEEHH8g+V9o8r4ZbI4nWvnsuMH2s4LD37vXCC5fYmV6NZGDA+o+VnD4I8uX7lsteLBefN8/oHg9UB5kjw+S/DWbnpyDfEr7PoLlvx4AqTv63AAu38BWN+XfYdZ9v6Dyf9HkldYlzTt/lAUafl8PgOhjSPKH3DzFfm0Q2CcN1v/0ADNQjFkUIS9Xe4j7303ufW1QeVD4MSHJ/1iw89Xig5z7PkHl+dH0G6XlaEzVEt8H8teFXzrg+GlY7wCXD1k1d7qbAwfYfSD5b2jyuxg/bpPgYIamPy1w/kMuH9L45yRoybDD5d+jybPtn9Thpy20DwAWB6ztAQAAAACAAACgAgAAhQAAAA==eF5jZGBg2JGidhBIMQgvMAXTyICRCHmPQ8pg8d5fxmTJ35BXBIuzhxqRJX/hlwxYvHG9AVnyH7aKg8W/CeuRJS9QKAQWf3VCiyz5hAgusPi9WnWy5B3m/j8AYl8yUiFLvkH4KVj+2HMFsuRnrNQGi8u/tsAqr8GuCRaf6WpOsjwALgpVHA==AQAAAACAAAD4AQAA1gAAAA==eF579///fyNH24MMQLAjRQ1Mg8BrqDgjmvgHoPia2dZgcY9DynDxz3jE1b5ZgsVvyCvCxb/gEP8NFF8QYAHmX/glgyEOpP6fRxJ/C+RLrTYDi7/fKg4X/wYVZwKyPyCJswPxVFZTsLhAoRBc/DtQ/RSgOAuaOMidfAnGBzmB7IQILrg4MxDzA8UZ0MR/AtV37DI8yApkO8z9fwDZHJA4A5o4yBxmUYODdUD5euGncHEOqPhvNHGgOfXybg7gcJuxUhtu7yeguh3P7MDiGuyaBMUBZd6hgQ==AQAAAACAAAD4AQAAZgAAAA==eF5jYEAGH+wZsALqix/1VVAs6b1gQ6x6XOJ74+QvXXp8HsOcyzuXBDouvk60+cZPVySlp54kWn3Nq8XPtKcdx1B/JXG/bKz0UQzxpyEzbZYuvU20+XUr8udknrlCtHpSxQF8qzqVAQAAAACAAAD4AQAAZgAAAA==eF5jYACBD/YMWAHtxY/6KiiW9F7YQ6x6XOJ74+QvXXp8HsOcyzuXBDouvk60+cZPVySlp54kWn3Nq8XPtKcdx1B/JXG/bKz0UQzxpyEzbZYuvU20+XUr8udknrlCtHpixQGJ6D6VAQAAAACAAAD4AQAAFAAAAA==eF5jYMAGPthjFR4VHzbiAH0sGNw=AQAAAACAAABIAwAAhwAAAA==eF51080JQjEUBeEUYQkpJCtBBFHxB7tJaa80N2dQJu/O8ltcyIG09t9ptL3mWX4Y+94Lp4s8dxbvhadxlefO4r3wNG9y3mXnXfa03eW5s3gvnB5y9rGzjz2Np5x97OxjT/MlZx87+9jT9pazj5197PSRs4+dfexpHuW8q/DlX3D/51+tySeVAQAAAACAAACgAgAAgQAAAA==eF51kskSglAQA99DQRRxwQUVFPH/P9KD6UvKmUtXupLb5PS7QizFpfmFWAU9ci2uzLNbBz1yI27Ms9sGPfJObM2z2wc98lE8mGfXBT3yWTyZZ3cJeuRevJpndwt65Id4N89uCHrkpziaZ/ey3iS+xTn9Pzz/lM1/RP6M3hdg+wUpAQAAAACAAACoAAAAPwAAAA==eF4txRECgDAAAMAgCIIgGARBEAwGgyAIgsEg2P/fEuxObhy6ybMXrw7evPvw6ejk7Mu3H78urv7c/APeCAOdAQAAAACAAAAVAAAACwAAAA==eF7j5MQCAAg0AL4=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_generate_values.py b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_generate_values.py
new file mode 100755
index 0000000000000000000000000000000000000000..db5b66ac85ffbc46d130465e9a32bde7dcca2e91
--- /dev/null
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_generate_values.py
@@ -0,0 +1,45 @@
+#!/usr/bin/env python
+from vtk import *
+
+
+def createStressArray(points):
+ sigma = vtkDoubleArray()
+ sigma.SetNumberOfComponents(4)
+ sigma.SetNumberOfTuples(points.GetNumberOfPoints())
+ for i in range(points.GetNumberOfPoints()):
+ y = points.GetPoint(i)[1]
+ sigma.SetTuple4(i, 2200 * 9.81 * 0.8 * (y - 100),
+ -0.5 * 1e6 + 2200 * 9.81 * (y - 100), 0, 0)
+ #sigma.SetTuple4(i, 4*i, 4*i + 1, 4*i + 2, 4*i + 3) # for debugging
+ sigma.SetName('nonequilibrium_stress')
+ return sigma
+
+
+def createStressArrayForCells(mesh):
+ cell_centers = vtkCellCenters()
+ cell_centers.SetInputData(mesh)
+ cell_centers.Update()
+ return createStressArray(cell_centers.GetOutput())
+
+
+def createStressArrayForNodes(mesh):
+ return createStressArray(mesh.GetPoints())
+
+
+def writeDataToFile(mesh, filename):
+ w = vtkXMLUnstructuredGridWriter()
+ w.SetFileName(filename)
+ w.SetInputData(mesh)
+ w.Write()
+
+def addStressArraysToMesh(mesh):
+ mesh.GetPointData().AddArray(createStressArrayForNodes(mesh))
+ mesh.GetCellData().AddArray(createStressArrayForCells(mesh))
+
+
+r = vtkXMLUnstructuredGridReader()
+r.SetFileName('soil_column.vtu')
+r.Update()
+m = r.GetOutput()
+addStressArraysToMesh(m)
+writeDataToFile(m, 'soil_column_nonequilibrium_sigma.vtu')
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_cell_properties_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_cell_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..11360b516b40c054704b13dedd6dc7a54eccc76b
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_cell_properties_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_cell_properties_val2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_cell_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..65a9cad4a55e0d92f0e9935d6ba6edd0caa3ad9c
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_cell_properties_val2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_msh_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_msh_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..17aa8c68a05004f51ebe7a65dd85927eb643073c
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_msh_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_msh_ele2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_msh_ele2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..a560ad6f88ee45e1d9f3dc6115d81379ab26bc9d
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_msh_ele2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_msh_ele_g2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_msh_ele_g2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..781a37c6d4d15c10e6ea236719a91f05976e83ac
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_msh_ele_g2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_msh_nod2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_msh_nod2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..0b17968d442e4df54351ddea1ade72fbfbe8627e
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_nonequilibrium_sigma_partitioned_msh_nod2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_q8_bottom.vtu b/Tests/Data/Mechanics/InitialStates/soil_column_q8_bottom.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..c68af3d102c880bb315b922f13990afaf485978d
--- /dev/null
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_q8_bottom.vtu
@@ -0,0 +1,33 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="5" NumberOfCells="2">
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="0" RangeMax="124">
+ KAAAAHwAAAAAAAAAAAAAAAAAAAACAAAAAAAAAAEAAAAAAAAAUwAAAAAAAAA=
+ </DataArray>
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="binary" RangeMin="0" RangeMax="1">
+ EAAAAAAAAAAAAAAAAQAAAAAAAAA=
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="10">
+ eAAAAAAAAAAAAB5AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACRAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABRAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAARAAAAAAAAAAAAAAAAAAAAAAA==
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="4">
+ MAAAAAEAAAAAAAAAAwAAAAAAAAAEAAAAAAAAAAMAAAAAAAAAAgAAAAAAAAAAAAAAAAAAAA==
+ </DataArray>
+ <DataArray type="Int64" Name="offsets" format="binary" RangeMin="3" RangeMax="6">
+ EAAAAAMAAAAAAAAABgAAAAAAAAA=
+ </DataArray>
+ <DataArray type="UInt8" Name="types" format="binary" RangeMin="21" RangeMax="21">
+ AgAAABUV
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_q8_left.vtu b/Tests/Data/Mechanics/InitialStates/soil_column_q8_left.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..fe15564751da543a537179d5d0e8d8b7921c2d7c
--- /dev/null
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_q8_left.vtu
@@ -0,0 +1,33 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="41" NumberOfCells="20">
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="0" RangeMax="82">
+ SAEAAAAAAAAAAAAAQAAAAAAAAAADAAAAAAAAAEEAAAAAAAAAQwAAAAAAAABEAAAAAAAAAD8AAAAAAAAADwAAAAAAAABCAAAAAAAAABIAAAAAAAAARgAAAAAAAABFAAAAAAAAABUAAAAAAAAARwAAAAAAAAAYAAAAAAAAAEgAAAAAAAAADAAAAAAAAAA8AAAAAAAAAAkAAAAAAAAAOQAAAAAAAABQAAAAAAAAAFIAAAAAAAAAJAAAAAAAAAAGAAAAAAAAADYAAAAAAAAATQAAAAAAAAAzAAAAAAAAAFEAAAAAAAAAIQAAAAAAAABKAAAAAAAAADAAAAAAAAAATwAAAAAAAAAtAAAAAAAAAB4AAAAAAAAATgAAAAAAAAAqAAAAAAAAACcAAAAAAAAATAAAAAAAAAAbAAAAAAAAAEsAAAAAAAAASQAAAAAAAAA=
+ </DataArray>
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="binary" RangeMin="0" RangeMax="38">
+ oAAAACYAAAAAAAAAJAAAAAAAAAAiAAAAAAAAACAAAAAAAAAAHgAAAAAAAAAcAAAAAAAAABoAAAAAAAAAGAAAAAAAAAAWAAAAAAAAABQAAAAAAAAAEgAAAAAAAAAQAAAAAAAAAA4AAAAAAAAADAAAAAAAAAAKAAAAAAAAAAgAAAAAAAAABgAAAAAAAAAEAAAAAAAAAAIAAAAAAAAAAAAAAAAAAAA=
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="100">
+ 2AMAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAHkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAFEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAKUAAAAAAAAAAAAAAAAAAAAAAAAAAAACANkAAAAAAAAAAAAAAAAAAAAAAAAAAAACAO0AAAAAAAAAAAAAAAAAAAAAAAAAAAAAABEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAOUAAAAAAAAAAAAAAAAAAAAAAAAAAAACAMUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPkAAAAAAAAAAAAAAAAAAAAAAAAAAAADAQkAAAAAAAAAAAAAAAAAAAAAAAAAAAABAQEAAAAAAAAAAAAAAAAAAAAAAAAAAAACAQUAAAAAAAAAAAAAAAAAAAAAAAAAAAABARUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAREAAAAAAAAAAAAAAAAAAAAAAAAAAAADAR0AAAAAAAAAAAAAAAAAAAAAAAAAAAAAANEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAWUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAALkAAAAAAAAAAAAAAAAAAAAAAAAAAAADAV0AAAAAAAAAAAAAAAAAAAAAAAAAAAADgVUAAAAAAAAAAAAAAAAAAAAAAAAAAAABgWEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAATkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAJEAAAAAAAAAAAAAAAAAAAAAAAAAAAACAVkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAgUkAAAAAAAAAAAAAAAAAAAAAAAAAAAABAVUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAgV0AAAAAAAAAAAAAAAAAAAAAAAAAAAACAS0AAAAAAAAAAAAAAAAAAAAAAAAAAAADATEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAVEAAAAAAAAAAAAAAAAAAAAAAAAAAAACgVEAAAAAAAAAAAAAAAAAAAAAAAAAAAADAUkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAASUAAAAAAAAAAAAAAAAAAAAAAAAAAAABgU0AAAAAAAAAAAAAAAAAAAAAAAAAAAACAUUAAAAAAAAAAAAAAAAAAAAAAAAAAAABAUEAAAAAAAAAAAAAAAAAAAAAAAAAAAADgUEAAAAAAAAAAAAAAAAAAAAAAAAAAAACARkAAAAAAAAAAAAAAAAAAAAAAAAAAAABAT0AAAAAAAAAAAAAAAAAAAAAAAAAAAABASkAAAAAAAAAAAA==
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="40">
+ 4AEAABEAAAAAAAAAEwAAAAAAAAAVAAAAAAAAABMAAAAAAAAAGAAAAAAAAAAbAAAAAAAAABgAAAAAAAAAGgAAAAAAAAAUAAAAAAAAABoAAAAAAAAAHgAAAAAAAAAfAAAAAAAAAB4AAAAAAAAAIAAAAAAAAAAiAAAAAAAAACAAAAAAAAAAIwAAAAAAAAAZAAAAAAAAACMAAAAAAAAAJAAAAAAAAAAlAAAAAAAAACQAAAAAAAAAFgAAAAAAAAAnAAAAAAAAABYAAAAAAAAAHAAAAAAAAAAdAAAAAAAAABwAAAAAAAAAIQAAAAAAAAAoAAAAAAAAACEAAAAAAAAAJgAAAAAAAAAPAAAAAAAAACYAAAAAAAAADgAAAAAAAAANAAAAAAAAAA4AAAAAAAAADAAAAAAAAAAKAAAAAAAAAAwAAAAAAAAACQAAAAAAAAALAAAAAAAAAAkAAAAAAAAABwAAAAAAAAAFAAAAAAAAAAcAAAAAAAAAEAAAAAAAAAAEAAAAAAAAABAAAAAAAAAAEgAAAAAAAAAIAAAAAAAAABIAAAAAAAAAFwAAAAAAAAADAAAAAAAAABcAAAAAAAAAAgAAAAAAAAABAAAAAAAAAAIAAAAAAAAAAAAAAAAAAAAGAAAAAAAAAA==
+ </DataArray>
+ <DataArray type="Int64" Name="offsets" format="binary" RangeMin="3" RangeMax="60">
+ oAAAAAMAAAAAAAAABgAAAAAAAAAJAAAAAAAAAAwAAAAAAAAADwAAAAAAAAASAAAAAAAAABUAAAAAAAAAGAAAAAAAAAAbAAAAAAAAAB4AAAAAAAAAIQAAAAAAAAAkAAAAAAAAACcAAAAAAAAAKgAAAAAAAAAtAAAAAAAAADAAAAAAAAAAMwAAAAAAAAA2AAAAAAAAADkAAAAAAAAAPAAAAAAAAAA=
+ </DataArray>
+ <DataArray type="UInt8" Name="types" format="binary" RangeMin="21" RangeMax="21">
+ FAAAABUVFRUVFRUVFRUVFRUVFRUVFRUV
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_q8_right.vtu b/Tests/Data/Mechanics/InitialStates/soil_column_q8_right.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..6b141282313d04049854cb1910ed9e9bad3b2d18
--- /dev/null
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_q8_right.vtu
@@ -0,0 +1,33 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="41" NumberOfCells="20">
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="2" RangeMax="164">
+ SAEAAKQAAAAAAAAAowAAAAAAAAChAAAAAAAAAKAAAAAAAAAAngAAAAAAAAAsAAAAAAAAAJ0AAAAAAAAAKQAAAAAAAAAmAAAAAAAAAJsAAAAAAAAAkQAAAAAAAAA4AAAAAAAAAJwAAAAAAAAACAAAAAAAAACSAAAAAAAAADsAAAAAAAAAnwAAAAAAAAALAAAAAAAAAJUAAAAAAAAALwAAAAAAAACTAAAAAAAAACAAAAAAAAAAPgAAAAAAAACiAAAAAAAAAA4AAAAAAAAAmAAAAAAAAACUAAAAAAAAABEAAAAAAAAAIwAAAAAAAAAUAAAAAAAAADIAAAAAAAAAlgAAAAAAAAAXAAAAAAAAAJcAAAAAAAAAGgAAAAAAAAAdAAAAAAAAADUAAAAAAAAABQAAAAAAAACZAAAAAAAAAAIAAAAAAAAAmgAAAAAAAAA=
+ </DataArray>
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="binary" RangeMin="1" RangeMax="39">
+ oAAAAAEAAAAAAAAAAwAAAAAAAAAFAAAAAAAAAAcAAAAAAAAACQAAAAAAAAALAAAAAAAAAA0AAAAAAAAADwAAAAAAAAARAAAAAAAAABMAAAAAAAAAFQAAAAAAAAAXAAAAAAAAABkAAAAAAAAAGwAAAAAAAAAdAAAAAAAAAB8AAAAAAAAAIQAAAAAAAAAjAAAAAAAAACUAAAAAAAAAJwAAAAAAAAA=
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="10" RangeMax="100.49875621">
+ 2AMAAAAAAAAAACRAAAAAAABgWEAAAAAAAAAAAAAAAAAAACRAAAAAAAAgV0AAAAAAAAAAAAAAAAAAACRAAAAAAACgVEAAAAAAAAAAAAAAAAAAACRAAAAAAABgU0AAAAAAAAAAAAAAAAAAACRAAAAAAADgUEAAAAAAAAAAAAAAAAAAACRAAAAAAACAUUAAAAAAAAAAAAAAAAAAACRAAAAAAABAT0AAAAAAAAAAAAAAAAAAACRAAAAAAABAUEAAAAAAAAAAAAAAAAAAACRAAAAAAAAATkAAAAAAAAAAAAAAAAAAACRAAAAAAABASkAAAAAAAAAAAAAAAAAAACRAAAAAAAAABEAAAAAAAAAAAAAAAAAAACRAAAAAAACAVkAAAAAAAAAAAAAAAAAAACRAAAAAAADATEAAAAAAAAAAAAAAAAAAACRAAAAAAAAAJEAAAAAAAAAAAAAAAAAAACRAAAAAAAAAHkAAAAAAAAAAAAAAAAAAACRAAAAAAADAV0AAAAAAAAAAAAAAAAAAACRAAAAAAAAgUkAAAAAAAAAAAAAAAAAAACRAAAAAAAAALkAAAAAAAAAAAAAAAAAAACRAAAAAAACANkAAAAAAAAAAAAAAAAAAACRAAAAAAADAUkAAAAAAAAAAAAAAAAAAACRAAAAAAAAAKUAAAAAAAAAAAAAAAAAAACRAAAAAAAAASUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAWUAAAAAAAAAAAAAAAAAAACRAAAAAAADgVUAAAAAAAAAAAAAAAAAAACRAAAAAAAAANEAAAAAAAAAAAAAAAAAAACRAAAAAAADAQkAAAAAAAAAAAAAAAAAAACRAAAAAAACAMUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAOUAAAAAAAAAAAAAAAAAAACRAAAAAAACAS0AAAAAAAAAAAAAAAAAAACRAAAAAAAAAPkAAAAAAAAAAAAAAAAAAACRAAAAAAAAAVEAAAAAAAAAAAAAAAAAAACRAAAAAAACAO0AAAAAAAAAAAAAAAAAAACRAAAAAAACAQUAAAAAAAAAAAAAAAAAAACRAAAAAAABAQEAAAAAAAAAAAAAAAAAAACRAAAAAAAAAREAAAAAAAAAAAAAAAAAAACRAAAAAAACARkAAAAAAAAAAAAAAAAAAACRAAAAAAABAVUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAFEAAAAAAAAAAAAAAAAAAACRAAAAAAABARUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACRAAAAAAADAR0AAAAAAAAAAAA==
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="40">
+ 4AEAACcAAAAAAAAAJQAAAAAAAAAKAAAAAAAAACUAAAAAAAAADQAAAAAAAAAOAAAAAAAAAA0AAAAAAAAAEQAAAAAAAAAUAAAAAAAAABEAAAAAAAAAGAAAAAAAAAAaAAAAAAAAABgAAAAAAAAAGwAAAAAAAAASAAAAAAAAABsAAAAAAAAAHQAAAAAAAAAfAAAAAAAAAB0AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACAAAAAAAAAAIgAAAAAAAAAZAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAmAAAAAAAAACMAAAAAAAAAFQAAAAAAAAAoAAAAAAAAABUAAAAAAAAAHAAAAAAAAAAJAAAAAAAAABwAAAAAAAAACAAAAAAAAAAMAAAAAAAAAAgAAAAAAAAABwAAAAAAAAAGAAAAAAAAAAcAAAAAAAAABQAAAAAAAAAEAAAAAAAAAAUAAAAAAAAAEwAAAAAAAAAQAAAAAAAAABMAAAAAAAAAHgAAAAAAAAADAAAAAAAAAB4AAAAAAAAAJAAAAAAAAAACAAAAAAAAACQAAAAAAAAACwAAAAAAAAAXAAAAAAAAAAsAAAAAAAAADwAAAAAAAAABAAAAAAAAAA8AAAAAAAAAFgAAAAAAAAAAAAAAAAAAAA==
+ </DataArray>
+ <DataArray type="Int64" Name="offsets" format="binary" RangeMin="3" RangeMax="60">
+ oAAAAAMAAAAAAAAABgAAAAAAAAAJAAAAAAAAAAwAAAAAAAAADwAAAAAAAAASAAAAAAAAABUAAAAAAAAAGAAAAAAAAAAbAAAAAAAAAB4AAAAAAAAAIQAAAAAAAAAkAAAAAAAAACcAAAAAAAAAKgAAAAAAAAAtAAAAAAAAADAAAAAAAAAAMwAAAAAAAAA2AAAAAAAAADkAAAAAAAAAPAAAAAAAAAA=
+ </DataArray>
+ <DataArray type="UInt8" Name="types" format="binary" RangeMin="21" RangeMax="21">
+ FAAAABUVFRUVFRUVFRUVFRUVFRUVFRUV
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_q8_top.vtu b/Tests/Data/Mechanics/InitialStates/soil_column_q8_top.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..7fe3fbac6429df6a1ae2786718fea3942ba21f19
--- /dev/null
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_q8_top.vtu
@@ -0,0 +1,33 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="5" NumberOfCells="2">
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="60" RangeMax="144">
+ KAAAAGcAAAAAAAAAPgAAAAAAAACQAAAAAAAAAD0AAAAAAAAAPAAAAAAAAAA=
+ </DataArray>
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="binary" RangeMin="38" RangeMax="39">
+ EAAAACcAAAAAAAAAJgAAAAAAAAA=
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="100" RangeMax="100.49875621">
+ eAAAAAAAAAAAAARAAAAAAAAAWUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAWUAAAAAAAAAAAAAAAAAAAB5AAAAAAAAAWUAAAAAAAAAAAAAAAAAAABRAAAAAAAAAWUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAWUAAAAAAAAAAAA==
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="4">
+ MAAAAAEAAAAAAAAAAwAAAAAAAAACAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAAAAAAAAAAAA==
+ </DataArray>
+ <DataArray type="Int64" Name="offsets" format="binary" RangeMin="3" RangeMax="6">
+ EAAAAAMAAAAAAAAABgAAAAAAAAA=
+ </DataArray>
+ <DataArray type="UInt8" Name="types" format="binary" RangeMin="21" RangeMax="21">
+ AgAAABUV
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_right.vtu b/Tests/Data/Mechanics/InitialStates/soil_column_right.vtu
index 6b141282313d04049854cb1910ed9e9bad3b2d18..5fbb21163362dfdb9a3aa2819be2bc91a102864f 100644
--- a/Tests/Data/Mechanics/InitialStates/soil_column_right.vtu
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_right.vtu
@@ -1,33 +1,24 @@
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
<UnstructuredGrid>
- <Piece NumberOfPoints="41" NumberOfCells="20">
+ <Piece NumberOfPoints="20" NumberOfCells="18" >
<PointData>
- <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="2" RangeMax="164">
- SAEAAKQAAAAAAAAAowAAAAAAAAChAAAAAAAAAKAAAAAAAAAAngAAAAAAAAAsAAAAAAAAAJ0AAAAAAAAAKQAAAAAAAAAmAAAAAAAAAJsAAAAAAAAAkQAAAAAAAAA4AAAAAAAAAJwAAAAAAAAACAAAAAAAAACSAAAAAAAAADsAAAAAAAAAnwAAAAAAAAALAAAAAAAAAJUAAAAAAAAALwAAAAAAAACTAAAAAAAAACAAAAAAAAAAPgAAAAAAAACiAAAAAAAAAA4AAAAAAAAAmAAAAAAAAACUAAAAAAAAABEAAAAAAAAAIwAAAAAAAAAUAAAAAAAAADIAAAAAAAAAlgAAAAAAAAAXAAAAAAAAAJcAAAAAAAAAGgAAAAAAAAAdAAAAAAAAADUAAAAAAAAABQAAAAAAAACZAAAAAAAAAAIAAAAAAAAAmgAAAAAAAAA=
- </DataArray>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="2" RangeMax="62" offset="0" />
</PointData>
<CellData>
- <DataArray type="UInt64" Name="bulk_element_ids" format="binary" RangeMin="1" RangeMax="39">
- oAAAAAEAAAAAAAAAAwAAAAAAAAAFAAAAAAAAAAcAAAAAAAAACQAAAAAAAAALAAAAAAAAAA0AAAAAAAAADwAAAAAAAAARAAAAAAAAABMAAAAAAAAAFQAAAAAAAAAXAAAAAAAAABkAAAAAAAAAGwAAAAAAAAAdAAAAAAAAAB8AAAAAAAAAIQAAAAAAAAAjAAAAAAAAACUAAAAAAAAAJwAAAAAAAAA=
- </DataArray>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="1" RangeMax="39" offset="220" />
</CellData>
<Points>
- <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="10" RangeMax="100.49875621">
- 2AMAAAAAAAAAACRAAAAAAABgWEAAAAAAAAAAAAAAAAAAACRAAAAAAAAgV0AAAAAAAAAAAAAAAAAAACRAAAAAAACgVEAAAAAAAAAAAAAAAAAAACRAAAAAAABgU0AAAAAAAAAAAAAAAAAAACRAAAAAAADgUEAAAAAAAAAAAAAAAAAAACRAAAAAAACAUUAAAAAAAAAAAAAAAAAAACRAAAAAAABAT0AAAAAAAAAAAAAAAAAAACRAAAAAAABAUEAAAAAAAAAAAAAAAAAAACRAAAAAAAAATkAAAAAAAAAAAAAAAAAAACRAAAAAAABASkAAAAAAAAAAAAAAAAAAACRAAAAAAAAABEAAAAAAAAAAAAAAAAAAACRAAAAAAACAVkAAAAAAAAAAAAAAAAAAACRAAAAAAADATEAAAAAAAAAAAAAAAAAAACRAAAAAAAAAJEAAAAAAAAAAAAAAAAAAACRAAAAAAAAAHkAAAAAAAAAAAAAAAAAAACRAAAAAAADAV0AAAAAAAAAAAAAAAAAAACRAAAAAAAAgUkAAAAAAAAAAAAAAAAAAACRAAAAAAAAALkAAAAAAAAAAAAAAAAAAACRAAAAAAACANkAAAAAAAAAAAAAAAAAAACRAAAAAAADAUkAAAAAAAAAAAAAAAAAAACRAAAAAAAAAKUAAAAAAAAAAAAAAAAAAACRAAAAAAAAASUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAWUAAAAAAAAAAAAAAAAAAACRAAAAAAADgVUAAAAAAAAAAAAAAAAAAACRAAAAAAAAANEAAAAAAAAAAAAAAAAAAACRAAAAAAADAQkAAAAAAAAAAAAAAAAAAACRAAAAAAACAMUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAOUAAAAAAAAAAAAAAAAAAACRAAAAAAACAS0AAAAAAAAAAAAAAAAAAACRAAAAAAAAAPkAAAAAAAAAAAAAAAAAAACRAAAAAAAAAVEAAAAAAAAAAAAAAAAAAACRAAAAAAACAO0AAAAAAAAAAAAAAAAAAACRAAAAAAACAQUAAAAAAAAAAAAAAAAAAACRAAAAAAABAQEAAAAAAAAAAAAAAAAAAACRAAAAAAAAAREAAAAAAAAAAAAAAAAAAACRAAAAAAACARkAAAAAAAAAAAAAAAAAAACRAAAAAAABAVUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAFEAAAAAAAAAAAAAAAAAAACRAAAAAAABARUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACRAAAAAAADAR0AAAAAAAAAAAA==
- </DataArray>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="10" RangeMax="100.49875621" offset="420" />
</Points>
<Cells>
- <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="40">
- 4AEAACcAAAAAAAAAJQAAAAAAAAAKAAAAAAAAACUAAAAAAAAADQAAAAAAAAAOAAAAAAAAAA0AAAAAAAAAEQAAAAAAAAAUAAAAAAAAABEAAAAAAAAAGAAAAAAAAAAaAAAAAAAAABgAAAAAAAAAGwAAAAAAAAASAAAAAAAAABsAAAAAAAAAHQAAAAAAAAAfAAAAAAAAAB0AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACAAAAAAAAAAIgAAAAAAAAAZAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAmAAAAAAAAACMAAAAAAAAAFQAAAAAAAAAoAAAAAAAAABUAAAAAAAAAHAAAAAAAAAAJAAAAAAAAABwAAAAAAAAACAAAAAAAAAAMAAAAAAAAAAgAAAAAAAAABwAAAAAAAAAGAAAAAAAAAAcAAAAAAAAABQAAAAAAAAAEAAAAAAAAAAUAAAAAAAAAEwAAAAAAAAAQAAAAAAAAABMAAAAAAAAAHgAAAAAAAAADAAAAAAAAAB4AAAAAAAAAJAAAAAAAAAACAAAAAAAAACQAAAAAAAAACwAAAAAAAAAXAAAAAAAAAAsAAAAAAAAADwAAAAAAAAABAAAAAAAAAA8AAAAAAAAAFgAAAAAAAAAAAAAAAAAAAA==
- </DataArray>
- <DataArray type="Int64" Name="offsets" format="binary" RangeMin="3" RangeMax="60">
- oAAAAAMAAAAAAAAABgAAAAAAAAAJAAAAAAAAAAwAAAAAAAAADwAAAAAAAAASAAAAAAAAABUAAAAAAAAAGAAAAAAAAAAbAAAAAAAAAB4AAAAAAAAAIQAAAAAAAAAkAAAAAAAAACcAAAAAAAAAKgAAAAAAAAAtAAAAAAAAADAAAAAAAAAAMwAAAAAAAAA2AAAAAAAAADkAAAAAAAAAPAAAAAAAAAA=
- </DataArray>
- <DataArray type="UInt8" Name="types" format="binary" RangeMin="21" RangeMax="21">
- FAAAABUVFRUVFRUVFRUVFRUVFRUVFRUV
- </DataArray>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1068" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="1460" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1660" />
</Cells>
</Piece>
</UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _oAAAAD4AAAAAAAAAAgAAAAAAAAAFAAAAAAAAABQAAAAAAAAACAAAAAAAAAALAAAAAAAAAA4AAAAAAAAAEQAAAAAAAAAXAAAAAAAAACYAAAAAAAAAKQAAAAAAAAAsAAAAAAAAABoAAAAAAAAANQAAAAAAAAAvAAAAAAAAAB0AAAAAAAAAOAAAAAAAAAAyAAAAAAAAACAAAAAAAAAAOwAAAAAAAAA=kAAAAAEAAAAAAAAAAwAAAAAAAAAFAAAAAAAAAAcAAAAAAAAACQAAAAAAAAALAAAAAAAAAA0AAAAAAAAADwAAAAAAAAARAAAAAAAAABMAAAAAAAAAGQAAAAAAAAAbAAAAAAAAAB0AAAAAAAAAHwAAAAAAAAAhAAAAAAAAACMAAAAAAAAAJQAAAAAAAAAnAAAAAAAAAA==4AEAAAAAAAAAACRAAAAAAAAAWUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACRAAAAAAAAAFEAAAAAAAAAAAAAAAAAAACRAAAAAAAAAPkAAAAAAAAAAAAAAAAAAACRAAAAAAAAAJEAAAAAAAAAAAAAAAAAAACRAAAAAAAAALkAAAAAAAAAAAAAAAAAAACRAAAAAAAAANEAAAAAAAAAAAAAAAAAAACRAAAAAAAAAOUAAAAAAAAAAAAAAAAAAACRAAAAAAACAQUAAAAAAAAAAAAAAAAAAACRAAAAAAAAATkAAAAAAAAAAAAAAAAAAACRAAAAAAABAUEAAAAAAAAAAAAAAAAAAACRAAAAAAACAUUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAREAAAAAAAAAAAAAAAAAAACRAAAAAAABAVUAAAAAAAAAAAAAAAAAAACRAAAAAAADAUkAAAAAAAAAAAAAAAAAAACRAAAAAAACARkAAAAAAAAAAAAAAAAAAACRAAAAAAACAVkAAAAAAAAAAAAAAAAAAACRAAAAAAAAAVEAAAAAAAAAAAAAAAAAAACRAAAAAAAAASUAAAAAAAAAAAAAAAAAAACRAAAAAAADAV0AAAAAAAAAAAA==IAEAAAEAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABgAAAAAAAAAGAAAAAAAAAAcAAAAAAAAABwAAAAAAAAADAAAAAAAAAAMAAAAAAAAACAAAAAAAAAAIAAAAAAAAAAwAAAAAAAAADAAAAAAAAAAPAAAAAAAAAA8AAAAAAAAAEgAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACgAAAAAAAAALAAAAAAAAAAsAAAAAAAAADgAAAAAAAAAOAAAAAAAAABEAAAAAAAAAEQAAAAAAAAANAAAAAAAAAA0AAAAAAAAAEAAAAAAAAAAQAAAAAAAAABMAAAAAAAAAEwAAAAAAAAAAAAAAAAAAAA==kAAAAAIAAAAAAAAABAAAAAAAAAAGAAAAAAAAAAgAAAAAAAAACgAAAAAAAAAMAAAAAAAAAA4AAAAAAAAAEAAAAAAAAAASAAAAAAAAABQAAAAAAAAAFgAAAAAAAAAYAAAAAAAAABoAAAAAAAAAHAAAAAAAAAAeAAAAAAAAACAAAAAAAAAAIgAAAAAAAAAkAAAAAAAAAA==EgAAAAMDAwMDAwMDAwMDAwMDAwMDAw==
+ </AppendedData>
</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_cell_properties_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_cell_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..a1edd9e70a232bad77bbe97dacfa5b9d3cf13646
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_cell_properties_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_cell_properties_val2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_cell_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..861e0326983d6f99351c339e5b9cf8a3f161d640
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_cell_properties_val2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_msh_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_msh_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..9698a3750c08a2be26d11503adcba4cb877b9a33
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_msh_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_msh_ele2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_msh_ele2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..342990120cc48e5744c3b68bac0f2f6be13a28a3
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_msh_ele2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_msh_ele_g2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_msh_ele_g2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_msh_nod2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_msh_nod2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..d774cc1c4477f3fb9d1ed523037bf5af1634d03a
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_msh_nod2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_node_properties_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_node_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..2001f1550ac676d8e1dc8c21797c25687ecb7d5e
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_node_properties_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_node_properties_val2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_node_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..f5e5c07824870f4d303eebd7f1027181ff485cde
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_right_partitioned_node_properties_val2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_top.vtu b/Tests/Data/Mechanics/InitialStates/soil_column_top.vtu
index 7fe3fbac6429df6a1ae2786718fea3942ba21f19..1904dbe33d73a723126deb88612e13c79197a33b 100644
--- a/Tests/Data/Mechanics/InitialStates/soil_column_top.vtu
+++ b/Tests/Data/Mechanics/InitialStates/soil_column_top.vtu
@@ -1,33 +1,24 @@
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
<UnstructuredGrid>
- <Piece NumberOfPoints="5" NumberOfCells="2">
+ <Piece NumberOfPoints="3" NumberOfCells="2" >
<PointData>
- <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="60" RangeMax="144">
- KAAAAGcAAAAAAAAAPgAAAAAAAACQAAAAAAAAAD0AAAAAAAAAPAAAAAAAAAA=
- </DataArray>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="60" RangeMax="62" offset="0" />
</PointData>
<CellData>
- <DataArray type="UInt64" Name="bulk_element_ids" format="binary" RangeMin="38" RangeMax="39">
- EAAAACcAAAAAAAAAJgAAAAAAAAA=
- </DataArray>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="38" RangeMax="39" offset="40" />
</CellData>
<Points>
- <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="100" RangeMax="100.49875621">
- eAAAAAAAAAAAAARAAAAAAAAAWUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAWUAAAAAAAAAAAAAAAAAAAB5AAAAAAAAAWUAAAAAAAAAAAAAAAAAAABRAAAAAAAAAWUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAWUAAAAAAAAAAAA==
- </DataArray>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="100" RangeMax="100.49875621" offset="68" />
</Points>
<Cells>
- <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="4">
- MAAAAAEAAAAAAAAAAwAAAAAAAAACAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAAAAAAAAAAAA==
- </DataArray>
- <DataArray type="Int64" Name="offsets" format="binary" RangeMin="3" RangeMax="6">
- EAAAAAMAAAAAAAAABgAAAAAAAAA=
- </DataArray>
- <DataArray type="UInt8" Name="types" format="binary" RangeMin="21" RangeMax="21">
- AgAAABUV
- </DataArray>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="172" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="220" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="248" />
</Cells>
</Piece>
</UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _GAAAADwAAAAAAAAAPgAAAAAAAAA9AAAAAAAAAA==EAAAACcAAAAAAAAAJgAAAAAAAAA=SAAAAAAAAAAAAAAAAAAAAAAAWUAAAAAAAAAAAAAAAAAAACRAAAAAAAAAWUAAAAAAAAAAAAAAAAAAABRAAAAAAAAAWUAAAAAAAAAAAA==IAAAAAEAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAAAAAAAAAAAEAAAAAIAAAAAAAAABAAAAAAAAAA=AgAAAAMD
+ </AppendedData>
</VTKFile>
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_cell_properties_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_cell_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..90902ab6bf1f45687b2a6d1fb9d99a0c8b2edb53
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_cell_properties_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_cell_properties_val2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_cell_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..cef10a7fcc98e8ee3b3a21bf54c522d830ec2d9a
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_cell_properties_val2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_msh_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_msh_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..9e80b79f4690a96b549f91b7af5a28e048625e7d
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_msh_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_msh_ele2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_msh_ele2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..7d0392c3f5170652bdda7bdf7255b1f8af19a804
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_msh_ele2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_msh_ele_g2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_msh_ele_g2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_msh_nod2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_msh_nod2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..0a2f718d36c9b3275aa63db485a2a7f9f03b4248
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_msh_nod2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_node_properties_cfg2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_node_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..9baa18973467d450c731266c5011b0254b36381f
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_node_properties_cfg2.bin differ
diff --git a/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_node_properties_val2.bin b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_node_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..76fe9e92ae7a4216a4d31a2e4f94d295a597bb28
Binary files /dev/null and b/Tests/Data/Mechanics/InitialStates/soil_column_top_partitioned_node_properties_val2.bin differ