rpc/fts/montecarlo/McMoveManager.h

#ifndef RPC_MC_MOVE_MANAGER_H
#define RPC_MC_MOVE_MANAGER_H
 
/*
* PSCF - Polymer Self-Consistent Field Theory
*
* Copyright 2016 - 2022, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
#include "McMove.h"                      // base class template parameter
#include <util/param/Manager.h>          // base class template
#include <util/containers/DArray.h>      // member template
 
namespace Util { class Random; }
 
namespace Pscf {
namespace Rpc {
 
   using namespace Util;
 
   template <int D> class System;
   template <int D> class McSimulator;
 
   template <int D>
   class McMoveManager : public Manager< McMove<D> >
   {
 
   public:
 
      McMoveManager(McSimulator<D>& simulator, System<D>& system);
 
      ~McMoveManager();
 
      virtual void readParameters(std::istream &in);
 
      void setup();
 
      McMove<D>& chooseMove();
 
      void output() const;
 
      double probability(int i) const;
      
      using Manager< McMove<D> >::size;
      
      void outputTimers(std::ostream& out) const;
      
      void clearTimers();
      
      bool needsCc();
      
      bool needsDc();
 
   protected:
      
      using Manager< McMove<D> >::setClassName;
   
   private:
 
      // Private data members
 
      DArray<double>  probabilities_;
      
      McSimulator<D>* simulatorPtr_;
 
      System<D>* systemPtr_;
 
      Random* randomPtr_;
 
      // Private member functions
 
      virtual Factory< McMove<D> >* newDefaultFactory() const;
 
   };
 
   // Inline functions
 
   /*
   * Return probability of move number i
   */
   template <int D>
   inline double McMoveManager<D>::probability(int i) const
   {
      assert(i >= 0);  
      assert(i < size());  
      return probabilities_[i];
   }
 
   #ifndef RPC_MC_MOVE_MANAGER_TPP
   // Suppress implicit instantiation
   extern template class McMoveManager<1>;
   extern template class McMoveManager<2>;
   extern template class McMoveManager<3>;
   #endif
 
}
}
#endif