ObservableEnsemble Class Reference

#include <Ensemble.hpp>

Collaboration diagram for ObservableEnsemble:

Public Types
using	MatrixType = ObjectOnGPU< Eigen::MatrixX< Scalar > >
using	EigenSolver = MatrixUtils::EigenSolver< MatrixType >
using	StateSpace = ManyBodySpinSpace
using	OperatorSpace = mBodyOperatorSpace< Scalar >
using	MatrixType = Eigen::MatrixX< Scalar >
using	EigenSolver = Eigen::SelfAdjointEigenSolver< MatrixType >
using	StateSpace = ManyBodySpinSpace
using	OperatorSpace = mBodyOperatorSpace< Scalar >
using	MatrixType = Eigen::MatrixX< Scalar >
using	EigenSolver = Eigen::SelfAdjointEigenSolver< MatrixType >
using	MatrixType = ObjectOnGPU< Eigen::MatrixX< Scalar > >
using	EigenSolver = MatrixUtils::EigenSolver< MatrixType >
using	StateSpace = ManyBodySpinSpace
using	OperatorSpace = mBodyOperatorSpace< Scalar >
using	MatrixType = Eigen::MatrixX< Scalar >
using	EigenSolver = Eigen::SelfAdjointEigenSolver< MatrixType >
using	StateSpace = ManyBodySpinSpace
using	OperatorSpace = mBodyOperatorSpace< Scalar >
using	MatrixType = ObjectOnGPU< Eigen::MatrixX< Scalar > >
using	EigenSolver = MatrixUtils::EigenSolver< MatrixType >
using	StateSpace = ManyBodySpinSpace
using	OperatorSpace = mBodyOperatorSpace< Scalar >
using	MatrixType = Eigen::MatrixX< Scalar >
using	EigenSolver = Eigen::SelfAdjointEigenSolver< MatrixType >
using	StateSpace = ManyBodySpinSpace
using	OperatorSpace = mBodyOperatorSpace< Scalar >

Public Member Functions
	ObservableEnsemble (size_t L, size_t m)
	ObservableEnsemble (size_t L, char **CL_argv)
	ObservableEnsemble ()=default
	ObservableEnsemble (const ObservableEnsemble &)=delete
ObservableEnsemble &	operator= (const ObservableEnsemble &)=delete
	ObservableEnsemble (ObservableEnsemble &&)=default
ObservableEnsemble &	operator= (ObservableEnsemble &&)=default
StateSpace &	stateSpace ()
OperatorSpace &	operatorSpace ()
void	reset ()
void	discard (size_t n)
MatrixType	sample ()
	ObservableEnsemble (size_t L, size_t m)
	ObservableEnsemble (size_t L, char **CL_argv)
	ObservableEnsemble ()=default
	ObservableEnsemble (const ObservableEnsemble &)=delete
ObservableEnsemble &	operator= (const ObservableEnsemble &)=delete
	ObservableEnsemble (ObservableEnsemble &&)=default
ObservableEnsemble &	operator= (ObservableEnsemble &&)=default
StateSpace &	stateSpace ()
OperatorSpace &	operatorSpace ()
void	reset ()
void	discard (size_t n)
MatrixType	sample ()
	ObservableEnsemble ()=default
	ObservableEnsemble (const ObservableEnsemble &)=delete
ObservableEnsemble &	operator= (const ObservableEnsemble &)=delete
	ObservableEnsemble (ObservableEnsemble &&)=default
ObservableEnsemble &	operator= (ObservableEnsemble &&)=default
void	reset ()
void	discard (size_t n)
MatrixType	sample ()
auto &	stateSpace ()
auto &	operatorSpace ()
	ObservableEnsemble (size_t L, size_t m)
	ObservableEnsemble (size_t L, char **CL_argv)
	ObservableEnsemble ()=default
	ObservableEnsemble (const ObservableEnsemble &)=delete
ObservableEnsemble &	operator= (const ObservableEnsemble &)=delete
	ObservableEnsemble (ObservableEnsemble &&)=default
ObservableEnsemble &	operator= (ObservableEnsemble &&)=default
StateSpace &	stateSpace ()
OperatorSpace &	operatorSpace ()
void	reset ()
void	discard (size_t n)
MatrixType	sample ()
	ObservableEnsemble (size_t L, size_t m)
	ObservableEnsemble (size_t L, char **CL_argv)
	ObservableEnsemble ()=default
	ObservableEnsemble (const ObservableEnsemble &)=delete
