ddMd/integrators/NphIntegrator.h

#ifndef DDMD_NPH_INTEGRATOR_H
#define DDMD_NPH_INTEGRATOR_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 "TwoStepIntegrator.h"      // base class

namespace DdMd
{

   class Simulation;
   using namespace Util;

   class NphIntegrator : public TwoStepIntegrator
   {

   public:

      NphIntegrator(Simulation& simulation);

      ~NphIntegrator();

      void readParameters(std::istream& in);

      void loadParameters(Serializable::IArchive& ar);

      void save(Serializable::OArchive& ar);

   protected:

      void initDynamicalState();

      void setup();

      virtual void integrateStep1();

      virtual void integrateStep2();

   private:

      Vector nu_;
      double  dt_;
      double W_;
      double P_target_;
      double T_kinetic_;
      Vector P_curr_diag_;
      LatticeSystem mode_;
      unsigned int ndof_;

      Vector exp_v_fac_;
      Vector sinhx_fac_v_;
      double mtk_term_2_;
      DArray<double> prefactors_;

   };

}
#endif