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

Fixing issue #11: consistent use of ncolour and nsp

Browse Source
This commit is contained in:
Alessandro Lupo 2022-11-23 18:45:01 +00:00
parent 88bdd4344b
commit 1ba429345b
4 changed files with 76 additions and 72 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -126,8 +126,8 @@ public:
{
if (isSp2n == true)
{
const int nSp = Nrepresentation/2;
Sp<nSp>::GaussianFundamentalLieAlgebraMatrix(pRNG, Pmu);
//const int nSp = Nrepresentation/2;
Sp<Nrepresentation>::GaussianFundamentalLieAlgebraMatrix(pRNG, Pmu);
} else
{
@ -148,18 +148,18 @@ public:
if (isSp2n == true)
{
P = Ta(P);
const int nsp = Nc / 2;
//const int nsp = Nc / 2;
Sp<nsp>::iSp2nMatrix<Complex> gen;
Sp<Nc>::iSp2nMatrix<Complex> gen;
auto Psum = P;
Psum = Zero();
for (int a = 0; a < Sp<nsp>::AlgebraDimension; a++)
for (int a = 0; a < Sp<Nrepresentation>::AlgebraDimension; a++)
{
Sp<nsp>::generator(a, gen);
Sp<Nrepresentation>::generator(a, gen);
auto coeff = 2. * trace(P * gen);
Psum += coeff * gen;
@ -223,8 +223,8 @@ public:
{
if (isSp2n == true)
{
const int nSp = Nrepresentation/2;
Sp<nSp>::HotConfiguration(pRNG, U);
//const int nSp = Nrepresentation/2;
Sp<Nrepresentation>::HotConfiguration(pRNG, U);
} else
{
Group::HotConfiguration(pRNG, U);
@ -235,8 +235,8 @@ public:
{
if (isSp2n == true)
{
const int nSp = Nrepresentation/2;
Sp<nSp>::TepidConfiguration(pRNG, U);
//const int nSp = Nrepresentation/2;
Sp<Nrepresentation>::TepidConfiguration(pRNG, U);
} else
{
Group::TepidConfiguration(pRNG, U);
@ -248,8 +248,8 @@ public:
{
if (isSp2n == true)
{
const int nSp = Nrepresentation/2;
Sp<nSp>::ColdConfiguration(pRNG, U);
//const int nSp = Nrepresentation/2;
Sp<Nrepresentation>::ColdConfiguration(pRNG, U);
} else
{
Group::ColdConfiguration(pRNG, U);

99
Grid/qcd/utils/Sp2n.h
View File

@ -5,21 +5,16 @@
NAMESPACE_BEGIN(Grid);
// Sp(2N)
// ncolour = N
// ncolour = 2N
// need to be careful with n and 2n
// I am defining the Sp class for Sp(2n) to be such that the template variable ncolour
// is n inside 2n and the typedef at the end of the file should eliminate possible confusion.
// the other routines, like projectOnSp2n, N will be the dimension of the actual number of colors for consistency with the sun routines
template <int ncolour>
class Sp {
public:
static const int Dimension = ncolour*2;
static const int AlgebraDimension = ncolour*(2*ncolour +1);
static int su2subgroups(void) { return (ncolour * (ncolour - 1)) / 2; }
static const int nnsp = ncolour;
static const int nsp = ncolour/2;
static const int Dimension = ncolour;
static const int AlgebraDimension = nsp*(2*nsp +1);
static int su2subgroups(void) { return (nsp * (nsp - 1)) / 2; }
template <typename vtype>
@ -65,20 +60,17 @@ public:
// there are 6 types named a,b,c,d and w,z
// abcd are N(N-1)/2 each while wz are N each
template <class cplx>
template <class cplx>
static void generator(int lieIndex, iSp2nMatrix<cplx> &ta) {
// map lie index into type of generators: diagonal, abcd type, wz type
int diagIndex;
int aIndex, bIndex, cIndex, dIndex;
int wIndex, zIndex; // a,b,c,d are N(N-1)/2 and w,z are N
int mod = ncolour * (ncolour-1) * 0.5;
int offdiag = 2*ncolour*ncolour; // number of generators not in the cartan subalgebra
int mod = nsp * (nsp-1) * 0.5;
int offdiag = 2*nsp*nsp; // number of generators not in the cartan subalgebra
int wmod = 4*mod;
int zmod = wmod+ncolour;
int zmod = wmod+nsp;
if (lieIndex >= offdiag) {
diagIndex = lieIndex - offdiag; // 0, ... ,N-1
//std::cout << GridLogMessage << "diag type " << std::endl;
@ -136,7 +128,7 @@ public:
RealD nrm = 1.0 / 2;
ta()()(diagIndex,diagIndex) = nrm;
ta()()(diagIndex+ncolour,diagIndex+ncolour) = -nrm;
ta()()(diagIndex+nsp,diagIndex+nsp) = -nrm;
}
template <class cplx>
@ -152,8 +144,8 @@ public:
su2SubGroupIndex(i1, i2, aIndex);
ta()()(i1,i2) = 1;
ta()()(i2,i1) = 1;
ta()()(i1+ncolour,i2+ncolour) = -1;
ta()()(i2+ncolour,i1+ncolour) = -1;
ta()()(i1+nsp,i2+nsp) = -1;
ta()()(i2+nsp,i1+nsp) = -1;
ta = ta * nrm;
}
@ -174,8 +166,8 @@ public:
ta()()(i1,i2) = i;
ta()()(i2,i1) = -i;
ta()()(i1+ncolour,i2+ncolour) = i;
ta()()(i2+ncolour,i1+ncolour) = -i;
ta()()(i1+nsp,i2+nsp) = i;
ta()()(i2+nsp,i1+nsp) = -i;
ta = ta * nrm;
}
@ -191,10 +183,10 @@ public:
RealD nrm = 1 / (2 * std::sqrt(2) );
su2SubGroupIndex(i1, i2, cIndex);
ta()()(i1,i2+ncolour) = 1;
ta()()(i2,i1+ncolour) = 1;
ta()()(i1+ncolour,i2) = 1;
ta()()(i2+ncolour,i1) = 1;
ta()()(i1,i2+nsp) = 1;
ta()()(i2,i1+nsp) = 1;
ta()()(i1+nsp,i2) = 1;
ta()()(i2+nsp,i1) = 1;
ta = ta * nrm;
}
@ -210,10 +202,10 @@ public:
RealD nrm = 1 / (2 * std::sqrt(2) );
su2SubGroupIndex(i1, i2, dIndex);
ta()()(i1,i2+ncolour) = i;
ta()()(i2,i1+ncolour) = i;
ta()()(i1+ncolour,i2) = -i;
ta()()(i2+ncolour,i1) = -i;
ta()()(i1,i2+nsp) = i;
ta()()(i2,i1+nsp) = i;
ta()()(i1+nsp,i2) = -i;
ta()()(i2+nsp,i1) = -i;
ta = ta * nrm;
}
@ -226,8 +218,8 @@ public:
ta = Zero();
RealD nrm = 1.0 / 2; //check
ta()()(wIndex,wIndex+ncolour) = 1;
ta()()(wIndex+ncolour,wIndex) = 1;
ta()()(wIndex,wIndex+nsp) = 1;
ta()()(wIndex+nsp,wIndex) = 1;
ta = ta * nrm;
}
@ -240,8 +232,8 @@ public:
ta = Zero();
RealD nrm = 1.0 / 2; //check
cplx i(0.0, 1.0);
ta()()(zIndex,zIndex+ncolour) = i;
ta()()(zIndex+ncolour,zIndex) = -i;
ta()()(zIndex,zIndex+nsp) = i;
ta()()(zIndex+nsp,zIndex) = -i;
ta = ta * nrm;
}
@ -251,11 +243,11 @@ public:
// Map a su2 subgroup number to the pair of rows that are non zero
////////////////////////////////////////////////////////////////////////
static void su2SubGroupIndex(int &i1, int &i2, int su2_index) {
assert((su2_index >= 0) && (su2_index < (ncolour * (ncolour - 1)) / 2));
assert((su2_index >= 0) && (su2_index < (nsp * (nsp - 1)) / 2));
int spare = su2_index;
for (i1 = 0; spare >= (ncolour - 1 - i1); i1++) {
spare = spare - (ncolour - 1 - i1); // remove the Nc-1-i1 terms
for (i1 = 0; spare >= (nsp - 1 - i1); i1++) {
spare = spare - (nsp - 1 - i1); // remove the Nc-1-i1 terms
}
i2 = i1 + 1 + spare;
}
@ -268,7 +260,7 @@ public:
for (int gen = 0; gen < AlgebraDimension; gen++) {
Matrix ta;
generator(gen, ta);
std::cout << GridLogMessage << "Nc (2n) = " << 2*ncolour << std::endl;
std::cout << GridLogMessage << "Nc = " << ncolour << std::endl;
std::cout << GridLogMessage << " t_" << gen << std::endl;
std::cout << GridLogMessage << ta << std::endl;
}
@ -464,10 +456,17 @@ public:
ColourMatrix Omega;
Omega = Zero();
for (int i = 0; i < ncolour; i++)
std::cout << GridLogMessage << "I am a ColourMatrix" << std::endl;
//for (int i = 0; i < ncolour; i++) wrong?!
//{
// Omega()()(i, 2*ncolour-1-i) = 1.;
// Omega()()(2*ncolour-1-i, i) = -1;
//}
for (int i = 0; i < nsp; i++)
{
Omega()()(i, 2*ncolour-1-i) = 1.;
Omega()()(2*ncolour-1-i, i) = -1;
Omega()()(i, nsp+i) = 1.;
Omega()()(nsp+i, i) = -1;
}
auto diff = Omega - (in * Omega * transpose(in) );
@ -498,6 +497,8 @@ public:
Omega = Zero();
identity = 1.;
std::cout << GridLogMessage << "I am a GaugeField " << std::endl;
U = PeekIndex<LorentzIndex>(in,1);
OmegaInvariance(U);
@ -516,10 +517,12 @@ public:
Omega = Zero();
identity = 1.;
for (int i = 0; i < ncolour; i++)
std::cout << GridLogMessage << "I am a LatticeColourMatrix " << std::endl;
for (int i = 0; i < nsp; i++)
{
Omega()()(i, ncolour+i) = 1.;
Omega()()(ncolour+i, i) = -1;
Omega()()(i, nsp+i) = 1.;
Omega()()(nsp+i, i) = -1;
}
std::cout << GridLogMessage << "Omega = " << Omega()() << std::endl;
@ -568,10 +571,10 @@ public:
}
}
typedef Sp<1> Sp2;
typedef Sp<2> Sp4;
typedef Sp<3> Sp6;
typedef Sp<4> Sp8;
typedef Sp<2> Sp2;
typedef Sp<4> Sp4;
typedef Sp<6> Sp6;
typedef Sp<8> Sp8;
NAMESPACE_END(Grid);
#endif

18
tests/sp2n/Test_Sp_start.cc
View File

@ -21,7 +21,7 @@ int main (int argc, char **argv)
double vol = Umu.Grid()->gSites();
const int nsp = Nc / 2;
const int nsp = Sp<Nc>::nsp;
identity = 1.;
Cidentity = 1.;
@ -34,8 +34,8 @@ int main (int argc, char **argv)
std::cout << GridLogMessage << "-------" << std::endl;
std::cout << GridLogMessage << "Checking Cold Configuration " << std::endl;
std::cout << GridLogMessage << std::endl;
Sp<nsp>::ColdConfiguration(pRNG,Umu);
Sp<nsp>::OmegaInvariance(Umu);
Sp<Nc>::ColdConfiguration(pRNG,Umu);
Sp<Nc>::OmegaInvariance(Umu);
std::cout << GridLogMessage << std::endl;
Umu = Umu - identity;
@ -47,8 +47,8 @@ int main (int argc, char **argv)
std::cout << GridLogMessage << "-------" << std::endl;
std::cout << GridLogMessage << "Checking Hot Configuration" << std::endl;
std::cout << GridLogMessage << std::endl;
Sp<nsp>::HotConfiguration(pRNG,Umu);
Sp<nsp>::OmegaInvariance(Umu);
Sp<Nc>::HotConfiguration(pRNG,Umu);
Sp<Nc>::OmegaInvariance(Umu);
std::cout << GridLogMessage << std::endl;
U = PeekIndex<LorentzIndex>(Umu,1);
std::cout << GridLogMessage << "Checking unitarity " << std::endl;
@ -62,8 +62,8 @@ int main (int argc, char **argv)
std::cout << GridLogMessage << "-------" << std::endl;
std::cout << GridLogMessage << "Checking Tepid Configuration" << std::endl;
std::cout << GridLogMessage << std::endl;
Sp<nsp>::TepidConfiguration(pRNG,Umu);
Sp<nsp>::OmegaInvariance(Umu);
Sp<Nc>::TepidConfiguration(pRNG,Umu);
Sp<Nc>::OmegaInvariance(Umu);
std::cout << GridLogMessage << std::endl;
U = PeekIndex<LorentzIndex>(Umu,1);
std::cout << GridLogMessage << "Checking unitarity " << std::endl;
@ -85,7 +85,7 @@ int main (int argc, char **argv)
Complex i(0., 1.);
Sp<nsp>::HotConfiguration(pRNG,Umu);
Sp<Nc>::HotConfiguration(pRNG,Umu);
U = PeekIndex<LorentzIndex>(Umu,0);
for (int c1 = 0; c1 < nsp; c1++) //check on W
{
@ -123,7 +123,7 @@ int main (int argc, char **argv)
std::cout << GridLogMessage << "Hot start ok " << std::endl;
Sp<nsp>::HotConfiguration(pRNG,Umu);
Sp<Nc>::HotConfiguration(pRNG,Umu);
U = PeekIndex<LorentzIndex>(Umu,0);
for (int c1 = 0; c1 < nsp; c1++) //check on W
{

7
tests/sp2n/Test_project_on_Sp.cc
View File

@ -20,7 +20,8 @@ int main (int argc, char **argv)
LatticeColourMatrixD aux(&Grid);
LatticeColourMatrixD identity(&Grid);
const int nsp = Nc / 2;
//const int nsp = Nc / 2;
const int nsp = Sp<Nc>::nsp;
identity = 1.0;
RealD epsilon = 0.01;
@ -47,7 +48,7 @@ int main (int argc, char **argv)
{
std::cout << GridLogMessage << "This matrix should not leave Omega invariant, expect a warning" << std::endl;
}
Sp<nsp>::OmegaInvariance(U);
Sp<Nc>::OmegaInvariance(U);
std::cout <<GridLogMessage << std::endl;
U = U + epsilon*identity;
@ -78,7 +79,7 @@ int main (int argc, char **argv)
// actual sp2n check
std::cout << GridLogMessage << "Checking invariance after projection "<< std::endl;
Sp<nsp>::OmegaInvariance(U);
Sp<Nc>::OmegaInvariance(U);
// checks on elements