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

Improved FlightRecorder

Browse Source
This commit is contained in:
Peter Boyle 2024-03-22 15:42:32 +00:00
parent a477c25e8c
commit d01e5fa838
1 changed files with 20 additions and 31 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

51
tests/Test_dwf_mixedcg_prec.cc
View File

@ -34,6 +34,7 @@ using namespace Grid;
#define HOST_NAME_MAX _POSIX_HOST_NAME_MAX #define HOST_NAME_MAX _POSIX_HOST_NAME_MAX
#endif #endif
NAMESPACE_BEGIN(Grid); NAMESPACE_BEGIN(Grid);
template<class Matrix,class Field> template<class Matrix,class Field>
class SchurDiagMooeeOperatorParanoid : public SchurOperatorBase<Field> { class SchurDiagMooeeOperatorParanoid : public SchurOperatorBase<Field> {
@ -143,14 +144,21 @@ int main (int argc, char ** argv)
time_t start = time(NULL); time_t start = time(NULL);
uint32_t csum, csumref; FlightRecorder::ContinueOnFail = 0;
csumref=0; FlightRecorder::PrintEntireLog = 0;
FlightRecorder::ChecksumComms = 1;
FlightRecorder::ChecksumCommsSend=0;
if(char *s=getenv("GRID_PRINT_ENTIRE_LOG")) FlightRecorder::PrintEntireLog = atoi(s);
if(char *s=getenv("GRID_CHECKSUM_RECV_BUF")) FlightRecorder::ChecksumComms = atoi(s);
if(char *s=getenv("GRID_CHECKSUM_SEND_BUF")) FlightRecorder::ChecksumCommsSend = atoi(s);
int iter=0; int iter=0;
do { do {
if ( iter == 0 ) { if ( iter == 0 ) {
SetGridNormLoggingMode(GridNormLoggingModeRecord); FlightRecorder::SetLoggingMode(FlightRecorder::LoggingModeRecord);
} else { } else {
SetGridNormLoggingMode(GridNormLoggingModeVerify); FlightRecorder::SetLoggingMode(FlightRecorder::LoggingModeVerify);
} }
std::cerr << "******************* SINGLE PRECISION SOLVE "<<iter<<std::endl; std::cerr << "******************* SINGLE PRECISION SOLVE "<<iter<<std::endl;
result_o = Zero(); result_o = Zero();
@ -162,31 +170,22 @@ int main (int argc, char ** argv)
flops+= CGsiteflops*FrbGrid->gSites()*iters; flops+= CGsiteflops*FrbGrid->gSites()*iters;
std::cout << " SinglePrecision iterations/sec "<< iters/(t2-t1)*1000.*1000.<<std::endl; std::cout << " SinglePrecision iterations/sec "<< iters/(t2-t1)*1000.*1000.<<std::endl;
std::cout << " SinglePrecision GF/s "<< flops/(t2-t1)/1000.<<std::endl; std::cout << " SinglePrecision GF/s "<< flops/(t2-t1)/1000.<<std::endl;
std::cout << " SinglePrecision error count "<< FlightRecorder::ErrorCount()<<std::endl;
csum = crc(result_o); assert(FlightRecorder::ErrorCount()==0);
if ( csumref == 0 ) { std::cout << " FlightRecorder is OK! "<<std::endl;
csumref = csum;
} else {
if ( csum != csumref ) {
std::cerr << host<<" FAILURE " <<iter <<" csum "<<std::hex<<csum<< " != "<<csumref <<std::dec<<std::endl;
assert(0);
} else {
std::cout << host <<" OK " <<iter <<" csum "<<std::hex<<csum<<std::dec<<" -- OK! "<<std::endl;
}
}
iter ++; iter ++;
} while (time(NULL) < (start + nsecs/2) ); } while (time(NULL) < (start + nsecs/2) );
std::cout << GridLogMessage << "::::::::::::: Starting double precision CG" << std::endl; std::cout << GridLogMessage << "::::::::::::: Starting double precision CG" << std::endl;
ConjugateGradient<LatticeFermionD> CG(1.0e-8,10000); ConjugateGradient<LatticeFermionD> CG(1.0e-8,10000);
csumref=0;
int i=0; int i=0;
do { do {
if ( i == 0 ) { if ( i == 0 ) {
SetGridNormLoggingMode(GridNormLoggingModeRecord); FlightRecorder::SetLoggingMode(FlightRecorder::LoggingModeRecord);
} else { } else {
SetGridNormLoggingMode(GridNormLoggingModeVerify); FlightRecorder::SetLoggingMode(FlightRecorder::LoggingModeVerify);
} }
std::cerr << "******************* DOUBLE PRECISION SOLVE "<<i<<std::endl; std::cerr << "******************* DOUBLE PRECISION SOLVE "<<i<<std::endl;
result_o_2 = Zero(); result_o_2 = Zero();
@ -199,19 +198,9 @@ int main (int argc, char ** argv)
std::cout << " DoublePrecision iterations/sec "<< iters/(t2-t1)*1000.*1000.<<std::endl; std::cout << " DoublePrecision iterations/sec "<< iters/(t2-t1)*1000.*1000.<<std::endl;
std::cout << " DoublePrecision GF/s "<< flops/(t2-t1)/1000.<<std::endl; std::cout << " DoublePrecision GF/s "<< flops/(t2-t1)/1000.<<std::endl;
std::cout << " DoublePrecision error count "<< FlightRecorder::ErrorCount()<<std::endl;
csum = crc(result_o); assert(FlightRecorder::ErrorCount()==0);
std::cout << " FlightRecorder is OK! "<<std::endl;
if ( csumref == 0 ) {
csumref = csum;
} else {
if ( csum != csumref ) {
std::cerr << i <<" csum "<<std::hex<<csum<< " != "<<csumref <<std::dec<<std::endl;
assert(0);
} else {
std::cout << i <<" csum "<<std::hex<<csum<<std::dec<<" -- OK! "<<std::endl;
}
}
i++; i++;
} while (time(NULL) < (start + nsecs) ); } while (time(NULL) < (start + nsecs) );