datan.cpp

#include <iostream>
#include <vector> 
#include <fstream>
#include <cmath> 
#include <algorithm>

//inFile.is_open() is True if inFile is open. However, when we extract the data of inFile, 
//if there is a white space or a line break, the method is_open() closes the file. To prevent
//this we must use a while loop with the inFile.good() method

//----------.dat files structure for the Metropolis algorithm-------------//
//For 2D XY model:   
// Ntherm, Nmeas, Nstep, L, Beta
// E, dE
// Cv, dCv
// chi, dchi
// V, dV 
// AV, dAV
// computing time (in seconds)



int L, Nmeas, Ntherm, Nstep;
double beta, E, Cv, M, chi, v, av, t; 
double dE, dCv, dM, dchi, dv, dav;

int main(){
    int alg;
    char Files[100];
    char name[100];
    int d;
    std::cout << "Dimension of the model (2 or 3): ";
    std::cin >> d;
    std::cout << "File with list of names: ";
    std::cin >> Files;
    std::cout << " " << std::endl;

    double values;
    std::vector<std::string> Names;
    std::ifstream inFile; 
    inFile.open(Files);
    int cont = 0; 
    if (inFile.is_open() ) {  
        while ( inFile.good() ) {
        inFile >> name; 
        Names.push_back(name);    
        cont += 1;
        } 
    } 
    inFile.close();

    std::ofstream Cvfile; std::ofstream Chifile; std::ofstream Efile; std::ofstream Timefile; 
    std::ofstream Mfile; std::ofstream Vfile; std::ofstream AVfile;  

    for(int i = 0; i<cont-1; i++){
        std::vector<double> DATA; 
        //We extract the data from the .dat files//
        std::ifstream Data;
        Data.open(Names[i]);
        if (Data.is_open() ) {  
            while ( Data.good() ) {
            Data >> values; 
            DATA.push_back(values);    
            } 
        } 
        Data.close();   
        Ntherm = (int) DATA[0]; Nmeas = (int) DATA[1]; Nstep = (int) DATA[2]; L = (int) DATA[3];
        beta = DATA[4]; E = DATA[5]; dE = DATA[6]; Cv = DATA[7]; dCv = DATA[8]; M = DATA[9]; dM=DATA[10];
        chi = DATA[11]; dchi = DATA[12]; 
        v = DATA[13]; dv = DATA[14]; av = DATA[15]; dav = DATA[16]; t = DATA[17];
        
        //------------------------------------//

        //Now we generate the files with all the information//
        char NameCv[100], Cv_str[100], NameChi[100], Chi_str[100], NameE[100], E_str[100], NameTime[100], Time_str[100];
        char NameV[100], V_str[100], NameAV[100], AV_str[100], NameM[100], M_str[100];

        if (i == 0){
            sprintf(NameE, "%dDXY_E_L%d_Meas%d.txt", d,L, Nmeas);
            sprintf(NameCv, "%dDXY_Cv_L%d_Meas%d.txt", d, L, Nmeas);
            sprintf(NameChi, "%dDXY_Chi_L%d_Meas%d.txt", d, L, Nmeas);
            sprintf(NameTime, "%dDXY_Time_L%d_Meas%d.txt", d, L, Nmeas);
            sprintf(NameV, "%dDXY_Vort_L%d_Meas%d.txt", d, L, Nmeas);
            sprintf(NameAV, "%dDXY_AVort_L%d_Meas%d.txt", d, L, Nmeas);
            sprintf(NameM, "%dDXY_M_L%d_Meas%d.txt", d, L, Nmeas);
            Cvfile.open(NameCv); Chifile.open(NameChi); Efile.open(NameE); Timefile.open(NameTime); 
            Vfile.open(NameV); AVfile.open(NameAV); Mfile.open(NameM); 
        }
        sprintf(E_str, "%-30.17g%-30.17g%-30.17g\n", beta, E, dE);
        sprintf(Cv_str, "%-30.17g%-30.17g%-30.17g\n", beta, Cv, dCv);
        sprintf(Chi_str, "%-30.17g%-30.17g%-30.17g\n", beta, chi, dchi);
        sprintf(Time_str, "%-30.17g%-30.17g\n", beta, t);
        sprintf(M_str, "%-30.17g%-30.17g%-30.17g\n", beta, M, dM);
        sprintf(V_str, "%-30.17g%-30.17g%-30.17g\n", beta, v, dv);
        sprintf(AV_str, "%-30.17g%-30.17g%-30.17g\n", beta, av, dav);
        Efile << E_str;
        Cvfile << Cv_str;
        Chifile << Chi_str;
        Mfile << M_str;
        Timefile << Time_str;
        Vfile << V_str;
        AVfile << AV_str;
        std::cout << "------------------------------------------------------------------------" << std::endl;
        std::cout << "| beta = " << beta << "  T = " << 1/beta << "  L = " << L << std::endl;
        std::cout << "| Ntherm = " << Ntherm << "  Nmeas = " << Nmeas << "  Nstep = " << Nstep << std::endl;
        std::cout << "| E = " << E << " +- " << dE  << std::endl;
        std::cout << "| Cv = " << Cv << " +- " << dCv << std::endl;
        std::cout << "| M = " << M << " +- " << dM << std::endl;
        std::cout << "| Chi = " << chi << " +- " << dchi << std::endl; 
        std::cout << "| V = " << v << " +- " << dv << std::endl; 
        std::cout << "| Av = " << av << " +- " << dav << std::endl; 
        std::cout << "| Time = " << t << " seconds" << std::endl; 
        std::cout << "------------------------------------------------------------------------" << std::endl;
    }
    Efile.close();
    Cvfile.close();
    Chifile.close();
    Timefile.close();
    Mfile.close();
    Vfile.close();
    AVfile.close();

    return 0;
}