mBodyOperatorSpace< ScalarType > Class Template Reference

#include <OperatorSpace.hpp>

Inheritance diagram for mBodyOperatorSpace< ScalarType >:

Collaboration diagram for mBodyOperatorSpace< ScalarType >:

Public Types
using	Scalar = ScalarType
using	Base = ManyBodyOperatorSpaceBase< mBodyOperatorSpace< ScalarType > >
using	BaseSpace = typename HilbertSpaceTraits< mBodyOperatorSpace< ScalarType > >::BaseSpace
using	LocalSpace = typename HilbertSpaceTraits< mBodyOperatorSpace< ScalarType > >::LocalSpace
Public Types inherited from ManyBodyOperatorSpaceBase< mBodyOperatorSpace< ScalarType > >
using	BaseSpace = typename HilbertSpaceTraits< mBodyOperatorSpace< ScalarType > >::BaseSpace
using	LocalSpace = typename HilbertSpaceTraits< mBodyOperatorSpace< ScalarType > >::LocalSpace
Public Types inherited from OperatorSpaceBase< Derived >
typedef Eigen::NumTraits< ScalarType >::Real	Real

Public Member Functions
__host__ __device__	mBodyOperatorSpace (int m, BaseSpace const &baseSpace)
	Construct an mBodyOperatorSpace object by copying a BaseSpace object.
__host__ __device__	mBodyOperatorSpace (int m, int sysSize, int dimLoc)
	Construct an mBodyOperatorSpace object by copying a BaseSpace object.
__host__ __device__	mBodyOperatorSpace ()=default
	Default constructor.
__host__ __device__	mBodyOperatorSpace (mBodyOperatorSpace const &other)=default
	Copy constructor.
__host__ __device__	mBodyOperatorSpace (mBodyOperatorSpace &&other)=default
	Move constructor.
__host__ __device__	~mBodyOperatorSpace ()=default
	Destructor.
Public Member Functions inherited from ManyBodyOperatorSpaceBase< mBodyOperatorSpace< ScalarType > >
__host__ __device__	ManyBodyOperatorSpaceBase (BaseSpace const &baseSpace, int systemSize, LocalSpace const &locSpace)
	Construct a ManyBodyOperatorSpaceBase object by copying.
__host__ __device__	ManyBodyOperatorSpaceBase (BaseSpace &&baseSpace, int systemSize, LocalSpace &&locSpace)
	Construct a ManyBodyOperatorSpaceBase object by moving.
__host__ __device__	ManyBodyOperatorSpaceBase ()=default
	Default constructor.
__host__ __device__	ManyBodyOperatorSpaceBase (ManyBodyOperatorSpaceBase const &other)=default
	Copy constructor.
__host__ __device__	ManyBodyOperatorSpaceBase (ManyBodyOperatorSpaceBase &&other)=default
	Move constructor.
__host__ __device__	~ManyBodyOperatorSpaceBase ()=default
	Destructor.
__host__ __device__ ManyBodyOperatorSpaceBase &	operator= (ManyBodyOperatorSpaceBase const &other)=default
	Copy assignment operator.
__host__ __device__ ManyBodyOperatorSpaceBase &	operator= (ManyBodyOperatorSpaceBase &&other)=default
	Move assignment operator.
Public Member Functions inherited from OperatorSpaceBase< Derived >
__host__ __device__	OperatorSpaceBase (BaseSpace const &baseSpace)
	Construct a OperatorSpaceBase object by copying a BaseSpace object.
__host__ __device__	OperatorSpaceBase (BaseSpace &&baseSpace)
	Construct a OperatorSpaceBase object by moving a BaseSpace object.
__host__ __device__	OperatorSpaceBase ()=default
	Default constructor.
__host__ __device__	OperatorSpaceBase (OperatorSpaceBase const &other)=default
	Copy constructor.
__host__ __device__	OperatorSpaceBase (OperatorSpaceBase &&other)=default
	Move constructor.
__host__ __device__	~OperatorSpaceBase ()=default
	Destructor.
__host__ __device__ OperatorSpaceBase &	operator= (OperatorSpaceBase const &other)=default
	Copy assignment operator.
__host__ __device__ OperatorSpaceBase &	operator= (OperatorSpaceBase &&other)=default
	Move assignment operator.
