TridiagonalHamiltonian.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_QUANTUMOPERATOR_TRIDIAGONALHAMILTONIAN_H_INCLUDED
#define CPPQEDCORE_QUANTUMOPERATOR_TRIDIAGONALHAMILTONIAN_H_INCLUDED

#include "TridiagonalHamiltonianFwd.h"

#include "Hamiltonian.h"

#include "Tridiagonal.h"

#include <list>


namespace quantumoperator {


namespace details {

using namespace structure;

template<int RANK, bool IS_TIME_DEPENDENT>
class TDH_Base 
  : public HamiltonianTimeDependenceDispatched<RANK,
                                               mpl::if_c<IS_TIME_DEPENDENT,
                                                         mpl::integral_c<TimeDependence,ONE_TIME>,
                                                         mpl::integral_c<TimeDependence, NO_TIME>
                                                         >::type::value
                                               >
{
protected:
  typedef quantumoperator::Tridiagonal<RANK> Tridiagonal;
  typedef std::list<Tridiagonal>             Tridiagonals;

  typedef typename Hamiltonian<RANK>::StateVectorLow StateVectorLow; 

  TDH_Base(const Tridiagonals& hOverIs) : hOverIs_(hOverIs) {}

  const Tridiagonals& getH_OverIs() const {return hOverIs_;}
        Tridiagonals& getH_OverIs()       {return const_cast<Tridiagonals&>(static_cast<const TDH_Base*>(this)->getH_OverIs());}

private:

  void addContribution_v(OneTime, const StateVectorLow&, StateVectorLow&) const;
  void addContribution_v( NoTime, const StateVectorLow&, StateVectorLow&) const;


  mutable Tridiagonals hOverIs_;

};


} // details




template<int RANK, bool IS_TIME_DEPENDENT>
class TridiagonalHamiltonian : public details::TDH_Base<RANK,IS_TIME_DEPENDENT>
{
private:
  typedef details::TDH_Base<RANK,IS_TIME_DEPENDENT> Base;

public:
  typedef typename Base::Tridiagonal  Tridiagonal ;
  typedef typename Base::Tridiagonals Tridiagonals;

  TridiagonalHamiltonian(const Tridiagonals& hOverIs) : Base(hOverIs) {} 
  TridiagonalHamiltonian(const Tridiagonal & hOverI ) : Base(Tridiagonals(1,hOverI)) {} 

};

#undef BASE_class

} // quantumoperator


#endif // CPPQEDCORE_QUANTUMOPERATOR_TRIDIAGONALHAMILTONIAN_H_INCLUDED