From 623e3ecb5cd7c8fc76be6f1d2051b8c6cf1fc39c Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Tue, 28 May 2024 12:59:04 +0200
Subject: [PATCH] [T/TRM] Increase accuracy of bedding vector

The original input is not close enough to unity.
---
 .../TaskCDECOVALEX2023/Decovalex-0-TRF.prj    |  6 ++++--
 .../Decovalex-0-TRF.prj                       |  6 ++++--
 .../BodyForce/square.prj                      |  6 ++++--
 .../TaskCDECOVALEX2023/Decovalex-0.prj        | 20 ++++++++-----------
 4 files changed, 20 insertions(+), 18 deletions(-)

diff --git a/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/Decovalex-0-TRF.prj b/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/Decovalex-0-TRF.prj
index 8865cb39463..140edaf9951 100644
--- a/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/Decovalex-0-TRF.prj
+++ b/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/Decovalex-0-TRF.prj
@@ -559,12 +559,14 @@
         <parameter>
             <name>b0</name>
             <type>Constant</type>
-            <values> 0.829037572555042 0.559192903470747</values>
+            <!-- 34 degrees bedding angle -->
+            <values> 0.829037572555041692 0.55919290347074683016</values>
         </parameter>
         <parameter>
             <name>b1</name>
             <type>Constant</type>
-            <values>-0.559192903470747 0.829037572555042</values>
+            <!-- 34 degrees bedding angle -->
+            <values>-0.55919290347074683016 0.829037572555041692</values>
         </parameter>
         <parameter>
             <name>ThermalConductivityClay</name>
diff --git a/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/WithPicardNonLinearSolverAndPETSc/Decovalex-0-TRF.prj b/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/WithPicardNonLinearSolverAndPETSc/Decovalex-0-TRF.prj
index 08853826e3d..92f99f0549d 100644
--- a/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/WithPicardNonLinearSolverAndPETSc/Decovalex-0-TRF.prj
+++ b/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/WithPicardNonLinearSolverAndPETSc/Decovalex-0-TRF.prj
@@ -559,12 +559,14 @@
         <parameter>
             <name>b0</name>
             <type>Constant</type>
-            <values> 0.829037572555042 0.559192903470747</values>
+            <!-- 34 degrees bedding angle -->
+            <values> 0.829037572555041692 0.55919290347074683016</values>
         </parameter>
         <parameter>
             <name>b1</name>
             <type>Constant</type>
-            <values>-0.559192903470747 0.829037572555042</values>
+            <!-- 34 degrees bedding angle -->
+            <values>-0.55919290347074683016 0.829037572555041692</values>
         </parameter>
         <parameter>
             <name>ThermalConductivityClay</name>
diff --git a/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj b/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj
index 10d9d7bb72a..b429403bc6f 100644
--- a/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj
+++ b/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj
@@ -194,12 +194,14 @@
         <parameter>
             <name>b0</name>
             <type>Constant</type>
-            <values> 0.829037572555042 0.559192903470747</values>
+            <!-- 34 degrees bedding angle -->
+            <values> 0.829037572555041692 0.55919290347074683016</values>
         </parameter>
         <parameter>
             <name>b1</name>
             <type>Constant</type>
-            <values>-0.559192903470747 0.829037572555042</values>
+            <!-- 34 degrees bedding angle -->
+            <values>-0.55919290347074683016 0.829037572555041692</values>
         </parameter>
         <parameter>
             <name>ThermalConductivityClay</name>
diff --git a/Tests/Data/ThermoRichardsMechanics/TaskCDECOVALEX2023/Decovalex-0.prj b/Tests/Data/ThermoRichardsMechanics/TaskCDECOVALEX2023/Decovalex-0.prj
index daff056884b..5420e1018eb 100644
--- a/Tests/Data/ThermoRichardsMechanics/TaskCDECOVALEX2023/Decovalex-0.prj
+++ b/Tests/Data/ThermoRichardsMechanics/TaskCDECOVALEX2023/Decovalex-0.prj
@@ -37,16 +37,11 @@
                 <displacement>displacement</displacement>
             </process_variables>
             <secondary_variables>
-                <secondary_variable internal_name="saturation"
-                                    output_name="saturation"/>
-                <secondary_variable internal_name="sigma"
-                                    output_name="sigma"/>
-                <secondary_variable internal_name="epsilon"
-                                    output_name="epsilon"/>
-                <secondary_variable internal_name="velocity"
-                                    output_name="velocity"/>
-                <secondary_variable internal_name="liquid_density"
-                                    output_name="liquid_density"/>
+                <secondary_variable internal_name="saturation" output_name="saturation"/>
+                <secondary_variable internal_name="sigma" output_name="sigma"/>
+                <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+                <secondary_variable internal_name="velocity" output_name="velocity"/>
+                <secondary_variable internal_name="liquid_density" output_name="liquid_density"/>
             </secondary_variables>
             <specific_body_force>0 0</specific_body_force>
             <initial_stress>Initial_stress</initial_stress>
@@ -552,14 +547,15 @@
         </output>
     </time_loop>
     <local_coordinate_system>
-        <basis_vector_0 implicit = "true" />
+        <basis_vector_0 implicit="true"/>
         <basis_vector_1>b1</basis_vector_1>
     </local_coordinate_system>
     <parameters>
         <parameter>
             <name>b1</name>
             <type>Constant</type>
-            <values>-0.559192903470747 0.829037572555042</values>
+            <!-- 34 degrees bedding angle -->
+            <values>-0.55919290347074683016 0.829037572555041692</values>
         </parameter>
         <parameter>
             <name>ThermalConductivityClay</name>
-- 
GitLab