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

Cleaner printing

Browse Source
This commit is contained in:
Peter Boyle
2025-04-04 16:46:57 -04:00
parent 938c47480f
commit 6fec3c15ca
1 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
tests/debug/Test_general_coarse_pvdagm.cc
View File

@@ -47,20 +47,20 @@ public:
void OpDir (const Field &in, Field &out,int dir,int disp) { assert(0); } void OpDir (const Field &in, Field &out,int dir,int disp) { assert(0); }
void OpDirAll (const Field &in, std::vector<Field> &out){ assert(0); }; void OpDirAll (const Field &in, std::vector<Field> &out){ assert(0); };
void Op (const Field &in, Field &out){ void Op (const Field &in, Field &out){
std::cout << "Op: PVdag M "<<std::endl; // std::cout << "Op: PVdag M "<<std::endl;
Field tmp(in.Grid()); Field tmp(in.Grid());
_Mat.M(in,tmp); _Mat.M(in,tmp);
_PV.Mdag(tmp,out); _PV.Mdag(tmp,out);
} }
void AdjOp (const Field &in, Field &out){ void AdjOp (const Field &in, Field &out){
std::cout << "AdjOp: Mdag PV "<<std::endl; // std::cout << "AdjOp: Mdag PV "<<std::endl;
Field tmp(in.Grid()); Field tmp(in.Grid());
_PV.M(in,tmp); _PV.M(in,tmp);
_Mat.Mdag(tmp,out); _Mat.Mdag(tmp,out);
} }
void HermOpAndNorm(const Field &in, Field &out,RealD &n1,RealD &n2){ assert(0); } void HermOpAndNorm(const Field &in, Field &out,RealD &n1,RealD &n2){ assert(0); }
void HermOp(const Field &in, Field &out){ void HermOp(const Field &in, Field &out){
std::cout << "HermOp: Mdag PV PVdag M"<<std::endl; // std::cout << "HermOp: Mdag PV PVdag M"<<std::endl;
Field tmp(in.Grid()); Field tmp(in.Grid());
// _Mat.M(in,tmp); // _Mat.M(in,tmp);
// _PV.Mdag(tmp,out); // _PV.Mdag(tmp,out);
@@ -83,14 +83,14 @@ public:
void OpDir (const Field &in, Field &out,int dir,int disp) { assert(0); } void OpDir (const Field &in, Field &out,int dir,int disp) { assert(0); }
void OpDirAll (const Field &in, std::vector<Field> &out){ assert(0); }; void OpDirAll (const Field &in, std::vector<Field> &out){ assert(0); };
void Op (const Field &in, Field &out){ void Op (const Field &in, Field &out){
std::cout << "Op: PVdag M "<<std::endl; // std::cout << "Op: PVdag M "<<std::endl;
Field tmp(in.Grid()); Field tmp(in.Grid());
_Mat.M(in,tmp); _Mat.M(in,tmp);
_PV.Mdag(tmp,out); _PV.Mdag(tmp,out);
out = out + shift * in; out = out + shift * in;
} }
void AdjOp (const Field &in, Field &out){ void AdjOp (const Field &in, Field &out){
std::cout << "AdjOp: Mdag PV "<<std::endl; // std::cout << "AdjOp: Mdag PV "<<std::endl;
Field tmp(in.Grid()); Field tmp(in.Grid());
_PV.M(tmp,out); _PV.M(tmp,out);
_Mat.Mdag(in,tmp); _Mat.Mdag(in,tmp);
@@ -98,7 +98,7 @@ public:
} }
void HermOpAndNorm(const Field &in, Field &out,RealD &n1,RealD &n2){ assert(0); } void HermOpAndNorm(const Field &in, Field &out,RealD &n1,RealD &n2){ assert(0); }
void HermOp(const Field &in, Field &out){ void HermOp(const Field &in, Field &out){
std::cout << "HermOp: Mdag PV PVdag M"<<std::endl; // std::cout << "HermOp: Mdag PV PVdag M"<<std::endl;
Field tmp(in.Grid()); Field tmp(in.Grid());
Op(in,tmp); Op(in,tmp);
AdjOp(tmp,out); AdjOp(tmp,out);