Add skeleton semi-analytic solver module

meta/FlexibleSUSY.m

@@ -5,7 +5,7 @@
              "Phases`", "LoopMasses`", "WriteOut`", "Constraint`",
              "ThresholdCorrections`", "ConvergenceTester`", "Utils`",
              "ThreeLoopSM`", "ThreeLoopMSSM`", "Observables`",
-             "GMuonMinus2`", "EffectiveCouplings`",
+             "GMuonMinus2`", "EffectiveCouplings`", "SemiAnalytic`",
              "FlexibleEFTHiggsMatching`", "FSMathLink`",
              "FlexibleTower`"}];
 

meta/SemiAnalytic.m

@@ -0,0 +1,60 @@
+
+BeginPackage["SemiAnalytic`", {"SARAH`", "CConversion`", "Constraint`", "Parameters`",
+                               "TextFormatting`", "WriteOut`"}];
+
+Begin["`Private`"];
+
+IsDimensionOne[par_] :=
+    Module[{dimOnePars},
+           dimOnePars = { SARAH`BetaTijk };
+           If[SARAH`SupersymmetricModel,
+              dimOnePars = Append[dimOnePars, SARAH`BetaMi];,
+              dimOnePars = Append[dimOnePars, SARAH`BetaMuij];
+             ];
+           dimOnePars = (Parameters`StripIndices[#[[1]]])& /@ (Join @@ dimOnePars);
+           MemberQ[dimOnePars, Parameters`StripIndices[par]]
+          ];
+
+IsDiracGauginoMass[par_] :=
+    Module[{diracMasses = {}},
+           If[SARAH`SupersymmetricModel,
+              diracMasses = Parameters`StripIndices[#[[1]]]& /@ SARAH`BetaDGi;
+             ];
+           MemberQ[diracMasses, Parameters`StripIndices[par]]
+          ];
+
+IsScalarMass[par_] :=
+    Module[{scalarMasses},
+           If[SARAH`SupersymmetricModel,
+              scalarMasses = Parameters`StripIndices[#[[1]]]& /@ SARAH`Betam2ij;,
+              scalarMasses = Parameters`StripIndices[#[[1]]]& /@ SARAH`BetaBij;
+             ];
+           MemberQ[scalarMasses, Parameters`StripIndices[par]]
+          ];
+
+IsSoftBilinear[par_] :=
+    Module[{softBilinears = {}},
+           If[SARAH`SupersymmetricModel,
+              softBilinears = Parameters`StripIndices[#[[1]]]& /@ SARAH`BetaBij;
+             ];
+           MemberQ[softBilinears, Parameters`StripIndices[par]]
+          ];
+
+IsSoftLinear[par_] :=
+    Module[{softLinears = {}},
+           If[SARAH`SupersymmetricModel,
+              softLinears = Parameters`StripIndices[#[[1]]]& /@ SARAH`BetaLSi;
+             ];
+           MemberQ[softLinears, Parameters`StripIndices[par]]
+          ];
+
+IsAllowedSemiAnalyticParameter[par_] :=
+    Or[IsDimensionOne[par],
+       IsDiracGauginoMass[par],
+       IsScalarMass[par],
+       IsSoftBilinear[par],
+       IsSoftLinear[par]];
+
+End[];
+
+EndPackage[];

meta/module.mk

@@ -68,6 +68,7 @@ META_SRC     := \
 		$(DIR)/ReadSLHA.m \
 		$(DIR)/RGIntegrator.m \
 		$(DIR)/SelfEnergies.m \
+		$(DIR)/SemiAnalytic.m \
 		$(DIR)/TestSuite.m \
 		$(DIR)/TextFormatting.m \
 		$(DIR)/ThreeLoopMSSM.m \