From 8170fb10d893039f4aa87a057a54558dc018d3a7 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Thu, 4 Apr 2024 11:34:39 +0200
Subject: [PATCH] [PL/TH2M] Default initialize scalars to NaN.

---
 ProcessLib/TH2M/ConstitutiveRelations/Base.h              | 4 ++--
 ProcessLib/TH2M/ConstitutiveRelations/FluidDensity.h      | 4 ++--
 ProcessLib/TH2M/ConstitutiveRelations/MassMoleFractions.h | 8 ++++----
 .../TH2M/ConstitutiveRelations/SolidThermalExpansion.h    | 4 ++--
 4 files changed, 10 insertions(+), 10 deletions(-)

diff --git a/ProcessLib/TH2M/ConstitutiveRelations/Base.h b/ProcessLib/TH2M/ConstitutiveRelations/Base.h
index d0f26028eca..a2bcea4a0dc 100644
--- a/ProcessLib/TH2M/ConstitutiveRelations/Base.h
+++ b/ProcessLib/TH2M/ConstitutiveRelations/Base.h
@@ -47,8 +47,8 @@ struct MediaData
 
 struct TemperatureData
 {
-    double T;
-    double T_prev;
+    double T = nan;
+    double T_prev = nan;
 };
 
 using ReferenceTemperatureData =
diff --git a/ProcessLib/TH2M/ConstitutiveRelations/FluidDensity.h b/ProcessLib/TH2M/ConstitutiveRelations/FluidDensity.h
index 067bc887273..22a0fa60d87 100644
--- a/ProcessLib/TH2M/ConstitutiveRelations/FluidDensity.h
+++ b/ProcessLib/TH2M/ConstitutiveRelations/FluidDensity.h
@@ -19,10 +19,10 @@ namespace ConstitutiveRelations
 struct FluidDensityData
 {
     // gas phase density
-    double rho_GR = 0.;
+    double rho_GR = nan;
 
     // liquid phase density
-    double rho_LR = 0.;
+    double rho_LR = nan;
 
     static auto reflect()
     {
diff --git a/ProcessLib/TH2M/ConstitutiveRelations/MassMoleFractions.h b/ProcessLib/TH2M/ConstitutiveRelations/MassMoleFractions.h
index a63b2514a47..a180d179d83 100644
--- a/ProcessLib/TH2M/ConstitutiveRelations/MassMoleFractions.h
+++ b/ProcessLib/TH2M/ConstitutiveRelations/MassMoleFractions.h
@@ -18,10 +18,10 @@ namespace ConstitutiveRelations
 {
 struct MassMoleFractionsData
 {
-    double xnCG = 0.;
-    double xmCG = 0.;
-    double xnWL = 0.;
-    double xmWL = 0.;
+    double xnCG = nan;
+    double xmCG = nan;
+    double xnWL = nan;
+    double xmWL = nan;
 
     static auto reflect()
     {
diff --git a/ProcessLib/TH2M/ConstitutiveRelations/SolidThermalExpansion.h b/ProcessLib/TH2M/ConstitutiveRelations/SolidThermalExpansion.h
index 9f225f30f7d..3eeab68c3d2 100644
--- a/ProcessLib/TH2M/ConstitutiveRelations/SolidThermalExpansion.h
+++ b/ProcessLib/TH2M/ConstitutiveRelations/SolidThermalExpansion.h
@@ -19,8 +19,8 @@ template <int DisplacementDim>
 struct SolidThermalExpansionData
 {
     KelvinVector<DisplacementDim> solid_linear_thermal_expansivity_vector;
-    double beta_T_SR;  /// Isotropic solid phase volumetric thermal expansion
-                       /// coefficient.
+    double beta_T_SR = nan;  /// Isotropic solid phase volumetric thermal
+                             /// expansion coefficient.
     double thermal_volume_strain = nan;
 };
 
-- 
GitLab