change calculate_amuon() interface: use const ref to model class

in order to guarantee that the model class will not be modified after the function call
FlexibleSUSY · Sep 13, 2015 · 1fb57f3 · 1fb57f3
1 parent d84a60d
commit 1fb57f3
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 2 deletions.
diff --git a/templates/g_muon_minus_2.cpp.in b/templates/g_muon_minus_2.cpp.in
@@ -527,8 +527,11 @@ static inline double OneLoopFunctionF2N(double x)
    return -3.0 / (y * y * y) * (1.0 - x * x + 2.0 * x * std::log(x));
 }
 
-double calculate_amuon(@ModelName@<Two_scale> &model)
+double calculate_amuon(const @ModelName@<Two_scale>& model_)
 {
+   // make copy since model is modified via call to calculate_MFe_pole()
+   @ModelName@<Two_scale> model(model_);
+
 #ifdef ENABLE_THREADS
 @GMuonMinus2_ThreadedCalculation@
 #else

diff --git a/templates/g_muon_minus_2.hpp.in b/templates/g_muon_minus_2.hpp.in
@@ -41,7 +41,7 @@ namespace @ModelName@_GMuonMinus2 {
 * @fn calculate_amuon
 * @brief Calculates GMM2/2.
 */
-double calculate_amuon(@ModelName@<Two_scale> &model);
+double calculate_amuon(const @ModelName@<Two_scale>& model);
 };
 } // namespace flexiblesusy