schr_nonpert_gaugebc.h

Go to the documentation of this file.

// -*- C++ -*-
/*! \file
 *  \brief Schroedinger BC - use for non-pertubative tuning of clover action
 */
 
#ifndef __schr_nonpert_gaugebc_h__
#define __schr_nonpert_gaugebc_h__
 
#include "actions/gauge/gaugebcs/schr_sf_gaugebc.h"
#include "actions/gauge/gaugebcs/schr_gaugebc_params.h"
 
namespace Chroma 
{ 
 
  /*! @ingroup gaugebcs */
  namespace SchrNonPertGaugeBCEnv 
  {
    extern const std::string name;
    bool registerAll();
  }
 
 
  //! Concrete class for Schroedinger BC - use for nonpertubative tuning
  /*! @ingroup gaugebcs
   *
   *  Schroedinger BC for gauge actions
   */
  class SchrNonPertGaugeBC : public SchrSFGaugeBC
  {
  public:
    //! Only full constructor
    SchrNonPertGaugeBC(const SchrGaugeBCParams& p);
 
    //! Destructor is automatic
    ~SchrNonPertGaugeBC() {}
 
    //! Decay direction
    int getDir() const {return param.decay_dir;}
 
    //! Mask which lattice sites have fixed gauge links
    const multi1d<LatticeBoolean>& lSFmask() const {return mask;}
 
    //! Fixed gauge links on only the lSFmask() sites
    const multi1d<LatticeColorMatrix>& SFBndFld() const {return fld;}
 
  protected:
    //! Maximum plaquette size. This is what knows about 1x1 plaq or 1x2 rect.
    /*! \return 1 for 1x1 plaq or 2 for 1x2 rect in decay_dir */
    int getMaxExtent() const {return param.loop_extent;}
 
    //! Multiplier on phases
    const Real& SchrPhiMult() const {return param.SchrPhiMult;}
 
    //! Get the angles on the boundaries
    const Phases_t& getPhases() const {return phases;}
 
    //! Initialize the phases
    void initPhases();
 
  private:
    // Hide default constuctor
    SchrNonPertGaugeBC() {}
    void operator=(const SchrNonPertGaugeBC&) {}
 
  private:
    SchrGaugeBCParams            param;
    Phases_t                     phases;
    multi1d<LatticeBoolean>      mask;
    multi1d<LatticeColorMatrix>  fld;
  };
  
}
 
#endif