ObservableEnsemble &	operator= (const ObservableEnsemble &)=delete
	ObservableEnsemble (ObservableEnsemble &&)=default
ObservableEnsemble &	operator= (ObservableEnsemble &&)=default
StateSpace &	stateSpace ()
OperatorSpace &	operatorSpace ()
void	reset ()
void	discard (size_t n)
MatrixType	sample ()
	ObservableEnsemble (size_t L, size_t m)
	ObservableEnsemble (size_t L, char **CL_argv)
	ObservableEnsemble ()=default
	ObservableEnsemble (const ObservableEnsemble &)=delete
ObservableEnsemble &	operator= (const ObservableEnsemble &)=delete
	ObservableEnsemble (ObservableEnsemble &&)=default
ObservableEnsemble &	operator= (ObservableEnsemble &&)=default
StateSpace &	stateSpace ()
OperatorSpace &	operatorSpace ()
void	reset ()
void	discard (size_t n)
MatrixType	sample ()
	ObservableEnsemble (size_t L, size_t m)
	ObservableEnsemble (size_t L, char **CL_argv)
	ObservableEnsemble ()=default
	ObservableEnsemble (const ObservableEnsemble &)=delete
ObservableEnsemble &	operator= (const ObservableEnsemble &)=delete
	ObservableEnsemble (ObservableEnsemble &&)=default
ObservableEnsemble &	operator= (ObservableEnsemble &&)=default
StateSpace &	stateSpace ()
OperatorSpace &	operatorSpace ()
void	reset ()
void	discard (size_t n)
MatrixType	sample ()

Static Public Member Functions
static std::string	dataname (size_t rep)
static std::string	usage (size_t nargs)
static constexpr size_t	nargs (void)
static std::string	dataname (size_t rep)
static std::string	usage (size_t nargs)
static constexpr size_t	nargs (void)
static std::string	dataname (size_t rep)
static std::string	usage (size_t nargs)
static constexpr size_t	nargs (void)
static std::string	dataname (size_t rep)
static std::string	usage (size_t nargs)
static constexpr size_t	nargs (void)
static std::string	dataname (size_t rep)
static std::string	usage (size_t nargs)
static constexpr size_t	nargs (void)
static std::string	dataname (size_t rep)
static std::string	usage (size_t nargs)
static constexpr size_t	nargs (void)
static std::string	dataname (size_t rep)
static std::string	usage (size_t nargs)
static constexpr size_t	nargs (void)

Static Public Attributes
static size_t	class_M = 0

Private Attributes
size_t	m_L
size_t	m_M
StateSpace	m_stateSpace
OperatorSpace	m_opSpace
GaussianRandomMatrixGenerator_GPU< OperatorSpace >	m_mbGRME
GaussianRandomMatrixGenerator< OperatorSpace >	m_mbGRME

Member Typedef Documentation

EigenSolver [1/7]

using ObservableEnsemble::EigenSolver = MatrixUtils::EigenSolver<MatrixType>

EigenSolver [2/7]

using ObservableEnsemble::EigenSolver = Eigen::SelfAdjointEigenSolver<MatrixType>

EigenSolver [3/7]

using ObservableEnsemble::EigenSolver = Eigen::SelfAdjointEigenSolver<MatrixType>

EigenSolver [4/7]

using ObservableEnsemble::EigenSolver = MatrixUtils::EigenSolver<MatrixType>

EigenSolver [5/7]

using ObservableEnsemble::EigenSolver = Eigen::SelfAdjointEigenSolver<MatrixType>

EigenSolver [6/7]

using ObservableEnsemble::EigenSolver = MatrixUtils::EigenSolver<MatrixType>

EigenSolver [7/7]

using ObservableEnsemble::EigenSolver = Eigen::SelfAdjointEigenSolver<MatrixType>

MatrixType [1/7]

using ObservableEnsemble::MatrixType = ObjectOnGPU< Eigen::MatrixX<Scalar> >

MatrixType [2/7]

using ObservableEnsemble::MatrixType = Eigen::MatrixX<Scalar>

MatrixType [3/7]

using ObservableEnsemble::MatrixType = Eigen::MatrixX<Scalar>

MatrixType [4/7]

using ObservableEnsemble::MatrixType = ObjectOnGPU< Eigen::MatrixX<Scalar> >

MatrixType [5/7]

using ObservableEnsemble::MatrixType = Eigen::MatrixX<Scalar>

MatrixType [6/7]

using ObservableEnsemble::MatrixType = ObjectOnGPU< Eigen::MatrixX<Scalar> >

