10 #ifndef EIGEN_SPARSEUTIL_H
11 #define EIGEN_SPARSEUTIL_H
16 #define EIGEN_DBG_SPARSE(X)
18 #define EIGEN_DBG_SPARSE(X) X
21 #define EIGEN_SPARSE_INHERIT_ASSIGNMENT_OPERATOR(Derived, Op) \
22 template<typename OtherDerived> \
23 EIGEN_STRONG_INLINE Derived& operator Op(const Eigen::SparseMatrixBase<OtherDerived>& other) \
25 return Base::operator Op(other.derived()); \
27 EIGEN_STRONG_INLINE Derived& operator Op(const Derived& other) \
29 return Base::operator Op(other); \
32 #define EIGEN_SPARSE_INHERIT_SCALAR_ASSIGNMENT_OPERATOR(Derived, Op) \
33 template<typename Other> \
34 EIGEN_STRONG_INLINE Derived& operator Op(const Other& scalar) \
36 return Base::operator Op(scalar); \
39 #define EIGEN_SPARSE_INHERIT_ASSIGNMENT_OPERATORS(Derived) \
40 EIGEN_SPARSE_INHERIT_ASSIGNMENT_OPERATOR(Derived, =) \
41 EIGEN_SPARSE_INHERIT_ASSIGNMENT_OPERATOR(Derived, +=) \
42 EIGEN_SPARSE_INHERIT_ASSIGNMENT_OPERATOR(Derived, -=) \
43 EIGEN_SPARSE_INHERIT_SCALAR_ASSIGNMENT_OPERATOR(Derived, *=) \
44 EIGEN_SPARSE_INHERIT_SCALAR_ASSIGNMENT_OPERATOR(Derived, /=)
46 #define _EIGEN_SPARSE_PUBLIC_INTERFACE(Derived, BaseClass) \
47 typedef BaseClass Base; \
48 typedef typename Eigen::internal::traits<Derived >::Scalar Scalar; \
49 typedef typename Eigen::NumTraits<Scalar>::Real RealScalar; \
50 typedef typename Eigen::internal::nested<Derived >::type Nested; \
51 typedef typename Eigen::internal::traits<Derived >::StorageKind StorageKind; \
52 typedef typename Eigen::internal::traits<Derived >::Index Index; \
53 enum { RowsAtCompileTime = Eigen::internal::traits<Derived >::RowsAtCompileTime, \
54 ColsAtCompileTime = Eigen::internal::traits<Derived >::ColsAtCompileTime, \
55 Flags = Eigen::internal::traits<Derived >::Flags, \
56 CoeffReadCost = Eigen::internal::traits<Derived >::CoeffReadCost, \
57 SizeAtCompileTime = Base::SizeAtCompileTime, \
58 IsVectorAtCompileTime = Base::IsVectorAtCompileTime }; \
59 using Base::derived; \
60 using Base::const_cast_derived;
62 #define EIGEN_SPARSE_PUBLIC_INTERFACE(Derived) \
63 _EIGEN_SPARSE_PUBLIC_INTERFACE(Derived, Eigen::SparseMatrixBase<Derived >)
65 const int CoherentAccessPattern = 0x1;
66 const int InnerRandomAccessPattern = 0x2 | CoherentAccessPattern;
67 const int OuterRandomAccessPattern = 0x4 | CoherentAccessPattern;
68 const int RandomAccessPattern = 0x8 | OuterRandomAccessPattern | InnerRandomAccessPattern;
70 template<
typename Derived>
class SparseMatrixBase;
71 template<
typename _Scalar,
int _Flags = 0,
typename _Index =
int>
class SparseMatrix;
72 template<
typename _Scalar,
int _Flags = 0,
typename _Index =
int>
class DynamicSparseMatrix;
73 template<
typename _Scalar,
int _Flags = 0,
typename _Index =
int>
class SparseVector;
74 template<
typename _Scalar,
int _Flags = 0,
typename _Index =
int>
class MappedSparseMatrix;
76 template<
typename MatrixType,
int Mode>
class SparseTriangularView;
78 template<
typename Lhs,
typename Rhs>
class SparseDiagonalProduct;
79 template<
typename MatrixType>
class SparseView;
81 template<
typename Lhs,
typename Rhs>
class SparseSparseProduct;
82 template<
typename Lhs,
typename Rhs>
class SparseTimeDenseProduct;
83 template<
typename Lhs,
typename Rhs>
class DenseTimeSparseProduct;
84 template<
typename Lhs,
typename Rhs,
bool Transpose>
class SparseDenseOuterProduct;
86 template<
typename Lhs,
typename Rhs>
struct SparseSparseProductReturnType;
87 template<typename Lhs, typename Rhs, int InnerSize = internal::traits<Lhs>::ColsAtCompileTime>
struct DenseSparseProductReturnType;
88 template<typename Lhs, typename Rhs, int InnerSize = internal::traits<Lhs>::ColsAtCompileTime>
struct SparseDenseProductReturnType;
89 template<
typename MatrixType,
int UpLo>
class SparseSymmetricPermutationProduct;
93 template<
typename T,
int Rows,
int Cols>
struct sparse_eval;
95 template<
typename T>
struct eval<T,
Sparse>
96 :
public sparse_eval<T, traits<T>::RowsAtCompileTime,traits<T>::ColsAtCompileTime>
99 template<
typename T,
int Cols>
struct sparse_eval<T,1,Cols> {
100 typedef typename traits<T>::Scalar _Scalar;
101 typedef typename traits<T>::Index _Index;
106 template<
typename T,
int Rows>
struct sparse_eval<T,Rows,1> {
107 typedef typename traits<T>::Scalar _Scalar;
108 typedef typename traits<T>::Index _Index;
110 typedef SparseVector<_Scalar, ColMajor, _Index> type;
113 template<
typename T,
int Rows,
int Cols>
struct sparse_eval {
114 typedef typename traits<T>::Scalar _Scalar;
115 typedef typename traits<T>::Index _Index;
118 typedef SparseMatrix<_Scalar, _Options, _Index> type;
121 template<
typename T>
struct sparse_eval<T,1,1> {
122 typedef typename traits<T>::Scalar _Scalar;
124 typedef Matrix<_Scalar, 1, 1> type;
127 template<
typename T>
struct plain_matrix_type<T,Sparse>
129 typedef typename traits<T>::Scalar _Scalar;
130 typedef typename traits<T>::Index _Index;
133 typedef SparseMatrix<_Scalar, _Options, _Index> type;
146 template<
typename Scalar,
typename Index=
unsigned int>
150 Triplet() : m_row(0), m_col(0), m_value(0) {}
152 Triplet(
const Index& i,
const Index& j,
const Scalar& v = Scalar(0))
153 : m_row(i), m_col(j), m_value(v)
157 const Index&
row()
const {
return m_row; }
160 const Index&
col()
const {
return m_col; }
163 const Scalar&
value()
const {
return m_value; }
171 #endif // EIGEN_SPARSEUTIL_H