QuantumTrajectory.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)
// -*- C++ -*-
#ifndef CPPQEDCORE_QUANTUMTRAJECTORY_QUANTUMTRAJECTORY_H_INCLUDED
#define CPPQEDCORE_QUANTUMTRAJECTORY_QUANTUMTRAJECTORY_H_INCLUDED

#include "QuantumSystem.h"
#include "Structure.h"

#include "Trajectory.h"


namespace quantumtrajectory {


template<int RANK>
inline double initialTimeStep(typename structure::QuantumSystem<RANK>::Ptr qs)
{
  return trajectory::initialTimeStep(qs->highestFrequency());
}


template<int RANK, typename BASE>
class QuantumTrajectory : public BASE
{
protected:
  typedef structure::QuantumSystemWrapper<RANK,true> QuantumSystemWrapper;

  template<typename... ArgumentPack>
  QuantumTrajectory(typename structure::QuantumSystem<RANK>::Ptr qs, bool isNoisy, ArgumentPack&&... argumentPack) 
    : BASE(std::forward<ArgumentPack>(argumentPack)...), t0_(0), qs_(qs, isNoisy) {};

  const QuantumSystemWrapper getQSW() const {return qs_;}

  cpputils::iarchive&  readStateMore_v(cpputils::iarchive &iar) override
    { BASE::readStateMore_v(iar); if (getQSW().getEx()) setT0(); return iar; }


  double getT0() const {return t0_;}
  
  void setT0(double t0) const {t0_=t0;}
  
  void setT0() const {t0_=BASE::getTime();}

  template<typename CONSTRUCT>
  void checkDimension(const CONSTRUCT& construct) const {if (construct!=*qs_.getQS()) throw DimensionalityMismatchException();}
  
private:
  mutable double t0_;
  const QuantumSystemWrapper qs_;

};



} // quantumtrajectory

#endif // CPPQEDCORE_QUANTUMTRAJECTORY_QUANTUMTRAJECTORY_H_INCLUDED