MCWF_TrajectoryLogger.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_MCWF_TRAJECTORYLOGGER_H_INCLUDED
#define CPPQEDCORE_QUANTUMTRAJECTORY_MCWF_TRAJECTORYLOGGER_H_INCLUDED

#include "MCWF_TrajectoryFwd.h"

#include "Archive.h"

#include "core_config.h"

#ifndef DO_NOT_USE_BOOST_SERIALIZATION
#include <boost/serialization/list.hpp>
#include <boost/serialization/utility.hpp>
#endif // DO_NOT_USE_BOOST_SERIALIZATION

#include <iosfwd>
#include <list>


namespace quantumtrajectory {

namespace ensemble {

using namespace mcwf;

typedef std::list<Logger> LoggerList;


std::ostream& displayLog(std::ostream&, const LoggerList&, size_t nBins, size_t nJumpsPerBin);
  
} // ensemble


namespace mcwf {


class Logger
{
public:
  typedef std::list<std::pair<double,size_t> > MCWF_Trajectory; 

  Logger(int logLevel, 
              bool isHamiltonian, 
              size_t nLindblads 
             );

  void step(); 

  void stepBack(double dp, double    dtDid, double newDtTry, double t, bool logControl); 
  void overshot(double dp, double oldDtTry, double newDtTry          , bool logControl); 

  void processNorm(double norm); 

  void jumpOccured(double t, size_t lindbladNo); 

  void logFailedSteps(size_t); 

  void hamiltonianCalled(); 

  std::ostream& onEnd(std::ostream&) const; 
  
private:
#ifndef DO_NOT_USE_BOOST_SERIALIZATION
  friend class boost::serialization::access;
  template<class Archive>
  void serialize(Archive& ar, const unsigned int) {ar & nSteps_ & nOvershot_ & nToleranceOvershot_ & nFailedSteps_ & nHamiltonianCalls_
                                                      & dpMaxOvershoot_ & dpToleranceMaxOvershoot_ & normMaxDeviation_
                                                      & traj_;}
#endif // DO_NOT_USE_BOOST_SERIALIZATION

  friend std::ostream& ensemble::displayLog(std::ostream&, const ensemble::LoggerList&, size_t, size_t);
  
  const int logLevel_;
  const bool isHamiltonian_;
  const size_t nLindblads_;

  size_t nSteps_, nOvershot_, nToleranceOvershot_, nFailedSteps_, nHamiltonianCalls_;
  double dpMaxOvershoot_, dpToleranceMaxOvershoot_, normMaxDeviation_;

  MCWF_Trajectory traj_;
  
};


} } // quantumtrajectory::mcwf


#endif // CPPQEDCORE_QUANTUMTRAJECTORY_MCWF_TRAJECTORYLOGGER_H_INCLUDED