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

Debug

Browse Source
This commit is contained in:
Peter Boyle 2023-05-11 12:22:24 -04:00
parent ccd21f96ff
commit 1b8a834beb
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
tests/Test_general_stencil.cc
View File

@ -115,6 +115,7 @@ int main(int argc, char ** argv)
if (SE->_permute & 0x2 ) { permute(check[i],tmp,1); tmp=check[i];} if (SE->_permute & 0x2 ) { permute(check[i],tmp,1); tmp=check[i];}
if (SE->_permute & 0x4 ) { permute(check[i],tmp,2); tmp=check[i];} if (SE->_permute & 0x4 ) { permute(check[i],tmp,2); tmp=check[i];}
if (SE->_permute & 0x8 ) { permute(check[i],tmp,3); tmp=check[i];} if (SE->_permute & 0x8 ) { permute(check[i],tmp,3); tmp=check[i];}
// std::cout<<GridLogMessage<<"stencil["<<i<<"] "<< check[i]<< " perm "<<(uint32_t)SE->_permute <<std::endl;
} }
Real nrmC = norm2(Check); Real nrmC = norm2(Check);
@ -138,18 +139,17 @@ int main(int argc, char ** argv)
ddiff = check -bar; ddiff = check -bar;
diff =norm2(ddiff); diff =norm2(ddiff);
if ( diff > 0){ if ( diff > 0){
std::cout <<"Coor (" << coor[0]<<","<<coor[1]<<","<<coor[2]<<","<<coor[3] std::cout <<"Diff at Coor (" << coor[0]<<","<<coor[1]<<","<<coor[2]<<","<<coor[3]
<<") " <<check<<" vs "<<bar<<std::endl; <<") stencil " <<check<<" vs cshift "<<bar<<std::endl;
} }
}}}} }}}}
if (nrm > 1.0e-4) { if (nrm > 1.0e-4) {
autoView( check , Check, CpuRead); autoView( check , Check, CpuRead);
autoView( bar , Bar, CpuRead); autoView( bar , Bar, CpuRead);
for(int i=0;i<check.size();i++){ for(int i=0;i<check.size();i++){
std::cout << i<<" Check "<<check[i]<< "\n"<<i<<" Bar "<<bar[i]<<std::endl; std::cout << i<<" ERROR Check \n"<<check[i]<< "\n"<<i<<" Bar \n"<<bar[i]<<std::endl;
} }
} }
if (nrm > 1.0e-4) exit(-1); if (nrm > 1.0e-4) exit(-1);