EigenExpValues_test.cu File Reference

Classes
class	MagmaController

Functions
int	main (int argc, char **argv)

Function Documentation

main()

int main	(	int	argc,
		char **	argv
	)

                                {
    if(argc != 2) {
        std::cerr << "Usage: 0.(This) 1.(dim)\n"
                  << "       argc=" << argc << std::endl;
        std::exit(EXIT_FAILURE);
    }
    int const dim = std::atoi(argv[1]);
    MagmaController::getInstance();
 
    std::mt19937                     mt;
    std::normal_distribution<double> Gaussian;
    using Scalar       = Complex_t<double>;
    auto randomComplex = [&]() { return Complex_t<double>(Gaussian(mt), Gaussian(mt)); };
 
    Eigen::MatrixX<Scalar> hMat1 = Eigen::MatrixX<Scalar>::NullaryExpr(dim, dim, randomComplex);
    hMat1                        = (hMat1 + hMat1.adjoint().eval()) / 2.0;
    Eigen::MatrixX<Scalar> hMat2 = Eigen::MatrixX<Scalar>::NullaryExpr(dim, dim, randomComplex);
    hMat2                        = (hMat2 + hMat2.adjoint().eval()) / 2.0;
    ObjectOnGPU< decltype(hMat1) > dMat1(hMat1);
    ObjectOnGPU< decltype(hMat2) > dMat2(hMat2);
    // std::cout << "hMat1:\n" << hMat1 << std::endl;
    // std::cout << "dMat1:\n" << dMat1 << std::endl;
    // std::cout << "hMat2:\n" << hMat2 << std::endl;
    // std::cout << "dMat2:\n" << dMat2 << std::endl;
 
    MatrixUtils::EigenSolver< decltype(dMat1) > dEigSolver(std::move(dMat1));
    auto expvals_d = EigenExpValues(dEigSolver.eigenvectors(), dMat2);
    std::cout << expvals_d << std::endl;
 
    // MatrixUtils::EigenSolver< decltype(hMat1) > hEigSolver(std::move(hMat1));
    // auto expvals_h = EigenExpValues(hEigSolver.eigenvectors(), hMat2);
    // std::cout << expvals_h << std::endl;
 
    // std::cout << "Max difference = " << (expvals_d - expvals_h).cwiseAbs().maxCoeff() << std::endl;
 
    return EXIT_SUCCESS;
}