gaugeConfig.cpp

Go to the documentation of this file.

#include "gaugeConfig.h"

const std::string GaugeConfig::class_name = "GaugeConfig";

//====================================================================
GaugeConfig::GaugeConfig(const string& type)
  : m_vl(CommonParameters::Vlevel()), m_fieldio()
{
  if (type == "Text") {
    m_fieldio = new FieldIO_Text(IO_Format::Gauge::ILDG);
  } else if (type == "Text_4x4x4x8") {
    m_fieldio = new FieldIO_Text_4x4x4x8(IO_Format::Gauge::ILDG);
  } else if (type == "Binary") {
    m_fieldio = new FieldIO_Binary(IO_Format::Gauge::ILDG);
  } else if (type == "Fortran_JLQCD") {
    m_fieldio = new FieldIO_Fortran(IO_Format::Gauge::JLQCD);
  } else if (type == "Fortran_ILDG") {
    m_fieldio = new FieldIO_Fortran(IO_Format::Gauge::ILDG);
  } else if (type == "ILDG") {
    m_fieldio = new FieldIO_LIME(IO_Format::Gauge::ILDG);
  } else if (type == "Binary_Parallel") {
    m_fieldio = new FieldIO_Binary_Parallel(IO_Format::Gauge::ILDG);
  } else if (type == "Binary_Distributed") {
    m_fieldio = new FieldIO_Binary_Distributed(IO_Format::Gauge::ILDG);
  } else if (type == "ILDG_Parallel") {
    m_fieldio = new FieldIO_LIME_MPIIO(IO_Format::Gauge::ILDG);
  } else {
    vout.crucial("%s: unsupported type \"%s\".\n", class_name.c_str(), type.c_str());
    exit(EXIT_FAILURE);
  }
}


//====================================================================
GaugeConfig::~GaugeConfig()
{
  if (m_fieldio) delete m_fieldio;
  m_fieldio = 0;
}


//====================================================================
void GaugeConfig::read_file(Field *U, const string& filename)
{
  if (!m_fieldio) {
    vout.crucial("%s::read_file(): FieldIO not set.\n", class_name.c_str());
    exit(EXIT_FAILURE);
  }

  m_fieldio->read_file(U, filename);
}


//====================================================================
void GaugeConfig::read_file(unique_ptr<Field_G>& U, const string& filename)
{
  if (!m_fieldio) {
    vout.crucial("%s::read_file(): FieldIO not set.\n", class_name.c_str());
    exit(EXIT_FAILURE);
  }

  m_fieldio->read_file((Field *)U.get(), filename);
}


//====================================================================
void GaugeConfig::write_file(Field *U, const string& filename)
{
  if (!m_fieldio) {
    vout.crucial("%s::write_file(): FieldIO not set.\n", class_name.c_str());
    exit(EXIT_FAILURE);
  }

  string filename_header = filename.substr(0, 9);

  if (filename_header != "NO_OUTPUT") {
    m_fieldio->write_file(U, filename);
  }
}


//====================================================================
void GaugeConfig::write_file(unique_ptr<Field_G>& U, const string& filename)
{
  if (!m_fieldio) {
    vout.crucial("%s::write_file(): FieldIO not set.\n", class_name.c_str());
    exit(EXIT_FAILURE);
  }

  string filename_header = filename.substr(0, 9);

  if (filename_header != "NO_OUTPUT") {
    m_fieldio->write_file((Field *)U.get(), filename);
  }
}


//====================================================================
//============================================================END=====