fts/brownian/ExplicitBdStep.h

#ifndef RP_EXPLICIT_BD_STEP_H
#define RP_EXPLICIT_BD_STEP_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 "BdStep.h"                    // base class
#include <util/containers/DArray.h>    // member
 
namespace Pscf {
namespace Rp {
 
   using namespace Util;

   template <int D, class T>
   class ExplicitBdStep : public T::BdStep
   {
 
   public:
 
      // Protected constructor and destructor (see below).
 
      // Not copyable or assignable
      ExplicitBdStep(ExplicitBdStep<D,T> const & ) = delete;
      ExplicitBdStep<D,T>& operator= (ExplicitBdStep<D,T> const &) = delete;

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

      void setup() override;

      bool step() override;
 
   protected:

      ExplicitBdStep(typename T::BdSimulator& simulator);
 
      ~ExplicitBdStep() = default;
 
      using BdStepT = typename T::BdStep;
      using BdStepT::system;
      using BdStepT::simulator;
      using BdStepT::vecRandom;
 
   private:
 
      // Local copy of w fields
      DArray< typename T::RField > w_;
 
      // Change in one component of wc
      typename T::RField dwc_;
 
      // Normal distributed random numbers
      typename T::RField gaussianField_;
 
      // Prefactor of -dc_ in deterministic drift term
      double mobility_;
 
   };
 
}
}
#endif