CUSP库使用实践

CUSP最新版本0.6,项目已经停止开发多年，比较老旧，在新版本CUDA下编译使用CUSP会有一些小错误，主要是头文件缺失或是缺少引用错误。

#include 
#include 
#include 

//compile: nvcc -O3 -I ....github.dircusplibrary cuda_ex4.cu 
//nvprof .a.exe

void print_version()
{
    int thrust_major=THRUST_MAJOR_VERSION;
    int thrust_minor=THRUST_MINOR_VERSION;

    int cusp_major=CUSP_MAJOR_VERSION;
    int cusp_minor=CUSP_MINOR_VERSION;
    std::cout<<"Thrust v"<
#include 

void test_sparse()
{
    cusp::csr_matrix csr_host(5,8,12);
    cusp::csr_matrix csr_device(csr_host);
    cusp::hyb_matrix hyb_device(csr_device);
}

#include
#include
#include
#include 

void test_cg()
{
    typedef int IndexType;
    typedef cusp::device_memory MemorySpace;
    typedef float ValueType;
    static const int ROWS=1000;
#if 0    
    cusp::hyb_matrix A;
#else 
    cusp::coo_matrix A;
#endif
    cusp::gallery::poisson5pt(A,ROWS,ROWS);
    cusp::array1d x(A.num_rows,0);
    cusp::array1d b(A.num_rows,1);
    cusp::monitor monitor(b,1000,1e-8,0,true);
#if 0
    cusp::identity_operatorM(A.num_rows,A.num_rows);
    cusp::krylov::cg(A,x,b,monitor,M);
    monitor.print();
#else 
    // setup preconditioner
    cusp::precond::aggregation::smoothed_aggregation M(A);
    // solve
    cusp::krylov::cg(A, x, b, monitor, M);
    // report status
    monitor.print();
    // print hierarchy information
    std::cout << "nPreconditioner statistics" << std::endl;
    M.print();
#endif

}

int main(int argc, char **argv)
{
    print_version();
    test_sparse();
    //test_cg();
    return(0);
}

• 快速入手