__host__ __device__ int	dim () const
__host__ __device__ BaseSpace const &	baseSpace () const
__host__ __device__ int	baseDim () const
template<typename RealType >
__host__ __device__ void	action (int &resBasisNum, Complex_t< RealType > &coeff, int opNum, int basisNum) const
template<typename RealType , class EigenDerived >
__host__ __device__ void	action (int &resBasisNum, Complex_t< RealType > &coeff, int opNum, int basisNum, Eigen::DenseBase< EigenDerived > &config) const
__host__ __device__ std::pair< int, Complex_t< Real > >	action (int opNum, int basisNum) const
__host__ void	basisOp (Eigen::SparseMatrix< ScalarType > &res, int opNum) const
__host__ Eigen::SparseMatrix< ScalarType >	basisOp (int opNum) const
Public Member Functions inherited from ManyBodySpaceBase< Derived >
__host__ __device__	ManyBodySpaceBase (int systemSize, LocSpace_t const &locSpace)
	Constructor1.
__host__ __device__	ManyBodySpaceBase (int systemSize, LocSpace_t &&locSpace)
	Constructor2.
template<typename... Args>
__host__ __device__	ManyBodySpaceBase (int systemSize, Args... args)
	Constructor3.
__host__ __device__	ManyBodySpaceBase ()=default
	Default constructor.
__host__ __device__	ManyBodySpaceBase (ManyBodySpaceBase const &other)
	Copy constructor.
__host__ __device__	ManyBodySpaceBase (ManyBodySpaceBase &&other)
	Move constructor.
__host__ __device__	~ManyBodySpaceBase ()=default
	Destructor.
__host__ __device__ ManyBodySpaceBase &	operator= (ManyBodySpaceBase const &other)=delete
	Copy assignment operator (deleted)
__host__ __device__ Derived &	operator= (ManyBodySpaceBase &&other)
	Move assignment operator.
__host__ __device__ bool	operator== (ManyBodySpaceBase const &other) const
	Equality operator.
__host__ __device__ LocSpace_t const &	locSpace () const
__host__ __device__ int	sysSize () const
__host__ __device__ int	dimLoc () const
__host__ __device__ int	locState (int state, int pos) const
__host__ Derived const &	printInString (int basisNum) const
	Print a basis state specified by the input (basisNum) in a string of integers.
__host__ __device__ Eigen::RowVectorXi	ordinalToConfig (int basisNum) const
template<class EigenDerived >
__host__ __device__ int	configToOrdinal (Eigen::DenseBase< EigenDerived > &config) const
__host__ __device__ int	reverse (int state) const
	Reverse the input state.
__host__ __device__ Vector_t const &	parityPair () const
__host__ __device__ int	parityPair (int state) const
__host__ __device__ void	computeParityPair () const
__host__ __device__ int	translate (int state, int trans) const
	Translate the input state to the left by one.
template<class EigenDerived >
__host__ __device__ int	translate (int state, int trans, Eigen::DenseBase< EigenDerived > &work) const
__host__ __device__ int	transEqDim () const
__host__ __device__ Vector_t const &	transEqClassRep () const
__host__ __device__ int	transEqClassRep (int id) const
__host__ __device__ int	transEqClassRep (int id, int trans) const
__host__ __device__ Vector_t const &	transPeriod () const
__host__ __device__ int	transPeriod (int id) const
__host__ __device__ Vector_t const &	stateToTransEqClass () const
__host__ __device__ int	stateToTransEqClass (int state) const
__host__ __device__ int	stateToTransPeriod (int state) const
__host__ __device__ void	computeTransEqClass () const
Public Member Functions inherited from HilbertSpace< Derived >
__host__ __device__ int	dim () const
__host__ __device__ bool	operator== (HilbertSpace const &other) const

Private Types
template<class T >
using	Vector_t = Eigen::VectorX< T >
template<class T >
using	Array_t = Eigen::ArrayX< T >

Private Attributes
int	m_mBody = 0
int	m_opConfDim = 0
IntegerComposition	m_ActingSites
BaseConverter< int >	m_opConfig
Array_t< Array_t< int > >	m_config

Operator overloads
friend	HilbertSpace< mBodyOperatorSpace >
friend	OperatorSpaceBase< mBodyOperatorSpace >
friend	ManyBodySpaceBase< mBodyOperatorSpace >
__host__ __device__ mBodyOperatorSpace &	operator= (mBodyOperatorSpace const &other)=default
	Copy assignment operator.
__host__ __device__ mBodyOperatorSpace &	operator= (mBodyOperatorSpace &&other)=default
	Move assignment operator.
