rpc/fts/montecarlo/RealMove.tpp

#ifndef RPC_REAL_MOVE_TPP
#define RPC_REAL_MOVE_TPP
 
/*
* 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 "RealMove.h"
#include "McMove.h" 
#include <util/param/ParamComposite.h>
#include <rpc/System.h>
#include <util/archives/Serializable_includes.h>
#include <util/random/Random.h>
 
 
namespace Pscf {
namespace Rpc {
 
   using namespace Util;
 
   /*
   * Constructor.
   */
   template <int D>
   RealMove<D>::RealMove(McSimulator<D>& simulator) 
    : McMove<D>(simulator),
      isAllocated_(false)
   { setClassName("RealMove"); }
 
   /*
   * Destructor, empty default implementation.
   */
   template <int D>
   RealMove<D>::~RealMove()
   {}
 
   /*
   * ReadParameters, empty default implementation.
   */
   template <int D>
   void RealMove<D>::readParameters(std::istream &in)
   {
 
      //Read the probability
      readProbability(in);
 
      // attampt move range [A, -A]
      read(in, "A", stepSize_);
   }
   
   template <int D>
   void RealMove<D>::setup()
   {  
      McMove<D>::setup();
      const int nMonomer = system().mixture().nMonomer();
      IntVec<D> const & dimensions = system().domain().mesh().dimensions();
      if (!isAllocated_){
         wFieldTmp_.allocate(nMonomer);
         for (int i = 0; i < nMonomer; ++i) {
            wFieldTmp_[i].allocate(dimensions);
         }
         isAllocated_ = true;
      }
   }
   
   
   
   /*
   * Attempt unconstrained move
   */
   template <int D>
   void RealMove<D>::attemptMove()
   {
      const int nMonomer = system().mixture().nMonomer();
      const int meshSize = system().domain().mesh().size();
      for (int i = 0; i < nMonomer; i++){
         for (int k = 0; k < meshSize; k++){
            //Random number generator
            double r = random().uniform(-stepSize_,stepSize_);
            wFieldTmp_[i][k] = system().w().rgrid()[i][k] + r;
         }
      }
      system().setWRGrid(wFieldTmp_);
   }
 
 
   /*
   * Trivial default implementation - do nothing
   */
   template <int D>
   void RealMove<D>::output()
   {}
   
   template<int D>
   void RealMove<D>::outputTimers(std::ostream& out)
   {
      out << "\n";
      out << "RealMove time contributions:\n";
      McMove<D>::outputTimers(out);
   }
 
}
}
#endif