MatrixType [7/7]

using ObservableEnsemble::MatrixType = Eigen::MatrixX<Scalar>

OperatorSpace [1/6]

using ObservableEnsemble::OperatorSpace = mBodyOperatorSpace<Scalar>

OperatorSpace [2/6]

using ObservableEnsemble::OperatorSpace = mBodyOperatorSpace<Scalar>

OperatorSpace [3/6]

using ObservableEnsemble::OperatorSpace = mBodyOperatorSpace<Scalar>

OperatorSpace [4/6]

using ObservableEnsemble::OperatorSpace = mBodyOperatorSpace<Scalar>

OperatorSpace [5/6]

using ObservableEnsemble::OperatorSpace = mBodyOperatorSpace<Scalar>

OperatorSpace [6/6]

using ObservableEnsemble::OperatorSpace = mBodyOperatorSpace<Scalar>

StateSpace [1/6]

using ObservableEnsemble::StateSpace = ManyBodySpinSpace

StateSpace [2/6]

using ObservableEnsemble::StateSpace = ManyBodySpinSpace

StateSpace [3/6]

using ObservableEnsemble::StateSpace = ManyBodySpinSpace

StateSpace [4/6]

using ObservableEnsemble::StateSpace = ManyBodySpinSpace

StateSpace [5/6]

using ObservableEnsemble::StateSpace = ManyBodySpinSpace

StateSpace [6/6]

using ObservableEnsemble::StateSpace = ManyBodySpinSpace

Constructor & Destructor Documentation

ObservableEnsemble() [1/33]

ObservableEnsemble::ObservableEnsemble ( size_t  L, size_t  m  )

inline

                                               : m_L(L), m_M(m) {
            assert(m_L >= 1);
            class_M      = m_M;
            m_stateSpace = ManyBodySpinSpace(m_L, g_dimLoc);
            m_opSpace    = mBodyOperatorSpace<Scalar>(m_M, m_stateSpace);
            m_mbGRME     = GaussianRandomMatrixGenerator_GPU<OperatorSpace>(m_opSpace, m_M);
        }

ObservableEnsemble() [2/33]

ObservableEnsemble::ObservableEnsemble ( size_t  L, char **  CL_argv  )

inline

            : ObservableEnsemble(L, std::atoi(CL_argv[0])) {}

ObservableEnsemble() [3/33]

ObservableEnsemble::ObservableEnsemble ( )

default

ObservableEnsemble() [4/33]

ObservableEnsemble::ObservableEnsemble ( const ObservableEnsemble &  )

delete

ObservableEnsemble() [5/33]

ObservableEnsemble::ObservableEnsemble ( ObservableEnsemble &&  )

default

ObservableEnsemble() [6/33]

ObservableEnsemble::ObservableEnsemble ( size_t  L, size_t  m  )

inline

                                               : m_L(L), m_M(m) {
            assert(m_L >= 1);
            class_M      = m_M;
            m_stateSpace = ManyBodySpinSpace(m_L, g_dimLoc);
            m_opSpace    = mBodyOperatorSpace<Scalar>(m_M, m_stateSpace);
            m_mbGRME     = GaussianRandomMatrixGenerator< decltype(m_opSpace) >(m_opSpace, m_M);
        }

ObservableEnsemble() [7/33]

ObservableEnsemble::ObservableEnsemble ( size_t  L, char **  CL_argv  )

inline

            : ObservableEnsemble(L, std::atoi(CL_argv[0])) {}

ObservableEnsemble() [8/33]

ObservableEnsemble::ObservableEnsemble ( )

default

ObservableEnsemble() [9/33]

ObservableEnsemble::ObservableEnsemble ( const ObservableEnsemble &  )

delete

ObservableEnsemble() [10/33]

ObservableEnsemble::ObservableEnsemble ( ObservableEnsemble &&  )

default

ObservableEnsemble() [11/33]

ObservableEnsemble::ObservableEnsemble ( )

default

ObservableEnsemble() [12/33]

ObservableEnsemble::ObservableEnsemble ( const ObservableEnsemble &  )

delete

ObservableEnsemble() [13/33]

ObservableEnsemble::ObservableEnsemble ( ObservableEnsemble &&  )

default

ObservableEnsemble() [14/33]

ObservableEnsemble::ObservableEnsemble ( size_t  L, size_t  m  )

