 #include <algorithm>
 #include <fstream>
+#include <cmath>
 
 #include <cusp/io/matrix_market.h>
 #include <cusp/csr_matrix.h>

	if (fileSol.length() > 0)
		cusp::io::write_matrix_market_file(x, fileSol);

	PrintStats(success, mySolver, mySpmv);

	// Calculate the actual residual and its norm.
	if (verbose) {
		Vector r(A.num_rows);
		VectorView r_view(r);
		mySpmv(x_view, r_view);


	cusp::io::read_matrix_market_file(A, fileMat);


	if (fileRhs.length() > 0)
		cusp::io::read_matrix_market_file(b, fileRhs);
	else