__host__ __device__ int	m () const
__host__ __device__ int	dim_impl () const
template<typename RealType >
__host__ __device__ void	action_impl (int &resBasisNum, Complex_t< RealType > &coeff, int opOrdinal, int basisNum) const
template<typename RealType , class Derived >
__host__ __device__ void	action_impl (int &resBasisNum, Complex_t< RealType > &coeff, int opOrdinal, int basisNum, Eigen::DenseBase< Derived > &config) const
__host__ __device__ int	locState_impl (int opOrdinal, int pos) const
__host__ __device__ int	translate_impl (int opOrdinal, int trans) const
template<class EigenDerived >
__host__ __device__ int	translate_impl (int opOrdinal, int trans, Eigen::DenseBase< EigenDerived > &workSpace) const
__host__ __device__ Eigen::RowVectorXi	ordinalToConfig_impl (int basisNum) const
template<class EigenDerived >
__host__ __device__ int	configToOrdinal_impl (Eigen::DenseBase< EigenDerived > &config) const

Additional Inherited Members
Protected Types inherited from ManyBodySpaceBase< Derived >
using	Vector_t = Eigen::VectorX< int >
Protected Attributes inherited from ManyBodySpaceBase< Derived >
int	m_systemSize = 0
LocSpace_t	m_locSpace
int	m_transEqDim = -1
Vector_t	m_transEqClassRep
Vector_t	m_transPeriod
Vector_t	m_stateToTransEqClass
Vector_t	m_parityPair

Member Typedef Documentation

Array_t

template<typename ScalarType >

template<class T >

using mBodyOperatorSpace< ScalarType >::Array_t = Eigen::ArrayX<T>

private

Base

template<typename ScalarType >

using mBodyOperatorSpace< ScalarType >::Base = ManyBodyOperatorSpaceBase< mBodyOperatorSpace<ScalarType> >

BaseSpace

template<typename ScalarType >

using mBodyOperatorSpace< ScalarType >::BaseSpace = typename HilbertSpaceTraits< mBodyOperatorSpace<ScalarType> >::BaseSpace

LocalSpace

template<typename ScalarType >

using mBodyOperatorSpace< ScalarType >::LocalSpace = typename HilbertSpaceTraits< mBodyOperatorSpace<ScalarType> >::LocalSpace

Scalar

template<typename ScalarType >

using mBodyOperatorSpace< ScalarType >::Scalar = ScalarType

Vector_t

template<typename ScalarType >

template<class T >

using mBodyOperatorSpace< ScalarType >::Vector_t = Eigen::VectorX<T>

private

Constructor & Destructor Documentation

mBodyOperatorSpace() [1/5]

template<typename ScalarType >

__host__ __device__ mBodyOperatorSpace< ScalarType >::mBodyOperatorSpace ( int  m, BaseSpace const &  baseSpace  )

inline

Construct an mBodyOperatorSpace object by copying a BaseSpace object.

            : Base(baseSpace, baseSpace.sysSize(), OperatorSpace<ScalarType>(baseSpace.locSpace())),
              m_mBody{m},
              m_opConfDim(powi(baseSpace.dimLoc() * baseSpace.dimLoc() - 1, m)),
              m_ActingSites(m, baseSpace.sysSize(), 1),
              m_opConfig(baseSpace.dimLoc() * baseSpace.dimLoc() - 1, m) {
            debug_constructor_printf(1);
#ifndef __CUDA_ARCH__
            if(m > baseSpace.sysSize()) {
                std::cerr << "ERROR: " << __PRETTY_FUNCTION__
                          << "\nMessage:\t m(" << m << ") > "
                          << "baseSpace.sysSize()(" << baseSpace.sysSize() << ")" << std::endl;
                std::exit(EXIT_FAILURE);
            }
            m_config.resize(get_max_threads());
            for(int j = 0; j != get_max_threads(); ++j) m_config[j].resize(this->sysSize());
#endif
        };

mBodyOperatorSpace() [2/5]

template<typename ScalarType >

__host__ __device__ mBodyOperatorSpace< ScalarType >::mBodyOperatorSpace ( int  m, int  sysSize, int  dimLoc  )

inline

Construct an mBodyOperatorSpace object by copying a BaseSpace object.

            : mBodyOperatorSpace(m, BaseSpace(sysSize, dimLoc)) {
            debug_constructor_printf(2);
        };

mBodyOperatorSpace() [3/5]

template<typename ScalarType >

__host__ __device__ mBodyOperatorSpace< ScalarType >::mBodyOperatorSpace ( )

default

Default constructor.

mBodyOperatorSpace() [4/5]

template<typename ScalarType >

__host__ __device__ mBodyOperatorSpace< ScalarType >::mBodyOperatorSpace ( mBodyOperatorSpace< ScalarType > const &  other )

default