inline

                                               : m_L(L), m_M(m) {
            assert(m_L >= 1);
            class_M      = m_M;
            m_stateSpace = ManyBodySpinSpace(m_L, g_dimLoc);
            m_opSpace    = mBodyOperatorSpace<Scalar>(m_M, m_stateSpace);
            m_mbGRME     = GaussianRandomMatrixGenerator_GPU<OperatorSpace>(m_opSpace, m_M);
        }

ObservableEnsemble() [15/33]

ObservableEnsemble::ObservableEnsemble ( size_t  L, char **  CL_argv  )

inline

            : ObservableEnsemble(L, std::atoi(CL_argv[0])) {}

ObservableEnsemble() [16/33]

ObservableEnsemble::ObservableEnsemble ( )

default

ObservableEnsemble() [17/33]

ObservableEnsemble::ObservableEnsemble ( const ObservableEnsemble &  )

delete

ObservableEnsemble() [18/33]

ObservableEnsemble::ObservableEnsemble ( ObservableEnsemble &&  )

default

ObservableEnsemble() [19/33]

ObservableEnsemble::ObservableEnsemble ( size_t  L, size_t  m  )

inline

                                               : m_L(L), m_M(m) {
            assert(m_L >= 1);
            class_M      = m_M;
            m_stateSpace = ManyBodySpinSpace(m_L, g_dimLoc);
            m_opSpace    = mBodyOperatorSpace<Scalar>(m_M, m_stateSpace);
            m_mbGRME     = GaussianRandomMatrixGenerator< decltype(m_opSpace) >(m_opSpace, m_M);
        }

ObservableEnsemble() [20/33]

ObservableEnsemble::ObservableEnsemble ( size_t  L, char **  CL_argv  )

inline

            : ObservableEnsemble(L, std::atoi(CL_argv[0])) {}

ObservableEnsemble() [21/33]

ObservableEnsemble::ObservableEnsemble ( )

default

ObservableEnsemble() [22/33]

ObservableEnsemble::ObservableEnsemble ( const ObservableEnsemble &  )

delete

ObservableEnsemble() [23/33]

ObservableEnsemble::ObservableEnsemble ( ObservableEnsemble &&  )

default

ObservableEnsemble() [24/33]

ObservableEnsemble::ObservableEnsemble ( size_t  L, size_t  m  )

inline

                                               : m_L(L), m_M(m) {
            assert(m_L >= 1);
            class_M      = m_M;
            m_stateSpace = ManyBodySpinSpace(m_L, g_dimLoc);
            m_opSpace    = mBodyOperatorSpace<Scalar>(m_M, m_stateSpace);
            m_mbGRME     = GaussianRandomMatrixGenerator_GPU<OperatorSpace>(m_opSpace, m_M);
        }

ObservableEnsemble() [25/33]

ObservableEnsemble::ObservableEnsemble ( size_t  L, char **  CL_argv  )

inline

            : ObservableEnsemble(L, std::atoi(CL_argv[0])) {}

ObservableEnsemble() [26/33]

ObservableEnsemble::ObservableEnsemble ( )

default

ObservableEnsemble() [27/33]

ObservableEnsemble::ObservableEnsemble ( const ObservableEnsemble &  )

delete

ObservableEnsemble() [28/33]

ObservableEnsemble::ObservableEnsemble ( ObservableEnsemble &&  )

default

ObservableEnsemble() [29/33]

ObservableEnsemble::ObservableEnsemble ( size_t  L, size_t  m  )

inline

                                               : m_L(L), m_M(m) {
            assert(m_L >= 1);
            class_M      = m_M;
            m_stateSpace = ManyBodySpinSpace(m_L, g_dimLoc);
            m_opSpace    = mBodyOperatorSpace<Scalar>(m_M, m_stateSpace);
            m_mbGRME     = GaussianRandomMatrixGenerator< decltype(m_opSpace) >(m_opSpace, m_M);
        }

ObservableEnsemble() [30/33]

ObservableEnsemble::ObservableEnsemble ( size_t  L, char **  CL_argv  )

inline

            : ObservableEnsemble(L, std::atoi(CL_argv[0])) {}

ObservableEnsemble() [31/33]

ObservableEnsemble::ObservableEnsemble ( )

default

ObservableEnsemble() [32/33]

ObservableEnsemble::ObservableEnsemble ( const ObservableEnsemble &  )

delete

ObservableEnsemble() [33/33]

ObservableEnsemble::ObservableEnsemble ( ObservableEnsemble &&  )

default

Member Function Documentation

dataname() [1/7]

static std::string ObservableEnsemble::dataname ( size_t  rep )

