TimeAveragingMCWF_Trajectory.h

Go to the documentation of this file.

// Copyright András Vukics 2006–2014. Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE.txt)
#ifndef   CPPQEDCORE_QUANTUMTRAJECTORY_TIMEAVERAGINGMCWF_TRAJECTORY_H_INCLUDED
#define   CPPQEDCORE_QUANTUMTRAJECTORY_TIMEAVERAGINGMCWF_TRAJECTORY_H_INCLUDED

#include "TimeAveragingMCWF_TrajectoryFwd.h"

#include "MCWF_Trajectory.h"


namespace quantumtrajectory {



template<int RANK>
class TimeAveragingMCWF_Trajectory : public MCWF_Trajectory<RANK>
{
private:
  typedef MCWF_Trajectory<RANK> Base;
  
public:
  typedef typename Base::StateVector    StateVector   ;
  typedef typename Base::StateVectorLow StateVectorLow;
  
  typedef typename Base::Averaged Averaged;
  
  typedef typename Averaged::Averages Averages;

  template<typename SYS>
  TimeAveragingMCWF_Trajectory(
                               StateVector& psi,
                               const SYS& sys,
                               const mcwf::Pars& p,
                               double relaxationTime, 
                               const StateVectorLow& scaleAbs=StateVectorLow()
                               )
    : Base(psi,sys,p,scaleAbs), relaxationTime_(relaxationTime), averages_(this->getQSW().template nAvr<structure::LA_Av>()), sum_(0), av_(this->getQSW().getAv())
    {
      averages_=0.;
    }

  const Averages getAverages() const {return averages_;}
  
private:
  std::ostream& display_v(std::ostream&, int) const override;
  
  std::ostream& logOnEnd_v(std::ostream& os) const override;
  
  cpputils::iarchive&  readStateMore_v(cpputils::iarchive& iar)       override {return Base:: readStateMore_v(iar) & averages_ & sum_;}
  cpputils::oarchive& writeStateMore_v(cpputils::oarchive& oar) const override {return Base::writeStateMore_v(oar) & averages_ & sum_;}
  
  const std::string trajectoryID_v() const override {return "TimeAveragingMCWF_Trajectory";}

  const double relaxationTime_;
  
  mutable Averages averages_;
  mutable long sum_;
  
  const typename Averaged::Ptr av_;

};


} // quantumtrajectory


#endif // CPPQEDCORE_QUANTUMTRAJECTORY_TIMEAVERAGINGMCWF_TRAJECTORY_H_INCLUDED