correcting counting of EWSB eqs and phases

meta/Constraint.m

@@ -491,7 +491,9 @@
     StringJoin[CheckPerturbativityForParameter[#,thresh]& /@ pars];
 
 GetNumberOfParameters[FlexibleSUSY`FSSolveEWSBFor[pars_List]] :=
-    Plus @@ (BetaFunction`CountNumberOfParameters[Parameters`GetType[#]]& /@ pars);
+    Length[pars];
+
+GetNumberOfParameters[{FlexibleSUSY`Phase[_], _}] := 1;
 
 GetNumberOfParameters[{par_[indices__] /; And @@ (IsIndex /@ {indices}), _}] :=
     BetaFunction`CountNumberOfParameters[CConversion`GetScalarElementType[Parameters`GetType[par]]];
@@ -502,6 +504,12 @@
 CalculateNumberOfParameters[settings_List] :=
     ToString[Plus @@ (GetNumberOfParameters /@ settings)];
 
+(* phases are equal if their real parts are equal *)
+CalculateDifference[{p:FlexibleSUSY`Phase[_], value_}, offset_Integer, vector_String, type_] :=
+    vector <> "(" <> ToString[offset] <> ") = Re(" <>
+    CConversion`ToValidCSymbolString[p] <> " - " <>
+    CConversion`RValueToCFormString[value] <> ");\n";
+
 CalculateDifference[{par_, Automatic}, offset_Integer, vector_String, type_] :=
     "calculate_" <> CConversion`ToValidCSymbolString[par] <> "_DRbar();\n" <>
     CalculateDifference[{par, "MODELPARAMETER(" <> CConversion`ToValidCSymbolString[par] <> ")"}, offset, vector, type];
@@ -553,6 +561,7 @@
 CalculateDifference[{par_, value_}, offset_Integer, vector_String] :=
     CalculateDifference[{par,value}, offset, vector, Parameters`GetType[par]];
 
+(* @todo: handle case of complex parameters and/or Re[Mu] etc. *)
 CalculateDifference[FlexibleSUSY`FSSolveEWSBFor[pars_List], offset_Integer, vector_String] :=
     Module[{i, currentOffset = offset, result = "", parStr},
            For[i = 1, i <= Length[pars], i++,