inlinestatic

                                              {
            return "m" + std::to_string(class_M) + "_Observable" + "_No" + std::to_string(rep);
        };

dataname() [2/7]

static std::string ObservableEnsemble::dataname ( size_t  rep )

inlinestatic

                                              {
            return "m" + std::to_string(class_M) + "_Observable" + "_No" + std::to_string(rep);
        };

dataname() [3/7]

static std::string ObservableEnsemble::dataname ( size_t  rep )

inlinestatic

{ return "Observable_No" + std::to_string(rep); };

dataname() [4/7]

static std::string ObservableEnsemble::dataname ( size_t  rep )

inlinestatic

                                              {
            return "m" + std::to_string(class_M) + "_Observable" + "_No" + std::to_string(rep);
        };

dataname() [5/7]

static std::string ObservableEnsemble::dataname ( size_t  rep )

inlinestatic

                                              {
            return "m" + std::to_string(class_M) + "_Observable" + "_No" + std::to_string(rep);
        };

dataname() [6/7]

static std::string ObservableEnsemble::dataname ( size_t  rep )

inlinestatic

                                              {
            return "m" + std::to_string(class_M) + "_Observable" + "_No" + std::to_string(rep);
        };

dataname() [7/7]

static std::string ObservableEnsemble::dataname ( size_t  rep )

inlinestatic

                                              {
            return "m" + std::to_string(class_M) + "_Observable" + "_No" + std::to_string(rep);
        };

discard() [1/7]

void ObservableEnsemble::discard ( size_t  n )

inline

{ m_mbGRME.discard(n); }

discard() [2/7]

void ObservableEnsemble::discard ( size_t  n )

inline

{ m_mbGRME.discard(n); }

discard() [3/7]

void ObservableEnsemble::discard ( size_t  n )

inline

                               {
            // Write codes here
        }

discard() [4/7]

void ObservableEnsemble::discard ( size_t  n )

inline

{ m_mbGRME.discard(n); }

discard() [5/7]

void ObservableEnsemble::discard ( size_t  n )

inline

{ m_mbGRME.discard(n); }

discard() [6/7]

void ObservableEnsemble::discard ( size_t  n )

inline

{ m_mbGRME.discard(n); }

discard() [7/7]

void ObservableEnsemble::discard ( size_t  n )

inline

{ m_mbGRME.discard(n); }

nargs() [1/7]

static constexpr size_t ObservableEnsemble::nargs ( void  )

inlinestaticconstexpr

{ return 1; };

nargs() [2/7]

static constexpr size_t ObservableEnsemble::nargs ( void  )

inlinestaticconstexpr

{ return 1; };

nargs() [3/7]

static constexpr size_t ObservableEnsemble::nargs ( void  )

inlinestaticconstexpr

{ return 1; };

nargs() [4/7]

static constexpr size_t ObservableEnsemble::nargs ( void  )

inlinestaticconstexpr

{ return 1; };

nargs() [5/7]

static constexpr size_t ObservableEnsemble::nargs ( void  )

inlinestaticconstexpr

{ return 1; };

nargs() [6/7]

static constexpr size_t ObservableEnsemble::nargs ( void  )

inlinestaticconstexpr

{ return 1; };

nargs() [7/7]

static constexpr size_t ObservableEnsemble::nargs ( void  )

inlinestaticconstexpr

{ return 1; };

operator=() [1/14]

ObservableEnsemble & ObservableEnsemble::operator= ( const ObservableEnsemble &  )

delete

operator=() [2/14]

ObservableEnsemble & ObservableEnsemble::operator= ( const ObservableEnsemble &  )

delete

operator=() [3/14]

ObservableEnsemble & ObservableEnsemble::operator= ( const ObservableEnsemble &  )

delete

operator=() [4/14]

ObservableEnsemble & ObservableEnsemble::operator= ( const ObservableEnsemble &  )

delete

operator=() [5/14]

ObservableEnsemble & ObservableEnsemble::operator= ( const ObservableEnsemble &  )

delete

operator=() [6/14]

ObservableEnsemble & ObservableEnsemble::operator= ( const ObservableEnsemble &  )

delete

operator=() [7/14]

ObservableEnsemble & ObservableEnsemble::operator= ( const ObservableEnsemble &  )

delete

operator=() [8/14]

ObservableEnsemble & ObservableEnsemble::operator= ( ObservableEnsemble &&  )

default

operator=() [9/14]

ObservableEnsemble & ObservableEnsemble::operator= ( ObservableEnsemble &&  )

default

operator=() [10/14]

