SparseCompressed_test.cpp File Reference

Functions
int	main ()

Function Documentation

main()

int main

(

void

)

           {
    Eigen::MatrixXcd dense = Eigen::MatrixXcd::Zero(4, 4);
    std::cout << dense << std::endl;
 
    SparseCompressed mat2(4, 4);
    std::cout << mat2 << std::endl;
 
    SparseCompressed mat(4, 4);
    std::cout << mat << std::endl;
 
    int                                    dim         = 6;
    int                                    maxNonZeros = std::rand() % (dim * dim);
    Eigen::SparseMatrix<Complex_t<double>> sparse(dim, dim);
    for(int j = 0; j != maxNonZeros; ++j) {
        sparse.coeffRef(std::rand() % dim, std::rand() % dim)
            = {(std::rand() + 1) / double(RAND_MAX + 1), (std::rand() + 1) / double(RAND_MAX + 1)};
    }
    sparse.makeCompressed();
    std::cout << sparse << std::endl;
    std::cout << "sparse.nonZeros() = " << sparse.nonZeros() << std::endl;
 
    mat.resize(dim, dim, sparse.nonZeros());
    mat.nonZeros() = sparse.nonZeros();
    for(int j = 0; j != dim + 1; ++j) mat.outerIndexPtr()[j] = sparse.outerIndexPtr()[j];
    for(int j = 0; j != mat.nonZeros(); ++j) {
        mat.innerIndexPtr()[j] = sparse.innerIndexPtr()[j];
        mat.valuePtr()[j]      = sparse.valuePtr()[j];
    }
    std::cout << mat << std::endl;
    std::cout << "mat.nonZeros() = " << sparse.nonZeros() << std::endl;
 
    mat2 = mat.adjoint();
    std::cout << mat << std::endl;
    std::cout << mat2 << std::endl;
 
    return EXIT_SUCCESS;
}