Liouvillean.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_STRUCTURE_LIOUVILLEAN_H_INCLUDED
#define CPPQEDCORE_STRUCTURE_LIOUVILLEAN_H_INCLUDED

#include "LiouvilleanFwd.h"

#include "LiouvilleanAveragedCommon.h"
#include "TimeDependence.h"

#include "LazyDensityOperator.h"
#include "StateVector.h"


namespace structure {


template<int RANK>
class Liouvillean : public quantumdata::Types<RANK,LiouvilleanAveragedCommonRanked<RANK> >
{
public:
  static const int N_RANK=RANK;

  typedef boost::shared_ptr<const Liouvillean> Ptr;

private:
  typedef quantumdata::Types<RANK,LiouvilleanAveragedCommonRanked<RANK> > Base;

public:
  typedef quantumdata::StateVector<RANK> StateVector;
  
  typedef typename Base::    StateVectorLow     StateVectorLow;
  typedef typename Base::DensityOperatorLow DensityOperatorLow;

  typedef typename Base::LazyDensityOperator LazyDensityOperator;

  typedef typename Base::DArray1D Rates; 

  virtual ~Liouvillean() {}
  

  const Rates rates(double t, const StateVector& psi) const {return Base::average(t,psi);}

  void actWithJ(double t,            
                StateVectorLow& psi, 
                size_t m             
               ) const {return actWithJ_v(t,psi,m);}

private:
  virtual void actWithJ_v(double, StateVectorLow&, size_t) const = 0;

};



template<int RANK, bool IS_TIME_DEPENDENT>
class LiouvilleanTimeDependenceDispatched : public Liouvillean<RANK>
{
public:
  typedef typename Liouvillean<RANK>::StateVectorLow      StateVectorLow     ;
  typedef typename Liouvillean<RANK>::LazyDensityOperator LazyDensityOperator;
  typedef typename Liouvillean<RANK>::Rates               Rates              ;
  
  typedef typename time::DispatcherIsTimeDependent<IS_TIME_DEPENDENT>::type Time;

private:
  virtual void         actWithJ_v(double t, StateVectorLow& psi, size_t lindbladNo) const final {actWithJ_v(Time(t),psi,lindbladNo);}   
  virtual const Rates   average_v(double t, const LazyDensityOperator&  matrix) const final {return rates_v(Time(t),matrix);} 

  virtual void        actWithJ_v(Time, StateVectorLow&, size_t   ) const = 0;
  virtual const Rates    rates_v(Time, const LazyDensityOperator&) const = 0;

};


} // structure

#endif // CPPQEDCORE_STRUCTURE_LIOUVILLEAN_H_INCLUDED