ObservableEnsemble & ObservableEnsemble::operator= ( ObservableEnsemble &&  )

default

operator=() [11/14]

ObservableEnsemble & ObservableEnsemble::operator= ( ObservableEnsemble &&  )

default

operator=() [12/14]

ObservableEnsemble & ObservableEnsemble::operator= ( ObservableEnsemble &&  )

default

operator=() [13/14]

ObservableEnsemble & ObservableEnsemble::operator= ( ObservableEnsemble &&  )

default

operator=() [14/14]

ObservableEnsemble & ObservableEnsemble::operator= ( ObservableEnsemble &&  )

default

operatorSpace() [1/7]

OperatorSpace & ObservableEnsemble::operatorSpace ( )

inline

{ return m_opSpace; }

operatorSpace() [2/7]

OperatorSpace & ObservableEnsemble::operatorSpace ( )

inline

{ return m_opSpace; }

operatorSpace() [3/7]

auto & ObservableEnsemble::operatorSpace ( )

inline

{ return; }

operatorSpace() [4/7]

OperatorSpace & ObservableEnsemble::operatorSpace ( )

inline

{ return m_opSpace; }

operatorSpace() [5/7]

OperatorSpace & ObservableEnsemble::operatorSpace ( )

inline

{ return m_opSpace; }

operatorSpace() [6/7]

OperatorSpace & ObservableEnsemble::operatorSpace ( )

inline

{ return m_opSpace; }

operatorSpace() [7/7]

OperatorSpace & ObservableEnsemble::operatorSpace ( )

inline

{ return m_opSpace; }

reset() [1/7]

void ObservableEnsemble::reset ( )

inline

{ m_mbGRME.reset(); }

reset() [2/7]

void ObservableEnsemble::reset ( )

inline

{ m_mbGRME.reset(); }

reset() [3/7]

void ObservableEnsemble::reset ( )

inline

                     {
            // Write codes here
        }

reset() [4/7]

void ObservableEnsemble::reset ( )

inline

{ m_mbGRME.reset(); }

reset() [5/7]

void ObservableEnsemble::reset ( )

inline

{ m_mbGRME.reset(); }

reset() [6/7]

void ObservableEnsemble::reset ( )

inline

{ m_mbGRME.reset(); }

reset() [7/7]

void ObservableEnsemble::reset ( )

inline

{ m_mbGRME.reset(); }

sample() [1/7]

MatrixType ObservableEnsemble::sample ( )

inline

                            {
            SubSpaceManager& subSpaces = SubSpaceManager::get_instance();
            auto&            subSpace  = subSpaces[m_L];
            int const        dim       = subSpace.dim();
            MatrixType       res(dim, dim);
            MatrixType       op = m_mbGRME();
            // std::cout << "observable sample:\n" << op << std::endl;
 
            constexpr void (*kernel)(
                Eigen::MatrixBase< std::remove_reference_t<decltype(*res.ptr())> >*,
                Eigen::MatrixBase< std::remove_reference_t<decltype(*op.ptr())> > const*,
                SubSpace<typename SubSpaceManager::TotalSpace, Scalar> const*)
                = &ObservableEnsemble_sample_kernel;
            struct cudaFuncAttributes attr;
            cuCHECK(cudaFuncGetAttributes(&attr, kernel));
            int const nThread = min(dim, (int)sqrt(attr.maxThreadsPerBlock));
            int const nBlock  = (int)(dim % nThread == 0 ? dim / nThread : dim / nThread + 1);
            debug_print(__func__ << ":\t subSpace.dim()=" << dim << ", nBlock=" << nBlock
                                 << ", nThread=" << nThread);
 
            auto dSubSpace = ObjectOnGPU<typename SubSpaceManager::SubHilbertSpace>(subSpaces[m_L]);
            kernel<<<dim3(nBlock, nBlock, 1), dim3(nThread, nThread, 1)>>>(res.ptr(), op.ptr(),
                                                                           dSubSpace.ptr());
            cuCHECK(cudaDeviceSynchronize());
 
            // std::cout << "observable res:\n" << res << std::endl;
            return res;
        }

sample() [2/7]

MatrixType ObservableEnsemble::sample ( )

inline

                            {
            SubSpaceManager& subSpaces    = SubSpaceManager::get_instance();
            auto&            subSpaceBase = subSpaces[m_L].basis();
            auto             op           = m_mbGRME();
            // std::cout << "observable sample:\n" << op << std::endl;
            MatrixType res = subSpaceBase.adjoint() * op * subSpaceBase;
 
            // std::cout << "observable res:\n" << res << std::endl;
            return res;
        }

