11 #ifndef EIGEN_MATRIXBASE_H
12 #define EIGEN_MATRIXBASE_H
52 #ifndef EIGEN_PARSED_BY_DOXYGEN
54 typedef typename internal::traits<Derived>::StorageKind StorageKind;
55 typedef typename internal::traits<Derived>::Index
Index;
56 typedef typename internal::traits<Derived>::Scalar Scalar;
57 typedef typename internal::packet_traits<Scalar>::type PacketScalar;
72 using Base::const_cast_derived;
78 using Base::lazyAssign;
80 using Base::operator+=;
81 using Base::operator-=;
82 using Base::operator*=;
83 using Base::operator/=;
85 typedef typename Base::CoeffReturnType CoeffReturnType;
87 typedef typename Base::RowXpr RowXpr;
88 typedef typename Base::ColXpr ColXpr;
89 #endif // not EIGEN_PARSED_BY_DOXYGEN
93 #ifndef EIGEN_PARSED_BY_DOXYGEN
97 #endif // not EIGEN_PARSED_BY_DOXYGEN
110 internal::traits<Derived>::RowsAtCompileTime,
111 internal::traits<Derived>::ColsAtCompileTime,
113 internal::traits<Derived>::MaxRowsAtCompileTime,
114 internal::traits<Derived>::MaxColsAtCompileTime
117 #ifndef EIGEN_PARSED_BY_DOXYGEN
121 typedef typename internal::conditional<NumTraits<Scalar>::IsComplex,
124 >::type AdjointReturnType;
131 internal::traits<Derived>::RowsAtCompileTime,
132 internal::traits<Derived>::ColsAtCompileTime> BasisReturnType;
133 #endif // not EIGEN_PARSED_BY_DOXYGEN
135 #define EIGEN_CURRENT_STORAGE_BASE_CLASS Eigen::MatrixBase
136 # include "../plugins/CommonCwiseUnaryOps.h"
137 # include "../plugins/CommonCwiseBinaryOps.h"
138 # include "../plugins/MatrixCwiseUnaryOps.h"
139 # include "../plugins/MatrixCwiseBinaryOps.h"
140 # ifdef EIGEN_MATRIXBASE_PLUGIN
141 # include EIGEN_MATRIXBASE_PLUGIN
143 #undef EIGEN_CURRENT_STORAGE_BASE_CLASS
153 template <
typename OtherDerived>
156 template <
typename OtherDerived>
159 template<
typename OtherDerived>
160 Derived&
operator=(
const ReturnByValue<OtherDerived>& other);
162 #ifndef EIGEN_PARSED_BY_DOXYGEN
163 template<
typename ProductDerived,
typename Lhs,
typename Rhs>
164 Derived& lazyAssign(
const ProductBase<ProductDerived, Lhs,Rhs>& other);
166 template<
typename MatrixPower,
typename Lhs,
typename Rhs>
167 Derived& lazyAssign(
const MatrixPowerProduct<MatrixPower, Lhs,Rhs>& other);
168 #endif // not EIGEN_PARSED_BY_DOXYGEN
170 template<
typename OtherDerived>
172 template<
typename OtherDerived>
175 template<
typename OtherDerived>
179 template<
typename OtherDerived>
183 template<
typename OtherDerived>
186 template<
typename OtherDerived>
189 template<
typename OtherDerived>
192 template<
typename DiagonalDerived>
193 const DiagonalProduct<Derived, DiagonalDerived, OnTheRight>
196 template<
typename OtherDerived>
197 typename internal::scalar_product_traits<typename internal::traits<Derived>::Scalar,
typename internal::traits<OtherDerived>::Scalar>::ReturnType
200 #ifdef EIGEN2_SUPPORT
201 template<
typename OtherDerived>
206 RealScalar
norm()
const;
213 const AdjointReturnType
adjoint()
const;
224 template<
int Index>
typename DiagonalIndexReturnType<Index>::Type
diagonal();
225 template<
int Index>
typename ConstDiagonalIndexReturnType<Index>::Type
diagonal()
const;
229 #if (defined _MSC_VER) && (_MSC_VER >= 1500) // 2008 or later
230 typename MatrixBase::template DiagonalIndexReturnType<DynamicIndex>::Type
diagonal(
Index index);
231 typename MatrixBase::template ConstDiagonalIndexReturnType<DynamicIndex>::Type
diagonal(
Index index)
const;
233 typename DiagonalIndexReturnType<DynamicIndex>::Type
diagonal(
Index index);
234 typename ConstDiagonalIndexReturnType<DynamicIndex>::Type
diagonal(
Index index)
const;
237 #ifdef EIGEN2_SUPPORT
238 template<
unsigned int Mode>
typename internal::eigen2_part_return_type<Derived, Mode>::type part();
239 template<
unsigned int Mode>
const typename internal::eigen2_part_return_type<Derived, Mode>::type part()
const;
243 template<
template<
typename T,
int N>
class U>
244 const DiagonalWrapper<ConstDiagonalReturnType> part()
const
246 #endif // EIGEN2_SUPPORT
248 template<
unsigned int Mode>
struct TriangularViewReturnType {
typedef TriangularView<Derived, Mode> Type; };
249 template<
unsigned int Mode>
struct ConstTriangularViewReturnType {
typedef const TriangularView<const Derived, Mode> Type; };
251 template<
unsigned int Mode>
typename TriangularViewReturnType<Mode>::Type
triangularView();
252 template<
unsigned int Mode>
typename ConstTriangularViewReturnType<Mode>::Type
triangularView()
const;
254 template<
unsigned int UpLo>
struct SelfAdjointViewReturnType {
typedef SelfAdjointView<Derived, UpLo> Type; };
255 template<
unsigned int UpLo>
struct ConstSelfAdjointViewReturnType {
typedef const SelfAdjointView<const Derived, UpLo> Type; };
257 template<
unsigned int UpLo>
typename SelfAdjointViewReturnType<UpLo>::Type selfadjointView();
258 template<
unsigned int UpLo>
typename ConstSelfAdjointViewReturnType<UpLo>::Type selfadjointView()
const;
260 const SparseView<Derived> sparseView(
const Scalar& m_reference = Scalar(0),
261 const typename NumTraits<Scalar>::Real& m_epsilon = NumTraits<Scalar>::dummy_precision())
const;
262 static const IdentityReturnType
Identity();
265 static const BasisReturnType
Unit(
Index i);
266 static const BasisReturnType
UnitX();
267 static const BasisReturnType
UnitY();
268 static const BasisReturnType
UnitZ();
269 static const BasisReturnType
UnitW();
271 const DiagonalWrapper<const Derived>
asDiagonal()
const;
272 const PermutationWrapper<const Derived> asPermutation()
const;
277 bool isIdentity(
const RealScalar& prec = NumTraits<Scalar>::dummy_precision())
const;
278 bool isDiagonal(
const RealScalar& prec = NumTraits<Scalar>::dummy_precision())
const;
280 bool isUpperTriangular(
const RealScalar& prec = NumTraits<Scalar>::dummy_precision())
const;
281 bool isLowerTriangular(
const RealScalar& prec = NumTraits<Scalar>::dummy_precision())
const;
283 template<
typename OtherDerived>
284 bool isOrthogonal(
const MatrixBase<OtherDerived>& other,
285 const RealScalar& prec = NumTraits<Scalar>::dummy_precision())
const;
286 bool isUnitary(
const RealScalar& prec = NumTraits<Scalar>::dummy_precision())
const;
292 template<
typename OtherDerived>
300 template<
typename OtherDerived>
308 template<
bool Enable>
inline typename internal::add_const_on_value_type<typename internal::conditional<Enable,ForceAlignedAccess<Derived>,Derived&>::type>::type
forceAlignedAccessIf()
const;
309 template<
bool Enable>
inline typename internal::conditional<Enable,ForceAlignedAccess<Derived>,Derived&>::type
forceAlignedAccessIf();
311 Scalar
trace()
const;
315 template<
int p> RealScalar
lpNorm()
const;
318 const MatrixBase<Derived>& matrix()
const {
return *
this; }
327 const FullPivLU<PlainObject>
fullPivLu()
const;
330 #if EIGEN2_SUPPORT_STAGE < STAGE20_RESOLVE_API_CONFLICTS
331 const LU<PlainObject>
lu()
const;
334 #ifdef EIGEN2_SUPPORT
335 const LU<PlainObject> eigen2_lu()
const;
338 #if EIGEN2_SUPPORT_STAGE > STAGE20_RESOLVE_API_CONFLICTS
339 const PartialPivLU<PlainObject>
lu()
const;
342 #ifdef EIGEN2_SUPPORT
343 template<
typename ResultType>
344 void computeInverse(MatrixBase<ResultType> *result)
const {
349 const internal::inverse_impl<Derived>
inverse()
const;
350 template<
typename ResultType>
355 const RealScalar& absDeterminantThreshold = NumTraits<Scalar>::dummy_precision()
357 template<
typename ResultType>
361 const RealScalar& absDeterminantThreshold = NumTraits<Scalar>::dummy_precision()
367 const LLT<PlainObject>
llt()
const;
368 const LDLT<PlainObject>
ldlt()
const;
376 #ifdef EIGEN2_SUPPORT
377 const QR<PlainObject> qr()
const;
385 JacobiSVD<PlainObject>
jacobiSvd(
unsigned int computationOptions = 0)
const;
387 #ifdef EIGEN2_SUPPORT
388 SVD<PlainObject> svd()
const;
393 #ifndef EIGEN_PARSED_BY_DOXYGEN
395 template<
typename OtherDerived>
struct cross_product_return_type {
396 typedef typename internal::scalar_product_traits<typename internal::traits<Derived>::Scalar,
typename internal::traits<OtherDerived>::Scalar>::ReturnType Scalar;
397 typedef Matrix<Scalar,MatrixBase::RowsAtCompileTime,MatrixBase::ColsAtCompileTime> type;
399 #endif // EIGEN_PARSED_BY_DOXYGEN
400 template<
typename OtherDerived>
401 typename cross_product_return_type<OtherDerived>::type
402 cross(
const MatrixBase<OtherDerived>& other)
const;
403 template<
typename OtherDerived>
408 #if EIGEN2_SUPPORT_STAGE > STAGE20_RESOLVE_API_CONFLICTS
409 ScalarMultipleReturnType
operator*(
const UniformScaling<Scalar>& s)
const;
412 typedef Homogeneous<Derived, HomogeneousReturnTypeDirection> HomogeneousReturnType;
419 typedef Block<
const Derived,
420 internal::traits<Derived>::ColsAtCompileTime==1 ? SizeMinusOne : 1,
421 internal::traits<Derived>::ColsAtCompileTime==1 ? 1 : SizeMinusOne> ConstStartMinusOne;
422 typedef CwiseUnaryOp<internal::scalar_quotient1_op<typename internal::traits<Derived>::Scalar>,
423 const ConstStartMinusOne > HNormalizedReturnType;
430 template<
typename EssentialPart>
432 Scalar& tau, RealScalar& beta)
const;
433 template<
typename EssentialPart>
437 template<
typename EssentialPart>
444 template<
typename OtherScalar>
446 template<
typename OtherScalar>
451 typedef typename internal::stem_function<Scalar>::type StemFunction;
452 const MatrixExponentialReturnValue<Derived> exp()
const;
453 const MatrixFunctionReturnValue<Derived> matrixFunction(StemFunction f)
const;
454 const MatrixFunctionReturnValue<Derived> cosh()
const;
455 const MatrixFunctionReturnValue<Derived> sinh()
const;
456 const MatrixFunctionReturnValue<Derived> cos()
const;
457 const MatrixFunctionReturnValue<Derived> sin()
const;
458 const MatrixSquareRootReturnValue<Derived> sqrt()
const;
459 const MatrixLogarithmReturnValue<Derived> log()
const;
460 const MatrixPowerReturnValue<Derived> pow(
const RealScalar& p)
const;
462 #ifdef EIGEN2_SUPPORT
463 template<
typename ProductDerived,
typename Lhs,
typename Rhs>
464 Derived&
operator+=(
const Flagged<ProductBase<ProductDerived, Lhs,Rhs>, 0,
467 template<
typename ProductDerived,
typename Lhs,
typename Rhs>
468 Derived&
operator-=(
const Flagged<ProductBase<ProductDerived, Lhs,Rhs>, 0,
473 template<
typename OtherDerived>
474 Derived& lazyAssign(
const Flagged<OtherDerived, 0, EvalBeforeAssigningBit>& other)
475 {
return lazyAssign(other._expression()); }
477 template<
unsigned int Added>
478 const Flagged<Derived, Added, 0> marked()
const;
479 const Flagged<Derived, 0, EvalBeforeAssigningBit> lazy()
const;
481 inline const Cwise<Derived> cwise()
const;
482 inline Cwise<Derived> cwise();
484 VectorBlock<Derived> start(
Index size);
485 const VectorBlock<const Derived> start(
Index size)
const;
486 VectorBlock<Derived> end(
Index size);
487 const VectorBlock<const Derived> end(
Index size)
const;
488 template<
int Size> VectorBlock<Derived,Size> start();
489 template<
int Size>
const VectorBlock<const Derived,Size> start()
const;
490 template<
int Size> VectorBlock<Derived,Size> end();
491 template<
int Size>
const VectorBlock<const Derived,Size> end()
const;
498 MatrixBase() : Base() {}
501 explicit MatrixBase(
int);
503 template<
typename OtherDerived>
explicit MatrixBase(
const MatrixBase<OtherDerived>&);
506 template<
typename OtherDerived> Derived&
operator+=(
const ArrayBase<OtherDerived>& )
507 {EIGEN_STATIC_ASSERT(std::ptrdiff_t(
sizeof(
typename OtherDerived::Scalar))==-1,YOU_CANNOT_MIX_ARRAYS_AND_MATRICES);
return *
this;}
509 template<
typename OtherDerived> Derived&
operator-=(
const ArrayBase<OtherDerived>& )
510 {EIGEN_STATIC_ASSERT(std::ptrdiff_t(
sizeof(
typename OtherDerived::Scalar))==-1,YOU_CANNOT_MIX_ARRAYS_AND_MATRICES);
return *
this;}
515 #endif // EIGEN_MATRIXBASE_H