clover_fermact_params_w.cc

Go to the documentation of this file.

/*! \file
 *  \brief Clover fermion action parameters
 */
 
#include "chromabase.h"
#include "actions/ferm/fermacts/clover_fermact_params_w.h"
 
#include "io/param_io.h"
 
namespace Chroma
{
 
  //! Default constructor
  CloverFermActParams::CloverFermActParams()
  {
    Mass = Real(0);
    u0   = Real(1);
    clovCoeffR = clovCoeffT = Real(0);
    max_norm=0;
    max_norm_usedP=false;
    twisted_m_usedP = false;
    twisted_m = Real(0);
 
  }
 
  //! Read parameters
  CloverFermActParams::CloverFermActParams(XMLReader& xml, const std::string& path)
  {
    XMLReader paramtop(xml, path);
    
    // Read the stuff for the action
    if (paramtop.count("Mass") != 0) 
    {
      read(paramtop, "Mass", Mass);
      if (paramtop.count("Kappa") != 0) 
      {
        QDPIO::cerr << "Error: found both a Kappa and a Mass tag" << std::endl;
        QDP_abort(1);
      }
    }
    else if (paramtop.count("Kappa") != 0)
    {
      Real Kappa;
      read(paramtop, "Kappa", Kappa);
      Mass = kappaToMass(Kappa);    // Convert Kappa to Mass
      QDPIO::cout << "Kappa is " << Kappa << "  Mass is " << Mass << std::endl << std::flush;
    }
    else
    {
      QDPIO::cerr << "Error: neither Mass or Kappa found" << std::endl;
      QDP_abort(1);
    }
 
    // Read optional u0
    if (paramtop.count("u0") != 0) 
      read(paramtop, "u0", u0);
    else {
      u0 = Real(1);
      QDPIO::cout << "u0 is " << u0 << std::endl << std::flush;
    }
    // Read optional anisoParam
    if (paramtop.count("AnisoParam") != 0) 
      read(paramtop, "AnisoParam", anisoParam);
 
    // If aniso, read all clover coeff
    if (anisoParam.anisoP)
    {
      read(paramtop, "clovCoeffR", clovCoeffR);
      read(paramtop, "clovCoeffT", clovCoeffT);
    }
    else
    {
      Real clovCoeff;
      read(paramtop, "clovCoeff", clovCoeff);
      clovCoeffR = clovCoeff;
      clovCoeffT = clovCoeff;
    }
 
    if( paramtop.count("MaxNorm") != 0 ){
      read(paramtop, "MaxNorm", max_norm);
      QDPIO::cout << "MaxNorm="<<max_norm<<std::endl;
      max_norm_usedP=true;
    }
    else { 
      max_norm_usedP=false;
      max_norm=Real(0);
    }
 
 
 
    if( paramtop.count("TwistedM") != 0 ) { 
      twisted_m_usedP = true;
      read(paramtop, "TwistedM", twisted_m);
    }
    else {                                     
      twisted_m_usedP = false;
    }
 
  }
 
  //! Read parameters
  void read(XMLReader& xml, const std::string& path, CloverFermActParams& param)
  {
    CloverFermActParams tmp(xml, path);
    param = tmp;
  }
 
 
  //! Write parameters
  void write(XMLWriter& xml, const std::string& path, const CloverFermActParams& param)
  {
    push(xml, path);
 
    write(xml, "Mass", param.Mass);
    write(xml, "u0", param.u0);
 
    if (param.anisoParam.anisoP)
    {
      write(xml, "clovCoeffR", param.clovCoeffR);
      write(xml, "clovCoeffT", param.clovCoeffT);
    }
    else
    {
      write(xml, "clovCoeff", param.clovCoeffR);
    }
 
    if (param.max_norm_usedP){
      write(xml, "MaxNorm", param.max_norm);
    }
 
 
    if (param.twisted_m_usedP == true ) { 
      write(xml, "TwistedM", param.twisted_m);
    }
 
    pop(xml);
  }
 
 
}