Update NMSSMCPV tests to use separate EWSB solver

FlexibleSUSY · Mar 1, 2017 · c24f20b · c24f20b
1 parent e354519
commit c24f20b
Show file tree

Hide file tree

Showing 2 changed files with 27 additions and 0 deletions.
diff --git a/test/test_NMSSMCPV_ewsb.cpp b/test/test_NMSSMCPV_ewsb.cpp
@@ -8,6 +8,7 @@
 #include "test_NMSSMCPV.hpp"
 #include "wrappers.hpp"
 #include "ew_input.hpp"
+#include "NMSSMCPV_two_scale_ewsb_solver.hpp"
 #include "NMSSMCPV_two_scale_model.hpp"
 
 using namespace flexiblesusy;
@@ -36,6 +37,10 @@ BOOST_AUTO_TEST_CASE( test_NMSSMCPV_ewsb_tree_level )
    const double precision = 1.0e-5;
    setup_NMSSMCPV(m, input);
 
+   NMSSMCPV_ewsb_solver<Two_scale> ewsb_solver;
+   m.set_ewsb_solver(
+      std::make_shared<NMSSMCPV_ewsb_solver<Two_scale> >(ewsb_solver));
+
    // initial guess
    m.set_mHu2(-Sqr(input.m0));
    m.set_mHd2(Sqr(input.m0));
@@ -76,6 +81,10 @@ BOOST_AUTO_TEST_CASE( test_NMSSMCPV_ewsb_one_loop )
    const double precision = 1.0e-10;
    setup_NMSSMCPV(m, input);
 
+   NMSSMCPV_ewsb_solver<Two_scale> ewsb_solver;
+   m.set_ewsb_solver(
+      std::make_shared<NMSSMCPV_ewsb_solver<Two_scale> >(ewsb_solver));
+
    // initial guess
    m.set_mHu2(-Sqr(input.m0));
    m.set_mHd2(Sqr(input.m0));
@@ -198,6 +207,10 @@ BOOST_AUTO_TEST_CASE( test_NMSSMCPV_one_loop_tadpoles_real_limit )
    NMSSMCPV<Two_scale> m;
    setup_NMSSMCPV(m, input);
 
+   NMSSMCPV_ewsb_solver<Two_scale> ewsb_solver;
+   m.set_ewsb_solver(
+      std::make_shared<NMSSMCPV_ewsb_solver<Two_scale> >(ewsb_solver));
+
    // initial guess
    m.set_mHu2(-Sqr(input.m0));
    m.set_mHd2(Sqr(input.m0));
@@ -234,6 +247,10 @@ BOOST_AUTO_TEST_CASE( test_NMSSMCPV_one_loop_tadpoles )
    const double precision = 1.0e-5;
    setup_NMSSMCPV(m, input);
 
+   NMSSMCPV_ewsb_solver<Two_scale> ewsb_solver;
+   m.set_ewsb_solver(
+      std::make_shared<NMSSMCPV_ewsb_solver<Two_scale> >(ewsb_solver));
+
    // initial guess
    m.set_mHu2(-Sqr(input.m0));
    m.set_mHd2(Sqr(input.m0));

diff --git a/test/test_NMSSMCPV_tree_level_spectrum.cpp b/test/test_NMSSMCPV_tree_level_spectrum.cpp
@@ -7,6 +7,7 @@
 #include "test.hpp"
 #include "test_NMSSMCPV.hpp"
 #include "NMSSM_two_scale_model.hpp"
+#include "NMSSMCPV_two_scale_ewsb_solver.hpp"
 #include "NMSSMCPV_two_scale_model.hpp"
 
 #define COMPARE_MASS(p,dev) TEST_CLOSE(a.get_##p(), b.get_##p(), dev);
@@ -96,6 +97,11 @@ BOOST_AUTO_TEST_CASE( test_NMSSMCPV_tree_level_spectrum )
    NMSSM<Two_scale> m2;
 
    setup_NMSSMCPV(m1, input);
+
+   NMSSMCPV_ewsb_solver<Two_scale> ewsb_solver;
+   m1.set_ewsb_solver(
+      std::make_shared<NMSSMCPV_ewsb_solver<Two_scale> >(ewsb_solver));
+
    m1.solve_ewsb_tree_level();
    copy_parameters(m1, m2);
    m1.calculate_DRbar_masses();
@@ -115,6 +121,10 @@ void check_goldstone_masses(const NMSSMCPV_input_parameters& input)
    NMSSMCPV<Two_scale> m;
    setup_NMSSMCPV(m, input);
 
+   NMSSMCPV_ewsb_solver<Two_scale> ewsb_solver;
+   m.set_ewsb_solver(
+      std::make_shared<NMSSMCPV_ewsb_solver<Two_scale> >(ewsb_solver));
+
    m.solve_ewsb_tree_level();
    m.calculate_DRbar_masses();
    m.reorder_DRbar_masses();