diff --git a/Eigen/src/Eigen2Support/CwiseOperators.h b/Eigen/src/Eigen2Support/CwiseOperators.h index ac1dc0a7e..8c53dd599 100644 --- a/Eigen/src/Eigen2Support/CwiseOperators.h +++ b/Eigen/src/Eigen2Support/CwiseOperators.h @@ -30,13 +30,7 @@ ***************************************************************************/ -/** \returns an expression of the coefficient-wise absolute value of \c *this - * - * Example: \include Cwise_abs.cpp - * Output: \verbinclude Cwise_abs.out - * - * \sa abs2() - */ +/** \deprecated ArrayBase::abs() */ template EIGEN_STRONG_INLINE const EIGEN_CWISE_UNOP_RETURN_TYPE(ei_scalar_abs_op) Cwise::abs() const @@ -44,13 +38,7 @@ Cwise::abs() const return _expression(); } -/** \returns an expression of the coefficient-wise squared absolute value of \c *this - * - * Example: \include Cwise_abs2.cpp - * Output: \verbinclude Cwise_abs2.out - * - * \sa abs(), square() - */ +/** \deprecated ArrayBase::abs2() */ template EIGEN_STRONG_INLINE const EIGEN_CWISE_UNOP_RETURN_TYPE(ei_scalar_abs2_op) Cwise::abs2() const @@ -58,13 +46,7 @@ Cwise::abs2() const return _expression(); } -/** \returns an expression of the coefficient-wise exponential of *this. - * - * Example: \include Cwise_exp.cpp - * Output: \verbinclude Cwise_exp.out - * - * \sa pow(), log(), sin(), cos() - */ +/** \deprecated ArrayBase::exp() */ template inline const EIGEN_CWISE_UNOP_RETURN_TYPE(ei_scalar_exp_op) Cwise::exp() const @@ -72,13 +54,7 @@ Cwise::exp() const return _expression(); } -/** \returns an expression of the coefficient-wise logarithm of *this. - * - * Example: \include Cwise_log.cpp - * Output: \verbinclude Cwise_log.out - * - * \sa exp() - */ +/** \deprecated ArrayBase::log() */ template inline const EIGEN_CWISE_UNOP_RETURN_TYPE(ei_scalar_log_op) Cwise::log() const @@ -86,13 +62,7 @@ Cwise::log() const return _expression(); } -/** \returns an expression of the Schur product (coefficient wise product) of *this and \a other - * - * Example: \include Cwise_product.cpp - * Output: \verbinclude Cwise_product.out - * - * \sa class CwiseBinaryOp, operator/(), square() - */ +/** \deprecated ArrayBase::operator*() */ template template EIGEN_STRONG_INLINE const EIGEN_CWISE_PRODUCT_RETURN_TYPE(ExpressionType,OtherDerived) @@ -101,13 +71,7 @@ Cwise::operator*(const MatrixBase &other) const return EIGEN_CWISE_PRODUCT_RETURN_TYPE(ExpressionType,OtherDerived)(_expression(), other.derived()); } -/** \returns an expression of the coefficient-wise quotient of *this and \a other - * - * Example: \include Cwise_quotient.cpp - * Output: \verbinclude Cwise_quotient.out - * - * \sa class CwiseBinaryOp, operator*(), inverse() - */ +/** \deprecated ArrayBase::operator/() */ template template EIGEN_STRONG_INLINE const EIGEN_CWISE_BINOP_RETURN_TYPE(ei_scalar_quotient_op) @@ -116,13 +80,7 @@ Cwise::operator/(const MatrixBase &other) const return EIGEN_CWISE_BINOP_RETURN_TYPE(ei_scalar_quotient_op)(_expression(), other.derived()); } -/** Replaces this expression by its coefficient-wise product with \a other. - * - * Example: \include Cwise_times_equal.cpp - * Output: \verbinclude Cwise_times_equal.out - * - * \sa operator*(), operator/=() - */ +/** \deprecated ArrayBase::operator*=() */ template template inline ExpressionType& Cwise::operator*=(const MatrixBase &other) @@ -130,13 +88,7 @@ inline ExpressionType& Cwise::operator*=(const MatrixBase template inline ExpressionType& Cwise::operator/=(const MatrixBase &other) @@ -144,13 +96,7 @@ inline ExpressionType& Cwise::operator/=(const MatrixBase template EIGEN_STRONG_INLINE const EIGEN_CWISE_BINOP_RETURN_TYPE(ei_scalar_min_op) @@ -159,13 +105,7 @@ Cwise::min(const MatrixBase &other) const return EIGEN_CWISE_BINOP_RETURN_TYPE(ei_scalar_min_op)(_expression(), other.derived()); } -/** \returns an expression of the coefficient-wise max of *this and \a other - * - * Example: \include Cwise_max.cpp - * Output: \verbinclude Cwise_max.out - * - * \sa class CwiseBinaryOp - */ +/** \deprecated ArrayBase::max() */ template template EIGEN_STRONG_INLINE const EIGEN_CWISE_BINOP_RETURN_TYPE(ei_scalar_max_op) @@ -180,15 +120,7 @@ Cwise::max(const MatrixBase &other) const // -- unary operators -- -/** \array_module - * - * \returns an expression of the coefficient-wise square root of *this. - * - * Example: \include Cwise_sqrt.cpp - * Output: \verbinclude Cwise_sqrt.out - * - * \sa pow(), square() - */ +/** \deprecated ArrayBase::sqrt() */ template inline const EIGEN_CWISE_UNOP_RETURN_TYPE(ei_scalar_sqrt_op) Cwise::sqrt() const @@ -196,15 +128,7 @@ Cwise::sqrt() const return _expression(); } -/** \array_module - * - * \returns an expression of the coefficient-wise cosine of *this. - * - * Example: \include Cwise_cos.cpp - * Output: \verbinclude Cwise_cos.out - * - * \sa sin(), exp(), EIGEN_FAST_MATH - */ +/** \deprecated ArrayBase::cos() */ template inline const EIGEN_CWISE_UNOP_RETURN_TYPE(ei_scalar_cos_op) Cwise::cos() const @@ -213,15 +137,7 @@ Cwise::cos() const } -/** \array_module - * - * \returns an expression of the coefficient-wise sine of *this. - * - * Example: \include Cwise_sin.cpp - * Output: \verbinclude Cwise_sin.out - * - * \sa cos(), exp(), EIGEN_FAST_MATH - */ +/** \deprecated ArrayBase::sin() */ template inline const EIGEN_CWISE_UNOP_RETURN_TYPE(ei_scalar_sin_op) Cwise::sin() const @@ -230,15 +146,7 @@ Cwise::sin() const } -/** \array_module - * - * \returns an expression of the coefficient-wise power of *this to the given exponent. - * - * Example: \include Cwise_pow.cpp - * Output: \verbinclude Cwise_pow.out - * - * \sa exp(), log() - */ +/** \deprecated ArrayBase::log() */ template inline const EIGEN_CWISE_UNOP_RETURN_TYPE(ei_scalar_pow_op) Cwise::pow(const Scalar& exponent) const @@ -247,15 +155,7 @@ Cwise::pow(const Scalar& exponent) const } -/** \array_module - * - * \returns an expression of the coefficient-wise inverse of *this. - * - * Example: \include Cwise_inverse.cpp - * Output: \verbinclude Cwise_inverse.out - * - * \sa operator/(), operator*() - */ +/** \deprecated ArrayBase::inverse() */ template inline const EIGEN_CWISE_UNOP_RETURN_TYPE(ei_scalar_inverse_op) Cwise::inverse() const @@ -263,15 +163,7 @@ Cwise::inverse() const return _expression(); } -/** \array_module - * - * \returns an expression of the coefficient-wise square of *this. - * - * Example: \include Cwise_square.cpp - * Output: \verbinclude Cwise_square.out - * - * \sa operator/(), operator*(), abs2() - */ +/** \deprecated ArrayBase::square() */ template inline const EIGEN_CWISE_UNOP_RETURN_TYPE(ei_scalar_square_op) Cwise::square() const @@ -279,15 +171,7 @@ Cwise::square() const return _expression(); } -/** \array_module - * - * \returns an expression of the coefficient-wise cube of *this. - * - * Example: \include Cwise_cube.cpp - * Output: \verbinclude Cwise_cube.out - * - * \sa square(), pow() - */ +/** \deprecated ArrayBase::cube() */ template inline const EIGEN_CWISE_UNOP_RETURN_TYPE(ei_scalar_cube_op) Cwise::cube() const @@ -298,15 +182,7 @@ Cwise::cube() const // -- binary operators -- -/** \array_module - * - * \returns an expression of the coefficient-wise \< operator of *this and \a other - * - * Example: \include Cwise_less.cpp - * Output: \verbinclude Cwise_less.out - * - * \sa MatrixBase::all(), MatrixBase::any(), operator>(), operator<=() - */ +/** \deprecated ArrayBase::operator<() */ template template inline const EIGEN_CWISE_BINOP_RETURN_TYPE(std::less) @@ -315,15 +191,7 @@ Cwise::operator<(const MatrixBase &other) const return EIGEN_CWISE_BINOP_RETURN_TYPE(std::less)(_expression(), other.derived()); } -/** \array_module - * - * \returns an expression of the coefficient-wise \<= operator of *this and \a other - * - * Example: \include Cwise_less_equal.cpp - * Output: \verbinclude Cwise_less_equal.out - * - * \sa MatrixBase::all(), MatrixBase::any(), operator>=(), operator<() - */ +/** \deprecated ArrayBase::<=() */ template template inline const EIGEN_CWISE_BINOP_RETURN_TYPE(std::less_equal) @@ -332,15 +200,7 @@ Cwise::operator<=(const MatrixBase &other) const return EIGEN_CWISE_BINOP_RETURN_TYPE(std::less_equal)(_expression(), other.derived()); } -/** \array_module - * - * \returns an expression of the coefficient-wise \> operator of *this and \a other - * - * Example: \include Cwise_greater.cpp - * Output: \verbinclude Cwise_greater.out - * - * \sa MatrixBase::all(), MatrixBase::any(), operator>=(), operator<() - */ +/** \deprecated ArrayBase::operator>() */ template template inline const EIGEN_CWISE_BINOP_RETURN_TYPE(std::greater) @@ -349,15 +209,7 @@ Cwise::operator>(const MatrixBase &other) const return EIGEN_CWISE_BINOP_RETURN_TYPE(std::greater)(_expression(), other.derived()); } -/** \array_module - * - * \returns an expression of the coefficient-wise \>= operator of *this and \a other - * - * Example: \include Cwise_greater_equal.cpp - * Output: \verbinclude Cwise_greater_equal.out - * - * \sa MatrixBase::all(), MatrixBase::any(), operator>(), operator<=() - */ +/** \deprecated ArrayBase::operator>=() */ template template inline const EIGEN_CWISE_BINOP_RETURN_TYPE(std::greater_equal) @@ -366,20 +218,7 @@ Cwise::operator>=(const MatrixBase &other) const return EIGEN_CWISE_BINOP_RETURN_TYPE(std::greater_equal)(_expression(), other.derived()); } -/** \array_module - * - * \returns an expression of the coefficient-wise == operator of *this and \a other - * - * \warning this performs an exact comparison, which is generally a bad idea with floating-point types. - * In order to check for equality between two vectors or matrices with floating-point coefficients, it is - * generally a far better idea to use a fuzzy comparison as provided by MatrixBase::isApprox() and - * MatrixBase::isMuchSmallerThan(). - * - * Example: \include Cwise_equal_equal.cpp - * Output: \verbinclude Cwise_equal_equal.out - * - * \sa MatrixBase::all(), MatrixBase::any(), MatrixBase::isApprox(), MatrixBase::isMuchSmallerThan() - */ +/** \deprecated ArrayBase::operator==() */ template template inline const EIGEN_CWISE_BINOP_RETURN_TYPE(std::equal_to) @@ -388,20 +227,7 @@ Cwise::operator==(const MatrixBase &other) const return EIGEN_CWISE_BINOP_RETURN_TYPE(std::equal_to)(_expression(), other.derived()); } -/** \array_module - * - * \returns an expression of the coefficient-wise != operator of *this and \a other - * - * \warning this performs an exact comparison, which is generally a bad idea with floating-point types. - * In order to check for equality between two vectors or matrices with floating-point coefficients, it is - * generally a far better idea to use a fuzzy comparison as provided by MatrixBase::isApprox() and - * MatrixBase::isMuchSmallerThan(). - * - * Example: \include Cwise_not_equal.cpp - * Output: \verbinclude Cwise_not_equal.out - * - * \sa MatrixBase::all(), MatrixBase::any(), MatrixBase::isApprox(), MatrixBase::isMuchSmallerThan() - */ +/** \deprecated ArrayBase::operator!=() */ template template inline const EIGEN_CWISE_BINOP_RETURN_TYPE(std::not_equal_to) @@ -412,12 +238,7 @@ Cwise::operator!=(const MatrixBase &other) const // comparisons to scalar value -/** \array_module - * - * \returns an expression of the coefficient-wise \< operator of *this and a scalar \a s - * - * \sa operator<(const MatrixBase &) const - */ +/** \deprecated ArrayBase::operator<(Scalar) */ template inline const EIGEN_CWISE_COMP_TO_SCALAR_RETURN_TYPE(std::less) Cwise::operator<(Scalar s) const @@ -426,12 +247,7 @@ Cwise::operator<(Scalar s) const typename ExpressionType::ConstantReturnType(_expression().rows(), _expression().cols(), s)); } -/** \array_module - * - * \returns an expression of the coefficient-wise \<= operator of *this and a scalar \a s - * - * \sa operator<=(const MatrixBase &) const - */ +/** \deprecated ArrayBase::operator<=(Scalar) */ template inline const EIGEN_CWISE_COMP_TO_SCALAR_RETURN_TYPE(std::less_equal) Cwise::operator<=(Scalar s) const @@ -440,12 +256,7 @@ Cwise::operator<=(Scalar s) const typename ExpressionType::ConstantReturnType(_expression().rows(), _expression().cols(), s)); } -/** \array_module - * - * \returns an expression of the coefficient-wise \> operator of *this and a scalar \a s - * - * \sa operator>(const MatrixBase &) const - */ +/** \deprecated ArrayBase::operator>(Scalar) */ template inline const EIGEN_CWISE_COMP_TO_SCALAR_RETURN_TYPE(std::greater) Cwise::operator>(Scalar s) const @@ -454,12 +265,7 @@ Cwise::operator>(Scalar s) const typename ExpressionType::ConstantReturnType(_expression().rows(), _expression().cols(), s)); } -/** \array_module - * - * \returns an expression of the coefficient-wise \>= operator of *this and a scalar \a s - * - * \sa operator>=(const MatrixBase &) const - */ +/** \deprecated ArrayBase::operator>=(Scalar) */ template inline const EIGEN_CWISE_COMP_TO_SCALAR_RETURN_TYPE(std::greater_equal) Cwise::operator>=(Scalar s) const @@ -468,17 +274,7 @@ Cwise::operator>=(Scalar s) const typename ExpressionType::ConstantReturnType(_expression().rows(), _expression().cols(), s)); } -/** \array_module - * - * \returns an expression of the coefficient-wise == operator of *this and a scalar \a s - * - * \warning this performs an exact comparison, which is generally a bad idea with floating-point types. - * In order to check for equality between two vectors or matrices with floating-point coefficients, it is - * generally a far better idea to use a fuzzy comparison as provided by MatrixBase::isApprox() and - * MatrixBase::isMuchSmallerThan(). - * - * \sa operator==(const MatrixBase &) const - */ +/** \deprecated ArrayBase::operator==(Scalar) */ template inline const EIGEN_CWISE_COMP_TO_SCALAR_RETURN_TYPE(std::equal_to) Cwise::operator==(Scalar s) const @@ -487,17 +283,7 @@ Cwise::operator==(Scalar s) const typename ExpressionType::ConstantReturnType(_expression().rows(), _expression().cols(), s)); } -/** \array_module - * - * \returns an expression of the coefficient-wise != operator of *this and a scalar \a s - * - * \warning this performs an exact comparison, which is generally a bad idea with floating-point types. - * In order to check for equality between two vectors or matrices with floating-point coefficients, it is - * generally a far better idea to use a fuzzy comparison as provided by MatrixBase::isApprox() and - * MatrixBase::isMuchSmallerThan(). - * - * \sa operator!=(const MatrixBase &) const - */ +/** \deprecated ArrayBase::operator!=(Scalar) */ template inline const EIGEN_CWISE_COMP_TO_SCALAR_RETURN_TYPE(std::not_equal_to) Cwise::operator!=(Scalar s) const @@ -508,15 +294,7 @@ Cwise::operator!=(Scalar s) const // scalar addition -/** \array_module - * - * \returns an expression of \c *this with each coeff incremented by the constant \a scalar - * - * Example: \include Cwise_plus.cpp - * Output: \verbinclude Cwise_plus.out - * - * \sa operator+=(), operator-() - */ +/** \deprecated ArrayBase::operator+(Scalar) */ template inline const typename Cwise::ScalarAddReturnType Cwise::operator+(const Scalar& scalar) const @@ -524,30 +302,14 @@ Cwise::operator+(const Scalar& scalar) const return typename Cwise::ScalarAddReturnType(m_matrix, ei_scalar_add_op(scalar)); } -/** \array_module - * - * Adds the given \a scalar to each coeff of this expression. - * - * Example: \include Cwise_plus_equal.cpp - * Output: \verbinclude Cwise_plus_equal.out - * - * \sa operator+(), operator-=() - */ +/** \deprecated ArrayBase::operator+=(Scalar) */ template inline ExpressionType& Cwise::operator+=(const Scalar& scalar) { return m_matrix.const_cast_derived() = *this + scalar; } -/** \array_module - * - * \returns an expression of \c *this with each coeff decremented by the constant \a scalar - * - * Example: \include Cwise_minus.cpp - * Output: \verbinclude Cwise_minus.out - * - * \sa operator+(), operator-=() - */ +/** \deprecated ArrayBase::operator-(Scalar) */ template inline const typename Cwise::ScalarAddReturnType Cwise::operator-(const Scalar& scalar) const @@ -555,16 +317,7 @@ Cwise::operator-(const Scalar& scalar) const return *this + (-scalar); } -/** \array_module - * - * Substracts the given \a scalar from each coeff of this expression. - * - * Example: \include Cwise_minus_equal.cpp - * Output: \verbinclude Cwise_minus_equal.out - * - * \sa operator+=(), operator-() - */ - +/** \deprecated ArrayBase::operator-=(Scalar) */ template inline ExpressionType& Cwise::operator-=(const Scalar& scalar) { diff --git a/Eigen/src/Geometry/Quaternion.h b/Eigen/src/Geometry/Quaternion.h index 24772089e..b8b41ebcf 100644 --- a/Eigen/src/Geometry/Quaternion.h +++ b/Eigen/src/Geometry/Quaternion.h @@ -136,26 +136,26 @@ public: template Scalar angularDistance(const QuaternionBase& other) const; - /** \returns an equivalent 3x3 rotation matrix */ + /** \returns an equivalent 3x3 rotation matrix */ Matrix3 toRotationMatrix() const; - /** \returns the quaternion which transform \a a into \a b through a rotation */ + /** \returns the quaternion which transform \a a into \a b through a rotation */ template Derived& setFromTwoVectors(const MatrixBase& a, const MatrixBase& b); template EIGEN_STRONG_INLINE Quaternion operator* (const QuaternionBase& q) const; template EIGEN_STRONG_INLINE Derived& operator*= (const QuaternionBase& q); - /** \returns the quaternion describing the inverse rotation */ + /** \returns the quaternion describing the inverse rotation */ Quaternion inverse() const; - /** \returns the conjugated quaternion */ + /** \returns the conjugated quaternion */ Quaternion conjugate() const; - /** \returns an interpolation for a constant motion between \a other and \c *this - * \a t in [0;1] - * see http://en.wikipedia.org/wiki/Slerp - */ + /** \returns an interpolation for a constant motion between \a other and \c *this + * \a t in [0;1] + * see http://en.wikipedia.org/wiki/Slerp + */ template Quaternion slerp(Scalar t, const QuaternionBase& other) const; /** \returns \c true if \c *this is approximately equal to \a other, within the precision diff --git a/Eigen/src/plugins/ArrayCwiseBinaryOps.h b/Eigen/src/plugins/ArrayCwiseBinaryOps.h index bb0e4dd22..dc17284e2 100644 --- a/Eigen/src/plugins/ArrayCwiseBinaryOps.h +++ b/Eigen/src/plugins/ArrayCwiseBinaryOps.h @@ -21,6 +21,24 @@ operator/(const EIGEN_CURRENT_STORAGE_BASE_CLASS &other) const return CwiseBinaryOp, Derived, OtherDerived>(derived(), other.derived()); } +/** \returns an expression of the coefficient-wise min of \c *this and \a other + * + * Example: \include Cwise_min.cpp + * Output: \verbinclude Cwise_min.out + * + * \sa max() + */ +EIGEN_MAKE_CWISE_BINARY_OP(min,ei_scalar_min_op) + +/** \returns an expression of the coefficient-wise max of \c *this and \a other + * + * Example: \include Cwise_max.cpp + * Output: \verbinclude Cwise_max.out + * + * \sa min() + */ +EIGEN_MAKE_CWISE_BINARY_OP(max,ei_scalar_max_op) + /** \returns an expression of the coefficient-wise \< operator of *this and \a other * * Example: \include Cwise_less.cpp diff --git a/doc/C01_QuickStartGuide.dox b/doc/C01_QuickStartGuide.dox index b342732ec..b346e2371 100644 --- a/doc/C01_QuickStartGuide.dox +++ b/doc/C01_QuickStartGuide.dox @@ -78,13 +78,20 @@ This slows compilation down but at least you don't have to worry anymore about i top \section TutorialCoreMatrixTypes Array, matrix and vector types -Eigen provides two kinds of dense objects: mathematical matrices and vectors which are both represented by the template class Matrix, and 1D and 2D arrays represented by the template class Array. While the former (Matrix) is specialized for the representation of mathematical objects, the latter (Array) represents a collection of scalar values arranged in a 1D or 2D fashion. In particular, all operations performed on arrays are coefficient wise. Conversion between the two worlds can be done using the MatrixBase::array() and ArrayBase::matrix() functions respectively without any overhead. See \ref TutorialCoreArithmeticOperators for further details. +Eigen provides two kinds of dense objects: mathematical matrices and vectors which are both represented by the template class Matrix, and 1D and 2D arrays represented by the template class Array. While the former (Matrix) is specialized for the representation of mathematical objects, the latter (Array) represents a collection of scalar values arranged in a 1D or 2D fashion. As a major difference, all operations performed on arrays are coefficient wise. Matrix and Array have a lot of similarities since they both inherits the DenseBase and DenseStorageBase classes. In the rest of this tutorial we will use the following symbols to emphasize the features which are specifics to a given kind of object: +\li \matrixworld for matrix/vector only features +\li \arrayworld for array only features -In most cases, you can simply use one of the \ref matrixtypedefs "convenience typedefs". +Note that conversion between the two worlds can be done using the MatrixBase::array() and ArrayBase::matrix() functions respectively without any overhead. -The template class Matrix, just like the class Array) take a number of template parameters, but for now it is enough to understand the 3 first ones (and the others can then be left unspecified): +In most cases, you can simply use one of the convenience typedefs for \ref matrixtypedefs "matrices" and \ref arraytypedefs "arrays". -\code Matrix \endcode +The template class Matrix (just like the class Array) take a number of template parameters, but for now it is enough to understand the 3 first ones (and the others can then be left unspecified): + +\code +Matrix +Array +\endcode \li \c Scalar is the scalar type, i.e. the type of the coefficients. That is, if you want a vector of floats, choose \c float here. \li \c RowsAtCompileTime and \c ColsAtCompileTime are the number of rows and columns of the matrix as known at compile-time. @@ -96,23 +103,36 @@ For dynamic-size, that is in order to left the number of rows or of columns unsp All combinations are allowed: you can have a matrix with a fixed number of rows and a dynamic number of columns, etc. The following are all valid: \code -Matrix // Dynamic number of columns -Matrix // Dynamic number of rows -Matrix // Fully dynamic -Matrix // Fully fixed +Matrix // Dynamic number of columns +Matrix // Dynamic number of rows +Matrix // Fully dynamic +Matrix // Fully fixed \endcode Fixed-size and partially-dynamic-size matrices may use all the same API calls as fully dynamic matrices, but the fixed dimension(s) must remain constant, or an assertion failure will occur. -top\section TutorialCoreCoefficients Coefficient access +Finally, note that the default typedefs for array containers is slighlty different as we have to distinghish between 1D and 2D arrays: +\code +ArrayXf // 1D dynamic array of floats +Array2i // 1D array of integers of size 2 +ArrayXXd // 2D fully dynamic array of doubles +Array44f // 2D array of floats of size 4x4 +\endcode -Eigen supports the following syntaxes for read and write coefficient access: + +top +\section TutorialCoreCoefficients Coefficient access + +Eigen supports the following syntaxes for read and write coefficient access of matrices, vectors and arrays: \code matrix(i,j); vector(i) vector[i] +\endcode +Vectors support also the following additional read-write accessors: +\code vector.x() // first coefficient vector.y() // second coefficient vector.z() // third coefficient @@ -121,9 +141,11 @@ vector.w() // fourth coefficient Notice that these coefficient access methods have assertions checking the ranges. So if you do a lot of coefficient access, these assertion can have an important cost. There are then two possibilities if you want avoid paying this cost: \li Either you can disable assertions altogether, by defining EIGEN_NO_DEBUG or NDEBUG. Notice that some IDEs like MS Visual Studio define NDEBUG automatically in "Release Mode". -\li Or you can disable the checks on a case-by-case basis by using the coeff() and coeffRef() methods: see MatrixBase::coeff(int,int) const, MatrixBase::coeffRef(int,int), etc. +\li Or you can disable the checks on a case-by-case basis by using the coeff() and coeffRef() methods: see DenseBase::coeff(int,int) const, DenseBase::coeffRef(int,int), etc. -top\section TutorialCoreMatrixInitialization Matrix and vector creation and initialization + +top +\section TutorialCoreMatrixInitialization Matrix and vector creation and initialization \subsection TutorialCtors Matrix constructors @@ -169,7 +191,8 @@ Vector4f w(1.2f, 3.4f, 5.6f, 7.8f); \endcode \subsection TutorialPredefMat Predefined Matrices -Eigen offers several static methods to create special matrix expressions, and non-static methods to assign these expressions to existing matrices: +Eigen offers several static methods to create special matrix expressions, and non-static methods to assign these expressions to existing matrices. +The following are @@ -180,13 +203,14 @@ Eigen offers several static methods to create special matrix expressions, and no - + + + + + + + + + \link DenseBase::block(int,int,int,int) (more) \endlink + \link DenseBase::block(int,int) (more) \endlink + \link DenseBase::corner(CornerType,int,int) (more) \endlink + \link DenseBase::corner(CornerType) (more) \endlink +\link DenseBase::minor() minor \endlink (read-write)
\code -Matrix3f x; +typedef {Matrix3f|Array33f} FixedXD; +FixedXD x; -x = Matrix3f::Zero(); -x = Matrix3f::Ones(); -x = Matrix3f::Constant(value); -x = Matrix3f::Identity(); -x = Matrix3f::Random(); +x = FixedXD::Zero(); +x = FixedXD::Ones(); +x = FixedXD::Constant(value); +x = FixedXD::Identity(); +x = FixedXD::Random(); x.setZero(); x.setOnes(); @@ -197,13 +221,14 @@ x.setRandom(); \code -MatrixXf x; +typedef {MatrixXf|ArrayXXf} Dynamic2D; +Dynamic2D x; -x = MatrixXf::Zero(rows, cols); -x = MatrixXf::Ones(rows, cols); -x = MatrixXf::Constant(rows, cols, value); -x = MatrixXf::Identity(rows, cols); -x = MatrixXf::Random(rows, cols); +x = Dynamic2D::Zero(rows, cols); +x = Dynamic2D::Ones(rows, cols); +x = Dynamic2D::Constant(rows, cols, value); +x = Dynamic2D::Identity(rows, cols); +x = Dynamic2D::Random(rows, cols); x.setZero(rows, cols); x.setOnes(rows, cols); @@ -214,13 +239,14 @@ x.setRandom(rows, cols); \code -VectorXf x; +typedef {VectorXf|ArrayXf} Dynamic1D; +Dynamic1D x; -x = VectorXf::Zero(size); -x = VectorXf::Ones(size); -x = VectorXf::Constant(size, value); -x = VectorXf::Identity(size); -x = VectorXf::Random(size); +x = Dynamic1D::Zero(size); +x = Dynamic1D::Ones(size); +x = Dynamic1D::Constant(size, value); +x = Dynamic1D::Identity(size); +x = Dynamic1D::Random(size); x.setZero(size); x.setOnes(size); @@ -231,7 +257,25 @@ x.setRandom(size);
\redstar the Random() and setRandom() functions require the inclusion of the Array module (\c \#include \c )
Basis vectors \link MatrixBase::Unit [details]\endlink
The following are for matrix only: \matrixworld
+\code +x = FixedXD::Identity(); +x.setIdentity(); +\endcode + +\code +x = Dynamic2D::Identity(rows, cols); +x.setIdentity(rows, cols); +\endcode + +
Basis vectors \matrixworld \link MatrixBase::Unit [details]\endlink
\code Vector3f::UnitX() // 1 0 0 Vector3f::UnitY() // 0 1 0 @@ -265,7 +309,7 @@ v = 6 6 6 \subsection TutorialCasting Casting -In Eigen, any matrices of same size and same scalar type are all naturally compatible. The scalar type can be explicitly casted to another one using the template MatrixBase::cast() function: +In Eigen, any matrices of same size and same scalar type are all naturally compatible. The scalar type can be explicitly casted to another one using the template DenseBase::cast() function: \code Matrix3d md(1,2,3); Matrix3f mf = md.cast(); @@ -280,6 +324,28 @@ res = a+b; // OK: res is resized to size 3x3 \endcode Of course, fixed-size matrices can't be resized. +An array object or expression can be directly assigned to a matrix, and vice versa: +\code +Matrix4f res; +Array44f a, b; +res = a * b; +\endcode +On the other hand, an array and a matrix expressions cannot be mixed in an expression, and one have to be converted to the other using the MatrixBase::array() \matrixworld and ArrayBase::matrix() \arrayworld functions respectively: +\code +Matrix4f m1, m2; +Array44f a1, a2; +m2 = a1 * m1.array(); // coeffwise product +a2 = a1.matrix() * m1; // matrix product +\endcode +Finally it is possible to declare a variable wrapping a matrix as an array object and vice versa: +\code +MatrixXf m1; +ArrayWrapper a1(m1); // a1 and m1 share the same coefficients +// now you can use a1 as an alias for m1.array() +ArrayXXf a2; +MatrixWrapper m2(a1); // a2 and m2 share the same coefficients +// ... +\endcode \subsection TutorialMap Map Any memory buffer can be mapped as an Eigen expression using the Map() static method: @@ -289,21 +355,21 @@ VectorXf::Map(&stlarray[0], stlarray.size()).squaredNorm(); \endcode Here VectorXf::Map returns an object of class Map, which behaves like a VectorXf except that it uses the existing array. You can write to this object, that will write to the existing array. You can also construct a named obtect to reuse it: \code -float array[rows*cols]; -Map m(array,rows,cols); +float data[rows*cols]; +Map m(data,rows,cols); m = othermatrix1 * othermatrix2; m.eigenvalues(); \endcode In the fixed-size case, no need to pass sizes: \code -float array[9]; -Map m(array); -Matrix3d::Map(array).setIdentity(); +float data[9]; +Map m(data); +Matrix3d::Map(data).setIdentity(); \endcode \subsection TutorialCommaInit Comma initializer -Eigen also offers a \ref MatrixBaseCommaInitRef "comma initializer syntax" which allows you to set all the coefficients of a matrix to specific values: +Eigen also offers a \ref MatrixBaseCommaInitRef "comma initializer syntax" which allows you to set all the coefficients of any dense objects (matrix, vector, array, block, etc.) to specific values: @@ -328,12 +394,12 @@ Eigen's comma initializer usually compiles to very optimized code without any ov +top +\section TutorialCoreArithmeticOperators Arithmetic Operators +In short, all arithmetic operators can be used right away as in the following example. Note however that for matrices and vectors arithmetic operators are only given their usual meaning from mathematics tradition while all array operators are performed coefficient wise. - -top\section TutorialCoreArithmeticOperators Arithmetic Operators - -In short, all arithmetic operators can be used right away as in the following example. Note however that arithmetic operators are only given their usual meaning from mathematics tradition. For other operations, such as taking the coefficient-wise product of two vectors, see the discussion of \link Cwise .cwise() \endlink below. Anyway, here is an example demonstrating basic arithmetic operators: +Here is an example demonstrating basic arithmetic operators: \code mat4 -= mat1*1.5 + mat2 * (mat3/4); \endcode @@ -342,7 +408,7 @@ a matrix addition ("+") and subtraction with assignment ("-=").
\include Tutorial_commainit_01.cpp
+
-matrix/vector product\code +matrix/vector product \matrixworld\code col2 = mat1 * col1; row2 = row1 * mat1; row1 *= mat1; mat3 = mat1 * mat2; mat3 *= mat1; \endcode @@ -357,107 +423,108 @@ scalar product\code mat3 = mat1 * s1; mat3 = s1 * mat1; mat3 *= s1; mat3 = mat1 / s1; mat3 /= s1;\endcode
+Other coefficient wise operators\code +mat1.cwiseProduct(mat2); mat1.cwiseQuotient(mat2); +mat1.cwiseMin(mat2); mat1.cwiseMax(mat2); +mat1.cwiseAbs2(); mat1.cwiseSqrt(); +mat1.cwiseAbs();\endcode +
-In Eigen, only traditional mathematical operators can be used right away. -But don't worry, thanks to the \link Cwise .cwise() \endlink operator prefix, -Eigen's matrices are also very powerful as a numerical container supporting -most common coefficient-wise operators. +In addition to the above operators, array objects supports all kind of coefficient wise operators which usually apply to scalar values. Recall that those operators can be used on matrices by converting them to arrays using the array() function (see \ref TutorialCasting Casting).
- - +
Coefficient wise \link Cwise::operator*() product \endlink\code mat3 = mat1.cwise() * mat2; \endcode +
Coefficient wise \link ArrayBase::operator*() product \arrayworld \endlink\code array3 = array1 * array2; \endcode
-Add a scalar to all coefficients \redstar\code -mat3 = mat1.cwise() + scalar; -mat3.array() += scalar; -mat3.array() -= scalar; +Add a scalar to all coefficients\code +array3 = array1 + scalar; +array3 += scalar; +array3 -= scalar; \endcode
-Coefficient wise \link Cwise::operator/() division \endlink \redstar\code -mat3 = mat1.array() / mat2.array(); \endcode +Coefficient wise \link ArrayBase::operator/() division \endlink \arrayworld\code +array3 = array1 / array2; \endcode
-Coefficient wise \link Cwise::inverse() reciprocal \endlink \redstar\code -mat3 = mat1.array().inverse(); \endcode +Coefficient wise \link ArrayBase::inverse() reciprocal \endlink \arrayworld\code +array3 = array1.inverse(); \endcode
-Coefficient wise comparisons \redstar \n +Coefficient wise comparisons \arrayworld \n (support all operators)\code -mat3 = mat1.array() < mat2.array(); -mat3 = mat1.array() <= mat2.array(); -mat3 = mat1.array() > mat2.array(); +array3 = array1 < array2; +array3 = array1 <= array2; +array3 = array1 > array2; etc. \endcode
-\b Trigo \redstar: \n -\link Cwise::sin sin \endlink, \link Cwise::cos cos \endlink\code -mat3 = mat1.array().sin(); +\b Trigo \arrayworld: \n +\link ArrayBase::sin sin \endlink, \link ArrayBase::cos cos \endlink\code +array3 = array1.sin(); etc. \endcode
-\b Power \redstar: \n \link Cwise::pow() pow \endlink, +\b Power \arrayworld: \n \link ArrayBase::pow() pow \endlink, \link ArrayBase::square square \endlink, \link ArrayBase::cube cube \endlink, \n \link ArrayBase::sqrt sqrt \endlink, \link ArrayBase::exp exp \endlink, \link ArrayBase::log log \endlink \code -mat3 = mat1.array().square(); -mat3 = mat1.array().pow(5); -mat3 = mat1.array().log(); +array3 = array1.square(); +array3 = array1.pow(5); +array3 = array1.log(); etc. \endcode
-\link Cwise::min min \endlink, \link Cwise::max max \endlink, \n -absolute value (\link Cwise::abs() abs \endlink, \link Cwise::abs2() abs2 \endlink) +\link ArrayBase::min min \endlink, \link ArrayBase::max max \endlink, \n +absolute value (\link ArrayBase::abs() abs \endlink, \link ArrayBase::abs2() abs2 \endlink \arrayworld) \code -mat3 = mat1.cwiseMin(mat2); -mat3 = mat1.cwiseMax(mat2); -mat3 = mat1.cwiseAbs(); -mat3 = mat1.cwiseAbs2(); +array3 = array1.min(array2); +array3 = array1.max(array2); +array3 = array1.abs(); +array3 = array1.abs2(); \endcode
-\redstar Those functions require the inclusion of the Array module (\c \#include \c ). -\b Side \b note: If you think that the \c .cwise() syntax is too verbose for your own taste and prefer to have non-conventional mathematical operators directly available, then feel free to extend MatrixBase as described \ref ExtendingMatrixBase "here". - -So far, we saw the notation \code mat1*mat2 \endcode for matrix product, and \code mat1.cwise()*mat2 \endcode for coefficient-wise product. What about other kinds of products, which in some other libraries also use arithmetic operators? In Eigen, they are accessed as follows -- note that here we are anticipating on further sections, for convenience. +So far, we saw the notation \code mat1*mat2 \endcode for matrix product, and \code array1*array2 \endcode for coefficient-wise product. What about other kinds of products, which in some other libraries also use arithmetic operators? In Eigen, they are accessed as follows -- note that here we are anticipating on further sections, for convenience. -
\link MatrixBase::dot() dot product \endlink (inner product)\code +
\link MatrixBase::dot() dot product \endlink (inner product) \matrixworld\code scalar = vec1.dot(vec2);\endcode
-outer product\code +outer product \matrixworld\code mat = vec1 * vec2.transpose();\endcode
-\link MatrixBase::cross() cross product \endlink\code +\link MatrixBase::cross() cross product \endlink \matrixworld\code #include vec3 = vec1.cross(vec2);\endcode
- -top\section TutorialCoreReductions Reductions +top +\section TutorialCoreReductions Reductions Eigen provides several reduction methods such as: -\link MatrixBase::minCoeff() minCoeff() \endlink, \link MatrixBase::maxCoeff() maxCoeff() \endlink, -\link MatrixBase::sum() sum() \endlink, \link MatrixBase::trace() trace() \endlink, -\link MatrixBase::norm() norm() \endlink, \link MatrixBase::squaredNorm() squaredNorm() \endlink, -\link MatrixBase::all() all() \endlink \redstar,and \link MatrixBase::any() any() \endlink \redstar. +\link DenseBase::minCoeff() minCoeff() \endlink, \link DenseBase::maxCoeff() maxCoeff() \endlink, +\link DenseBase::sum() sum() \endlink, \link MatrixBase::trace() trace() \endlink \matrixworld, +\link MatrixBase::norm() norm() \endlink \matrixworld, \link MatrixBase::squaredNorm() squaredNorm() \endlink \matrixworld, +\link DenseBase::all() all() \endlink \redstar,and \link DenseBase::any() any() \endlink \redstar. All reduction operations can be done matrix-wise, -\link MatrixBase::colwise() column-wise \endlink \redstar or -\link MatrixBase::rowwise() row-wise \endlink \redstar. Usage example: +\link DenseBase::colwise() column-wise \endlink \redstar or +\link DenseBase::rowwise() row-wise \endlink \redstar. Usage example:
\code 5 3 1 @@ -472,7 +539,7 @@ mat = 2 7 8 \endcode
-Also note that maxCoeff and minCoeff can takes optional arguments returning the coordinates of the respective min/max coeff: \link MatrixBase::maxCoeff(int*,int*) const maxCoeff(int* i, int* j) \endlink, \link MatrixBase::minCoeff(int*,int*) const minCoeff(int* i, int* j) \endlink. +Also note that maxCoeff and minCoeff can takes optional arguments returning the coordinates of the respective min/max coeff: \link DenseBase::maxCoeff(int*,int*) const maxCoeff(int* i, int* j) \endlink, \link DenseBase::minCoeff(int*,int*) const minCoeff(int* i, int* j) \endlink. \b Side \b note: The all() and any() functions are especially useful in combination with coeff-wise comparison operators. @@ -482,8 +549,8 @@ Also note that maxCoeff and minCoeff can takes optional arguments returning the top\section TutorialCoreMatrixBlocks Matrix blocks -Read-write access to a \link MatrixBase::col(int) column \endlink -or a \link MatrixBase::row(int) row \endlink of a matrix: +Read-write access to a \link DenseBase::col(int) column \endlink +or a \link DenseBase::row(int) row \endlink of a matrix (or array): \code mat1.row(i) = mat2.col(j); mat1.col(j1).swap(mat1.col(j2)); @@ -505,34 +572,34 @@ Read-write access to sub-vectors: Read-write access to sub-matrices:
\code mat1.block(i,j,rows,cols)\endcode - \link MatrixBase::block(int,int,int,int) (more) \endlink \code mat1.block(i,j)\endcode - \link MatrixBase::block(int,int) (more) \endlink the \c rows x \c cols sub-matrix \n starting from position (\c i,\c j)
\code mat1.corner(TopLeft,rows,cols) mat1.corner(TopRight,rows,cols) mat1.corner(BottomLeft,rows,cols) mat1.corner(BottomRight,rows,cols)\endcode - \link MatrixBase::corner(CornerType,int,int) (more) \endlink \code mat1.corner(TopLeft) mat1.corner(TopRight) mat1.corner(BottomLeft) mat1.corner(BottomRight)\endcode - \link MatrixBase::corner(CornerType) (more) \endlink the \c rows x \c cols sub-matrix \n taken in one of the four corners
\code mat4x4.minor(i,j) = mat3x3; mat3x3 = mat4x4.minor(i,j);\endcode -\link MatrixBase::minor() minor \endlink (read-write)
-top\section TutorialCoreDiagonalMatrices Diagonal matrices +top\section TutorialCoreDiagonalMatrices Diagonal matrices \matrixworld
@@ -549,23 +616,29 @@ mat3 = mat1 * vec2.asDiagonal();\endcode
-top\section TutorialCoreTransposeAdjoint Transpose and Adjoint operations + + +top +\section TutorialCoreTransposeAdjoint Transpose and Adjoint operations
-\link MatrixBase::transpose() transposition \endlink (read-write)\code +\link DenseBase::transpose() transposition \endlink (read-write)\code mat3 = mat1.transpose() * mat2; mat3.transpose() = mat1 * mat2.transpose(); \endcode
-\link MatrixBase::adjoint() adjoint \endlink (read only)\n\code +\link MatrixBase::adjoint() adjoint \endlink (read only) \matrixworld\n\code mat3 = mat1.adjoint() * mat2; \endcode
-top\section TutorialCoreDotNorm Dot-product, vector norm, normalization + + +top +\section TutorialCoreDotNorm Dot-product, vector norm, normalization \matrixworld
@@ -586,7 +659,10 @@ vec1.normalize();\endcode
-top\section TutorialCoreTriangularMatrix Dealing with triangular matrices + + +top +\section TutorialCoreTriangularMatrix Dealing with triangular matrices \matrixworld Currently, Eigen does not provide any explicit triangular matrix, with storage class. Instead, we can reference a triangular part of a square matrix or expression to perform special treatment on it. @@ -629,7 +705,9 @@ m1.adjoint().triangularView().solveInPlace(m2)\endcode -top\section TutorialCoreSelfadjointMatrix Dealing with symmetric/selfadjoint matrices + +top +\section TutorialCoreSelfadjointMatrix Dealing with symmetric/selfadjoint matrices \matrixworld Just as for triangular matrix, you can reference any triangular part of a square matrix to see it a selfadjoint matrix to perform special and optimized operations. Again the opposite triangular is never referenced and can be @@ -673,7 +751,8 @@ m1.selfadjointView().ldlt().solveInPlace(m2); -top\section TutorialCoreSpecialTopics Special Topics +top +\section TutorialCoreSpecialTopics Special Topics \ref TopicLazyEvaluation "Lazy Evaluation and Aliasing": Thanks to expression templates, Eigen is able to apply lazy evaluation wherever that is beneficial. diff --git a/doc/Doxyfile.in b/doc/Doxyfile.in index 64bb7bd01..b8957250d 100644 --- a/doc/Doxyfile.in +++ b/doc/Doxyfile.in @@ -208,6 +208,8 @@ ALIASES = "only_for_vectors=This is only for vectors (either row- "svd_module=This is defined in the %SVD module. \code #include \endcode" \ "label=\bug" \ "redstar=*" \ + "matrixworld=*" \ + "arrayworld=*" \ "nonstableyet=\warning This is not considered to be part of the stable public API yet. Changes may happen in future releases. See \ref Experimental \"Experimental parts of Eigen\"" \ "note_about_arbitrary_choice_of_solution=If there exists more than one solution, this method will arbitrarily choose one." \ "note_about_using_kernel_to_study_multiple_solutions=If you need a complete analysis of the space of solutions, take the one solution obtained by this method and add to it elements of the kernel, as determined by kernel()." \ diff --git a/doc/snippets/Cwise_abs.cpp b/doc/snippets/Cwise_abs.cpp index 7f94dbdc2..0aeec3a40 100644 --- a/doc/snippets/Cwise_abs.cpp +++ b/doc/snippets/Cwise_abs.cpp @@ -1,2 +1,2 @@ -Vector3d v(1,-2,-3); -cout << v.cwise().abs() << endl; +Array3d v(1,-2,-3); +cout << v.abs() << endl; diff --git a/doc/snippets/Cwise_abs2.cpp b/doc/snippets/Cwise_abs2.cpp index 3b338af3e..2c4f9b344 100644 --- a/doc/snippets/Cwise_abs2.cpp +++ b/doc/snippets/Cwise_abs2.cpp @@ -1,2 +1,2 @@ -Vector3d v(1,-2,-3); -cout << v.cwise().abs2() << endl; +Array3d v(1,-2,-3); +cout << v.abs2() << endl; diff --git a/doc/snippets/Cwise_cos.cpp b/doc/snippets/Cwise_cos.cpp index 4eb230a80..f589f0724 100644 --- a/doc/snippets/Cwise_cos.cpp +++ b/doc/snippets/Cwise_cos.cpp @@ -1,2 +1,2 @@ -Vector3d v(M_PI, M_PI/2, M_PI/3); -cout << v.cwise().cos() << endl; +Array3d v(M_PI, M_PI/2, M_PI/3); +cout << v.cos() << endl; diff --git a/doc/snippets/Cwise_cube.cpp b/doc/snippets/Cwise_cube.cpp index fd404365a..85e41dcce 100644 --- a/doc/snippets/Cwise_cube.cpp +++ b/doc/snippets/Cwise_cube.cpp @@ -1,2 +1,2 @@ -Vector3d v(2,3,4); -cout << v.cwise().cube() << endl; +Array3d v(2,3,4); +cout << v.cube() << endl; diff --git a/doc/snippets/Cwise_equal_equal.cpp b/doc/snippets/Cwise_equal_equal.cpp index ea49120f0..0ba96f660 100644 --- a/doc/snippets/Cwise_equal_equal.cpp +++ b/doc/snippets/Cwise_equal_equal.cpp @@ -1,2 +1,2 @@ -Vector3d v(1,2,3), w(3,2,1); -cout << (v.cwise()==w) << endl; +Array3d v(1,2,3), w(3,2,1); +cout << (v==w) << endl; diff --git a/doc/snippets/Cwise_exp.cpp b/doc/snippets/Cwise_exp.cpp index c22fd4cb7..db23618f9 100644 --- a/doc/snippets/Cwise_exp.cpp +++ b/doc/snippets/Cwise_exp.cpp @@ -1,2 +1,2 @@ -Vector3d v(1,2,3); -cout << v.cwise().exp() << endl; +Array3d v(1,2,3); +cout << v.exp() << endl; diff --git a/doc/snippets/Cwise_greater.cpp b/doc/snippets/Cwise_greater.cpp index 54641e1ce..40ad0296d 100644 --- a/doc/snippets/Cwise_greater.cpp +++ b/doc/snippets/Cwise_greater.cpp @@ -1,2 +1,2 @@ -Vector3d v(1,2,3), w(3,2,1); -cout << (v.cwise()>w) << endl; +Array3d v(1,2,3), w(3,2,1); +cout << (v>w) << endl; diff --git a/doc/snippets/Cwise_greater_equal.cpp b/doc/snippets/Cwise_greater_equal.cpp index 744a9ae1e..6a08f8948 100644 --- a/doc/snippets/Cwise_greater_equal.cpp +++ b/doc/snippets/Cwise_greater_equal.cpp @@ -1,2 +1,2 @@ -Vector3d v(1,2,3), w(3,2,1); -cout << (v.cwise()>=w) << endl; +Array3d v(1,2,3), w(3,2,1); +cout << (v>=w) << endl; diff --git a/doc/snippets/Cwise_inverse.cpp b/doc/snippets/Cwise_inverse.cpp index 7f6eee85d..3967a7ecf 100644 --- a/doc/snippets/Cwise_inverse.cpp +++ b/doc/snippets/Cwise_inverse.cpp @@ -1,2 +1,2 @@ -Vector3d v(2,3,4); -cout << v.cwise().inverse() << endl; +Array3d v(2,3,4); +cout << v.inverse() << endl; diff --git a/doc/snippets/Cwise_less.cpp b/doc/snippets/Cwise_less.cpp index 008e3d9e6..cafd3b6e0 100644 --- a/doc/snippets/Cwise_less.cpp +++ b/doc/snippets/Cwise_less.cpp @@ -1,2 +1,2 @@ -Vector3d v(1,2,3), w(3,2,1); -cout << (v.cwise()p0).all()) << endl; + << ((boxMin.array()p0.array()).all()) << endl; cout << "Is (" << p1.transpose() << ") inside the box: " - << ((boxMin.cwise()p1).all()) << endl; + << ((boxMin.array()p1.array()).all()) << endl; diff --git a/doc/snippets/MatrixBase_cwise.cpp b/doc/snippets/MatrixBase_array.cpp similarity index 54% rename from doc/snippets/MatrixBase_cwise.cpp rename to doc/snippets/MatrixBase_array.cpp index f246a1a5f..f215086db 100644 --- a/doc/snippets/MatrixBase_cwise.cpp +++ b/doc/snippets/MatrixBase_array.cpp @@ -1,4 +1,4 @@ Vector3d v(1,2,3); -v.cwise() += 3; -v.cwise() -= 2; +v.array() += 3; +v.array() -= 2; cout << v << endl; diff --git a/doc/snippets/MatrixBase_array_const.cpp b/doc/snippets/MatrixBase_array_const.cpp new file mode 100644 index 000000000..cd3b26a7c --- /dev/null +++ b/doc/snippets/MatrixBase_array_const.cpp @@ -0,0 +1,4 @@ +Vector3d v(-1,2,-3); +cout << "the absolute values:" << endl << v.array().abs() << endl; +cout << "the absolute values plus one:" << endl << v.array().abs()+1 << endl; +cout << "sum of the squares: " << v.array().square().sum() << endl; diff --git a/doc/snippets/MatrixBase_colwise.cpp b/doc/snippets/MatrixBase_colwise.cpp index c8b380cee..a048beffa 100644 --- a/doc/snippets/MatrixBase_colwise.cpp +++ b/doc/snippets/MatrixBase_colwise.cpp @@ -2,4 +2,4 @@ Matrix3d m = Matrix3d::Random(); cout << "Here is the matrix m:" << endl << m << endl; cout << "Here is the sum of each column:" << endl << m.colwise().sum() << endl; cout << "Here is the maximum absolute value of each column:" - << endl << m.cwise().abs().colwise().maxCoeff() << endl; + << endl << m.cwiseAbs().colwise().maxCoeff() << endl; diff --git a/doc/snippets/MatrixBase_cwise_const.cpp b/doc/snippets/MatrixBase_cwise_const.cpp deleted file mode 100644 index 2b7e481c8..000000000 --- a/doc/snippets/MatrixBase_cwise_const.cpp +++ /dev/null @@ -1,4 +0,0 @@ -Vector3d v(-1,2,-3); -cout << "the absolute values:" << endl << v.cwise().abs() << endl; -cout << "the absolute values plus one:" << endl << v.cwise().abs().cwise()+1 << endl; -cout << "sum of the squares: " << v.cwise().square().sum() << endl; diff --git a/doc/snippets/MatrixBase_lazy.cpp b/doc/snippets/MatrixBase_lazy.cpp deleted file mode 100644 index 0177b1e51..000000000 --- a/doc/snippets/MatrixBase_lazy.cpp +++ /dev/null @@ -1,5 +0,0 @@ -Matrix2d m; m << 1,2,3,4; -cout << (m*m).lazy().row(0) << endl; - // this computes only one row of the product. By contrast, - // if we did "(m*m).row(0);" then m*m would first be evaluated into - // a temporary, because the Product expression has the EvalBeforeNestingBit. diff --git a/doc/snippets/MatrixBase_noalias.cpp b/doc/snippets/MatrixBase_noalias.cpp new file mode 100644 index 000000000..3b54a79a6 --- /dev/null +++ b/doc/snippets/MatrixBase_noalias.cpp @@ -0,0 +1,3 @@ +Matrix2d a, b, c; a << 1,2,3,4; b << 5,6,7,8; +c.noalias() = a * b; // this computes the product directly to c +cout << c << endl; diff --git a/doc/snippets/MatrixBase_rowwise.cpp b/doc/snippets/MatrixBase_rowwise.cpp index 220546540..ae93964ea 100644 --- a/doc/snippets/MatrixBase_rowwise.cpp +++ b/doc/snippets/MatrixBase_rowwise.cpp @@ -2,4 +2,4 @@ Matrix3d m = Matrix3d::Random(); cout << "Here is the matrix m:" << endl << m << endl; cout << "Here is the sum of each row:" << endl << m.rowwise().sum() << endl; cout << "Here is the maximum absolute value of each row:" - << endl << m.cwise().abs().rowwise().maxCoeff() << endl; + << endl << m.cwiseAbs().rowwise().maxCoeff() << endl; diff --git a/doc/snippets/MatrixBase_select.cpp b/doc/snippets/MatrixBase_select.cpp index 54893be1e..ae5477f02 100644 --- a/doc/snippets/MatrixBase_select.cpp +++ b/doc/snippets/MatrixBase_select.cpp @@ -1,6 +1,6 @@ MatrixXi m(3, 3); -m << 1, 2, 3, - 4, 5, 6, +m << 1, 2, 3, + 4, 5, 6, 7, 8, 9; -m = (m.cwise() >= 5).select(-m, m); +m = (m.array() >= 5).select(-m, m); cout << m << endl; diff --git a/doc/snippets/PartialRedux_count.cpp b/doc/snippets/PartialRedux_count.cpp index 914a5dee3..c7b3097e4 100644 --- a/doc/snippets/PartialRedux_count.cpp +++ b/doc/snippets/PartialRedux_count.cpp @@ -1,3 +1,3 @@ Matrix3d m = Matrix3d::Random(); cout << "Here is the matrix m:" << endl << m << endl; -cout << "Here is the count of elements larger or equal than 0.5 of each row:" << endl << (m.cwise() >= 0.5).rowwise().count() << endl; +cout << "Here is the count of elements larger or equal than 0.5 of each row:" << endl << (m.array() >= 0.5).rowwise().count() << endl; diff --git a/doc/unsupported_modules.dox b/doc/unsupported_modules.dox index f606aca0a..503ecb8b3 100644 --- a/doc/unsupported_modules.dox +++ b/doc/unsupported_modules.dox @@ -4,8 +4,8 @@ namespace Eigen { /** \defgroup Unsupported_modules Unsupported modules * * The unsupported modules are contributions from various users. They are - * provided "as is", without any support. Nevertheless, they are subject to be - * included in Eigen in the future. + * provided "as is", without any support. Nevertheless, some of them are + * subject to be included in Eigen in the future. */ // please list here all unsupported modules