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

minor stylistic changes

Browse Source
This commit is contained in:
Alessandro Lupo 2023-04-06 12:08:17 +01:00
parent 4ea29b8f0f
commit 178376f24b
4 changed files with 15 additions and 34 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
Grid/qcd/action/gauge/GaugeImplTypes.h
View File

@ -133,19 +133,6 @@ public:
static inline Field projectForce(Field &P)
{
/* // this works
iScalar<iScalar<iMatrix<Complex, Nc> > > gen;
auto Psum = P;
Psum = Zero();
for (int a = 0; a < Group::AlgebraDimension; a++)
{
Group::generator(a, gen);
auto coeff = 2. * trace(P * gen);
Psum += coeff * gen;
return Psum;}*/
//this doesnt
Field ret(P.Grid());
Group::taProj(P, ret);
return ret;
@ -159,7 +146,7 @@ public:
autoView(P_v,P,AcceleratorRead);
accelerator_for(ss, P.Grid()->oSites(),1,{
for (int mu = 0; mu < Nd; mu++) {
U_v[ss](mu) = Exponentiate(P_v[ss](mu), ep, Nexp) * U_v[ss](mu); //
U_v[ss](mu) = Exponentiate(P_v[ss](mu), ep, Nexp) * U_v[ss](mu);
Group::ProjectOnGaugeGroup(U_v[ss](mu));
}
});
@ -206,9 +193,9 @@ typedef GaugeImplTypes<vComplex, Nc> GimplTypesR;
typedef GaugeImplTypes<vComplexF, Nc> GimplTypesF;
typedef GaugeImplTypes<vComplexD, Nc> GimplTypesD;
typedef GaugeImplTypes<vComplex, Nc, 12, Sp<Nc>> SymplGimplTypesR;
typedef GaugeImplTypes<vComplexF, Nc, 12, Sp<Nc>> SymplGimplTypesF;
typedef GaugeImplTypes<vComplexD, Nc, 12, Sp<Nc>> SymplGimplTypesD;
typedef GaugeImplTypes<vComplex, Nc, 12, Sp<Nc> > SymplGimplTypesR;
typedef GaugeImplTypes<vComplexF, Nc, 12, Sp<Nc> > SymplGimplTypesF;
typedef GaugeImplTypes<vComplexD, Nc, 12, Sp<Nc> > SymplGimplTypesD;
typedef GaugeImplTypes<vComplex, SU<Nc>::AdjointDimension> GimplAdjointTypesR;
typedef GaugeImplTypes<vComplexF, SU<Nc>::AdjointDimension> GimplAdjointTypesF;

10
tests/core/Test_reunitarise.cc
View File

@ -123,14 +123,15 @@ int main (int argc, char ** argv)
std::cout << "Determinant defect before projection " <<norm2(detU)<<std::endl;
tmp = U*adj(U) - ident;
std::cout << "Unitarity check before projection " << norm2(tmp)<<std::endl;
#if Nc==3
#if Nc==3
ProjectSU3(U);
detU= Determinant(U) ;
detU= detU -1.0;
std::cout << "Determinant ProjectSU3 defect " <<norm2(detU)<<std::endl;
tmp = U*adj(U) - ident;
std::cout << "Unitarity check after projection " << norm2(tmp)<<std::endl;
#endif
#endif
ProjectSUn(UU);
detUU= Determinant(UU);
detUU= detUU -1.0;
@ -140,8 +141,3 @@ int main (int argc, char ** argv)
Grid_finalize();
}

12
tests/sp2n/Test_project_on_Sp.cc
View File

@ -31,6 +31,7 @@ int main (int argc, char **argv)
identity = 1.0;
RealD epsilon = 0.01;
RealD Delta = 666.;
Complex i(0., 1.);
RealD u = 0.;
double vol = Umu.Grid()->gSites();
@ -118,7 +119,7 @@ int main (int argc, char **argv)
std::cout << GridLogMessage << "Testing ProjectOnGaugeGroup" << std::endl;
U = U + 2932.111*identity;
U = U + Delta*identity;
std::cout << GridLogMessage << "Apply ProjectOnGaugeGroup to deformed matrix" << std::endl;
Sp<Nc>::ProjectOnGaugeGroup(U);
aux = U*adj(U) - identity;
@ -162,7 +163,7 @@ int main (int argc, char **argv)
}
std::cout << GridLogMessage << "Testing ProjectGn" << std::endl;
U = U + 2932.111*identity;
U = U + Delta*identity;
std::cout << GridLogMessage << "Apply ProjectGn to deformed matrix" << std::endl;
Sp<Nc>::ProjectGn(U);
aux = U*adj(U) - identity;
@ -214,7 +215,7 @@ int main (int argc, char **argv)
std::cout << GridLogMessage << "Testing SpTa" << std::endl;
U = PeekIndex<LorentzIndex>(Umu,1);
U = U + 666.*identity;
U = U + Delta*identity;
std::cout << GridLogMessage << "Matrix deformed " << std::endl;
std::cout << GridLogMessage << "Apply SpTa to deformed matrix" << std::endl;
U = SpTa(U);
@ -279,7 +280,7 @@ int main (int argc, char **argv)
// test taProj
std::cout << GridLogMessage << "Testing taProj" << std::endl;
U = U + 666.*identity;
U = U + Delta*identity;
std::cout << GridLogMessage << "Matrix deformed " << std::endl;
std::cout << GridLogMessage << "Apply taProj to deformed matrix" << std::endl;
Sp<Nc>::taProj(U, Up);
@ -334,9 +335,6 @@ int main (int argc, char **argv)
}
}
Grid_finalize();
}