 			}
 		}
 	}
-	std::cout << "\nQT stuff:\n";
+	/*std::cout << "\nQT stuff:\n";
 	std::cout << "Original size: " << QTHost.values.size() << std::endl;
 	std::cout << "Actual size  : " << idx << std::endl;
 	NavierStokesSolver<memoryType>::QT = QTHost;

	cusp::io::write_matrix_market_file(NavierStokesSolver<memoryType>::Q, "Q.mtx");
	std::cout << "Wrote Q to file." << std::endl;
	cusp::io::write_matrix_market_file(NavierStokesSolver<memoryType>::QT, "QT.mtx");
	std::cout << "Wrote QT to file." << std::endl;
}

template <typename memoryType>

	CHost.sort_by_row_and_column();
	CHost.values[0] += CHost.values[0];
	NavierStokesSolver<memoryType>::C = CHost;
	//cusp::io::write_matrix_market_file(NavierStokesSolver<memoryType>::C, "C-generateC.mtx");
	cusp::blas::scal(NavierStokesSolver<memoryType>::C.values, dt);
}
*/

template class DFImprovedSolver<host_memory>;
template class DFImprovedSolver<device_memory>;