integrator_Leapfrog.h

Go to the documentation of this file.

#ifndef INTEGRATOR_LEAPFROG_INCLUDED
#define INTEGRATOR_LEAPFROG_INCLUDED

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

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


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

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

 public:
  Integrator_Leapfrog(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;
  }

  ~Integrator_Leapfrog()
  {
    // std::cout << "byebye from integrator_leapfrog.\n";
  }

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

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