rpg/fts/montecarlo/ForceBiasMove.h

#ifndef RPG_FORCE_BIAS_MOVE_H
#define RPG_FORCE_BIAS_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/containers/DArray.h> 
 
namespace Pscf {
namespace Rpg {
 
   using namespace Util;
   using namespace Pscf::Prdc;
   using namespace Pscf::Prdc::Cuda;
 
   template <int D>
   class ForceBiasMove : public McMove<D>
   {
 
   public:
 
      ForceBiasMove(McSimulator<D>& simulator);
 
      ~ForceBiasMove();
 
      void readParameters(std::istream &in);
 
      void output();
 
      void setup();
 
      bool move();
 
      void outputTimers(std::ostream& out);
      
      bool needsDc();
 
      // 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;
      using McMove<D>::incrementNAttempt;
      using McMove<D>::incrementNAccept;
      using McMove<D>::incrementNFail;
 
      using McMove<D>::attemptMoveTimer_;
      using McMove<D>::compressorTimer_;
      using McMove<D>::componentTimer_;
      using McMove<D>::hamiltonianTimer_;
      using McMove<D>::decisionTimer_;
      using McMove<D>::totalTimer_;
 
   private:
 
      DArray< RField<D> > w_;
 
      DArray< RField<D> > dc_;
 
      DArray<RField<D> >  dwc_;
 
      RField<D> gaussianField_;
      
      RField<D> biasField_;
      
      double mobility_;
      
   };
   class ForceBiasMove : public McMove<D> {…};
   
   // Public inline methods
 
   /*
   * Return whether dc fields need to be saved for ForceBiasMove.
   */
   template <int D>
   inline bool ForceBiasMove<D>::needsDc()
   {  return true; }
   inline bool ForceBiasMove<D>::needsDc() {…}
 
   #ifndef RPG_FORCE_BIAS_MOVE_TPP
   // Suppress implicit instantiation
   extern template class ForceBiasMove<1>;
   extern template class ForceBiasMove<2>;
   extern template class ForceBiasMove<3>;
   #endif
 
}
}
#endif