rpc/fts/montecarlo/ForceBiasMove.h

#ifndef RPC_FORCE_BIAS_MOVE_H
#define RPC_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/cpu/RField.h>                 // member
#include <util/containers/DArray.h>          // member
 
namespace Pscf {
namespace Rpc {
 
   using namespace Util;
   using namespace Pscf::Prdc::Cpu;
 
   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>::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_;
 
      double mobility_;
 
   };
   
   // Public inline methods
 
   /*
   * Return whether dc fields need to be saved for ForceBiasMove.
   */
   template <int D>
   inline bool ForceBiasMove<D>::needsDc()
   {  return true; }
 
   #ifndef RPC_FORCE_BIAS_MOVE_TPP
   // Suppress implicit instantiation
   extern template class ForceBiasMove<1>;
   extern template class ForceBiasMove<2>;
   extern template class ForceBiasMove<3>;
   #endif
 
}
}
#endif