testing MSSM 3L soft-breaking beta functions

some tests still fail due to missing MatMul[] head. Note: The JJ result misses the Tr1[1] term.
FlexibleSUSY · Dec 1, 2015 · 09ca19f · 09ca19f
1 parent cb4a527
commit 09ca19f
Showing 1 changed file with 18 additions and 1 deletion.
diff --git a/test/test_CMSSM_3loop_beta.m b/test/test_CMSSM_3loop_beta.m
@@ -30,7 +30,13 @@
 TestBetaEquality[lst_, c_, loop_] :=
     Module[{sa, sm},
            sm = ThreeLoopMSSM`BetaMSSM[c][[loop]] /. UniformTraces[];
-           sa = FindBetaFunction[lst, c][[loop]] /. a_[i1,i2] :> a /. UniformTraces[];
+           sa = FindBetaFunction[lst, c][[loop]] /. {
+               Kronecker[_,_] :> 1,
+               a_[i1,i2] :> a,
+               conj[a_] :> a,
+               Tr1[1] -> 0 (* why is this term missing in the JJ result? *)
+           } /. UniformTraces[];
+           Print["Difference: ", Expand[sm] - Expand[sa]]
            TestEquality[Expand[sm], Expand[sa]]
           ];
 
@@ -41,6 +47,17 @@
     TestBetaEquality[SARAH`BetaYijk , SARAH`UpYukawa           , l];
     TestBetaEquality[SARAH`BetaYijk , SARAH`DownYukawa         , l];
     TestBetaEquality[SARAH`BetaYijk , SARAH`ElectronYukawa     , l];
+    TestBetaEquality[SARAH`BetaMi   , MassB                    , l];
+    TestBetaEquality[SARAH`BetaMi   , MassWB                   , l];
+    TestBetaEquality[SARAH`BetaMi   , MassG                    , l];
+    TestBetaEquality[SARAH`BetaTijk , SARAH`TrilinearUp        , l];
+    TestBetaEquality[SARAH`BetaTijk , SARAH`TrilinearDown      , l];
+    TestBetaEquality[SARAH`BetaTijk , SARAH`TrilinearLepton    , l];
+    TestBetaEquality[SARAH`Betam2ij , SARAH`SoftSquark         , l];
+    TestBetaEquality[SARAH`Betam2ij , SARAH`SoftUp             , l];
+    TestBetaEquality[SARAH`Betam2ij , SARAH`SoftDown           , l];
+    TestBetaEquality[SARAH`Betam2ij , SARAH`SoftLeftLepton     , l];
+    TestBetaEquality[SARAH`Betam2ij , SARAH`SoftRightLepton    , l];
    ];
 
 PrintTestSummary[];