schr_triv_gaugebc.h

Go to the documentation of this file.

// -*- C++ -*-
/*! \file
 *  \brief Schroedinger functional trivial gauge boundary conditions
 */
 
#ifndef __schr_triv_gaugebc_h__
#define __schr_triv_gaugebc_h__
 
#include "actions/gauge/gaugebcs/schr_sf_gaugebc.h"
#include "actions/gauge/gaugebcs/schr_gaugebc_params.h"
 
namespace Chroma 
{ 
 
  /*! @ingroup gaugebcs */
  namespace SchrTrivialGaugeBCEnv 
  {
    extern const std::string name;
    bool registerAll();
  }
 
 
  //! Concrete class for gauge action boundary conditions with Schroedinger BC
  /*! @ingroup gaugebcs
   *
   *  Schroedinger BC for gauge actions
   */
  class SchrTrivialGaugeBC : public SchrSFGaugeBC
  {
  public:
    //! Only full constructor
    SchrTrivialGaugeBC(const SchrGaugeBCParams& p);
 
    //! Destructor is automatic
    ~SchrTrivialGaugeBC() {}
 
    //! 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 */
    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
    SchrTrivialGaugeBC() {}
    void operator=(const SchrTrivialGaugeBC&) {}
 
  private:
    SchrGaugeBCParams            param;
    Phases_t                     phases;
    multi1d<LatticeBoolean>      mask;
    multi1d<LatticeColorMatrix>  fld;
  };
  
}
 
#endif