Skip to content

Commit

Permalink
Update EWSB test cases
Browse files Browse the repository at this point in the history
  • Loading branch information
Dylan Harries committed Jan 23, 2017
1 parent 406137c commit d95c385
Showing 1 changed file with 11 additions and 9 deletions.
20 changes: 11 additions & 9 deletions test/test_EWSB.m
Expand Up @@ -342,34 +342,36 @@

Print["testing EWSB substitutions ..."];

Parameters`SetModelParameters[{\[Mu], B[\[Mu]], mHd2, mHu2}];

subEwsbEqs = {
\[Mu]^2 + x^2 + x y + z + 5,
B[\[Mu]] - x^2 + x y + z + 5
B[\[Mu]] - x^2 + x y + z + 5
};

ewsbSubs = {
{\[Mu], Sign[\[Mu]] Sqrt[MuSqr]}
Rule[\[Mu], Sign[\[Mu]] Sqrt[MuSqr]]
};

subEwsbOutputParameters = { MuSqr, B[\[Mu]] };

Parameters`SetRealParameters[subEwsbOutputParameters];

{subSolution, subFreePhases} = EWSB`FindSolutionAndFreePhases[subEwsbEqs, subEwsbOutputParameters, ewsbSubs];
{subSolution, subFreePhases} = EWSB`FindSolutionAndFreePhases[subEwsbEqs, subEwsbOutputParameters, "", ewsbSubs];

TestEquality[subFreePhases, {}];
TestEquality[Sort[Rule[#[[1]],Expand[#[[2]]]]& /@ subSolution],
Sort[Rule[#[[1]],Expand[#[[2]]]]& /@ {MuSqr -> -(x^2 + x y + z + 5) / Sign[\[Mu]]^2,
B[\[Mu]] -> x^2 - x y - z - 5}]];

ewsbSubs = {
{\[Mu], Sign[\[Mu]] Sqrt[MuSqr]},
{B[\[Mu]], BMu0}
Rule[\[Mu], Sign[\[Mu]] Sqrt[MuSqr]],
Rule[B[\[Mu]], BMu0]
};

subEwsbOutputParameters = { MuSqr, BMu0 };

{subSolution, subFreePhases} = EWSB`FindSolutionAndFreePhases[subEwsbEqs, subEwsbOutputParameters, ewsbSubs];
{subSolution, subFreePhases} = EWSB`FindSolutionAndFreePhases[subEwsbEqs, subEwsbOutputParameters, "", ewsbSubs];

TestEquality[subFreePhases, {}];
TestEquality[Sort[Rule[#[[1]],Expand[#[[2]]]]& /@ subSolution],
Expand All @@ -384,15 +386,15 @@
};

ewsbSubs = {
{mHd2, m0^2 + m12^2},
{mHu2, m0^2 + AzeroSqr}
Rule[mHd2, m0^2 + m12^2],
Rule[mHu2, m0^2 + AzeroSqr]
};

subEwsbOutputParameters = {m12, AzeroSqr};

Parameters`SetRealParameters[subEwsbOutputParameters];

{subSolution, subFreePhases} = EWSB`FindSolutionAndFreePhases[subEwsbEqs, subEwsbOutputParameters, ewsbSubs];
{subSolution, subFreePhases} = EWSB`FindSolutionAndFreePhases[subEwsbEqs, subEwsbOutputParameters, "", ewsbSubs];

TestEquality[subFreePhases, {Sign[m12]}];
TestEquality[Sort[Rule[#[[1]],Expand[#[[2]]]]& /@ subSolution],
Expand Down

0 comments on commit d95c385

Please sign in to comment.