sample() [3/7]

MatrixType ObservableEnsemble::sample ( )

inline

                            {
            MatrixType res;
            // Write codes to construct an observable here
            return res;
        }

sample() [4/7]

MatrixType ObservableEnsemble::sample ( )

inline

                            {
            SubSpaceManager& subSpaces = SubSpaceManager::get_instance();
            auto&            subSpace  = subSpaces[m_L];
            int const        dim       = subSpace.dim();
            MatrixType       res(dim, dim);
            MatrixType       op = m_mbGRME();
            // std::cout << "observable sample:\n" << op << std::endl;
 
            constexpr void (*kernel)(
                Eigen::MatrixBase< std::remove_reference_t<decltype(*res.ptr())> >*,
                Eigen::MatrixBase< std::remove_reference_t<decltype(*op.ptr())> > const*,
                SubSpace<typename SubSpaceManager::TotalSpace, Scalar> const*)
                = &ObservableEnsemble_sample_kernel;
            struct cudaFuncAttributes attr;
            cuCHECK(cudaFuncGetAttributes(&attr, kernel));
            int const nThread = min(dim, (int)sqrt(attr.maxThreadsPerBlock));
            int const nBlock  = (int)(dim % nThread == 0 ? dim / nThread : dim / nThread + 1);
            debug_print(__func__ << ":\t subSpace.dim()=" << dim << ", nBlock=" << nBlock
                                 << ", nThread=" << nThread);
 
            auto dSubSpace = ObjectOnGPU<typename SubSpaceManager::SubHilbertSpace>(subSpaces[m_L]);
            kernel<<<dim3(nBlock, nBlock, 1), dim3(nThread, nThread, 1)>>>(res.ptr(), op.ptr(),
                                                                           dSubSpace.ptr());
            cuCHECK(cudaDeviceSynchronize());
 
            // std::cout << "observable res:\n" << res << std::endl;
            return res;
        }

sample() [5/7]

MatrixType ObservableEnsemble::sample ( )

inline

                            {
            SubSpaceManager& subSpaces    = SubSpaceManager::get_instance();
            auto&            subSpaceBase = subSpaces[m_L].basis();
            auto             op           = m_mbGRME();
            // std::cout << "observable sample:\n" << op << std::endl;
            MatrixType res = subSpaceBase.adjoint() * op * subSpaceBase;
 
            // std::cout << "observable res:\n" << res << std::endl;
            return res;
        }

sample() [6/7]

MatrixType ObservableEnsemble::sample ( )

inline

                            {
            SubSpaceManager& subSpaces = SubSpaceManager::get_instance();
            auto&            subSpace  = subSpaces[m_L];
            int const        dim       = subSpace.dim();
            MatrixType       res(dim, dim);
            MatrixType       op = m_mbGRME();
            // std::cout << "observable sample:\n" << op << std::endl;
 
            constexpr void (*kernel)(
                Eigen::MatrixBase< std::remove_reference_t<decltype(*res.ptr())> >*,
                Eigen::MatrixBase< std::remove_reference_t<decltype(*op.ptr())> > const*,
                SubSpace<typename SubSpaceManager::TotalSpace, Scalar> const*)
                = &ObservableEnsemble_sample_kernel;
            struct cudaFuncAttributes attr;
            cuCHECK(cudaFuncGetAttributes(&attr, kernel));
            int const nThread = min(dim, (int)sqrt(attr.maxThreadsPerBlock));
            int const nBlock  = (int)(dim % nThread == 0 ? dim / nThread : dim / nThread + 1);
            debug_print(__func__ << ":\t subSpace.dim()=" << dim << ", nBlock=" << nBlock
                                 << ", nThread=" << nThread);
 
            auto dSubSpace = ObjectOnGPU<typename SubSpaceManager::SubHilbertSpace>(subSpaces[m_L]);
            kernel<<<dim3(nBlock, nBlock, 1), dim3(nThread, nThread, 1)>>>(res.ptr(), op.ptr(),
                                                                           dSubSpace.ptr());
            cuCHECK(cudaDeviceSynchronize());
 
            // std::cout << "observable res:\n" << res << std::endl;
            return res;
        }

sample() [7/7]

MatrixType ObservableEnsemble::sample ( )

