Implement clear and initialize consistenly for low-scale constraints

FlexibleSUSY · Jan 29, 2017 · b6d7953 · b6d7953
1 parent 3551ec3
commit b6d7953
Showing 1 changed file with 7 additions and 5 deletions.
diff --git a/templates/semi_analytic_low_scale_constraint.cpp.in b/templates/semi_analytic_low_scale_constraint.cpp.in
@@ -102,6 +102,9 @@ void @ModelName@_low_scale_constraint<Semi_analytic>::clear()
    initial_scale_guess = 0.;
    model = nullptr;
    qedqcd = softsusy::QedQcd();
+   ckm.setIdentity();
+   pmns.setIdentity();
+   neutrinoDRbar.setZero();
    MWDRbar = 0.;
    MZDRbar = 0.;
    AlphaS = 0.;
@@ -121,6 +124,9 @@ void @ModelName@_low_scale_constraint<Semi_analytic>::initialize()
 @scaleGuess@
    scale = initial_scale_guess;
 
+   ckm = qedqcd.get_complex_ckm();
+   pmns = qedqcd.get_complex_pmns();
+   neutrinoDRbar = Eigen::Matrix<double,3,3>::Zero();
    MWDRbar = 0.;
    MZDRbar = 0.;
    AlphaS = 0.;
@@ -129,8 +135,6 @@ void @ModelName@_low_scale_constraint<Semi_analytic>::initialize()
    new_g1 = 0.;
    new_g2 = 0.;
    new_g3 = 0.;
-   ckm = qedqcd.get_complex_ckm();
-   pmns = qedqcd.get_complex_pmns();
    self_energy_w_at_mw = 0.;
 }
 
@@ -294,8 +298,6 @@ void @ModelName@_low_scale_constraint<Semi_analytic>::calculate_@ElectronYukawa@
 
 void @ModelName@_low_scale_constraint<Semi_analytic>::calculate_MNeutrino_DRbar()
 {
-   check_model_ptr();
-
    neutrinoDRbar.setZero();
    neutrinoDRbar(0,0) = qedqcd.displayNeutrinoPoleMass(1);
    neutrinoDRbar(1,1) = qedqcd.displayNeutrinoPoleMass(2);
@@ -306,7 +308,7 @@ void @ModelName@_low_scale_constraint<Semi_analytic>::calculate_MNeutrino_DRbar(
  * Recalculates the W boson pole mass using the new gauge couplings.
  *
  * @param mw_pole current W pole mass
- * @param new W pole mass
+ * @return new W pole mass
  */
 double @ModelName@_low_scale_constraint<Semi_analytic>::recalculate_mw_pole(double mw_pole)
 {