Sigma.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_SIGMA_H_INCLUDED
#define   CPPQEDCORE_QUANTUMOPERATOR_SIGMA_H_INCLUDED

#include "SigmaFwd.h"

#include "Types.h"


namespace quantumoperator {


template<int L, int R, typename OTHER>
const DirectProduct<L,R,OTHER,true >
operator*(const Sigma<L,R>&, const OTHER&);



template<int L, int R, typename OTHER>
const DirectProduct<L,R,OTHER,false>
operator*(const OTHER&, const Sigma<L,R>&);

template<int L1, int R1, int L2, int R2>
const DirectProduct<L1,R1,Sigma<L2,R2>,true>
operator*(const Sigma<L1,R1>&, const Sigma<L2,R2>&);



template<int L, int R>
class Sigma
{
public:
  static const int N_RANK=1;

  typedef quantumdata::Types<1>::StateVectorLow StateVectorLow;

  void
  apply(const StateVectorLow& psi, StateVectorLow& dpsidt) const
  {
    dpsidt(L)+=psi(R);
  }  

  const Sigma<R,L> dagger() const {return Sigma<R,L>();}

};



template<int RANK, bool IS_HEAD>
const typename quantumdata::Types<RANK-1>::StateVectorLow
partialProject(const typename quantumdata::Types<RANK>::StateVectorLow& psi, int n);

} // quantumoperator

#endif // CPPQEDCORE_QUANTUMOPERATOR_SIGMA_H_INCLUDED