portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-11-04 14:04:32 +00:00
Code Releases Activity

ClangFormat

Browse Source
This commit is contained in:
Julian Lenz
2023-05-03 16:15:12 +01:00
parent 6a7bdca53b
commit 6e62f4f616
1 changed files with 25 additions and 28 deletions
Download Patch File Download Diff File Expand all files Collapse all files

53
tests/sp2n/Test_project_on_Sp.cc
View File

@@ -111,36 +111,33 @@ bool has_correct_algebra_block_structure(const T& U) {
const int nsp = Nc / 2; const int nsp = Nc / 2;
Complex i(0., 1.); Complex i(0., 1.);
std::cout << GridLogMessage << "Checking the structure is " << std::endl; std::cout << GridLogMessage << "Checking the structure is " << std::endl;
std::cout << GridLogMessage << "U = ( W X ) " << std::endl; std::cout << GridLogMessage << "U = ( W X ) " << std::endl;
std::cout << GridLogMessage << " ( -X^* W^* ) " << std::endl; std::cout << GridLogMessage << " ( -X^* W^* ) " << std::endl;
for (int c1 = 0; c1 < nsp; c1++) //check on W for (int c1 = 0; c1 < nsp; c1++) // check on W
{ {
for (int c2 = 0; c2 < nsp; c2++) for (int c2 = 0; c2 < nsp; c2++) {
{ auto W = PeekIndex<ColourIndex>(U, c1, c2);
auto W = PeekIndex<ColourIndex>(U,c1,c2); auto Wstar = PeekIndex<ColourIndex>(U, c1 + nsp, c2 + nsp);
auto Wstar = PeekIndex<ColourIndex>(U,c1+nsp,c2+nsp); auto Ww = conjugate(Wstar);
auto Ww = conjugate( Wstar ); auto amizero = sum(W - Ww);
auto amizero = sum(W - Ww); auto amizeroo = TensorRemove(amizero);
auto amizeroo = TensorRemove(amizero); assert(amizeroo.real() < 10e-6);
assert( amizeroo.real() < 10e-6 ); amizeroo *= i;
amizeroo *= i; assert(amizeroo.real() < 10e-6);
assert( amizeroo.real() < 10e-6 );
}
} }
for (int c1 = 0; c1 < nsp ; c1++) }
{ for (int c1 = 0; c1 < nsp; c1++) {
for (int c2 = 0; c2 < nsp; c2++) for (int c2 = 0; c2 < nsp; c2++) {
{ auto X = PeekIndex<ColourIndex>(U, c1, c2 + nsp);
auto X = PeekIndex<ColourIndex>(U,c1,c2+nsp); auto minusXstar = PeekIndex<ColourIndex>(U, c1 + nsp, c2);
auto minusXstar = PeekIndex<ColourIndex>(U,c1+nsp,c2); auto minusXx = conjugate(minusXstar);
auto minusXx = conjugate(minusXstar); auto amizero = sum(X + minusXx);
auto amizero = sum (X + minusXx); auto amizeroo = TensorRemove(amizero);
auto amizeroo = TensorRemove(amizero); assert(amizeroo.real() < 10e-6);
assert( amizeroo.real() < 10e-6 ); amizeroo *= i;
amizeroo *= i; assert(amizeroo.real() < 10e-6);
assert( amizeroo.real() < 10e-6 );
}
} }
}
return true; return true;
} }