action_G_Plaq_SF.cpp

Go to the documentation of this file.

 
#include "action_G_Plaq_SF.h"
 
#ifdef USE_FACTORY_AUTOREGISTER
namespace {
  bool init = Action_G_Plaq_SF::register_factory();
}
#endif
 
const std::string Action_G_Plaq_SF::class_name = "Action_G_Plaq_SF";
 
//====================================================================
void Action_G_Plaq_SF::set_parameters(const Parameters& params)
{
  std::string vlevel;
  if (!params.fetch_string("verbose_level", vlevel)) {
    m_vl = vout.set_verbose_level(vlevel);
  }
 
  //- fetch and check input parameters
  double              beta, ct0, ct1, ct2;
  std::vector<double> phi, phipr;
 
  int err = 0;
  err += params.fetch_double("beta", beta);
  err += params.fetch_double("ct0", ct0);
  err += params.fetch_double("ct1", ct1);
  err += params.fetch_double("ct2", ct2);
  err += params.fetch_double_vector("phi", phi);
  err += params.fetch_double_vector("phipr", phipr);
 
  if (err) {
    vout.crucial(m_vl, "Error ar %s: input parameter not found.\n", class_name.c_str());
    exit(EXIT_FAILURE);
  }
 
 
  const double gg = 6.0 / beta;
  const double ct = ct0 + ct1 * gg + ct2 * gg * gg;
 
  set_parameters(beta, &phi[0], &phipr[0], ct);
 
  //- post-process
  m_force_G->set_parameters(params);
}
 
 
//====================================================================
void Action_G_Plaq_SF::get_parameters(Parameters& params) const
{
  params.set_double("beta", m_beta);
  // params.set_double("ct0", m_ct0);
  // params.set_double("ct1", m_ct1);
  // params.set_double("ct2", m_ct2);
  params.set_double_vector("phi", m_phi);
  params.set_double_vector("phipr", m_phipr);
 
  params.set_double("ct", m_ct);
 
  params.set_string("verbose_level", vout.get_verbose_level(m_vl));
}
 
 
//====================================================================
 
void Action_G_Plaq_SF::set_parameters(const double beta,
                                      double *phi, double *phipr,
                                      const double ct)
{
  //- print input parameters
  vout.general(m_vl, "%s:\n", class_name.c_str());
  vout.general(m_vl, "  beta  = %12.6f\n", beta);
  vout.general(m_vl, "  phi1  = %12.6f\n", phi[0]);
  vout.general(m_vl, "  phi2  = %12.6f\n", phi[1]);
  vout.general(m_vl, "  phi3  = %12.6f\n", phi[2]);
  vout.general(m_vl, "  phipr1= %12.6f\n", phipr[0]);
  vout.general(m_vl, "  phipr2= %12.6f\n", phipr[1]);
  vout.general(m_vl, "  phipr3= %12.6f\n", phipr[2]);
  vout.general(m_vl, "  ct    = %12.6f\n", ct);
 
  //- range check
  // NB. beta,phi,phipr,ct = 0 is allowed.
 
  //- store values
  m_beta = beta;
 
  // m_phi   = phi;
  // m_phipr = phipr;
  // m_ct    = ct;
  //  m_phi = std::vector<double>(phi, phi+3);
  //  m_phipr = std::vector<double>(phipr, phi+3);
  m_phi.resize(3);
  m_phipr.resize(3);
  for (int i = 0; i < 3; ++i) {
    m_phi[i]   = phi[i];
    m_phipr[i] = phipr[i];
  }
  //  m_phipr = std::vector<double>(phipr, phi+3);
  m_ct = ct;
 
  //- post-process
  m_staple.set_parameters(m_phi, m_phipr);
}
 
 
//====================================================================
double Action_G_Plaq_SF::langevin(RandomNumbers *rand)
{
  const double H_U = calcH();
 
  return H_U;
}
 
 
//====================================================================
 
double Action_G_Plaq_SF::calcH()
{
  const int Nc = CommonParameters::Nc();
 
  const double plaq = m_staple.plaquette_ct(*m_U, m_ct);
 
  const double H_U = -(1.0 / Nc) * m_beta * plaq;
 
  vout.general(m_vl, "H_Gplaq    = %18.8f\n", H_U);
 
  return H_U;
}
 
 
//====================================================================
 
void Action_G_Plaq_SF::force(Field& force)
{
  //- check of argument type
  assert(force.nin() == m_U->nin());
  assert(force.nvol() == m_U->nvol());
  assert(force.nex() == m_U->nex());
 
  vout.general(m_vl, "  %s: %s\n", class_name.c_str(), m_label.c_str());
 
  force.set(0.0);
 
  m_force_G->force_core(force, m_U);
}
 
 
//====================================================================
//============================================================END=====