More precise log

josephmure · Apr 4, 2024 · 2cd2163 · 2cd2163
1 parent 30d024a
commit 2cd2163
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/lib/src/Uncertainty/Distribution/ClaytonCopula.cxx b/lib/src/Uncertainty/Distribution/ClaytonCopula.cxx
@@ -392,7 +392,9 @@ Point ClaytonCopula::computeQuantile(const Scalar prob,
   // Independent case
   if (theta_ == 0.0) return Point(2, std::sqrt(q));
   // General case
-  LOGWARN(OSS() << "In ClaytonCopula::computeQuantile: q = " << q << " and return value - 1 = " << std::exp((M_LN2 - log1p(std::pow(q, -theta_))) / theta_) - 1.0);
+  LOGWARN(OSS() << "In ClaytonCopula::computeQuantile: q - 1 = " << q - 1.0 << " and theta_ = " << theta_);
+  LOGWARN(OSS() << "M_LN2 - log1p(std::pow(q, -theta_)) = " << M_LN2 - log1p(std::pow(q, -theta_)));
+  LOGWARN(OSS() << "Return value - 1 = " << std::exp((M_LN2 - log1p(std::pow(q, -theta_))) / theta_) - 1.0);
   return Point(2, std::exp((M_LN2 - log1p(std::pow(q, -theta_))) / theta_));
 }