63template<
typename MatrixType,
int _DiagIndex>
class Diagonal
64 :
public internal::dense_xpr_base< Diagonal<MatrixType,_DiagIndex> >::type
68 enum { DiagIndex = _DiagIndex };
69 typedef typename internal::dense_xpr_base<Diagonal>::type Base;
70 EIGEN_DENSE_PUBLIC_INTERFACE(
Diagonal)
73 explicit inline Diagonal(MatrixType& matrix,
Index a_index = DiagIndex) : m_matrix(matrix), m_index(a_index)
75 eigen_assert( a_index <= m_matrix.cols() && -a_index <= m_matrix.rows() );
78 EIGEN_INHERIT_ASSIGNMENT_OPERATORS(
Diagonal)
81 inline Index rows()
const
83 return m_index.value()<0 ? numext::mini<Index>(m_matrix.cols(),m_matrix.rows()+m_index.value())
84 : numext::mini<Index>(m_matrix.rows(),m_matrix.cols()-m_index.value());
87 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
88 inline Index cols()
const EIGEN_NOEXCEPT {
return 1; }
90 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
91 inline Index innerStride()
const EIGEN_NOEXCEPT {
92 return m_matrix.outerStride() + 1;
95 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
96 inline Index outerStride()
const EIGEN_NOEXCEPT {
return 0; }
98 typedef typename internal::conditional<
99 internal::is_lvalue<MatrixType>::value,
102 >::type ScalarWithConstIfNotLvalue;
105 inline ScalarWithConstIfNotLvalue* data() {
return &(m_matrix.coeffRef(rowOffset(), colOffset())); }
107 inline const Scalar* data()
const {
return &(m_matrix.coeffRef(rowOffset(), colOffset())); }
112 EIGEN_STATIC_ASSERT_LVALUE(MatrixType)
113 return m_matrix.coeffRef(row+rowOffset(), row+colOffset());
117 inline const Scalar& coeffRef(
Index row,
Index)
const
119 return m_matrix.coeffRef(row+rowOffset(), row+colOffset());
123 inline CoeffReturnType coeff(
Index row,
Index)
const
125 return m_matrix.coeff(row+rowOffset(), row+colOffset());
129 inline Scalar& coeffRef(
Index idx)
131 EIGEN_STATIC_ASSERT_LVALUE(MatrixType)
132 return m_matrix.coeffRef(idx+rowOffset(), idx+colOffset());
136 inline const Scalar& coeffRef(
Index idx)
const
138 return m_matrix.coeffRef(idx+rowOffset(), idx+colOffset());
142 inline CoeffReturnType coeff(
Index idx)
const
144 return m_matrix.coeff(idx+rowOffset(), idx+colOffset());
148 inline const typename internal::remove_all<typename MatrixType::Nested>::type&
149 nestedExpression()
const
155 inline Index index()
const
157 return m_index.value();
161 typename internal::ref_selector<MatrixType>::non_const_type m_matrix;
162 const internal::variable_if_dynamicindex<Index, DiagIndex> m_index;
166 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE EIGEN_CONSTEXPR
167 Index absDiagIndex()
const EIGEN_NOEXCEPT {
return m_index.value()>0 ? m_index.value() : -m_index.value(); }
168 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE EIGEN_CONSTEXPR
169 Index rowOffset()
const EIGEN_NOEXCEPT {
return m_index.value()>0 ? 0 : -m_index.value(); }
170 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE EIGEN_CONSTEXPR
171 Index colOffset()
const EIGEN_NOEXCEPT {
return m_index.value()>0 ? m_index.value() : 0; }
173 template<
int LoadMode>
typename MatrixType::PacketReturnType packet(
Index)
const;
174 template<
int LoadMode>
typename MatrixType::PacketReturnType packet(
Index,
Index)
const;