schr_nonpert_gaugebc.cc

Go to the documentation of this file.

/*! \file
 *  \brief Schroedinger BC - use for non-pertubative tuning of clover action
 */
 
#include "actions/gauge/gaugebcs/schr_nonpert_gaugebc.h"
#include "actions/gauge/gaugebcs/gaugebc_factory.h"
 
namespace Chroma 
{
 
  namespace SchrNonPertGaugeBCEnv 
  { 
    //! Callback function to register with the factory
    GaugeBC< multi1d<LatticeColorMatrix>, multi1d<LatticeColorMatrix> >* createGaugeBC(XMLReader& xml, 
                                                                                       const std::string& path)
    {
      return new SchrNonPertGaugeBC(SchrGaugeBCParams(xml, path));
    }
 
    const std::string name = "SCHROEDINGER_NONPERT_GAUGEBC";
 
    //! Local registration flag
    static bool registered = false;
 
    //! Register all the factories
    bool registerAll() 
    {
      bool success = true; 
      if (! registered)
      {
        success &= TheGaugeBCFactory::Instance().registerObject(name, createGaugeBC);
        registered = true;
      }
      return success;
    }
  }
 
 
  // Only full constructor
  SchrNonPertGaugeBC::SchrNonPertGaugeBC(const SchrGaugeBCParams& p) : 
    param(p)
  {
    // Initialize the phases
    initPhases();
 
    // Initialize the boundary fields
    initBnd(fld, mask);
  }
 
 
  // Initialize the phases
  void SchrNonPertGaugeBC::initPhases()
  {
    phases.lower.resize(Nc);
    phases.upper.resize(Nc);
 
    Real ftmp = Chroma::twopi * 0.5 * SchrPhiMult();
 
    switch (Nc)
    {
    case 2:                     /*  SF nonperturbative improvement bc's */
      ftmp /= Real(4);
 
      /* lower boundary */
      phases.lower[0] = -ftmp;
      phases.lower[1] =  ftmp;
 
      /* upper boundary */
      phases.upper[0] = -Real(3) * ftmp;
      phases.upper[1] =  Real(3) * ftmp;
      break;
 
    case 3:                     /*  SF nonperturbative improvement bc's */
      ftmp /= Real(6);
 
      /* lower boundary */
      phases.lower[0] = -ftmp;
      phases.lower[1] = 0;
      phases.lower[2] = ftmp;
 
      /* upper boundary */
      phases.upper[0] = -Real(5) * ftmp;
      phases.upper[1] = Real(2) * ftmp;
      phases.upper[2] = Real(3) * ftmp;
      break;
 
    default:
      QDPIO::cerr << __func__ << ": unsupport number of colors" << std::endl;
      QDP_abort(1);
    }
  }
 
 
}