26 #ifndef EIGEN_TRANSPOSE_H
27 #define EIGEN_TRANSPOSE_H
46 template<
typename MatrixType>
47 struct traits<Transpose<MatrixType> > :
traits<MatrixType>
49 typedef typename MatrixType::Scalar Scalar;
50 typedef typename nested<MatrixType>::type MatrixTypeNested;
51 typedef typename remove_reference<MatrixTypeNested>::type MatrixTypeNestedPlain;
55 RowsAtCompileTime = MatrixType::ColsAtCompileTime,
56 ColsAtCompileTime = MatrixType::RowsAtCompileTime,
57 MaxRowsAtCompileTime = MatrixType::MaxColsAtCompileTime,
58 MaxColsAtCompileTime = MatrixType::MaxRowsAtCompileTime,
59 FlagsLvalueBit = is_lvalue<MatrixType>::value ?
LvalueBit : 0,
61 Flags1 = Flags0 | FlagsLvalueBit,
63 CoeffReadCost = MatrixTypeNestedPlain::CoeffReadCost,
64 InnerStrideAtCompileTime = inner_stride_at_compile_time<MatrixType>::ret,
65 OuterStrideAtCompileTime = outer_stride_at_compile_time<MatrixType>::ret
70 template<
typename MatrixType,
typename StorageKind>
class TransposeImpl;
73 :
public TransposeImpl<MatrixType,typename internal::traits<MatrixType>::StorageKind>
80 inline
Transpose(MatrixType& matrix) : m_matrix(matrix) {}
84 inline Index rows()
const {
return m_matrix.cols(); }
85 inline Index
cols()
const {
return m_matrix.rows(); }
88 const typename internal::remove_all<typename MatrixType::Nested>::type&
92 typename internal::remove_all<typename MatrixType::Nested>::type&
101 template<typename MatrixType, bool HasDirectAccess = has_direct_access<MatrixType>::ret>
102 struct TransposeImpl_base
104 typedef typename dense_xpr_base<Transpose<MatrixType> >::type type;
107 template<
typename MatrixType>
108 struct TransposeImpl_base<MatrixType, false>
110 typedef typename dense_xpr_base<Transpose<MatrixType> >::type type;
116 :
public internal::TransposeImpl_base<MatrixType>::type
120 typedef typename internal::TransposeImpl_base<MatrixType>::type
Base;
123 inline Index innerStride()
const {
return derived().nestedExpression().innerStride(); }
124 inline Index
outerStride()
const {
return derived().nestedExpression().outerStride(); }
126 typedef typename internal::conditional<
127 internal::is_lvalue<MatrixType>::value,
133 inline const Scalar*
data()
const {
return derived().nestedExpression().data(); }
138 return derived().nestedExpression().const_cast_derived().coeffRef(col, row);
144 return derived().nestedExpression().const_cast_derived().coeffRef(index);
147 inline const Scalar& coeffRef(Index
row, Index
col)
const
149 return derived().nestedExpression().coeffRef(col, row);
152 inline const Scalar& coeffRef(Index index)
const
154 return derived().nestedExpression().coeffRef(index);
157 inline CoeffReturnType coeff(Index
row, Index
col)
const
159 return derived().nestedExpression().coeff(col, row);
162 inline CoeffReturnType coeff(Index index)
const
164 return derived().nestedExpression().coeff(index);
167 template<
int LoadMode>
168 inline const PacketScalar packet(Index
row, Index
col)
const
170 return derived().nestedExpression().template packet<LoadMode>(
col,
row);
173 template<
int LoadMode>
174 inline void writePacket(Index
row, Index
col,
const PacketScalar& x)
176 derived().nestedExpression().const_cast_derived().template writePacket<LoadMode>(
col,
row, x);
179 template<
int LoadMode>
180 inline const PacketScalar packet(Index index)
const
182 return derived().nestedExpression().template packet<LoadMode>(index);
185 template<
int LoadMode>
186 inline void writePacket(Index index,
const PacketScalar& x)
188 derived().nestedExpression().const_cast_derived().template writePacket<LoadMode>(index, x);
211 template<
typename Derived>
212 inline Transpose<Derived>
223 template<
typename Derived>
249 template<
typename Derived>
253 return this->transpose();
263 template<
typename MatrixType,
264 bool IsSquare = (MatrixType::RowsAtCompileTime == MatrixType::ColsAtCompileTime) && MatrixType::RowsAtCompileTime!=
Dynamic>
265 struct inplace_transpose_selector;
267 template<
typename MatrixType>
268 struct inplace_transpose_selector<MatrixType,true> {
269 static void run(MatrixType& m) {
270 m.template triangularView<StrictlyUpper>().swap(m.transpose());
274 template<
typename MatrixType>
275 struct inplace_transpose_selector<MatrixType,false> {
276 static void run(MatrixType& m) {
277 if (m.rows()==m.cols())
278 m.template triangularView<StrictlyUpper>().swap(m.transpose());
280 m = m.transpose().eval();
304 template<
typename Derived>
307 internal::inplace_transpose_selector<Derived>::run(derived());
332 template<
typename Derived>
335 derived() = adjoint().
eval();
338 #ifndef EIGEN_NO_DEBUG
344 template<
typename BinOp,
typename NestedXpr,
typename Rhs>
345 struct blas_traits<SelfCwiseBinaryOp<BinOp,NestedXpr,Rhs> >
346 : blas_traits<NestedXpr>
349 static inline const XprType extract(
const XprType& x) {
return x; }
352 template<
bool DestIsTransposed,
typename OtherDerived>
353 struct check_transpose_aliasing_compile_time_selector
355 enum { ret =
bool(blas_traits<OtherDerived>::IsTransposed) != DestIsTransposed };
358 template<
bool DestIsTransposed,
typename BinOp,
typename DerivedA,
typename DerivedB>
359 struct check_transpose_aliasing_compile_time_selector<DestIsTransposed,CwiseBinaryOp<BinOp,DerivedA,DerivedB> >
361 enum { ret =
bool(blas_traits<DerivedA>::IsTransposed) != DestIsTransposed
362 ||
bool(blas_traits<DerivedB>::IsTransposed) != DestIsTransposed
366 template<
typename Scalar,
bool DestIsTransposed,
typename OtherDerived>
367 struct check_transpose_aliasing_run_time_selector
369 static bool run(
const Scalar* dest,
const OtherDerived& src)
371 return (
bool(blas_traits<OtherDerived>::IsTransposed) != DestIsTransposed) && (dest!=0 && dest==(Scalar*)
extract_data(src));
375 template<
typename Scalar,
bool DestIsTransposed,
typename BinOp,
typename DerivedA,
typename DerivedB>
376 struct check_transpose_aliasing_run_time_selector<Scalar,DestIsTransposed,CwiseBinaryOp<BinOp,DerivedA,DerivedB> >
378 static bool run(
const Scalar* dest,
const CwiseBinaryOp<BinOp,DerivedA,DerivedB>& src)
380 return ((blas_traits<DerivedA>::IsTransposed != DestIsTransposed) && (dest!=0 && dest==(Scalar*)
extract_data(src.lhs())))
381 || ((blas_traits<DerivedB>::IsTransposed != DestIsTransposed) && (dest!=0 && dest==(Scalar*)
extract_data(src.rhs())));
391 template<
typename Derived,
typename OtherDerived,
392 bool MightHaveTransposeAliasing
393 = check_transpose_aliasing_compile_time_selector
394 <blas_traits<Derived>::IsTransposed,OtherDerived>::ret
396 struct checkTransposeAliasing_impl
398 static void run(
const Derived& dst,
const OtherDerived& other)
400 eigen_assert((!check_transpose_aliasing_run_time_selector
401 <
typename Derived::Scalar,blas_traits<Derived>::IsTransposed,OtherDerived>
403 &&
"aliasing detected during tranposition, use transposeInPlace() "
404 "or evaluate the rhs into a temporary using .eval()");
409 template<
typename Derived,
typename OtherDerived>
410 struct checkTransposeAliasing_impl<Derived, OtherDerived, false>
412 static void run(
const Derived&,
const OtherDerived&)
419 template<
typename Derived>
420 template<
typename OtherDerived>
423 internal::checkTransposeAliasing_impl<Derived, OtherDerived>::run(derived(), other);
429 #endif // EIGEN_TRANSPOSE_H