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

Multiply TVs in Wilson MG with G5 instead of G5R5

Browse Source
This commit is contained in:
Daniel Richtmann 2017-12-22 13:07:56 +01:00
parent df152648d6
commit 39558cce52
No known key found for this signature in database
GPG Key ID: B33C490AF0772057
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
tests/solver/Test_wilson_ddalphaamg.cc
View File

@ -820,7 +820,7 @@ int main (int argc, char ** argv)
testVectorAnalyzer(HermOp, Aggregates.subspace); testVectorAnalyzer(HermOp, Aggregates.subspace);
for(int n=0;n<nb;n++){ for(int n=0;n<nb;n++){
Aggregates.subspace[n+nb] = g5 * Aggregates.subspace[n]; // multiply with g5 normally instead of G5R5 Aggregates.subspace[n+nb] = g5 * Aggregates.subspace[n]; // multiply with g5 normally instead of G5R5 since this specific to DWF
std::cout<<GridLogMessage<<n<<" subspace "<<norm2(Aggregates.subspace[n+nb])<<" "<<norm2(Aggregates.subspace[n]) <<std::endl; std::cout<<GridLogMessage<<n<<" subspace "<<norm2(Aggregates.subspace[n+nb])<<" "<<norm2(Aggregates.subspace[n]) <<std::endl;
} }
for(int n=0;n<nbasis;n++){ for(int n=0;n<nbasis;n++){