fts/analyzer/FourthOrderParameter.h

#ifndef RP_FOURTH_ORDER_PARAMETER_H
#define RP_FOURTH_ORDER_PARAMETER_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 <pscf/math/IntVec.h>                     // member
 
// Forward declaration
namespace Util {
   template <typename T> class Array;
}
 
namespace Pscf {
namespace Rp {
 
   using namespace Util;

   template <int D, class T>
   class FourthOrderParameter : public T::AverageAnalyzer
   {
 
   public:

      void setup() override;
 
   protected:

      typename T::RField prefactor_;

      int  kSize_;

      FourthOrderParameter(typename T::Simulator& simulator, 
                           typename T::System& system);

      ~FourthOrderParameter() = default;

      double compute() override;

      void computePrefactor(Array<double>& prefactor);
 
      // Alias for base class
      using AverageAnalyzerT = typename T::AverageAnalyzer;
 
      // Inherited protected member functions (selected).
      using AverageAnalyzerT::simulator;
      using AverageAnalyzerT::system;
 
   private:

      typename T::RFieldDft wK_;

      typename T::RField psi_;

      IntVec<D> kMeshDimensions_;

      bool isInitialized_;

      virtual void computePrefactor() = 0;
 
      using FFTT = typename T::FFT;
 
   };
 
}
}
#endif