fts/brownian/PredCorrBdStep.h

#ifndef RP_PRED_CORR_BD_STEP_H
#define RP_PRED_CORR_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 <util/containers/DArray.h>        // member
 
namespace Pscf {
namespace Rp {
 
   using namespace Util;

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

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

      void setup() override;

      bool step() override;
 
   protected:

      PredCorrBdStep(typename T::BdSimulator& simulator);
 
      ~PredCorrBdStep() = default;
 
      // Alias for base class.
      using BdStepT = typename T::BdStep;
 
      // Protected inherited member functions
      using BdStepT::system;
      using BdStepT::simulator;
 
   private:
 
      using RFieldT = typename T::RField;
 
      // Predicted values of fields (monomer fields)
      DArray< RFieldT > wp_;
 
      // Corrected (final) values of fields (monomer fields)
      DArray< RFieldT > wf_;
 
      // Initial deterministic forces (eigenvector components)
      DArray< RFieldT > dci_;
 
      // Random displacement components (eigenvector components)
      DArray< RFieldT > eta_;
 
      // Change in one component of wc
      RFieldT dwc_;
 
      // Change in pressure field component
      RFieldT dwp_;
 
      // Prefactor of -dc_ in deterministic drift term
      double mobility_;
 
   };
 
}
}
#endif