rpg/fts/montecarlo/McMove.h

#ifndef RPG_MC_MOVE_H
#define RPG_MC_MOVE_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/ParamComposite.h>    // base class
#include <util/misc/Timer.h>              // member
 
// Forward declarations
namespace Util { 
   class Random; 
}
namespace Pscf { 
   class CudaRandom; 
   namespace Rpg { 
      template <int D> class System;
      template <int D> class McSimulator;
   }
}
 
namespace Pscf {
namespace Rpg {
 
   using namespace Util;

   template <int D>
   class McMove : public ParamComposite
   {
 
   public:

      McMove(McSimulator<D>& simulator);

      virtual ~McMove();

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

      void setProbability(double probability);

      virtual void setup();

      virtual bool move();
      
      virtual bool needsCc()
      {  return false; }
      
      virtual bool needsDc()
      { return false; }
      
      virtual void outputTimers(std::ostream& out);
      
      virtual void clearTimers();
 
      // Accessor Functions

      double probability() const;

      long nAttempt() const;

      long nAccept() const;
      
      long nFail() const;

      virtual void output();
 
   protected:

      void incrementNAttempt();

      void incrementNAccept();
      
      void incrementNFail();

      System<D>& system();

      McSimulator<D>& simulator();

      Random& random();
      
      CudaRandom& cudaRandom();
      
      void readProbability(std::istream& in);

      virtual void attemptMove() 
      {};

      Timer attemptMoveTimer_;
      Timer compressorTimer_;
      Timer componentTimer_;
      Timer hamiltonianTimer_;
      Timer decisionTimer_;
      Timer totalTimer_;
 
   private:

      McSimulator<D>* simulatorPtr_;

      System<D>* systemPtr_;

      Random  *randomPtr_;
      
      CudaRandom  *cudaRandomPtr_;

      double  probability_;

      long  nAttempt_;

      long  nAccept_;
      
      long  nFail_;
      
   };
 
   // Public inline methods
 
   /*
   * Return number of moves that have been attempted.
   */
   template <int D>
   inline long McMove<D>::nAttempt() const
   {  return nAttempt_; }
 
   /*
   * Return number of moves that have been accepted.
   */
   template <int D>
   inline long McMove<D>::nAccept() const
   {  return nAccept_; }
   
   /*
   * Return number of moves that fail to converge.
   */
   template <int D>
   inline long McMove<D>::nFail() const
   {  return nFail_; }
 
   // Protected inline methods
 
   /*
   * Increment the number of attempted moves.
   */
   template <int D>
   inline void McMove<D>::incrementNAttempt()
   {  ++nAttempt_; }
 
   /*
   * Increment the number of accepted moves.
   */
   template <int D>
   inline void McMove<D>::incrementNAccept()
   {  ++nAccept_; }
   
   /*
   * Increment the number of fail moves.
   */
   template <int D>
   inline void McMove<D>::incrementNFail()
   {  ++nFail_; }
 
   /*
   * Get parent System object.
   */
   template <int D>
   inline System<D>& McMove<D>::system()
   {  return *systemPtr_; }
 
   /*
   * Get parent McSimulator object.
   */
   template <int D>
   inline McSimulator<D>& McMove<D>::simulator()
   {  return *simulatorPtr_; }
 
   /*
   * Get Random number generator.
   */
   template <int D>
   inline Random& McMove<D>::random()
   {  return *randomPtr_; }
 
   /*
   * Get CudaRandom number generator.
   */
   template <int D>
   inline CudaRandom& McMove<D>::cudaRandom()
   {  return *cudaRandomPtr_; }
   
   /*
   * Get the probability.
   */
   template <int D>
   inline double McMove<D>::probability() const
   {  return probability_; }
 
   /*
   * Set the probability.
   */
   template <int D>
   inline void McMove<D>::setProbability(double probability)
   {  probability_ = probability; }
 
   // Explicit instantiation declarations
   extern template class McMove<1>;
   extern template class McMove<2>;
   extern template class McMove<3>;
 
}
}
#endif