aniso_sym_gaugeact_params.cc

Go to the documentation of this file.

/*! \file
 *  \brief Anisotropic gaugeact useful for spectrum from hep-lat/9911003
 *
 *  Tree-level LW with tapole improvement, missing 1x2 in time, also including
 *  2-plaq term. Taken from Morningstar-Peardon, hep-lat/9911003
 */
 
#include "chromabase.h"
 
#include "actions/gauge/gaugeacts/aniso_sym_gaugeact_params.h"
 
namespace Chroma
{
 
  
  AnisoSymGaugeActParams::AnisoSymGaugeActParams(XMLReader& xml_in, const std::string& path) 
  {
    XMLReader paramtop(xml_in, path);
 
    try 
    {
      read(paramtop, "beta", beta);
      read(paramtop, "u_s", u_s);
      read(paramtop, "u_t", u_t);
      read(paramtop, "AnisoParam", aniso);
      if( paramtop.count("ZeroEnergy") == 1 ) {
        read(paramtop, "ZeroEnergy", sub_zero);
        use_subtraction = true;
      }
      else { 
        sub_zero = Real(0);
        use_subtraction = false;
      }
    }
    catch( const std::string& e ) { 
      QDPIO::cerr << "Error reading XML: " <<  e << std::endl;
      QDP_abort(1);
    }
  }
 
 
  void read(XMLReader& xml, const std::string& path, AnisoSymGaugeActParams& p) 
  {
    AnisoSymGaugeActParams tmp(xml, path);
    p=tmp;
  }
 
  void write(XMLWriter& xml, const std::string& path, const AnisoSymGaugeActParams& param) 
  {
    push(xml, path);
 
    write(xml, "beta", param.beta);
    write(xml, "u_s", param.u_s);
    write(xml, "u_t", param.u_t);
    if( param.use_subtraction ) { 
      write(xml, "ZeroEnergy", param.sub_zero);
    }
    write(xml, "AnisoParam", param.aniso);
 
    pop(xml);
  }
 
  AnisoSymSpatialGaugeActParams::AnisoSymSpatialGaugeActParams(XMLReader& xml_in, const std::string& path) 
  {
    XMLReader paramtop(xml_in, path);
 
    try 
    {
      read(paramtop, "beta", beta);
      read(paramtop, "u_s", u_s);
      if( paramtop.count("ZeroEnergy") == 1 ) {
        read(paramtop, "ZeroEnergy", sub_zero);
        use_subtraction = true;
      }
      else { 
        sub_zero = Real(0);
        use_subtraction = false;
      }
      read(paramtop, "AnisoParam", aniso);
    }
    catch( const std::string& e ) { 
      QDPIO::cerr << "Error reading XML: " <<  e << std::endl;
      QDP_abort(1);
    }
  }
 
 
  void read(XMLReader& xml, const std::string& path, AnisoSymSpatialGaugeActParams& p) 
  {
    AnisoSymSpatialGaugeActParams tmp(xml, path);
    p=tmp;
  }
 
  void write(XMLWriter& xml, const std::string& path, const AnisoSymSpatialGaugeActParams& param) 
  {
    push(xml, path);
 
    write(xml, "beta", param.beta);
    write(xml, "u_s", param.u_s);
    if( param.use_subtraction ) { 
      write(xml, "ZeroEnergy", param.sub_zero);
    }
    write(xml, "AnisoParam", param.aniso);
 
    pop(xml);
  }
 
}