integrator_Omelyan.h

Go to the documentation of this file.

#ifndef INTEGRATOR_OMELYAN_INCLUDED
#define INTEGRATOR_OMELYAN_INCLUDED

#include "action.h"
#include "director.h"
#include "randomNumbers.h"
#include "integrator.h"

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


//- parameters class
class Parameters_Integrator_Omelyan : virtual public Parameters
{
 public:
  Parameters_Integrator_Omelyan();
};
//- end

class Integrator_Omelyan : public Integrator
{
 private:
  double                  m_Estep;  // step size of each molecular dynamical evolution
  int                     m_Nstep;  // number of steps
  int                     m_level;  // level number
  double                  m_lambda; // Omelyan's optimized parameter
  Integrator              *m_integ_next;
  std::valarray<Action *> m_action; // actions

 public:
  Integrator_Omelyan(std::valarray<Action *> action,
                     Integrator *integ_next)
  {
    m_action.resize(action.size());
    for (int i = 0; i < action.size(); ++i) {
      m_action[i] = action[i];
    }
    m_integ_next = integ_next;
    m_Estep      = 0.0;
    m_Nstep      = 0;
    m_lambda     = 0.0;
  }

  ~Integrator_Omelyan() {}

  void set_parameters(const Parameters& params);
  void set_parameters(int level, double Estep, int Nstep,
                      double lambda_Omelyan);

  void evolve(Field_G& iP, Field_G& U);
};
#endif