langevin_Momentum.h

Go to the documentation of this file.

#ifndef LANGEVIN_MOMENTUM_INCLUDED
#define LANGEVIN_MOMENTUM_INCLUDED

#include "commonParameters.h"
#include "parameters.h"
#include "field_G.h"
#include "randomNumbers.h"
#include "generatorSet_Mat_SU_N.h"

#include "bridgeIO.h"
using Bridge::vout;

class RandomNumbers;


class Langevin_Momentum
{
 private:
  RandomNumbers        *m_rand;
  Bridge::VerboseLevel m_vl;

 public:
  Langevin_Momentum(RandomNumbers *rand) :
    m_vl(CommonParameters::Vlevel())
  {
    m_rand = rand;
  }

  ~Langevin_Momentum() {}

 private:
  // non-copyable
  Langevin_Momentum(const Langevin_Momentum&);
  Langevin_Momentum& operator=(const Langevin_Momentum&);

 public:
  double set_iP(Field_G& iP);

 private:
  double set_iP_SU3(Field_G& iP);

  double set_iP_general_SU_N(Field_G& iP);

  double set_iP_SU3_alt(Field_G& iP);
};
#endif