rpg/fts/analyzer/FourthOrderParameter.h

#ifndef RPG_FOURTH_ORDER_PARAMETER_H
#define RPG_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 "AverageAnalyzer.h"                      // base class
 
#include <prdc/cuda/RField.h>                     // member
#include <prdc/cuda/RFieldDft.h>                  // member
#include <pscf/math/IntVec.h>                     // member
 
 
namespace Pscf {
namespace Rpg {
 
   template <int D> class System;
   template <int D> class Simulator;
 
   using namespace Util;
   using namespace Pscf::Prdc;
   using namespace Pscf::Prdc::Cuda;

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

      FourthOrderParameter(Simulator<D>& simulator, System<D>& system);

      virtual ~FourthOrderParameter();
   
      virtual void setup();
      
      virtual double compute();
      
      virtual void outputValue(int step, double value);
            
      void computePrefactor();
      
      using AverageAnalyzer<D>::readParameters;
      using AverageAnalyzer<D>::nSamplePerOutput;
      using AverageAnalyzer<D>::setup;
      using AverageAnalyzer<D>::sample;
      using AverageAnalyzer<D>::output;
   
   protected:
 
      using AverageAnalyzer<D>::simulator;
      using AverageAnalyzer<D>::system;
      using AverageAnalyzer<D>::outputFile_;
      using ParamComposite::setClassName;
      
   private:
      
      int  kSize_;
      
      IntVec<D> kMeshDimensions_;
      
      bool isInitialized_;
      
      RFieldDft<D> wK_;
      
      RField<D> wc0_;
      
      RField<D> prefactor_;
      
      double FourthOrderParameter_;
 
   };
   
   #ifndef RPG_FOURTH_ORDER_PARAMETER_TPP
   // Suppress implicit instantiation
   extern template class FourthOrderParameter<1>;
   extern template class FourthOrderParameter<2>;
   extern template class FourthOrderParameter<3>;
   #endif
 
}
}
#endif