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

Time stamp

Browse Source
This commit is contained in:
Peter Boyle 2023-11-24 12:56:45 -05:00
parent 3e448435d3
commit 59abaeb5cd
1 changed files with 15 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
tests/debug/Test_general_coarse.cc
View File

@ -247,17 +247,31 @@ int main (int argc, char ** argv)
CoarseVector rh_phi(CoarseMrhs); CoarseVector rh_phi(CoarseMrhs);
CoarseVector rh_res(CoarseMrhs); CoarseVector rh_res(CoarseMrhs);
random(rh_CRNG,rh_phi); random(rh_CRNG,rh_phi);
mrhs.M(rh_phi,rh_res);
mrhs.M(rh_phi,rh_res);
const int ncall=100;
RealD t0=-usecond();
for(int i=0;i<ncall;i++){
mrhs.M(rh_phi,rh_res);
}
t0+=usecond();
RealD t1=0;
for(int r=0;r<nrhs;r++){ for(int r=0;r<nrhs;r++){
ExtractSlice(phi,rh_phi,r,0); ExtractSlice(phi,rh_phi,r,0);
ExtractSlice(chi,rh_res,r,0); ExtractSlice(chi,rh_res,r,0);
LittleDiracOp.M(phi,Aphi); LittleDiracOp.M(phi,Aphi);
t1-=usecond();
for(int i=0;i<ncall;i++){
LittleDiracOp.M(phi,Aphi);
}
t1+=usecond();
std::cout << r << " mrhs " << norm2(chi)<<std::endl; std::cout << r << " mrhs " << norm2(chi)<<std::endl;
std::cout << r << " srhs " << norm2(Aphi)<<std::endl; std::cout << r << " srhs " << norm2(Aphi)<<std::endl;
chi=chi-Aphi; chi=chi-Aphi;
std::cout << r << " diff " << norm2(chi)<<std::endl; std::cout << r << " diff " << norm2(chi)<<std::endl;
} }
std::cout << nrhs<< " mrhs " << t0/ncall/nrhs <<" us"<<std::endl;
std::cout << nrhs<< " srhs " << t1/ncall/nrhs <<" us"<<std::endl;
} }
Grid_finalize(); Grid_finalize();
return 0; return 0;