     cusp::array1d<ValueType,cusp::device_memory> y(A.num_rows, 0);
     
     // warm up
-    A.multiply(x, y);
+    cusp::multiply(A, x, y);
 
     // benchmark SpMV
     timer t;

        std::cout << "\n----------- benchmarking CG -----------\n";
        // compute residual
        cusp::array1d<ValueType, cusp::host_memory> r(A.num_rows,0);
        A.multiply(x,r);
        cusp::blas::axpy(b, r, ValueType(-1.0));

        ValueType residual_norm = cusp::blas::nrm2(r);

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

#include <iostream>
#include <string>
