outputterm.cpp

Go to the documentation of this file.

/******************************************************************************
* 
* MIT License
* 
* Copyright (c) 2018 Giovanni Pederiva
* 
* Permission is hereby granted, free of charge, to any person obtaining a copy
* of this software and associated documentation files (the "Software"), to deal
* in the Software without restriction, including without limitation the rights
* to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
* copies of the Software, and to permit persons to whom the Software is
* furnished to do so, subject to the following conditions:
* 
* The above copyright notice and this permission notice shall be included in 
* all copies or substantial portions of the Software.
* 
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
* AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
* OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
* SOFTWARE.
******************************************************************************/

#include <cstdio>
#include "InputOutput/outputterm.h"
#include "Math/lattice.h"
#include "ParallelTools/parallel.h"
#include "Observables/observable.h"

namespace LatticeIO {

    void OutputTerm::printInitialConditions(){
        if(Parallel::rank() == 0){
            printf("Total Lattice Size: %2i %2i %2i %2i\n", Parallel::latticeFullSize()[0],Parallel::latticeFullSize()[1],
                                                            Parallel::latticeFullSize()[2], Parallel::latticeFullSize()[3]);
            printf("Sublattice Size:    %2i %2i %2i %2i\n\n", Parallel::latticeSubSize()[0],Parallel::latticeSubSize()[1],
                                                              Parallel::latticeSubSize()[2],Parallel::latticeSubSize()[3]);
            printf("Number of Processors: \t\t%i\n", Parallel::numProcs());
            printf("Number of Active Processors: \t%i\n\n", Parallel::activeProcs());
        }
    }

    void OutputTerm::printThermStep(int step, std::vector<Observable*>& obsList, double acceptRatio){
        if(Parallel::rank() == 0){
            printf("Thermalization Step %4i\n", step);
            for(auto& obs : obsList)
                printf("\t%s\t = %lf \n", obs->getName(), obs->value());
            printf("\tAccept. Ratio = %lf\n", acceptRatio);
        }
    }

    void OutputTerm::printGenStep(int confNum, std::vector<Observable*>& obsList, double acceptRatio){
        if(Parallel::rank() == 0){
            printf("Configuration Number %4i\n", confNum);
            for(auto& obs : obsList)
                printf("\t%s\t = %lf \n", obs->getName(), obs->value());
            printf("\tAccept. Ratio = %lf\n", acceptRatio);
        }
    }

    void OutputTerm::writeObs(int confNum, std::vector<Observable*>& obsList){
        if(Parallel::rank() == 0){
            printf("Configuration Number %4i\n", confNum);
            //or(auto& obs : obsList)
            //    printf("\t%s\t = %lf \n", obs->getName(), obs->value());
            printf("\t%s\t = %g \n", "Plaq.", obsList[0]->plaq);
            printf("\t%s\t = %g \n", "Topc.", obsList[0]->topc);
            printf("\t%s\t = %g \n", "Energ.", obsList[0]->energy);
            }

    }

    void OutputTerm::writeFlowObs(double flowTime, std::vector<Observable*>& obsList){
        if(Parallel::rank() == 0){
            printf("Flow Time %2.2lf\n", flowTime);
            //for(auto& obs : obsList)
            //    printf("\t%s\t = %lf \n", obs->getName(), obs->value());
            printf("\t%s\t = %g \n", "Plaq.", obsList[0]->plaq);
            printf("\t%s\t = %g \n", "Topc.", obsList[0]->topc);
            printf("\t%s\t = %g \n", "Energ.", obsList[0]->energy);
        }
    }

} // end LatticeIO