diffusion.cu

#include <cusp/gallery/diffusion.h>
#include <cusp/coo_matrix.h>
#include <cusp/print.h>
 
int main(void)
{
    cusp::coo_matrix<int, float, cusp::device_memory> A;
 
    // create a matrix for a diffusion problem on a 4x4 grid
    cusp::gallery::diffusion<cusp::gallery::FD>(A, 4, 4);
 
    // print matrix
    cusp::print(A);
 
    return 0;
}