fts/analyzer/MaxOrderParameter.h

#ifndef RP_MAX_ORDER_PARAMETER_H
#define RP_MAX_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 "AverageAnalyzer.h"                      // base class
#include <pscf/math/IntVec.h>                     // member
#include <util/containers/DArray.h>               // member
#include <iostream>
 
namespace Pscf {
namespace Rp {
 
   using namespace Util;

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

      void setup() override;
 
   protected:
 
   protected:

      typename T::RField psi_;

      double maxPsi_;

      IntVec<D> Gmax_;

      int  kSize_;

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

      ~MaxOrderParameter() = default;

      void outputValue(int step, double value) override;

      void computePsi();

      void findMaximum(Array<typename T::Real> const & psi);
 
      // Alias for base class.
      using AverageAnalyzerT = typename T::AverageAnalyzer;
 
      // Inherited protected member variables (selected).
      using AverageAnalyzerT::simulator;
      using AverageAnalyzerT::system;
 
   private:

      typename T::RFieldDft wK_;

      IntVec<D> meshDimensions_;

      IntVec<D> kMeshDimensions_;
 
      using FFTT = typename T::FFT;
   };
 
}
}
#endif