fts/analyzer/AverageAnalyzer.h

#ifndef RP_AVERAGE_ANALYZER_H
#define RP_AVERAGE_ANALYZER_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/accumulators/Average.h>           // member
#include <iostream>
 
namespace Pscf {
namespace Rp {
 
   using namespace Util;

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

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

      virtual void setup();

      virtual void sample(long iStep);

      virtual void output();

      int nSamplePerOutput() const;
 
   protected:

      std::ofstream outputFile_;

      Average accumulator_;
 
      // Protected member functions

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

      ~AverageAnalyzer() = default;

      virtual double compute() = 0;

      virtual void outputValue(int step, double value);
 
   private:

      int nSamplePerOutput_;
 
      using AnalyzerT = typename T::Analyzer;
 
   };
 
   // Inline functions
 
   // Get nSamplePerOutput.
   template <int D, class T> inline 
   int AverageAnalyzer<D,T>::nSamplePerOutput() const
   {  return nSamplePerOutput_; }
 
}
}
#endif