fts/montecarlo/ShiftMove.h

#ifndef RP_SHIFT_MOVE_H
#define RP_SHIFT_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 "McMove.h"                    // base class
#include <util/containers/DArray.h>    // member
#include <pscf/math/IntVec.h>          // template with defaults
 
// Forward declaration
namespace Util {
   template <typename Data> class Array;
}
 
namespace Pscf {
namespace Rp {
 
   using namespace Util;

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

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

      void setup() override;
 
   protected:

      ShiftMove(typename T::McSimulator& simulator);

      ~ShiftMove() = default;

      void attemptMove() override;
 
      virtual void shiftFields(IntVec<D> const & shift) = 0;
 
      void shiftField(Array<double>& out, 
                      Array<double> const & in,
                      IntVec<D> shift, 
                      IntVec<D> dimensions) const;

      mutable DArray< typename T::RField > w_;
 
      // Alias for McMove base class.
      using McMoveT = typename T::McMove;
 
      // Inherited protected member functions (selected).
      using McMoveT::system;
      using McMoveT::simulator;
 
   private:

      int maxShift_;

      bool isAllocated_;
 
   };
 
}
}
#endif