builder_Integrator.h

Go to the documentation of this file.

#ifndef BUILDER_INTEGRATOR_INCLUDED
#define BUILDER_INTEGRATOR_INCLUDED

#include <cassert>

#include "integrator_UpdateU.h"

#include "integrator_Leapfrog.h"
#include "integrator_Omelyan.h"

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


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

class Builder_Integrator
{
 protected:
  Bridge::VerboseLevel m_vl;

 private:
  double                      m_Estep; 
  int                         m_Nprec; 
  double                      m_lambda_Omelyan;
  std::string                 m_str_integrator_type;
  std::valarray<Action *>     m_action;
  std::valarray<Director *>   m_director;
  int                         m_Nlevel;  
  std::valarray<int>          m_Naction; 
  std::valarray<int>          m_Nstep;   
  std::valarray<Integrator *> m_integs;  

 public:
  Builder_Integrator(std::valarray<Action *> action,
                     std::valarray<Director *> director)
    : m_vl(CommonParameters::Vlevel())
  {
    m_action.resize(action.size());
    for (int i = 0; i < action.size(); ++i) {
      m_action[i] = action[i];
    }
    m_director.resize(director.size());
    for (int i = 0; i < director.size(); ++i) {
      m_director[i] = director[i];
    }
    m_Estep          = 0.0;
    m_Nprec          = 0;
    m_lambda_Omelyan = 0.0;
  }

  Builder_Integrator(std::valarray<Action *> action)
    : m_vl(CommonParameters::Vlevel())
  {
    m_action.resize(action.size());
    for (int i = 0; i < action.size(); ++i) {
      m_action[i] = action[i];
    }
    m_Estep          = 0.0;
    m_Nprec          = 0;
    m_lambda_Omelyan = 0.0;
  }

  ~Builder_Integrator()
  {
    tidyup();
  }

  void set_parameters(const Parameters& params);
  void set_parameters(std::string str_integrator_type,
                      double Estep, int Nlevel,
                      const std::valarray<int>& Naction,
                      const std::valarray<int>& Nstep,
                      const int Nprec,
                      const double lambda_Omelyan);

  void set_parameter_verboselevel(const Bridge::VerboseLevel vl) { m_vl = vl; }

  Integrator *build();

  Integrator *build_leapfrog();
  Integrator *build_omelyan();

  void tidyup();
};
#endif