EigenExpValues.cuh File Reference

Functions
template<class Matrix1 , class Matrix2 >
auto	EigenExpValues (ObjectOnGPU< Matrix1 > const &eigenvectors, ObjectOnGPU< Matrix2 > const &observable) -> Eigen::VectorX< typename Matrix1::RealScalar >

Function Documentation

EigenExpValues()

template<class Matrix1 , class Matrix2 >

auto EigenExpValues ( ObjectOnGPU< Matrix1 > const &  eigenvectors, ObjectOnGPU< Matrix2 > const &  observable  ) -> Eigen::VectorX<typename Matrix1::RealScalar>

inline

                                                {
    auto const                                   dim  = eigenvectors.cols();
    auto                                         tmp1 = observable * eigenvectors;
    Eigen::VectorX<typename Matrix1::RealScalar> expvals(dim);
 
    magma_int_t dev = 0;
    magma_getdevice(&dev);
    std::vector<magma_queue_t> queue(omp_get_max_threads(), NULL);
    for(int j = 0; j < queue.size(); ++j) magma_queue_create(dev, &queue[j]);
#pragma omp parallel for
    for(int j = 0; j < dim; ++j) {
        expvals(j)
            = real(magma_dotc_wrapper(dim, eigenvectors.data() + j * eigenvectors.LD(), 1,
                                      tmp1.data() + j * tmp1.LD(), 1, queue[omp_get_thread_num()]));
    }
    for(int j = 0; j < queue.size(); ++j) magma_queue_destroy(queue[j]);
 
    return expvals;
}