Ensemble copy.hpp File Reference

Go to the source code of this file.

Classes
class	RandomMatrixEnsemble< Scalar_t >

Functions
template<class Derived1 , class Derived2 , class TotalSpace , typename Scalar >
__global__ void	constructHamiltonian_kernel (Eigen::MatrixBase< Derived1 > *dTotMat, Eigen::MatrixBase< Derived2 > const *dLocMat, SubSpace< TotalSpace, Scalar > const *subSpace)

Function Documentation

constructHamiltonian_kernel()

template<class Derived1 , class Derived2 , class TotalSpace , typename Scalar >

__global__ void constructHamiltonian_kernel	(	Eigen::MatrixBase< Derived1 > *	dTotMat,
		Eigen::MatrixBase< Derived2 > const *	dLocMat,
		SubSpace< TotalSpace, Scalar > const *	subSpace
	)

                                                                                          {
    int const idx = blockIdx.x * blockDim.x + threadIdx.x;
    int const idy = blockIdx.y * blockDim.y + threadIdx.y;
    if(idx >= (*dTotMat).rows() || idy >= (*dTotMat).cols()) return;
    if(idx < idy) return;
 
    int const dimLoc = (*dLocMat).cols();
 
    (*dTotMat)(idx, idy) = {0, 0};
    magma_index_t innerX, innerY, locIdX, locIdY;
    for(auto valuePosX = (*subSpace).basis().outerIndexPtr()[idx];
        valuePosX != (*subSpace).basis().outerIndexPtr()[idx + 1]; ++valuePosX) {
        innerX          = (*subSpace).basis().innerIndexPtr()[valuePosX];
        auto conjCoeffX = conj((*subSpace).basis().valuePtr()[valuePosX]);
 
        for(auto valuePosY = (*subSpace).basis().outerIndexPtr()[idy];
            valuePosY != (*subSpace).basis().outerIndexPtr()[idy + 1]; ++valuePosY) {
            innerY      = (*subSpace).basis().innerIndexPtr()[valuePosY];
            auto CoeffY = (*subSpace).basis().valuePtr()[valuePosY];
 
            // if((innerX - innerY) / dimLoc != 0) continue; // WRONG FORMULA used until 2020/10/07
            if(innerX / dimLoc != innerY / dimLoc) continue;
            locIdX = innerX % dimLoc;
            locIdY = innerY % dimLoc;
            (*dTotMat)(idx, idy) += conjCoeffX * (*dLocMat)(locIdX, locIdY) * CoeffY;
        }
    }
 
    (*dTotMat)(idy, idx) = conj((*dTotMat)(idx, idy));
    // printf("LDT=%d, id=(%d,%d), val=%+f%+f*i\t %s is yet to be tested.\n", LDT, idx, idy,
    //        real(dTotMat[idx + LDT * idy]), imag(dTotMat[idx + LDT * idy]), __func__);
    // printf("\t id=(%d,%d)\n", idx, idy);
    return;
}