IntraBondTensorAutoCorr.h

#ifndef MCMD_INTRA_BOND_TENSOR_AUTO_CORR_H
#define MCMD_INTRA_BOND_TENSOR_AUTO_CORR_H

/*
* Simpatico - Simulation Package for Polymeric and Molecular Liquids
*
* Copyright 2010 - 2017, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/

#include <mcMd/analyzers/SystemAnalyzer.h>    // base class template
#include <util/accumulators/AutoCorrArray.h>  // member template 
#include <util/space/Tensor.h>                // member template parameter
#include <util/containers/DArray.h>           // member template
#include <util/archives/serialize.h>          // used in method template
#include <util/global.h>                      // used in method template


namespace Simp {
   class Species;
}

namespace McMd
{

   using namespace Util;
   using namespace Simp;

   template <class SystemType>
   class IntraBondTensorAutoCorr : public SystemAnalyzer<SystemType>
   {
   
   public:
  
      IntraBondTensorAutoCorr(SystemType &system);

      virtual ~IntraBondTensorAutoCorr();
  
      virtual void readParameters(std::istream& in);
   
      virtual void loadParameters(Serializable::IArchive &ar);
   
      virtual void save(Serializable::OArchive &ar);
  
      template <class Archive>
      void serialize(Archive& ar, const unsigned int version);

      virtual void setup();
  
      virtual void sample(long iStep);
   
      virtual void output();

   protected:

      using SystemAnalyzer<SystemType>::read;
      using SystemAnalyzer<SystemType>::readOutputFileName;
      using SystemAnalyzer<SystemType>::readInterval;
      using SystemAnalyzer<SystemType>::loadParameter;
      using SystemAnalyzer<SystemType>::isAtInterval;
      using SystemAnalyzer<SystemType>::writeParam;
      using SystemAnalyzer<SystemType>::system;
      using SystemAnalyzer<SystemType>::outputFileName;
      using SystemAnalyzer<SystemType>::fileMaster;

   private:

      std::ofstream outputFile_;

      AutoCorrArray<Tensor, double>  accumulator_;
 
      DArray<Tensor> data_;
 
      Species* speciesPtr_;
   
      int speciesId_;
   
      int nMolecule_;

      int nBond_;
 
      int capacity_;
   
      bool isInitialized_;

   };

   /*
   * Serialize to/from an archive. 
   */
   template <class SystemType>
   template <class Archive>
   void IntraBondTensorAutoCorr<SystemType>::serialize(Archive& ar, 
                                                       const unsigned int version)
   {
      Analyzer::serialize(ar, version);
      ar & speciesId_;
      ar & capacity_;

      ar & nBond_;
      ar & nMolecule_;
      ar & accumulator_;
   }

}
#endif