Fix numerical threshold for ReduceVEV and corrected output (#162)

* Implemented Dilogarithm/Spencer's funtion and unit tests for it. * Espinosa-Konstandin - Example A * Espinosa-Konstandin - Example B * Correct typo * Fix wrong sign in numerical derivative of analytical solution when d2Vdl < 0 and alpha = 3 (O(4)-symmetric solution) * Fix recursive step that assumed l0 = 0. * Protect against infinite recursion in "BounceActionInt::ExactSolutionLin()" * Espinosa-Konstandin - Example D * Include "gsl_sf_expint" in tests * Espinosa-Konstandin - Example C * Espinosa-Konstandin - Example E * Espinosa-Konstandin - Example 2D * adjust output column names * correct output text Pf -> I(T) * add numerical threshold for ReduceVEV * add functional include * write false vacuum fraction to logginglevel * Version bump * Apply suggestions from code review * restore file --------- Co-authored-by: João Viana <[email protected]> Co-authored-by: Philipp Basler <[email protected]>
phbasler · Sep 4, 2024 · 7302d69 · 7302d69
1 parent 2b49797
commit 7302d69
Show file tree

Hide file tree

Showing 4 changed files with 5 additions and 5 deletions.
diff --git a/CMakeLists.txt b/CMakeLists.txt
@@ -6,7 +6,7 @@
 cmake_minimum_required(VERSION 3.23)
 project(
   BSMPT
-  VERSION 3.0.6
+  VERSION 3.0.7
   LANGUAGES C CXX
   DESCRIPTION
     "BSMPT - Beyond the Standard Model Phase Transitions : A C++ package for the computation of the EWPT in BSM models"

diff --git a/README.md b/README.md
@@ -4,7 +4,7 @@ SPDX-FileCopyrightText: 2021 Philipp Basler, Margarete Mühlleitner and Jonas M
 SPDX-License-Identifier: GPL-3.0-or-later
 -->
 
-Program: BSMPT version 3.0.6
+Program: BSMPT version 3.0.7
 
 Released by: Philipp Basler, Lisa Biermann, Margarete Mühlleitner, Jonas Müller, Rui Santos and João Viana
 

diff --git a/src/bounce_solution/bounce_solution.cpp b/src/bounce_solution/bounce_solution.cpp
@@ -744,7 +744,7 @@ double BounceSolution::CalcTempAtFalseVacFraction(const double &false_vac_frac)
   }
 
   Logger::Write(LoggingLevel::BounceDetailed,
-                "T ( Pf = " + std::to_string(int_at_false_vac_frac) +
+                "T ( Pf = " + std::to_string(std::exp(-int_at_false_vac_frac)) +
                     " ) is in interval [ " + std::to_string(T_down) + ", " +
                     std::to_string(T_up) + " ]");
 
@@ -769,7 +769,7 @@ double BounceSolution::CalcTempAtFalseVacFraction(const double &false_vac_frac)
 
       Logger::Write(LoggingLevel::BounceDetailed,
                     "Pf ( T = " + std::to_string(T_middle) +
-                        " ) = " + std::to_string(IatT));
+                        " ) = " + std::to_string(std::exp(-IatT)));
 
       if (IatT < int_at_false_vac_frac)
       {

diff --git a/src/minimum_tracer/minimum_tracer.cpp b/src/minimum_tracer/minimum_tracer.cpp
@@ -711,7 +711,7 @@ void MinimumTracer::ReduceVEV(std::vector<double> &vev)
     // Feed buffer
     for (size_t i = 0; i < vev.size(); i++)
       BinaryNumber.append(std::to_string(vev[i] * GroupElement(i, i) >=
-                                         0)); // Heaviside function
+                                         1e-5)); // Heaviside function
 
     // Converto to decimal and compare
     if (std::strtol(BinaryNumber.c_str(), &ptr, 2) > MaximumMeasure)