Copy constructor.

mBodyOperatorSpace() [5/5]

template<typename ScalarType >

__host__ __device__ mBodyOperatorSpace< ScalarType >::mBodyOperatorSpace ( mBodyOperatorSpace< ScalarType > &&  other )

default

Move constructor.

~mBodyOperatorSpace()

template<typename ScalarType >

__host__ __device__ mBodyOperatorSpace< ScalarType >::~mBodyOperatorSpace ( )

default

Destructor.

Member Function Documentation

action_impl() [1/2]

template<typename ScalarType >

template<typename RealType >

__host__ __device__ void mBodyOperatorSpace< ScalarType >::action_impl ( int &  resBasisNum, Complex_t< RealType > &  coeff, int  opOrdinal, int  basisNum  ) const

inlineprivate

                                                                                {
#ifndef __CUDA_ARCH__
            action_impl(resBasisNum, coeff, opOrdinal, basisNum, m_config[get_thread_num()]);
#else
            Eigen::ArrayXi config(this->sysSize());
            action_impl(resBasisNum, coeff, opOrdinal, basisNum, config);
#endif
        };

action_impl() [2/2]

template<typename ScalarType >

template<typename RealType , class Derived >

__host__ __device__ void mBodyOperatorSpace< ScalarType >::action_impl ( int &  resBasisNum, Complex_t< RealType > &  coeff, int  opOrdinal, int  basisNum, Eigen::DenseBase< Derived > &  config  ) const

private

                                              {
// debug_print(__PRETTY_FUNCTION__);
#ifndef NDEBUG
    if(workSpace.size() < this->sysSize())
        debug_printf("%s\n\tworkSpace.size()=%d\n", __PRETTY_FUNCTION__, (int)workSpace.size());
#endif
    assert(workSpace.size() >= this->sysSize());
    resBasisNum = basisNum;
    coeff       = {1.0, 0};
    int                 resLocBasisNum;
    Complex_t<RealType> resLocCoeff;
 
    int posConfNum = opOrdinal / m_opConfDim;
    int opConfNum  = opOrdinal % m_opConfDim;
    this->m_ActingSites.ordinalToPart(workSpace, posConfNum);
    int locOpNum, locBasisNum;
 
    int index = 0;
    int base  = 1;
    for(int pos = 0; pos != this->sysSize(); ++pos, base *= this->baseSpace().dimLoc()) {
        locBasisNum = this->baseSpace().locState(basisNum, pos);
        if(workSpace[pos] == 0) {
            this->locSpace().action(resLocBasisNum, resLocCoeff, 0, locBasisNum);
            coeff *= resLocCoeff;
            continue;
        }
        locOpNum = 1 + this->m_opConfig.digit(opConfNum, index++);
#ifndef NDEBUG
        if(locOpNum >= 4) {
            printf("Error: mBodyOperatorSpace::%s Out of bounds locOpNum = %d\n\n", __func__,
                   locOpNum);
    #ifndef __CUDA_ARCH__
            std::cout << "locOpNum = " << locOpNum << std::endl;
            this->m_opConfig.status();
            std::exit(EXIT_FAILURE);
    #endif
        }
#endif
 
        this->locSpace().action(resLocBasisNum, resLocCoeff, locOpNum, locBasisNum);
        if(abs(resLocCoeff) < 1.0e-8) {
            coeff = {0, 0};
            return;
        }
        coeff *= resLocCoeff;
        resBasisNum += (resLocBasisNum - locBasisNum) * base;
    }
    return;
}

configToOrdinal_impl()

template<typename ScalarType >

template<class EigenDerived >

__host__ __device__ int mBodyOperatorSpace< ScalarType >::configToOrdinal_impl ( Eigen::DenseBase< EigenDerived > &  config ) const

private

dim_impl()

template<typename ScalarType >

__host__ __device__ int mBodyOperatorSpace< ScalarType >::dim_impl ( ) const

inlineprivate

{ return m_ActingSites.dim() * m_opConfDim; }

locState_impl()

template<typename ScalarType >

__host__ __device__ int mBodyOperatorSpace< ScalarType >::locState_impl ( int  opOrdinal, int  pos  ) const

inlineprivate

                                                                            {
            int posConfNum = opOrdinal / m_opConfDim;
            if(this->m_ActingSites.locNumber(posConfNum, pos) == 0)
                return 0;
            else {
                int index = -1;
                for(int l = 0; l <= (pos % this->m_systemSize); ++l) {
                    index += this->m_ActingSites.locNumber(posConfNum, l);
                }
                int opConfNum = opOrdinal % m_opConfDim;
                return this->m_opConfig.digit(opConfNum, index) + 1;
            }
        }

