11 #ifndef EIGEN_SPARSEVIEW_H
12 #define EIGEN_SPARSEVIEW_H
18 template<
typename MatrixType>
19 struct traits<SparseView<MatrixType> > : traits<MatrixType>
22 typedef Sparse StorageKind;
24 Flags = int(traits<MatrixType>::Flags) & (
RowMajorBit)
30 template<
typename MatrixType>
31 class SparseView :
public SparseMatrixBase<SparseView<MatrixType> >
33 typedef typename MatrixType::Nested MatrixTypeNested;
34 typedef typename internal::remove_all<MatrixTypeNested>::type _MatrixTypeNested;
36 EIGEN_SPARSE_PUBLIC_INTERFACE(SparseView)
38 SparseView(const MatrixType& mat, const Scalar& m_reference = Scalar(0),
39 typename NumTraits<Scalar>::Real m_epsilon = NumTraits<Scalar>::dummy_precision()) :
40 m_matrix(mat), m_reference(m_reference), m_epsilon(m_epsilon) {}
44 inline Index rows()
const {
return m_matrix.rows(); }
45 inline Index cols()
const {
return m_matrix.cols(); }
47 inline Index innerSize()
const {
return m_matrix.innerSize(); }
48 inline Index outerSize()
const {
return m_matrix.outerSize(); }
51 MatrixTypeNested m_matrix;
53 typename NumTraits<Scalar>::Real m_epsilon;
56 template<
typename MatrixType>
57 class SparseView<MatrixType>::InnerIterator :
public _MatrixTypeNested::InnerIterator
60 typedef typename _MatrixTypeNested::InnerIterator IterBase;
61 InnerIterator(
const SparseView& view, Index outer) :
62 IterBase(view.m_matrix, outer), m_view(view)
67 EIGEN_STRONG_INLINE InnerIterator& operator++()
69 IterBase::operator++();
74 using IterBase::value;
77 const SparseView& m_view;
80 void incrementToNonZero()
82 while((
bool(*
this)) && internal::isMuchSmallerThan(value(), m_view.m_reference, m_view.m_epsilon))
84 IterBase::operator++();
89 template<
typename Derived>
90 const SparseView<Derived> MatrixBase<Derived>::sparseView(
const Scalar& m_reference,
91 typename NumTraits<Scalar>::Real m_epsilon)
const
93 return SparseView<Derived>(derived(), m_reference, m_epsilon);