portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-09 23:45:36 +00:00
Code Releases Activity

first attempt at U3 projection

Browse Source
This commit is contained in:
david clarke 2024-01-22 02:49:40 -07:00
parent 981c93d67a
commit f5b3d582b0
2 changed files with 42 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

46
Grid/qcd/smearing/HISQSmearing.h
View File

@ -83,7 +83,6 @@ struct SmearingParameters{
/*! @brief create fat links from link variables */
template<class Gimpl>
class Smear_HISQ : public Gimpl {
// TODO: this needs to be renamed, becaues the Naik guy is not part of the fat smear
private:
GridCartesian* const _grid;
@ -289,16 +288,16 @@ public:
std::vector<LF> V(Nd, u_smr.Grid());
std::vector<LF> Vnaik(Nd, u_naik.Grid());
for (int mu = 0; mu < Nd; mu++) {
U[mu] = PeekIndex<LorentzIndex>(u_thin, mu);
V[mu] = PeekIndex<LorentzIndex>(u_smr, mu);
U[mu] = PeekIndex<LorentzIndex>(u_thin, mu);
V[mu] = PeekIndex<LorentzIndex>(u_smr, mu);
}
for(int mu=0;mu<Nd;mu++) {
// Naik
Vnaik[mu] = lt.c_naik*Gimpl::CovShiftForward(U[mu],mu,
Gimpl::CovShiftForward(U[mu],mu,
Gimpl::CovShiftIdentityForward(U[mu],mu)));
Vnaik[mu] = lt.c_naik*Gimpl::CovShiftForward(U[mu],mu,
Gimpl::CovShiftForward(U[mu],mu,
Gimpl::CovShiftIdentityForward(U[mu],mu)));
// LePage
for (int nu_h=1;nu_h<Nd;nu_h++) {
@ -325,6 +324,41 @@ public:
}
};
// Intent: OUT--u_proj
// IN--u_mu
void projectU3(GF& u_proj, GF& u_mu) const {
LF V, Q, sqrtQinv;
Real c1, c2, g0, g1, g2, S, R, theta, u, v, w, den, f0, f1, f2;
// Follow MILC 10.1103/PhysRevD.82.074501, eqs (B2-B3) and (C1-C8)
for (int mu = 0; mu < Nd; mu++) {
V = PeekIndex<LorentzIndex>(u_mu, mu);
Q = adj(V[mu])*V[mu];
c1 = trace(Q*Q)/2.; // SU(N) matrices are traceless, so c0=0.
c2 = trace(Q*Q*Q)/3.;
S = c1/3.;
R = c2/2.;
theta = std::acos(R/std::pow(S,1.5));
g0 = 2.*std::sqrt(S)*std::cos(theta/3.-2*M_PI/3.);
g1 = 2.*std::sqrt(S)*std::cos(theta/3. );
g2 = 2.*std::sqrt(S)*std::cos(theta/3.+2*M_PI/3.);
// if (fabs(Q.determinant()/(g0*g1*g2)-1.0) > 1e-5) {}
u = std::sqrt(g0) + std::sqrt(g1) + std::sqrt(g2);
v = std::sqrt(g0*g1) + std::sqrt(g0*g2) + std::sqrt(g1*g2);
w = std::sqrt(g0*g1*g2);
den = w*(u*v-w);
f0 = (-w*(u*u+v)+u*v*v)/den;
f1 = (-w-u*u*u+2*u*v)/den;
f2 = u/den;
sqrtQinv = f0 + f1*Q + f2*Q*Q;
PokeIndex<LorentzIndex>(u_proj, V*sqrtQinv, mu);
}
};
// void derivative(const GaugeField& Gauge) const {
// };
};

3
tests/smearing/Test_fatLinks.cc
View File

@ -61,7 +61,7 @@ struct ConfParameters: Serializable {
void testSmear(GridCartesian& GRID, LatticeGaugeFieldD Umu, LatticeGaugeFieldD Usmr, LatticeGaugeFieldD Unaik,
LatticeGaugeFieldD Ucontrol, Real c1, Real cnaik, Real c3, Real c5, Real c7, Real clp) {
Smear_HISQ<PeriodicGimplD> hisq_fat(&GRID,c1,cnaik,c3,c5,c7,clp);
LatticeGaugeFieldD diff(&GRID);
LatticeGaugeFieldD diff(&GRID), Uproj(&GRID);
hisq_fat.smear(Usmr, Unaik, Umu);
if (cnaik < 1e-30) { // Testing anything but Naik term
diff = Ucontrol-Usmr;
@ -79,6 +79,7 @@ void testSmear(GridCartesian& GRID, LatticeGaugeFieldD Umu, LatticeGaugeFieldD U
} else {
Grid_error(" |Umu-Unaik|/|Umu| = ",absDiff);
}
hisq_fat.projectU3(Uproj,Usmr);
// NerscIO::writeConfiguration(Unaik,"nersc.l8t4b3360.naik");
}
}