Skip to content

Commit

Permalink
Revert "Attempt to deduce mass dimension of input parameters in simpl…
Browse files Browse the repository at this point in the history
…e cases"

This reverts commit fcc8967.
  • Loading branch information
Dylan Harries committed Feb 2, 2017
1 parent 544cdf9 commit a9cae3f
Show file tree
Hide file tree
Showing 2 changed files with 2 additions and 40 deletions.
3 changes: 0 additions & 3 deletions meta/FlexibleSUSY.m
Expand Up @@ -2795,9 +2795,6 @@ corresponding tadpole is real or imaginary (only in models with CP
Quit[1];
];

SemiAnalytic`FindAdditionalDimensionlessParameters[{FlexibleSUSY`LowScaleInput,
FlexibleSUSY`SUSYScaleInput,
FlexibleSUSY`HighScaleInput}];
semiAnalyticBCs = SemiAnalytic`SelectSemiAnalyticConstraint[{FlexibleSUSY`LowScaleInput,
FlexibleSUSY`SUSYScaleInput,
FlexibleSUSY`HighScaleInput}];
Expand Down
39 changes: 2 additions & 37 deletions meta/SemiAnalytic.m
Expand Up @@ -17,14 +17,10 @@
SetSemiAnalyticParameters::usage="";
GetSemiAnalyticParameters::usage="";
GetBoundaryValueParameters::usage="";
AddDimensionlessParameter::usage="Adds a parameter to be treated as
dimensionless when constructing the semi-analytic solutions.";

IsAllowedSemiAnalyticParameter::usage="";
IsSemiAnalyticParameter::usage="";

FindAdditionalDimensionlessParameters::usage="Attempts to deduce
the mass dimensions of user-defined parameters.";
GetSemiAnalyticSolutions::usage="Constructs the semi-analytic
solutions implied by the given list of boundary conditions.";
ExpandSemiAnalyticSolutions::usage="Expands the given solutions
Expand Down Expand Up @@ -66,7 +62,6 @@
Begin["`Private`"];

allSemiAnalyticParameters = {};
additionalDimensionlessPars = {};

GetName[SemiAnalyticSolution[name_, basis_List]] := name;

Expand All @@ -80,7 +75,7 @@

IsDimensionZero[par_] :=
Module[{dimZeroPars},
dimZeroPars = Join[Parameters`GetModelParametersWithMassDimension[0], additionalDimensionlessPars];
dimZeroPars = Parameters`GetModelParametersWithMassDimension[0];
MemberQ[Parameters`StripIndices[#]& /@ dimZeroPars, Parameters`StripIndices[par]]
];

Expand Down Expand Up @@ -145,8 +140,6 @@

GetSemiAnalyticParameters[] := allSemiAnalyticParameters;

AddDimensionlessParameter[par_ /; !IsDimensionZero[par]] := AppendTo[additionalDimensionlessPars, par];

IsSemiAnalyticSetting[setting_] :=
Intersection[Constraint`FindFixedParametersFromConstraint[{setting}],
allSemiAnalyticParameters] =!= {};
Expand Down Expand Up @@ -361,8 +354,7 @@
Module[{i, dimZeroPars, dimensions, expr,
monomials = {}, coeff, term, basis = {}},
If[pars =!= {},
dimZeroPars = Join[Parameters`GetModelParametersWithMassDimension[0],
additionalDimensionlessPars];
dimZeroPars = Parameters`GetModelParametersWithMassDimension[0];
dimensions = Parameters`GetParameterDimensions[#]& /@ pars;
expr = Plus @@ (SumOverElements[#]& /@ pars);
expr = Expand[expr /. subs];
Expand Down Expand Up @@ -453,33 +445,6 @@
GetLinearSystemSolutions[softLinears, boundaryConditionSubs, extraTerms]
];

SelectDimensionlessSettings[settings_List] :=
Module[{noMacros, dimZeroPars},
noMacros = Cases[settings, {_, value_ /; value =!= Automatic}];
dimZeroPars = Parameters`GetModelParametersWithMassDimension[0];
Select[noMacros, MemberQ[dimZeroPars, Parameters`StripIndices[#[[1]]]]&]
];

FindAdditionalDimensionlessParameters[constraints_List] :=
Module[{i, dimZeroSettings = {}, knownDimensionPars, unknownDimensionPars,
guessedDimZeroPars},
For[i = 1, i <= Length[constraints], i++,
dimZeroSettings = Join[dimZeroSettings, SelectDimensionlessSettings[constraints[[i]]]];
];
knownDimensionPars = Flatten[Table[Parameters`GetModelParametersWithMassDimension[i], {i, 0, 3}]];
unknownDimensionPars = Complement[Parameters`FindAllParameters[#[[2]]& /@ dimZeroSettings],
knownDimensionPars];
(* be conservative for now and consider only simplest case: settings of the form parameter = input *)
guessedDimZeroPars = DeleteDuplicates[Join[
Cases[dimZeroSettings, {_, v_ /; MemberQ[unknownDimensionPars, v]} :> v],
Cases[dimZeroSettings, {_, v_ CConversion`UNITMATRIX[n_Integer] /; MemberQ[unknownDimensionPars, v]} :> v],
Cases[dimZeroSettings, {_, v_ CConversion`UNITMATRIXCOMPLEX[n_Integer] /; MemberQ[unknownDimensionPars, v]} :> v]
]];
Print["Warning: the following parameters will be assumed"];
Print[" to be dimensionless: ", InputForm[guessedDimZeroPars]];
AddDimensionlessParameter /@ guessedDimZeroPars;
];

GetSemiAnalyticSolutions[settings_List] :=
Module[{boundaryValSubs, dimOnePars, dimOneBCs, dimOneSolns = {},
diracPars, diracBCs, diracSolns = {}, scalarMasses, scalarMassesBCs,
Expand Down

0 comments on commit a9cae3f

Please sign in to comment.