Ensemble.cuh File Reference

Classes
class	SubSpaceManager
class	Ensemble::internal::ScalableOperator< MatrixType >
class	HamiltonianEnsemble
class	ObservableEnsemble

Namespaces
namespace	Ensemble
namespace	Ensemble::internal

Typedefs
using	RealScalar = Real_t
using	Scalar = Complex_t< RealScalar >

Functions
template<class Derived1 , class Derived2 , class TotalSpace >
__global__ void	ObservableEnsemble_sample_kernel (Eigen::MatrixBase< Derived1 > *resPtr, Eigen::MatrixBase< Derived2 > const *opTotPtr, SubSpace< TotalSpace, typename Derived1::Scalar > const *subSpacePtr)

Variables
int	g_dimLoc = 2
int	g_momentum = 0

Typedef Documentation

RealScalar

using RealScalar = Real_t

Scalar

using Scalar = Complex_t<RealScalar>

Function Documentation

ObservableEnsemble_sample_kernel()

template<class Derived1 , class Derived2 , class TotalSpace >

__global__ void ObservableEnsemble_sample_kernel	(	Eigen::MatrixBase< Derived1 > *	resPtr,
		Eigen::MatrixBase< Derived2 > const *	opTotPtr,
		SubSpace< TotalSpace, typename Derived1::Scalar > const *	subSpacePtr
	)

                                                                        {
    int const idx   = blockIdx.x * blockDim.x + threadIdx.x;
    int const idy   = blockIdx.y * blockDim.y + threadIdx.y;
    auto&     basis = subSpacePtr->basis();
    if(idx >= resPtr->rows() || idy >= resPtr->cols()) return;
    if(idx > idy) return;
 
    (*resPtr)(idx, idy) = {0, 0};
    int innerX, innerY;
    for(auto valuePosX = basis.outerIndexPtr()[idx]; valuePosX != basis.outerIndexPtr()[idx + 1];
        ++valuePosX) {
        innerX          = basis.innerIndexPtr()[valuePosX];
        auto conjCoeffX = conj(basis.valuePtr()[valuePosX]);
 
        for(auto valuePosY = basis.outerIndexPtr()[idy];
            valuePosY != basis.outerIndexPtr()[idy + 1]; ++valuePosY) {
            innerY      = basis.innerIndexPtr()[valuePosY];
            auto CoeffY = basis.valuePtr()[valuePosY];
 
            (*resPtr)(idx, idy) += conjCoeffX * (*opTotPtr)(innerX, innerY) * CoeffY;
        }
    }
 
    if(idx == idy) (*resPtr)(idy, idx) = real((*resPtr)(idy, idx));
 
    (*resPtr)(idy, idx) = conj((*resPtr)(idx, idy));
}

Variable Documentation

g_dimLoc

int g_dimLoc = 2

g_momentum

int g_momentum = 0