fts/montecarlo/ForceBiasMove.h

#ifndef RP_FORCE_BIAS_MOVE_H
#define RP_FORCE_BIAS_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/containers/DArray.h>          // member
#include <iostream>
 
namespace Pscf {
namespace Rp {
 
   using namespace Util;

   template <int D, class T>
   class ForceBiasMove : public T::McMove
   {
 
   public:
 
      // Protected constructor and destructor (see below).

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

      void setup() override;

      bool move() override;

      void output() override;

      bool needsDc() override;
 
   protected:

      ForceBiasMove(typename T::McSimulator& simulator);

      ~ForceBiasMove() = default;

      using McMoveT = typename T::McMove;
 
      // Protected inherited member functions (selected).
      using McMoveT::system;
      using McMoveT::simulator;
 
   private:

      using RFieldT = typename T::RField;

      DArray< RFieldT > w_;

      DArray< RFieldT > dc_;

      DArray<RFieldT >  dwc_;

      RFieldT eta_;

      RFieldT biasField_;

      double mobility_;

      virtual
      void computeForceBias(RFieldT& result,
                            RFieldT const & di,
                            RFieldT const & df,
                            RFieldT const & dwc,
                            double mobility) = 0;
 
   };
 
   // Public inline methods
 
   /*
   * Specify if dc fields need to be saved.
   */
   template <int D, class T>
   inline bool ForceBiasMove<D,T>::needsDc()
   {  return true; }
 
}
}
#endif