 // where to perform the computation
 typedef cusp::device_memory MemorySpace;
 
-// which floating point type to use
+// which doubleing point type to use
 typedef float ValueType;
 
 // Reads an array from a Matrix Market sparse matrix file

  printf("Compute time %f ms\n", elapsed/1000.0);
  printf("Total time %f s\n", elapsed_total/(1000.0 * 1000.0));

  return 0;
}

  cusp::array1d<ValueType, MemorySpace> x(A.num_rows, 0);

  // Set stopping criteria (iteration count, relative error)
  cusp::verbose_monitor<ValueType> monitor(bb, 100 * A.num_rows, 1e-20);

  // set preconditioner (identity)
  cusp::identity_operator<ValueType, MemorySpace> M(A.num_rows, A.num_rows);


  printf("Starting cg...\n");
  gettimeofday(&ts_compute, NULL);
