From 1c9a47d4836a13950dbaaa6c173d6ce2e4484c44 Mon Sep 17 00:00:00 2001
From: Christoph Lehmann <christoph.lehmann@ufz.de>
Date: Tue, 9 Aug 2022 08:21:58 +0200
Subject: [PATCH] [T] Use function mimicking MaterialIDs to test computed
 coordinates

---
 ProcessLib/ThermoHydroMechanics/Tests.cmake    | 18 ++++++++++++++++++
 .../square_1e2_function.xml                    | 12 ++++++++++++
 2 files changed, 30 insertions(+)
 create mode 100644 Tests/Data/ThermoHydroMechanics/Linear/Beam_sealed_bimaterial/square_1e2_function.xml

diff --git a/ProcessLib/ThermoHydroMechanics/Tests.cmake b/ProcessLib/ThermoHydroMechanics/Tests.cmake
index 44dc1602d3e..a92655671d5 100644
--- a/ProcessLib/ThermoHydroMechanics/Tests.cmake
+++ b/ProcessLib/ThermoHydroMechanics/Tests.cmake
@@ -39,6 +39,24 @@ AddTest(
     expected_square_1e2_ts_10_t_100.000000.vtu square_1e2_ts_10_t_100.000000.vtu sigma sigma 1e-8 1e-8
 )
 
+# Same as above, but with function instead of group based parameter for Young's modulus
+AddTest(
+NAME ThermoHydroMechanics_square_1e2_sealed_bimaterial_function
+PATH ThermoHydroMechanics/Linear/Beam_sealed_bimaterial
+RUNTIME 5
+EXECUTABLE ogs
+EXECUTABLE_ARGS square_1e2_function.xml
+WRAPPER time
+TESTER vtkdiff
+REQUIREMENTS NOT OGS_USE_MPI
+DIFF_DATA
+expected_square_1e2_ts_10_t_100.000000.vtu square_1e2_function_ts_10_t_100.000000.vtu displacement displacement 1e-8 1e-8
+expected_square_1e2_ts_10_t_100.000000.vtu square_1e2_function_ts_10_t_100.000000.vtu pressure pressure 1e-8 1e-8
+expected_square_1e2_ts_10_t_100.000000.vtu square_1e2_function_ts_10_t_100.000000.vtu temperature temperature 1e-8 1e-8
+expected_square_1e2_ts_10_t_100.000000.vtu square_1e2_function_ts_10_t_100.000000.vtu epsilon epsilon 1e-8 1e-8
+expected_square_1e2_ts_10_t_100.000000.vtu square_1e2_function_ts_10_t_100.000000.vtu sigma sigma 1e-8 1e-8
+)
+
 # ThermoHydroMechanics; Small deformation, linear poroelastic, unsealed, bimaterial
 AddTest(
     NAME ThermoHydroMechanics_square_1e2_unsealed_bimaterial
diff --git a/Tests/Data/ThermoHydroMechanics/Linear/Beam_sealed_bimaterial/square_1e2_function.xml b/Tests/Data/ThermoHydroMechanics/Linear/Beam_sealed_bimaterial/square_1e2_function.xml
new file mode 100644
index 00000000000..be328afd51d
--- /dev/null
+++ b/Tests/Data/ThermoHydroMechanics/Linear/Beam_sealed_bimaterial/square_1e2_function.xml
@@ -0,0 +1,12 @@
+<?xml version='1.0' encoding='ISO-8859-1'?>
+<OpenGeoSysProjectDiff base_file="square_1e2.prj">
+    <remove sel="/*/parameters/parameter/name[text()=&quot;E&quot;]/.." />
+    <add sel="/*/parameters">
+        <parameter>
+            <name>E</name>
+            <type>Function</type>
+            <expression>if (x &lt; 2) 19; else if (x &lt; 2.1) 0.95; else 9.5;</expression>
+        </parameter>
+    </add>
+    <replace sel="/*/time_loop/output/prefix/text()">square_1e2_function</replace>
+</OpenGeoSysProjectDiff>
-- 
GitLab