fix: consistent naming for all Mathematica functions: start with Capi…

…tal letter Signed-off-by: Franz R. Sattler <[email protected]>
satfra · Dec 17, 2024 · 610a5e7 · 610a5e7
1 parent 66f2c2c
commit 610a5e7
Show file tree

Hide file tree

Showing 15 changed files with 1,723 additions and 1,400 deletions.
diff --git a/Examples/FourFermi/Four-Fermion.m b/Examples/FourFermi/Four-Fermion.m
@@ -180,7 +180,7 @@
 ex1=Map[{ReleaseHold[#]}&,ex1];
 Return[ex1]
 ]
-Finalize[expr_]:=simplifyAllMomenta[q,expandScalarProductsFiniteT[expr]]
+Finalize[expr_]:=SimplifyAllMomenta[q,ExpandScalarProductsFiniteT[expr]]
 
 
 kernelLambda1=<|"Path"->"qbqqbq","Name"->"lambda1","Type"->"Quadratureq0","Angles"->0,"d"->4,"AD"->False,"ctype"->"double","Device"->"CPU"|>;

diff --git a/Examples/FourFermi/Four-Fermion.nb b/Examples/FourFermi/Four-Fermion.nb
diff --git a/Examples/ONfiniteT/ON.nb b/Examples/ONfiniteT/ON.nb
diff --git a/Examples/ONfiniteT/flows/V/V.kernel b/Examples/ONfiniteT/flows/V/V.kernel
@@ -4,11 +4,14 @@ public:
   template <typename T1, typename T2, typename T3, typename T4, typename T5>
   static __forceinline__ __host__ __device__ auto kernel(const double q, const double k, const T1 N, const T2 T, const T3 rhoPhi, const T4 m2Pi, const T5 m2Sigma)
   {
+    // optimization definitions
+    const double DiFfRGREPL1 = RB(powr<2>(k), powr<2>(q));
+
     // flow
     return 0.25 *
-           (CothFiniteT((sqrt(m2Sigma + powr<2>(q) + RB(powr<2>(k), powr<2>(q)))), T) * (sqrt(m2Pi + powr<2>(q) + RB(powr<2>(k), powr<2>(q)))) +
-            (-1. + N) * CothFiniteT((sqrt(m2Pi + powr<2>(q) + RB(powr<2>(k), powr<2>(q)))), T) * (sqrt(m2Sigma + powr<2>(q) + RB(powr<2>(k), powr<2>(q))))) *
-           powr<-1>(sqrt((m2Pi + powr<2>(q) + RB(powr<2>(k), powr<2>(q))) * (m2Sigma + powr<2>(q) + RB(powr<2>(k), powr<2>(q))))) * RBdot(powr<2>(k), powr<2>(q));
+           (CothFiniteT((sqrt(DiFfRGREPL1 + m2Sigma + powr<2>(q))), T) * (sqrt(DiFfRGREPL1 + m2Pi + powr<2>(q))) +
+            (-1. + N) * CothFiniteT((sqrt(DiFfRGREPL1 + m2Pi + powr<2>(q))), T) * (sqrt(DiFfRGREPL1 + m2Sigma + powr<2>(q)))) *
+           powr<-1>(sqrt((DiFfRGREPL1 + m2Pi + powr<2>(q)) * (DiFfRGREPL1 + m2Sigma + powr<2>(q)))) * RBdot(powr<2>(k), powr<2>(q));
   }
 
   template <typename T1, typename T2, typename T3, typename T4, typename T5>

diff --git a/Examples/QuarkMesonLPAprime/QuarkMesonLPAprime.m b/Examples/QuarkMesonLPAprime/QuarkMesonLPAprime.m
@@ -129,11 +129,11 @@
 };
 
 sanity=(\[CapitalGamma]\[CapitalPi]\[CapitalPi][{-p,a1,p,a2}]+R\[CapitalPi]\[CapitalPi][{-p,a1,p,a2}])G\[CapitalPi]\[CapitalPi][{p,a2,-p,a1}]//.PreTraceRulesMesonic//FormTrace;
-sanity=sanity//.PostTraceRulesMesonic//expandScalarProductsFiniteT//FullSimplify;
+sanity=sanity//.PostTraceRulesMesonic//ExpandScalarProductsFiniteT//FullSimplify;
 Print[sanity," should be -1+\!\(\*SuperscriptBox[\(Nf\), \(2\)]\)"]
 
 sanity=(\[CapitalGamma]\[Sigma]\[Sigma][{-p,p}]+R\[Sigma]\[Sigma][{-p,p}])G\[Sigma]\[Sigma][{p,-p}]//.PreTraceRulesMesonic//FormTrace;
-sanity=sanity//.PostTraceRulesMesonic//expandScalarProductsFiniteT//FullSimplify;
+sanity=sanity//.PostTraceRulesMesonic//ExpandScalarProductsFiniteT//FullSimplify;
 Print[sanity," should be 1"]
 
 
@@ -306,14 +306,14 @@
 RQdot[__],RPhidot[__]
 };
 QCDSimp[expr_]:=Module[{ret},
-ret=expr//UseLorentzLinearity//expandScalarProductsFiniteT//simplifyAllMomenta[q,#]&;
+ret=expr//UseLorentzLinearity//ExpandScalarProductsFiniteT//SimplifyAllMomenta[q,#]&;
 ret=ret/.Map[(Head[#][a___]:>FullSimplify[Head[#][a]])&,otherList\[Union]dotsList];
 ret=Collect[ret,couplingsList\[Union]dotsList,QuickSimplify];
 ret=Collect[ret,couplingsList\[Union]dotsList,Simplify];
 Return[ret];
 ];
 QCDQuickSimp[expr_]:=Module[{ret},
-ret=expr//UseLorentzLinearity//expandScalarProductsFiniteT//simplifyAllMomenta[q,#]&;
+ret=expr//UseLorentzLinearity//ExpandScalarProductsFiniteT//SimplifyAllMomenta[q,#]&;
 ret=ret/.Map[(Head[#][a___]:>FullSimplify[Head[#][a]])&,otherList\[Union]dotsList];
 ret=Collect[ret,couplingsList\[Union]dotsList,QuickSimplify];
 Return[ret];
@@ -426,7 +426,7 @@
 $Assumptions=$Assumptions&&A0>=0&& gAqbq1>=0&&muq>=0&&q>p>0&&q>0&&k>0&&-1<cos<1&&-d1V-q^2-2 d2V rhoPhi-RB[k^2,q^2]<0&&-d1V-q^2-RB[k^2,q^2]<0&&-m2Quark-q^2-2 q RF[k^2,q^2]-RF[k^2,q^2]^2<0;
 
 ProjectedV=DiagramsV//.PreTraceRules;
-VLoop=expandScalarProductsFiniteT[Map[FormTrace[PreTrace[#]]&,ProjectedV]]//Simplify;
+VLoop=ExpandScalarProductsFiniteT[Map[FormTrace[PreTrace[#]]&,ProjectedV]]//Simplify;
 VLoop=VLoop//.PostTraceRules//.hPhi->Sqrt[(m2Quark 2Nf)/(2rhoPhi)]//.2 d2V rhoPhi->m2Sigma-d1V//QCDSimp;
 VLoop=(Map[MatsubaraSum[#,q0,T]&,VLoop]//Expand)//.Tanh[a__]:>Tanh[FullSimplify[a]]//SafeFiniteTFunctions[#,T]&//QCDSimp;
 

diff --git a/Examples/QuarkMesonLPAprime/QuarkMesonLPAprime.nb b/Examples/QuarkMesonLPAprime/QuarkMesonLPAprime.nb
diff --git a/Examples/YangMills/Full/Yang-Mills.m b/Examples/YangMills/Full/Yang-Mills.m
@@ -231,14 +231,14 @@
 dtZA[__],ZA[__],etac[__],RA[__],Rc[__],RAdot[__],Rcdot[__]
 };
 QCDSimp[expr_]:=Module[{ret},
-ret=expr//.Nc->3//.qf->q//UseLorentzLinearity//expandScalarProducts//simplifyAllMomenta[q,#]&;
+ret=expr//.Nc->3//.qf->q//UseLorentzLinearity//ExpandScalarProducts//SimplifyAllMomenta[q,#]&;
 ret=ret/.Map[(Head[#][a___]:>FullSimplify[Head[#][a]])&,couplingsList\[Union]otherList];
 ret=Collect[ret,couplingsList,QuickSimplify];
 ret=Collect[ret,couplingsList,Simplify];
 Return[ret]
 ];
 QCDQuickSimp[expr_]:=Module[{ret},
-ret=expr//.Nc->3//.qf->q//UseLorentzLinearity//expandScalarProducts//simplifyAllMomenta[q,#]&;
+ret=expr//.Nc->3//.qf->q//UseLorentzLinearity//ExpandScalarProducts//SimplifyAllMomenta[q,#]&;
 ret=ret/.Map[(Head[#][a___]:>FullSimplify[Head[#][a]])&,couplingsList\[Union]otherList];
 ret=Collect[ret,couplingsList,QuickSimplify];
 Return[ret]
@@ -335,11 +335,11 @@
 ProjectorZA=TBGetProjector["AA",1,{-p,v1,a1},{p,v2,a2}];
 
 (*Sanity Check*)
-sanity=(FormTrace[ProjectorZA ((\[CapitalGamma]AA[{-p,v2,a2,p,v1,a1}]-\[CapitalGamma]AA[{-q,v2,a2,q,v1,a1}])/sp[p,p])//.PreTraceRules]//.PostTraceRules//expandScalarProducts)//.{cos[p,q]->1,q->0}//Simplify;
+sanity=(FormTrace[ProjectorZA ((\[CapitalGamma]AA[{-p,v2,a2,p,v1,a1}]-\[CapitalGamma]AA[{-q,v2,a2,q,v1,a1}])/sp[p,p])//.PreTraceRules]//.PostTraceRules//ExpandScalarProducts)//.{cos[p,q]->1,q->0}//Simplify;
 Print["Projection check is ", sanity, ", should be ZA[p]"]
 
 (*Sanity Check*)
-sanity=simplifyAllMomenta[q,FormTrace[ProjectorZA \[CapitalGamma]AA[{-p,v2,a2,p,v1,a1}]//.PreTraceRules]//.PostTraceRules//expandScalarProducts//Simplify]//.p->0;
+sanity=SimplifyAllMomenta[q,FormTrace[ProjectorZA \[CapitalGamma]AA[{-p,v2,a2,p,v1,a1}]//.PreTraceRules]//.PostTraceRules//ExpandScalarProducts//Simplify]//.p->0;
 Print["Projection check is ", sanity, ", should be m2A"]
 
 
@@ -425,7 +425,7 @@
 (*Projection*)
 ProjectorZAcbc=TBGetProjector["Acbc",1,{p1,v1,a1},{-p2-p1,a3},{p2,a2}]//UseLorentzLinearity//FullSimplify;
 (*Sanity check*)
-sanity=Simplify[ProjectorZAcbc \[CapitalGamma]Acbc[{p1,v1,a1,p2,a2,-p1-p2,a3}]//.PreTraceRules//FormTrace]//.PostTraceRules//transf3PTo//expandScalarProducts//FullSimplify;
+sanity=Simplify[ProjectorZAcbc \[CapitalGamma]Acbc[{p1,v1,a1,p2,a2,-p1-p2,a3}]//.PreTraceRules//FormTrace]//.PostTraceRules//transf3PTo//ExpandScalarProducts//FullSimplify;
 Print["Projection check is ", sanity, ", should be ZAcbc[S0,S1,SPhi]"]