Add function testing for linearly independent equations

FlexibleSUSY · Dec 7, 2016 · 9cb28c6 · 9cb28c6
1 parent 7d2ea17
commit 9cb28c6
Showing 1 changed file with 14 additions and 0 deletions.
diff --git a/meta/EWSB.m b/meta/EWSB.m
@@ -1,6 +1,9 @@
 
 BeginPackage["EWSB`", {"SARAH`", "TextFormatting`", "CConversion`", "Parameters`", "TreeMasses`", "WriteOut`", "Utils`"}];
 
+GetLinearlyIndependentEqs::"Removes linearly dependent EWSB equations
+from a list of equations";
+
 FindSolutionAndFreePhases::usage="Finds solution to the EWSB and free
 phases / signs."
 
@@ -70,6 +73,17 @@
 CheckInEquations[parameter_, statement_, equations_List] :=
     And @@ (statement[parameter,#]& /@ equations);
 
+GetLinearlyIndependentEqs[eqs_List, parameters_List, substitutions_List:{}] :=
+    Module[{eqsToSolve, indepEqsToSolve, eqsToKeep},
+           If[substitutions =!= {},
+              eqsToSolve = eqs /. (Rule[#[[1]], #[[2]]]& /@ substitutions);,
+              eqsToSolve = eqs;
+             ];
+           indepEqsToSolve = Parameters`FilterOutLinearDependentEqs[eqsToSolve, parameters];
+           eqsToKeep = Position[eqsToSolve, p_ /; MemberQ[indepEqsToSolve, p]];
+           Extract[eqs, eqsToKeep]
+          ];
+
 CreateEWSBEqPrototype[higgs_] :=
     Module[{result = "", i, ctype},
            ctype = CConversion`CreateCType[CConversion`ScalarType[CConversion`realScalarCType]];