lw_tree_gaugeact.cc

Go to the documentation of this file.

/*! \file
 *  \brief Tree-level tadpole-improved Luscher-Weisz gauge action
 */
 
#include "chromabase.h"
#include "actions/gauge/gaugeacts/lw_tree_gaugeact.h"
#include "actions/gauge/gaugeacts/gaugeact_factory.h"
#include "actions/gauge/gaugestates/gauge_createstate_aggregate.h"
 
namespace Chroma
{
 
  namespace LWTreeGaugeActEnv 
  { 
    GaugeAction< multi1d<LatticeColorMatrix>, 
                 multi1d<LatticeColorMatrix> >* createGaugeAct(XMLReader& xml, 
                                                               const std::string& path) 
    {
      return new LWTreeGaugeAct(CreateGaugeStateEnv::reader(xml, path), 
                                LWTreeGaugeActParams(xml, path));
    }
 
    const std::string name = "LW_TREE_GAUGEACT";
 
    //! Local registration flag
    static bool registered = false;
 
    //! Register all the factories
    bool registerAll() 
    {
      bool success = true; 
      if (! registered)
      {
        success &= TheGaugeActFactory::Instance().registerObject(name, createGaugeAct);
        registered = true;
      }
      return success;
    }
  }
 
 
  LWTreeGaugeActParams::LWTreeGaugeActParams(XMLReader& xml_in, const std::string& path) {
    XMLReader paramtop(xml_in, path);
 
    try {
      read(paramtop, "./beta", beta);
      read(paramtop, "./u0", u0);
 
      //  Read optional anisoParam.
      if (paramtop.count("AnisoParam") != 0) 
        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, LWTreeGaugeActParams& p) {
    LWTreeGaugeActParams tmp(xml, path);
    p=tmp;
  }
 
 
  // Private initializer
  void
  LWTreeGaugeAct::init(Handle< CreateGaugeState<P,Q> > cgs)
  {
    START_CODE();
 
    // Fold in normalizations and create action
    // NOTE: the 5/3 is folded into beta, hence divided out of c1
    Real c0 = param.beta;
    plaq = new PlaqGaugeAct(cgs,c0,param.aniso);
 
    Real c1 = -c0 * (1/(20*param.u0*param.u0));
    rect = new RectGaugeAct(cgs,c1);
    
    END_CODE();
  } 
 
}