rpg/fts/montecarlo/RealMove.h

#ifndef RPG_REAL_MOVE_H
#define RPG_REAL_MOVE_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
#include <prdc/cuda/RField.h>
 
#include <util/param/ParamComposite.h>
#include <util/containers/DArray.h>
#include <util/global.h>
 
namespace Pscf {
namespace Rpg {
 
   using namespace Util;
   using namespace Pscf::Prdc;
   using namespace Pscf::Prdc::Cuda;
 
   template <int D>
   class RealMove : public McMove<D>
   {
 
   public:
 
      RealMove(McSimulator<D>& simulator);
 
      ~RealMove();
 
      void readParameters(std::istream &in);
      
      void output();
      
      void setup();
      
      void outputTimers(std::ostream& out);
      
      // Inherited public member function
      using McMove<D>::move;
      using McMove<D>::readProbability;
      using McMove<D>::clearTimers;
      using ParamComposite::read;
      using ParamComposite::setClassName;
 
   protected:
      
      using McMove<D>::system;
      using McMove<D>::simulator;
      using McMove<D>::random;
      using McMove<D>::cudaRandom;
 
      void attemptMove();
 
   private:
            
      // Change in one component of wc
      RField<D> dwc_;
      
      // New field values
      DArray< RField<D> > w_;
      
      // Normal-distributed random fields
      RField<D> gaussianField_;
      
      // The standard deviation of the Gaussian distribution
      double sigma_;
      
      // Has the variable been allocated?
      bool isAllocated_;
      
   };
      
   #ifndef RPG_REAL_MOVE_TPP
   // Suppress implicit instantiation
   extern template class RealMove<1>;
   extern template class RealMove<2>;
   extern template class RealMove<3>;
   #endif
 
}
}
#endif