int main()
{
Stopwatch stopwatch;
HANDLE hTimer = NULL;
// Create an unnamed waitable timer.
hTimer = CreateWaitableTimer(NULL, TRUE, NULL);
if (NULL == hTimer) {
printf("CreateWaitableTimer failed (%d)\n", GetLastError());
return 1;
}
// With Sleep
stopwatch.Start();
DoWait(hTimer, 1000);
printf("Wait 1000ms, Elapsed time : %d\n", stopwatch.ElapsedTime());
// With Sleep
stopwatch.Start();
DoWait(hTimer, 1234);
Sleep(100);
printf("Wait 1234ms with Sleep(100), Elapsed time : %d\n", stopwatch.ElapsedTime());
return 0;
}
int main(int argc, char* argv[])
{
if (argc != 6) {
cout << "usage: ./a.out ";
cout << "<coupling.arg> ";
cout << "<evo.arg> ";
cout << "<lattice.arg> ";
cout << "<plaquette.arg> ";
cout << "<retherm.arg> ";
cout << endl;
exit(EXIT_FAILURE);
}
// Get/set parameters
CouplingArg coupling_arg( argv[1] );
double beta = coupling_arg.beta;;
EvoArg evo_arg( argv[2] );
int n_conf = evo_arg.n_conf;
int n_chkpt = evo_arg.n_chkpt;
string cfg_file_stem = evo_arg.file_stem;
LatticeArg lattice_arg( argv[3] );
int sites[4];
int sitesX[4];
for (int mu=0; mu<4; ++mu) {
sites[mu] = lattice_arg.sites[mu];
if (sites[mu]%2 ==0) {
sitesX[mu] = sites[mu]/2;
} else {
cout << "Numer of lattice sites must be even!" << endl;
exit(EXIT_FAILURE);
}
}
enum lat_state state = lattice_arg.state;;
PlaquetteArg plaquette_arg( argv[4] );
int n_plaq = plaquette_arg.n_plaq;
string plaq_file_stem = plaquette_arg.file_stem;
RethermArg retherm_arg( argv[5]);
int retherm_sweeps = retherm_arg.sweeps;
std::cout << "Initializing the lattice" << "..." << std::endl;
Lattice latticeX(sitesX);
Lattice lattice(sites);
std::cout << "Initializing the monster..." << endl;;
PlaquetteMonster monster(sites, n_plaq);
vector<double> **plaq;
vector<double> ***plaqs;
plaq = new vector<double> * [n_plaq];
plaqs = new vector<double> ** [n_plaq];
for (int i=0; i<n_plaq; ++i) {
plaq[i] = new vector<double> [n_plaq];
plaqs[i] = new vector<double> * [n_plaq];
for (int j=0; j<n_plaq; ++j) {
plaqs[i][j]= new vector<double> [retherm_sweeps];
}
}
std::stringstream ss;
ofstream file;
cout << setprecision(15);
Stopwatch timer;
int conf = 0;
do {
timer.Restart();
if ( conf%n_chkpt==0 ) {
std::cout << "\t";
ss.str(std::string());
ss << cfg_file_stem << "-" << conf << ".bin";
if (!latticeX.Read(ss.str()) ) {
std::cout << "Failed to find lattice : " << std::endl;
std::cout << ss.str() << std::endl;
exit(EXIT_FAILURE);
};
std::cout << "\t";
std::cout << "Computing observables..." << std::endl;
monster.Initialize(latticeX);
for (int m=0; m<n_plaq/2; ++m)
for (int n=0; n<n_plaq/2; ++n) {
plaq[m][n].push_back( 1.0-monster.MeanPlaquette(m+1,n+1) );
}
if (retherm_sweeps>0) {
std::cout << "\t";
std::cout << "Refining lattice..." << std::endl;
if ( !InterpolatedTwoCellBoundaryRefine(latticeX, lattice) ) {
std::cout << "Failed to refine lattice!" << std::endl;
exit(EXIT_FAILURE);
}
for(int cool_sweeps=0; cool_sweeps<50; ++cool_sweeps) {
lattice.CoolUpdate(Subsets::two_cell_bulk_checker_board);
}
std::cout << "\t";
std::cout << "Performing rethermalization sweeps and computing observables..." << std::endl;
for (int h=0; h<retherm_sweeps; ++h) {
monster.Initialize(lattice);
for (int m=0; m<n_plaq; ++m)
for (int n=0; n<n_plaq; ++n) {
plaqs[m][n][h].push_back( 1.0-monster.MeanPlaquette(m+1,n+1) );
}
lattice.HeatBathUpdate(beta);
}
}
}
std::cout << "Total time for config " << conf << " is : " << timer.ElapsedTime() << std::endl;
++conf;
} while (conf < n_conf);
//// Write observables to a file ////
for (int m=0; m<n_plaq; ++m)
for (int n=0; n<n_plaq; ++n) {
ss.str(std::string());
ss << plaq_file_stem << "-" << m+1 << "_" << n+1 << ".dat";
std::cout << "Writing observables to file : ";
std::cout << ss.str() << std::endl;
file.open(ss.str().c_str());
file << setprecision(15);
for(int j=0; j<plaq[m][n].size(); ++j) {
file << plaq[m][n][j] << std::endl;
}
file.close();
}
for (int m=0; m<n_plaq; ++m)
for (int n=0; n<n_plaq; ++n)
for (int h=0; h<retherm_sweeps; ++h) {
ss.str(std::string());
ss << plaq_file_stem << "-" << m+1 << "_" << n+1 << "-" << h << ".dat";
std::cout << "Writing observables to file : ";
std::cout << ss.str() << std::endl;
file.open(ss.str().c_str());
file << setprecision(15);
for(int j=0; j<plaqs[m][n][h].size(); ++j) {
file << plaqs[m][n][h][j] << std::endl;
}
file.close();
}
//// DONE! ////
for (int i=0; i<n_plaq; ++i) {
for (int j=0; j<n_plaq; ++j) { delete [] plaqs[i][j]; }
delete [] plaq[i];
delete [] plaqs[i];
}
delete [] plaq;
delete [] plaqs;
cout << "Congratulations on a job... done!" << endl;
exit(EXIT_SUCCESS);
}
