     const size_t num_iterations = 500;
 
     for(size_t i = 0; i < num_iterations; i++)
-        A.multiply(x, y);
+        cusp::multiply(A, x, y);
 
     float time = t.seconds_elapsed() / num_iterations;
     float GFLOPs = (time == 0) ? 0 : (2 * A.num_entries / time) / 1e9;

#include <cusp/csr_matrix.h>
#include <cusp/io/matrix_market.h>
#include <cusp/blas.h>

#include <iostream>
#include <string>

    cusp::array1d<ValueType,cusp::device_memory> y(A.num_rows, 0);
    
    // warm up
    cusp::multiply(A, x, y);

    // benchmark SpMV
    timer t;
