From 0cf78bb9af30b753b204b8d0a177b665ae51a67e Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Mon, 2 Sep 2024 13:48:10 +0200
Subject: [PATCH] [T/MPL] Change temperature eval point due to limit
function is defined only on the range [273, +inf)
---
Tests/MaterialLib/TestMPLWaterVapourDensity.cpp | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/Tests/MaterialLib/TestMPLWaterVapourDensity.cpp b/Tests/MaterialLib/TestMPLWaterVapourDensity.cpp
index 16bde0ae8a8..88830c020ec 100644
--- a/Tests/MaterialLib/TestMPLWaterVapourDensity.cpp
+++ b/Tests/MaterialLib/TestMPLWaterVapourDensity.cpp
@@ -44,7 +44,9 @@ TEST(MaterialPropertyLib, WaterVapourDensity)
// The derivative of the water vapour with respect of temperature
{
- std::array const temperatures = {273.0, 293.0, 393.0, 420.0, 500.0};
+ // The first evaluation point 273.0001 is chosen to avoid lower limit
+ // (T=273K) of the function in when computing numerical derivative.
+ std::array const temperatures = {273.0001, 293.0, 393.0, 420.0, 500.0};
std::array const rho_vw_expected = {4.875989e-03, 1.692871e-02,
1.276865, 2.882635, 1.920386e+01};
--
GitLab