11 #ifndef EIGEN_MAPBASE_H
12 #define EIGEN_MAPBASE_H
14 #define EIGEN_STATIC_ASSERT_INDEX_BASED_ACCESS(Derived) \
15 EIGEN_STATIC_ASSERT((int(internal::traits<Derived>::Flags) & LinearAccessBit) || Derived::IsVectorAtCompileTime, \
16 YOU_ARE_TRYING_TO_USE_AN_INDEX_BASED_ACCESSOR_ON_AN_EXPRESSION_THAT_DOES_NOT_SUPPORT_THAT)
28 :
public internal::dense_xpr_base<Derived>::type
32 typedef typename internal::dense_xpr_base<Derived>::type Base;
34 RowsAtCompileTime = internal::traits<Derived>::RowsAtCompileTime,
35 ColsAtCompileTime = internal::traits<Derived>::ColsAtCompileTime,
36 SizeAtCompileTime = Base::SizeAtCompileTime
39 typedef typename internal::traits<Derived>::StorageKind StorageKind;
40 typedef typename internal::traits<Derived>::Index Index;
41 typedef typename internal::traits<Derived>::Scalar Scalar;
42 typedef typename internal::packet_traits<Scalar>::type PacketScalar;
43 typedef typename NumTraits<Scalar>::Real RealScalar;
44 typedef typename internal::conditional<
45 bool(internal::is_lvalue<Derived>::value),
54 using Base::MaxRowsAtCompileTime;
55 using Base::MaxColsAtCompileTime;
56 using Base::MaxSizeAtCompileTime;
57 using Base::IsVectorAtCompileTime;
59 using Base::IsRowMajor;
66 using Base::lazyAssign;
69 using Base::innerStride;
70 using Base::outerStride;
71 using Base::rowStride;
72 using Base::colStride;
75 using Base::operator=;
77 typedef typename Base::CoeffReturnType CoeffReturnType;
79 inline Index rows()
const {
return m_rows.value(); }
80 inline Index cols()
const {
return m_cols.value(); }
88 inline const Scalar* data()
const {
return m_data; }
90 inline const Scalar& coeff(Index row, Index col)
const
92 return m_data[col * colStride() + row * rowStride()];
95 inline const Scalar& coeff(Index index)
const
97 EIGEN_STATIC_ASSERT_INDEX_BASED_ACCESS(Derived)
98 return m_data[index * innerStride()];
101 inline const Scalar& coeffRef(Index row, Index col)
const
103 return this->m_data[col * colStride() + row * rowStride()];
106 inline const Scalar& coeffRef(Index index)
const
108 EIGEN_STATIC_ASSERT_INDEX_BASED_ACCESS(Derived)
109 return this->m_data[index * innerStride()];
112 template<
int LoadMode>
113 inline PacketScalar packet(Index row, Index col)
const
115 return internal::ploadt<PacketScalar, LoadMode>
116 (m_data + (col * colStride() + row * rowStride()));
119 template<
int LoadMode>
120 inline PacketScalar packet(Index index)
const
122 EIGEN_STATIC_ASSERT_INDEX_BASED_ACCESS(Derived)
123 return internal::ploadt<PacketScalar, LoadMode>(m_data + index * innerStride());
126 inline MapBase(PointerType data) : m_data(data), m_rows(RowsAtCompileTime), m_cols(ColsAtCompileTime)
128 EIGEN_STATIC_ASSERT_FIXED_SIZE(Derived)
132 inline MapBase(PointerType data, Index size)
134 m_rows(RowsAtCompileTime ==
Dynamic ? size : Index(RowsAtCompileTime)),
135 m_cols(ColsAtCompileTime ==
Dynamic ? size : Index(ColsAtCompileTime))
137 EIGEN_STATIC_ASSERT_VECTOR_ONLY(Derived)
138 eigen_assert(size >= 0);
139 eigen_assert(data == 0 || SizeAtCompileTime ==
Dynamic || SizeAtCompileTime == size);
143 inline MapBase(PointerType data, Index rows, Index cols)
144 : m_data(data), m_rows(rows), m_cols(cols)
146 eigen_assert( (data == 0)
147 || ( rows >= 0 && (RowsAtCompileTime ==
Dynamic || RowsAtCompileTime == rows)
148 && cols >= 0 && (ColsAtCompileTime ==
Dynamic || ColsAtCompileTime == cols)));
154 void checkSanity()
const
156 EIGEN_STATIC_ASSERT(EIGEN_IMPLIES(internal::traits<Derived>::Flags&
PacketAccessBit,
157 internal::inner_stride_at_compile_time<Derived>::ret==1),
158 PACKET_ACCESS_REQUIRES_TO_HAVE_INNER_STRIDE_FIXED_TO_1);
159 eigen_assert(EIGEN_IMPLIES(internal::traits<Derived>::Flags&
AlignedBit, (
size_t(m_data) % 16) == 0)
160 &&
"data is not aligned");
164 const internal::variable_if_dynamic<Index, RowsAtCompileTime> m_rows;
165 const internal::variable_if_dynamic<Index, ColsAtCompileTime> m_cols;
169 :
public MapBase<Derived, ReadOnlyAccessors>
173 typedef MapBase<Derived, ReadOnlyAccessors> Base;
175 typedef typename Base::Scalar Scalar;
176 typedef typename Base::PacketScalar PacketScalar;
177 typedef typename Base::Index Index;
178 typedef typename Base::PointerType PointerType;
185 using Base::coeffRef;
187 using Base::innerStride;
188 using Base::outerStride;
189 using Base::rowStride;
190 using Base::colStride;
192 typedef typename internal::conditional<
193 internal::is_lvalue<Derived>::value,
196 >::type ScalarWithConstIfNotLvalue;
198 inline const Scalar* data()
const {
return this->m_data; }
199 inline ScalarWithConstIfNotLvalue* data() {
return this->m_data; }
201 inline ScalarWithConstIfNotLvalue& coeffRef(Index row, Index col)
203 return this->m_data[col * colStride() + row * rowStride()];
206 inline ScalarWithConstIfNotLvalue& coeffRef(Index index)
208 EIGEN_STATIC_ASSERT_INDEX_BASED_ACCESS(Derived)
209 return this->m_data[index * innerStride()];
212 template<
int StoreMode>
213 inline
void writePacket(Index row, Index col, const PacketScalar& x)
215 internal::pstoret<Scalar, PacketScalar, StoreMode>
216 (this->m_data + (col * colStride() + row * rowStride()), x);
219 template<
int StoreMode>
220 inline void writePacket(Index index,
const PacketScalar& x)
222 EIGEN_STATIC_ASSERT_INDEX_BASED_ACCESS(Derived)
223 internal::pstoret<Scalar, PacketScalar, StoreMode>
224 (this->m_data + index * innerStride(), x);
227 explicit inline MapBase(PointerType data) : Base(data) {}
228 inline MapBase(PointerType data, Index size) : Base(data, size) {}
229 inline MapBase(PointerType data, Index rows, Index cols) : Base(data, rows, cols) {}
231 Derived& operator=(
const MapBase& other)
233 Base::Base::operator=(other);
237 using Base::Base::operator=;
242 #endif // EIGEN_MAPBASE_H