diff --git a/CMakeLists.txt b/CMakeLists.txt
index 69312602c..c4e14e316 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -6,7 +6,7 @@
 cmake_minimum_required(VERSION 3.23)
 project(
   BSMPT
-  VERSION 3.0.5
+  VERSION 3.0.8
   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
index 53b9077d4..8b55c5fff 100644
--- 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.5
+Program: BSMPT version 3.0.8
 
 Released by: Philipp Basler, Lisa Biermann, Margarete Mühlleitner, Jonas Müller, Rui Santos and João Viana
 
diff --git a/include/BSMPT/utility/const_velocity_spline.h b/include/BSMPT/utility/const_velocity_spline.h
index c2add41f7..2f4429f80 100644
--- a/include/BSMPT/utility/const_velocity_spline.h
+++ b/include/BSMPT/utility/const_velocity_spline.h
@@ -189,12 +189,12 @@ class cvspline
    * @param l
    * @return std::vector<double>
    */
-  std::vector<double> operator()(double l);
+  std::vector<double> operator()(double l) const;
   /**
    * @brief print the current knots of the spline
    *
    */
-  void print_path();
+  void print_path() const;
   /**
    * @brief save the knots of the splien into a file
    *
diff --git a/src/bounce_solution/bounce_solution.cpp b/src/bounce_solution/bounce_solution.cpp
index 78772fcc1..a1ab130eb 100644
--- a/src/bounce_solution/bounce_solution.cpp
+++ b/src/bounce_solution/bounce_solution.cpp
@@ -756,7 +756,7 @@ double BounceSolution::CalcTempAtFalseVacFraction(const double &false_vac_frac)
 
     while ((std::abs(T_up / T_down - 1) >
             RelativeTemperatureInCalcTempAtFalseVacFraction *
-                MarginOfCalcTempAtFalseVacFractionBeforeFailure) or
+                MarginOfCalcTempAtFalseVacFractionBeforeFailure) and
            (not almost_the_same(
                int_at_false_vac_frac,
                IatT,
diff --git a/src/utility/const_velocity_spline.cpp b/src/utility/const_velocity_spline.cpp
index c10b36f42..907e39d55 100644
--- a/src/utility/const_velocity_spline.cpp
+++ b/src/utility/const_velocity_spline.cpp
@@ -266,7 +266,7 @@ std::vector<double> cvspline::d2l(double l)
 
   return rr;
 }
-std::vector<double> cvspline::operator()(double l)
+std::vector<double> cvspline::operator()(double l) const
 {
   // Calculates the VEV at the position l in spline length,
   l = l_to_x(l); // Convert from spline length to linear length.
@@ -278,7 +278,7 @@ std::vector<double> cvspline::operator()(double l)
   return r;
 }
 
-void cvspline::print_path()
+void cvspline::print_path() const
 {
   int wid = 15;
 
diff --git a/standalone/TunnelingPath.cpp b/standalone/TunnelingPath.cpp
index 1954a8627..552aeb127 100644
--- a/standalone/TunnelingPath.cpp
+++ b/standalone/TunnelingPath.cpp
@@ -77,10 +77,10 @@ int main()
     bounce.CalculatePercolationTemp();
     if (bounce.GetPercolationTemp() == -1) continue;
     double errorTtoTp = 1e100;
-    BounceActionInt *ClosestBounceActionInt;
+    BounceActionInt const *ClosestBounceActionInt;
     std::cout << "Found a transitions with Tp =\t"
               << bounce.GetPercolationTemp() << " GeV.\n";
-    for (auto BAInt : bounce.SolutionList)
+    for (const auto &BAInt : bounce.SolutionList)
     {
       if (abs(BAInt.T - bounce.GetPercolationTemp()) < errorTtoTp)
       {