static std::ptrdiff_t num_nonzeros (matrix_type& cm) { return cm.nnz(); // Joerg, this isn't very intuitive :o( // return cm.non_zeros(); }