SrednickiAnsatz.cpp File Reference

Functions
int	main (int argc, char **argv)

Function Documentation

main()

int main	(	int	argc,
		char **	argv
	)

                                {
    if(argc != 4) {
        std::cerr << "Usage: 1.This 2.Directory 3.Nmin 4.Nmax" << std::endl;
        std::exit(EX_USAGE);
    }
    std::string const InputDir( std::string(argv[1])+"/ProcessedData");
    int const Nmin = std::atoi(argv[2]);
    int const Nmax = std::atoi(argv[3]);
 
    std::vector<std::string> dirList;
  {
    std::string str;
    int dataNo;
    for(const std::filesystem::directory_entry &entry : std::filesystem::directory_iterator(InputDir)) {
      if ( entry.is_directory() ) {
        str = entry.path().string();
        dataNo = std::atoi(str.substr(str.find("_No")+3).c_str());
        debug_print("# dir=" << str << ", dataNo=" << dataNo);
        if(dataNo>= dirList.size()) dirList.resize(dataNo+1);
        dirList.at(dataNo) = str;
      }
    }
  }
    // #pragma omp parallel for schedule(dynamic, 10)
  // for(size_t sample = 0;sample < dirList.size(); ++sample) {
    //
    //  #pragma omp critical
    //  std::cout << "# Dir: " << dirList[sample] << std::endl;
    //
    // }
    #pragma omp parallel for schedule(dynamic, 10)
  for(size_t sample = 0;sample < dirList.size(); ++sample) {
        #pragma omp critical
        std::cout << "# Dir(" << sample << "): " << dirList[sample] << std::endl;
 
        std::filesystem::create_directories(std::regex_replace(dirList[sample], std::regex("/LocalRandomMatrix"), "/LocalRandomMatrix/test"));
        int Ndata, itemp;
        double Energy, Stddev, dtemp;
        std::ifstream InFile;
        for(int n = Nmin;n <= Nmax; ++n) {
            std::stringstream ss("");
            ss << "/SrednickiAnsatz_Diagonal_N" << n << ".txt";
            std::string InFileName(dirList[sample] + ss.str());
            if( InFile.is_open() ) { InFile.close(); InFile.clear(); }
            InFile.open(InFileName);
                if( !InFile.is_open() ) {
                    std::cerr << "Error: Can't open file " << InFileName << std::endl;
                    continue;
                }
            std::string OutFName(std::regex_replace(InFileName, std::regex("/LocalRandomMatrix"), "/LocalRandomMatrix/test"));
            std::ofstream OutFile(OutFName); checkIsFileOpen(OutFile, OutFName);
            std::string str;
            OutFile << std::scientific;
            for(getline(InFile, str); !InFile.eof(); getline(InFile, str)) {
                if( str[0] == '#' || str[0] == '\n' ) {
                    OutFile << str << "\n";
                    continue;
                }
                ss.str(str); ss.clear(std::stringstream::goodbit);
                ss >> itemp >> Energy >> dtemp >> dtemp >> dtemp >> dtemp >> Ndata;
                Stddev = 0;
                for(int j = 0;j < Ndata; ++j) {
                    dtemp = Gaussian();
                    Stddev += dtemp*dtemp;
                }
                Stddev = sqrt( Stddev/(double)Ndata )/(double)std::log(Ndata);
                dtemp = NAN;
                OutFile << n      << " "
                        << Energy << " "
                                << dtemp  << " "
                                << Stddev << " "
                                << dtemp  << " "
                                << dtemp  << " "
                                << Ndata  << "\n";
            }
            OutFile.close();
        }
    }
    return 0;
}