integrator_UpdateU.h

Go to the documentation of this file.

#ifndef INTEGRATOR_UPDATEU_INCLUDED
#define INTEGRATOR_UPDATEU_INCLUDED

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

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


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

class Integrator_UpdateU : public Integrator
{
 public:
  static const std::string class_name;

  static const int default_nprec = 8;

 private:
  int m_Nprec;                          // precision of exponential series

  std::vector<Director *>   m_director; // directors
  std::vector<Integrator *> m_integs;   // list of update_p to send notification

 public:

  Integrator_UpdateU(const std::vector<Director *>& director = std::vector<Director *>())
    : m_Nprec(default_nprec), m_director(director), m_integs()
  {
  }

  ~Integrator_UpdateU() {}

  void set_parameters(const Parameters& params);
  void set_parameters(int Nprec);

  void set_parameter_nprec(const int Nprec);


  void evolve(const double step_size, Field_G& iP, Field_G& U);


  // cache management
  void invalidate_cache()
  {
    // lowest level. do nothing.
  }

  void append_notify(Integrator *const integ)
  {
    m_integs.push_back(integ);
  }

  void notify_update();
};
#endif