fts/montecarlo/McMoveManager.h

#ifndef RP_MC_MOVE_MANAGER_H
#define RP_MC_MOVE_MANAGER_H
 
/*
* PSCF - Polymer Self-Consistent Field
*
* Copyright 2015 - 2025, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
#include <util/param/Manager.h>          // base class template
#include <util/containers/DArray.h>      // member
 
// Forward declaration
namespace Util { class Random; }
 
namespace Pscf {
namespace Rp {
 
   using namespace Util;

   template <int D, class T>
   class McMoveManager : public Manager< typename T::McMove >
   {
 
   public:
 
      // Protected constructor and destructor (see below).

      virtual void readParameters(std::istream &in);

      void setup();

      typename T::McMove& chooseMove();

      void output() const;

      double probability(int i) const;

      void outputTimers(std::ostream& out) const;

      void clearTimers();

      bool needsCc();

      bool needsDc();
 
   protected:

      McMoveManager(typename T::McSimulator& simulator, 
                    typename T::System& system);

      ~McMoveManager() = default;
 
   private:
 
      using McMoveT = typename T::McMove;
      using Base = Manager< McMoveT >;
 
      // Private data members

      DArray<double>  probabilities_;

      typename T::McSimulator* simulatorPtr_;

      typename T::System* systemPtr_;

      Random* randomPtr_;
 
      // Private member functions

      virtual Factory< McMoveT >* newDefaultFactory() const;
 
   };
 
   // Inline functions
 
   /*
   * Return probability of move number i.
   */
   template <int D, class T> inline 
   double McMoveManager<D,T>::probability(int i) const
   {
      assert(i >= 0);
      assert(i < Base::size());
      return probabilities_[i];
   }
 
}
}
#endif