11 #ifndef EIGEN_XPRHELPER_H
12 #define EIGEN_XPRHELPER_H
17 #if (defined __GNUG__) && !((__GNUC__==4) && (__GNUC_MINOR__==3))
18 #define EIGEN_EMPTY_STRUCT_CTOR(X) \
19 EIGEN_STRONG_INLINE X() {} \
20 EIGEN_STRONG_INLINE X(const X& ) {}
22 #define EIGEN_EMPTY_STRUCT_CTOR(X)
27 typedef EIGEN_DEFAULT_DENSE_INDEX_TYPE DenseIndex;
32 class no_assignment_operator
35 no_assignment_operator& operator=(
const no_assignment_operator&);
39 template<
typename I1,
typename I2>
40 struct promote_index_type
42 typedef typename conditional<(sizeof(I1)<sizeof(I2)), I2, I1>::type type;
49 template<
typename T,
int Value>
class variable_if_dynamic
52 EIGEN_EMPTY_STRUCT_CTOR(variable_if_dynamic)
53 explicit variable_if_dynamic(T v) { EIGEN_ONLY_USED_FOR_DEBUG(v); assert(v == T(Value)); }
54 static T value() {
return T(Value); }
58 template<
typename T>
class variable_if_dynamic<T,
Dynamic>
61 variable_if_dynamic() { assert(
false); }
63 explicit variable_if_dynamic(T value) : m_value(value) {}
64 T value()
const {
return m_value; }
65 void setValue(T value) { m_value = value; }
68 template<
typename T>
struct functor_traits
77 template<
typename T>
struct packet_traits;
79 template<
typename T>
struct unpacket_traits
85 template<
typename _Scalar,
int _Rows,
int _Cols,
86 int _Options = AutoAlign |
89 : EIGEN_DEFAULT_MATRIX_STORAGE_ORDER_OPTION ),
92 >
class make_proper_matrix_type
95 IsColVector = _Cols==1 && _Rows!=1,
96 IsRowVector = _Rows==1 && _Cols!=1,
102 typedef Matrix<_Scalar, _Rows, _Cols, Options, _MaxRows, _MaxCols> type;
105 template<
typename Scalar,
int Rows,
int Cols,
int Options,
int MaxRows,
int MaxCols>
106 class compute_matrix_flags
114 ((Options&DontAlign)==0)
116 #if EIGEN_ALIGN_STATICALLY
117 ((!is_dynamic_size_storage) && (((MaxCols*MaxRows*
int(
sizeof(Scalar))) % 16) == 0))
125 is_dynamic_size_storage
132 packet_access_bit = packet_traits<Scalar>::Vectorizable && aligned_bit ?
PacketAccessBit : 0
139 template<
int _Rows,
int _Cols>
struct size_at_compile_time
148 template<typename T, typename StorageKind = typename traits<T>::StorageKind>
struct plain_matrix_type;
149 template<
typename T,
typename BaseClassType>
struct plain_matrix_type_dense;
150 template<
typename T>
struct plain_matrix_type<T,Dense>
152 typedef typename plain_matrix_type_dense<T,typename traits<T>::XprKind>::type type;
155 template<
typename T>
struct plain_matrix_type_dense<T,MatrixXpr>
157 typedef Matrix<typename traits<T>::Scalar,
158 traits<T>::RowsAtCompileTime,
159 traits<T>::ColsAtCompileTime,
161 traits<T>::MaxRowsAtCompileTime,
162 traits<T>::MaxColsAtCompileTime
166 template<
typename T>
struct plain_matrix_type_dense<T,ArrayXpr>
168 typedef Array<typename traits<T>::Scalar,
169 traits<T>::RowsAtCompileTime,
170 traits<T>::ColsAtCompileTime,
172 traits<T>::MaxRowsAtCompileTime,
173 traits<T>::MaxColsAtCompileTime
181 template<typename T, typename StorageKind = typename traits<T>::StorageKind>
struct eval;
183 template<
typename T>
struct eval<T,Dense>
185 typedef typename plain_matrix_type<T>::type type;
197 template<
typename _Scalar,
int _Rows,
int _Cols,
int _Options,
int _MaxRows,
int _MaxCols>
198 struct eval<Matrix<_Scalar, _Rows, _Cols, _Options, _MaxRows, _MaxCols>, Dense>
200 typedef const Matrix<_Scalar, _Rows, _Cols, _Options, _MaxRows, _MaxCols>& type;
203 template<
typename _Scalar,
int _Rows,
int _Cols,
int _Options,
int _MaxRows,
int _MaxCols>
204 struct eval<Array<_Scalar, _Rows, _Cols, _Options, _MaxRows, _MaxCols>, Dense>
206 typedef const Array<_Scalar, _Rows, _Cols, _Options, _MaxRows, _MaxCols>& type;
213 template<
typename T>
struct plain_matrix_type_column_major
215 enum { Rows = traits<T>::RowsAtCompileTime,
216 Cols = traits<T>::ColsAtCompileTime,
217 MaxRows = traits<T>::MaxRowsAtCompileTime,
218 MaxCols = traits<T>::MaxColsAtCompileTime
220 typedef Matrix<typename traits<T>::Scalar,
231 template<
typename T>
struct plain_matrix_type_row_major
233 enum { Rows = traits<T>::RowsAtCompileTime,
234 Cols = traits<T>::ColsAtCompileTime,
235 MaxRows = traits<T>::MaxRowsAtCompileTime,
236 MaxCols = traits<T>::MaxColsAtCompileTime
238 typedef Matrix<typename traits<T>::Scalar,
248 template<
typename T>
struct must_nest_by_value {
enum { ret =
false }; };
253 template <
typename T>
256 typedef typename conditional<
257 bool(traits<T>::Flags & NestByRefBit),
264 template<
typename T1,
typename T2>
265 struct transfer_constness
267 typedef typename conditional<
268 bool(internal::is_const<T1>::value),
269 typename internal::add_const_on_value_type<T2>::type,
294 template<typename T, int n=1, typename PlainObject = typename eval<T>::type>
struct nested
302 DynamicAsInteger = 10000,
303 ScalarReadCost = NumTraits<typename traits<T>::Scalar>::ReadCost,
304 ScalarReadCostAsInteger = ScalarReadCost ==
Dynamic ? int(DynamicAsInteger) : int(ScalarReadCost),
305 CoeffReadCost = traits<T>::CoeffReadCost,
306 CoeffReadCostAsInteger = CoeffReadCost ==
Dynamic ? int(DynamicAsInteger) : int(CoeffReadCost),
307 NAsInteger = n ==
Dynamic ? int(DynamicAsInteger) : n,
308 CostEvalAsInteger = (NAsInteger+1) * ScalarReadCostAsInteger + CoeffReadCostAsInteger,
309 CostNoEvalAsInteger = NAsInteger * CoeffReadCostAsInteger
312 typedef typename conditional<
314 int(CostEvalAsInteger) < int(CostNoEvalAsInteger)
317 typename ref_selector<T>::type
322 T* const_cast_ptr(
const T* ptr)
324 return const_cast<T*
>(ptr);
327 template<typename Derived, typename XprKind = typename traits<Derived>::XprKind>
328 struct dense_xpr_base
333 template<
typename Derived>
334 struct dense_xpr_base<Derived, MatrixXpr>
336 typedef MatrixBase<Derived> type;
339 template<
typename Derived>
340 struct dense_xpr_base<Derived, ArrayXpr>
342 typedef ArrayBase<Derived> type;
347 template<
typename Derived,
typename Scalar,
typename OtherScalar,
348 bool EnableIt = !is_same<Scalar,OtherScalar>::value >
349 struct special_scalar_op_base :
public DenseCoeffsBase<Derived>
356 template<
typename Derived,
typename Scalar,
typename OtherScalar>
357 struct special_scalar_op_base<Derived,Scalar,OtherScalar,true> :
public DenseCoeffsBase<Derived>
359 const CwiseUnaryOp<scalar_multiple2_op<Scalar,OtherScalar>, Derived>
360 operator*(
const OtherScalar& scalar)
const
362 return CwiseUnaryOp<scalar_multiple2_op<Scalar,OtherScalar>, Derived>
363 (*
static_cast<const Derived*
>(
this), scalar_multiple2_op<Scalar,OtherScalar>(scalar));
366 inline friend const CwiseUnaryOp<scalar_multiple2_op<Scalar,OtherScalar>, Derived>
367 operator*(
const OtherScalar& scalar,
const Derived& matrix)
368 {
return static_cast<const special_scalar_op_base&
>(matrix).
operator*(scalar); }
371 template<
typename XprType,
typename CastType>
struct cast_return_type
373 typedef typename XprType::Scalar CurrentScalarType;
374 typedef typename remove_all<CastType>::type _CastType;
375 typedef typename _CastType::Scalar NewScalarType;
376 typedef typename conditional<is_same<CurrentScalarType,NewScalarType>::value,
377 const XprType&,CastType>::type type;
380 template <
typename A,
typename B>
struct promote_storage_type;
382 template <
typename A>
struct promote_storage_type<A,A>
390 template<
typename ExpressionType,
typename Scalar =
typename ExpressionType::Scalar>
391 struct plain_row_type
393 typedef Matrix<Scalar, 1, ExpressionType::ColsAtCompileTime,
394 ExpressionType::PlainObject::Options |
RowMajor, 1, ExpressionType::MaxColsAtCompileTime> MatrixRowType;
395 typedef Array<Scalar, 1, ExpressionType::ColsAtCompileTime,
396 ExpressionType::PlainObject::Options |
RowMajor, 1, ExpressionType::MaxColsAtCompileTime> ArrayRowType;
398 typedef typename conditional<
399 is_same< typename traits<ExpressionType>::XprKind, MatrixXpr >::value,
405 template<
typename ExpressionType,
typename Scalar =
typename ExpressionType::Scalar>
406 struct plain_col_type
408 typedef Matrix<Scalar, ExpressionType::RowsAtCompileTime, 1,
409 ExpressionType::PlainObject::Options & ~
RowMajor, ExpressionType::MaxRowsAtCompileTime, 1> MatrixColType;
410 typedef Array<Scalar, ExpressionType::RowsAtCompileTime, 1,
411 ExpressionType::PlainObject::Options & ~
RowMajor, ExpressionType::MaxRowsAtCompileTime, 1> ArrayColType;
413 typedef typename conditional<
414 is_same< typename traits<ExpressionType>::XprKind, MatrixXpr >::value,
420 template<
typename ExpressionType,
typename Scalar =
typename ExpressionType::Scalar>
421 struct plain_diag_type
423 enum { diag_size = EIGEN_SIZE_MIN_PREFER_DYNAMIC(ExpressionType::RowsAtCompileTime, ExpressionType::ColsAtCompileTime),
424 max_diag_size = EIGEN_SIZE_MIN_PREFER_FIXED(ExpressionType::MaxRowsAtCompileTime, ExpressionType::MaxColsAtCompileTime)
426 typedef Matrix<Scalar, diag_size, 1, ExpressionType::PlainObject::Options & ~RowMajor, max_diag_size, 1> MatrixDiagType;
427 typedef Array<Scalar, diag_size, 1, ExpressionType::PlainObject::Options & ~RowMajor, max_diag_size, 1> ArrayDiagType;
429 typedef typename conditional<
430 is_same< typename traits<ExpressionType>::XprKind, MatrixXpr >::value,
436 template<
typename ExpressionType>
439 enum { value = !bool(is_const<ExpressionType>::value) &&
440 bool(traits<ExpressionType>::Flags &
LvalueBit) };
447 #endif // EIGEN_XPRHELPER_H