m()

template<typename ScalarType >

__host__ __device__ int mBodyOperatorSpace< ScalarType >::m ( ) const

inline

{ return m_mBody; }

operator=() [1/2]

template<typename ScalarType >

__host__ __device__ mBodyOperatorSpace & mBodyOperatorSpace< ScalarType >::operator= ( mBodyOperatorSpace< ScalarType > &&  other )

default

Move assignment operator.

Parameters

other

operator=() [2/2]

template<typename ScalarType >

__host__ __device__ mBodyOperatorSpace & mBodyOperatorSpace< ScalarType >::operator= ( mBodyOperatorSpace< ScalarType > const &  other )

default

Copy assignment operator.

Parameters

other

ordinalToConfig_impl()

template<typename ScalarType >

__host__ __device__ Eigen::RowVectorXi mBodyOperatorSpace< ScalarType >::ordinalToConfig_impl ( int  basisNum ) const

private

                        {
    int posConfNum = basisNum / m_opConfDim;
    int opConfNum  = basisNum % m_opConfDim;
 
    Eigen::RowVectorXi config(this->sysSize());
    this->m_ActingSites.ordinalToPart(config, posConfNum);
 
    int index = 0;
    for(int pos = 0; pos != this->sysSize(); ++pos) {
        if(config[pos] == 0) continue;
        config[pos] = 1 + this->m_opConfig.digit(opConfNum, index++);
    }
    return config;
}

translate_impl() [1/2]

template<typename ScalarType >

__host__ __device__ int mBodyOperatorSpace< ScalarType >::translate_impl ( int  opOrdinal, int  trans  ) const

inlineprivate

                                                                               {
            Eigen::ArrayXi work(this->sysSize());
            return this->translate(opOrdinal, trans, work);
        };

translate_impl() [2/2]

template<typename ScalarType >

template<class EigenDerived >

__host__ __device__ int mBodyOperatorSpace< ScalarType >::translate_impl ( int  opOrdinal, int  trans, Eigen::DenseBase< EigenDerived > &  workSpace  ) const

private

                                                                             {
    int posConfNum = opOrdinal / m_opConfDim;
    int opConfNum  = opOrdinal % m_opConfDim;
#ifndef __CUDA_ARCH__
    workSpace.resize(this->sysSize());
#else
    assert(workSpace.size() >= this->sysSize());
#endif
#ifndef NDEBUG
    if(workSpace.size() < this->sysSize())
        debug_printf("%s\n\tworkSpace.size()=%d\n", __PRETTY_FUNCTION__, (int)workSpace.size());
#endif
    m_ActingSites.ordinalToPart(workSpace, posConfNum);
    int opTrans = 0;
    for(int l = 0; l != trans; ++l) opTrans += workSpace[l];
 
    posConfNum = m_ActingSites.translate(posConfNum, trans, workSpace);
    opConfNum  = this->m_opConfig.shiftDigits(opConfNum, opTrans);
    return opConfNum + posConfNum * m_opConfDim;
}

Member Data Documentation

HilbertSpace< mBodyOperatorSpace >

template<typename ScalarType >

friend mBodyOperatorSpace< ScalarType >::HilbertSpace< mBodyOperatorSpace >

private

m_ActingSites

template<typename ScalarType >

IntegerComposition mBodyOperatorSpace< ScalarType >::m_ActingSites

private

m_config

template<typename ScalarType >

Array_t< Array_t<int> > mBodyOperatorSpace< ScalarType >::m_config

mutableprivate

m_mBody

template<typename ScalarType >

int mBodyOperatorSpace< ScalarType >::m_mBody = 0

private

m_opConfDim

template<typename ScalarType >

int mBodyOperatorSpace< ScalarType >::m_opConfDim = 0

private

m_opConfig

template<typename ScalarType >

BaseConverter<int> mBodyOperatorSpace< ScalarType >::m_opConfig

private

ManyBodySpaceBase< mBodyOperatorSpace >

template<typename ScalarType >

friend mBodyOperatorSpace< ScalarType >::ManyBodySpaceBase< mBodyOperatorSpace >

private

OperatorSpaceBase< mBodyOperatorSpace >

template<typename ScalarType >

friend mBodyOperatorSpace< ScalarType >::OperatorSpaceBase< mBodyOperatorSpace >

private

---

The documentation for this class was generated from the following file:

• /Users/shoki/GitHub/Locality/Headers/HilbertSpace/OperatorSpace.hpp