statmech_gpu.cuh File Reference

Functions
template<typename TypeComplex >
magma_int_t	EigenValues (std::vector< double > &eigenvalue, magma_int_t const dim, matrix_gpu< TypeComplex > &dHamiltonian)
template<typename TypeComplex >
double	SpectralRange (double &OpMin, magma_int_t const dim, matrix_gpu< TypeComplex > &dHamiltonian)
template<typename TypeComplex >
magma_int_t	Diagonalize (magma_int_t const dim, matrix_gpu< TypeComplex > &dHamiltonian, std::vector< double > &eigenvalue_h)
template<typename TypeComplex >
magma_int_t	EigenMatrixElements (std::vector< double > &eigenvalue_h, magma_int_t const dim, matrix_gpu< TypeComplex > &dHamiltonian, matrix_gpu< TypeComplex > &dOperator, GPUconfig const &GPUconf)
template<typename TypeComplex >
magma_int_t	EigenExpValue (std::vector< double > &EXPvalue_h, std::vector< double > &eigenvalue_h, magma_int_t const dim, matrix_gpu< TypeComplex > &dHamiltonian, matrix_gpu< TypeComplex > &dOperator, GPUconfig const &GPUconf)

Function Documentation

Diagonalize()

template<typename TypeComplex >

magma_int_t Diagonalize	(	magma_int_t const	dim,
		matrix_gpu< TypeComplex > &	dHamiltonian,
		std::vector< double > &	eigenvalue_h
	)

EigenExpValue()

template<typename TypeComplex >

magma_int_t EigenExpValue	(	std::vector< double > &	EXPvalue_h,
		std::vector< double > &	eigenvalue_h,
		magma_int_t const	dim,
		matrix_gpu< TypeComplex > &	dHamiltonian,
		matrix_gpu< TypeComplex > &	dOperator,
		GPUconfig const &	GPUconf
	)

                                                                                        {
    debug_print("# " << __func__ << "_gpu: dim=" << dim);
    magma_int_t info = EigenMatrixElements(eigenvalue_h, dim, dHamiltonian, dOperator, GPUconf);
 
    std::vector<TypeComplex> vector(dim);
    magma_getvector(dim, sizeof(TypeComplex), dOperator.ptr(), dOperator.LD() + 1, &*vector.begin(),
                    1, GPUconf.queue());
 
#pragma omp parallel for
    for(size_t i = 0; i < dim; ++i) {
        EXPvalue_h[i] = (double)real(vector[i]);
        if(isnan(EXPvalue_h[i])) {
#pragma omp critical
            std::cerr << "Error(" << __func__ << "): EXPvalue_h[" << i << "] = NaN." << std::endl;
        }
    }
    return info;
}

EigenMatrixElements()

template<typename TypeComplex >

magma_int_t EigenMatrixElements	(	std::vector< double > &	eigenvalue_h,
		magma_int_t const	dim,
		matrix_gpu< TypeComplex > &	dHamiltonian,
		matrix_gpu< TypeComplex > &	dOperator,
		GPUconfig const &	GPUconf
	)

                                                                                              {
    debug_print("# " << __func__ << "_gpu: dim=" << dim);
    magma_int_t const LD   = dHamiltonian.LD();
    magma_int_t const LDT  = magma_roundup(dim, 32);
    magma_int_t       info = Diagonalize(dim, dHamiltonian, eigenvalue_h);
    if(info != 0) return info;
 
    matrix_gpu<TypeComplex> ztemp(LDT, dim);
    matrixProduct_hemm(MagmaLeft, MagmaLower, dim, dim, dOperator, dHamiltonian, ztemp, GPUconf);
    matrixProduct_gemm(MagmaConjTrans, MagmaNoTrans, dim, dim, dim, dHamiltonian, ztemp, dOperator,
                       GPUconf);
    return info;
}

EigenValues()

template<typename TypeComplex >

magma_int_t EigenValues	(	std::vector< double > &	eigenvalue,
		magma_int_t const	dim,
		matrix_gpu< TypeComplex > &	dHamiltonian
	)

SpectralRange()

template<typename TypeComplex >

double SpectralRange	(	double &	OpMin,
		magma_int_t const	dim,
		matrix_gpu< TypeComplex > &	dHamiltonian
	)