inline

                            {
            SubSpaceManager& subSpaces    = SubSpaceManager::get_instance();
            auto&            subSpaceBase = subSpaces[m_L].basis();
            auto             op           = m_mbGRME();
            // std::cout << "observable sample:\n" << op << std::endl;
            MatrixType res = subSpaceBase.adjoint() * op * subSpaceBase;
 
            // std::cout << "observable res:\n" << res << std::endl;
            return res;
        }

stateSpace() [1/7]

StateSpace & ObservableEnsemble::stateSpace ( )

inline

{ return m_stateSpace; }

stateSpace() [2/7]

StateSpace & ObservableEnsemble::stateSpace ( )

inline

{ return m_stateSpace; }

stateSpace() [3/7]

auto & ObservableEnsemble::stateSpace ( )

inline

{ return; }

stateSpace() [4/7]

StateSpace & ObservableEnsemble::stateSpace ( )

inline

{ return m_stateSpace; }

stateSpace() [5/7]

StateSpace & ObservableEnsemble::stateSpace ( )

inline

{ return m_stateSpace; }

stateSpace() [6/7]

StateSpace & ObservableEnsemble::stateSpace ( )

inline

{ return m_stateSpace; }

stateSpace() [7/7]

StateSpace & ObservableEnsemble::stateSpace ( )

inline

{ return m_stateSpace; }

usage() [1/7]

static std::string ObservableEnsemble::usage ( size_t  nargs )

inlinestatic

                                             {
            std::string usage = "ObservableEnsemble Usage: ";
            usage += std::to_string(nargs + 0) + ".(m) ";
            return usage;
        };

usage() [2/7]

static std::string ObservableEnsemble::usage ( size_t  nargs )

inlinestatic

                                             {
            std::string usage = "ObservableEnsemble Usage: ";
            usage += std::to_string(nargs + 0) + ".(m) ";
            return usage;
        };

usage() [3/7]

static std::string ObservableEnsemble::usage ( size_t  nargs )

inlinestatic

                                             {
            std::string usage = "ObservableEnsemble Usage: ";
            usage += std::to_string(nargs + 0) + ".(...) ";
            // Add codes here
            return usage;
        };

usage() [4/7]

static std::string ObservableEnsemble::usage ( size_t  nargs )

inlinestatic

                                             {
            std::string usage = "ObservableEnsemble Usage: ";
            usage += std::to_string(nargs + 0) + ".(m) ";
            return usage;
        };

usage() [5/7]

static std::string ObservableEnsemble::usage ( size_t  nargs )

inlinestatic

                                             {
            std::string usage = "ObservableEnsemble Usage: ";
            usage += std::to_string(nargs + 0) + ".(m) ";
            return usage;
        };

usage() [6/7]

static std::string ObservableEnsemble::usage ( size_t  nargs )

inlinestatic

                                             {
            std::string usage = "ObservableEnsemble Usage: ";
            usage += std::to_string(nargs + 0) + ".(m) ";
            return usage;
        };

usage() [7/7]

static std::string ObservableEnsemble::usage ( size_t  nargs )

inlinestatic

                                             {
            std::string usage = "ObservableEnsemble Usage: ";
            usage += std::to_string(nargs + 0) + ".(m) ";
            return usage;
        };

Member Data Documentation

class_M

size_t ObservableEnsemble::class_M = 0

static

m_L

size_t ObservableEnsemble::m_L

private

m_M

size_t ObservableEnsemble::m_M

private

m_mbGRME [1/2]

GaussianRandomMatrixGenerator< OperatorSpace > ObservableEnsemble::m_mbGRME

private

m_mbGRME [2/2]

GaussianRandomMatrixGenerator<OperatorSpace> ObservableEnsemble::m_mbGRME

private

m_opSpace

OperatorSpace ObservableEnsemble::m_opSpace

private

m_stateSpace

StateSpace ObservableEnsemble::m_stateSpace

private

---

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

• /Users/shoki/GitHub/Locality/PBC_TI/CRMT/Ensemble.cuh
• /Users/shoki/GitHub/Locality/PBC_TI/CRMT/Ensemble.hpp
• /Users/shoki/GitHub/Locality/PBC_TI/Framework/Ensemble.template.hpp
• /Users/shoki/GitHub/Locality/PBC_TI/ShortRange_mBody/Ensemble.cuh
• /Users/shoki/GitHub/Locality/PBC_TI/ShortRange_mBody/Ensemble.hpp
• /Users/shoki/GitHub/Locality/PBC_TI/transverseIsing/Ensemble.cuh
• /Users/shoki/GitHub/Locality/PBC_TI/transverseIsing/Ensemble.hpp