rpg/fts/analyzer/FourthOrderParameter.h

#ifndef RPG_FOURTH_ORDER_PARAMETER_H
#define RPG_FOURTH_ORDER_PARAMETER_H
 
/*
* PSCF - Polymer Self-Consistent Field Theory
*
* Copyright 2016 - 2022, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
#include "AverageAnalyzer.h"                      // base class
 
#include <util/containers/DArray.h>               // member
#include <util/accumulators/Average.h>            // member
#include <prdc/cuda/RField.h>                     // member
#include <prdc/cuda/RFieldDft.h>                  // member
#include <map>                                    // member
 
#include <string>
#include <iostream>
 
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