1 #ifndef VIENNACL_LINALG_DETAIL_AMG_AMG_DEBUG_HPP
2 #define VIENNACL_LINALG_DETAIL_AMG_AMG_DEBUG_HPP
31 #define VIENNACL_AMG_MATRIXTYPE boost::numeric::ublas::matrix<ScalarType>
33 #define VIENNACL_AMG_MATRIXTYPE MatrixType
45 #ifdef VIENNACL_AMG_DEBUG
46 template<
typename MatrixT>
49 typedef typename MatrixT::value_type
ScalarType;
50 typedef typename VIENNACL_AMG_MATRIXTYPE::iterator1 InternalRowIterator;
51 typedef typename VIENNACL_AMG_MATRIXTYPE::iterator2 InternalColIterator;
55 for (InternalRowIterator row_iter = mat2.begin1(); row_iter != mat2.end1(); ++row_iter)
57 for (InternalColIterator col_iter = row_iter.begin(); col_iter != row_iter.end(); ++col_iter)
59 std::cout << *col_iter <<
" ";
61 std::cout << std::endl;
63 std::cout << std::endl;
66 template<
typename VectorT>
69 for (
typename VectorT::const_iterator iter = vec.begin(); iter != vec.end(); ++iter)
71 std::cout << *iter <<
" ";
73 std::cout << std::endl;
76 template<
typename MatrixT>
79 template<
typename VectorT>
A reader and writer for the matrix market format is implemented here.
#define VIENNACL_AMG_MATRIXTYPE
void printvector(VectorT const &)
void printmatrix(MatrixT &, int)