Bugfix: extend IncreaseIndexLiterals[] function to nested functions

The following example call of IncreaseIndexLiterals[] now works: In[]:= IncreaseIndexLiterals[Hpm[{2}][PL]] Out[]= Hpm[{3}][PL] Thanks to Markus Bach.
FlexibleSUSY · Feb 16, 2017 · 2bb7259 · 2bb7259
1 parent 3827dba
commit 2bb7259
Showing 1 changed file with 3 additions and 4 deletions.
diff --git a/meta/Parameters.m b/meta/Parameters.m
@@ -1344,12 +1344,11 @@
                                           allOutputParameters]];
 
 IncreaseIndexLiterals[expr_, num_Integer, heads_List] :=
-    Module[{indexedSymbols, rules, decrExpr, allHeads},
+    Module[{indexedSymbols, rules, allHeads},
            allHeads = Join[heads /. FlexibleSUSY`M -> Identity, {SARAH`Delta, SARAH`ThetaStep}];
-           indexedSymbols = Cases[{expr}, s_[__] /; MemberQ[allHeads, s], Infinity];
+           indexedSymbols = Extract[{expr}, Position[{expr}, s_[__] /; MemberQ[allHeads, s], Infinity]];
            rules = Rule[#, IncreaseIndices[#,num]] & /@ indexedSymbols;
-           decrExpr = expr /. rules;
-           Return[decrExpr]
+           expr /. rules
           ];
 
 DecreaseIndexLiterals[expr_] :=