make sure the loop factor is factored out

FlexibleSUSY · Dec 2, 2015 · a024111 · a024111
1 parent f675a3f
commit a024111
Showing 1 changed file with 12 additions and 1 deletion.
diff --git a/meta/BetaFunction.m b/meta/BetaFunction.m
@@ -106,9 +106,20 @@
              ]
           ];
 
+FactorOutLoopFactor[expr_] :=
+    Module[{i, prefactors = {CConversion`oneOver16PiSqr, CConversion`twoLoop,
+                             CConversion`threeLoop, CConversion`oneOver16PiSqr^4}},
+           For[i = 1, i <= Length[prefactors], i++,
+               If[Coefficient[expr, prefactors[[i]]] =!= 0,
+                  Return[Simplify[prefactors[[i]] Expand[expr / prefactors[[i]]]]];
+                 ];
+              ];
+           expr
+          ];
+
 (* split expression into sub-expressions of given maximum size *)
 SplitExpression[expr_, size_Integer] :=
-    Factor /@ (Plus @@@ Utils`SplitList[ToList[expr, Plus], size]);
+    FactorOutLoopFactor /@ (Plus @@@ Utils`SplitList[ToList[expr, Plus], size]);
 
 NeedToSplitExpression[expr_, threshold_Integer] :=
     Length[ToList[expr, Plus]] > threshold;