Please, help us to better know about our user community by answering the following short survey: https://forms.gle/wpyrxWi18ox9Z5ae9
Eigen  3.4.0
 
Loading...
Searching...
No Matches
VectorwiseOp.h
1// This file is part of Eigen, a lightweight C++ template library
2// for linear algebra.
3//
4// Copyright (C) 2008-2019 Gael Guennebaud <gael.guennebaud@inria.fr>
5// Copyright (C) 2006-2008 Benoit Jacob <jacob.benoit.1@gmail.com>
6//
7// This Source Code Form is subject to the terms of the Mozilla
8// Public License v. 2.0. If a copy of the MPL was not distributed
9// with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
10
11#ifndef EIGEN_PARTIAL_REDUX_H
12#define EIGEN_PARTIAL_REDUX_H
13
14namespace Eigen {
15
32template< typename MatrixType, typename MemberOp, int Direction>
33class PartialReduxExpr;
34
35namespace internal {
36template<typename MatrixType, typename MemberOp, int Direction>
37struct traits<PartialReduxExpr<MatrixType, MemberOp, Direction> >
38 : traits<MatrixType>
39{
40 typedef typename MemberOp::result_type Scalar;
41 typedef typename traits<MatrixType>::StorageKind StorageKind;
42 typedef typename traits<MatrixType>::XprKind XprKind;
43 typedef typename MatrixType::Scalar InputScalar;
44 enum {
45 RowsAtCompileTime = Direction==Vertical ? 1 : MatrixType::RowsAtCompileTime,
46 ColsAtCompileTime = Direction==Horizontal ? 1 : MatrixType::ColsAtCompileTime,
47 MaxRowsAtCompileTime = Direction==Vertical ? 1 : MatrixType::MaxRowsAtCompileTime,
48 MaxColsAtCompileTime = Direction==Horizontal ? 1 : MatrixType::MaxColsAtCompileTime,
49 Flags = RowsAtCompileTime == 1 ? RowMajorBit : 0,
50 TraversalSize = Direction==Vertical ? MatrixType::RowsAtCompileTime : MatrixType::ColsAtCompileTime
51 };
52};
53}
54
55template< typename MatrixType, typename MemberOp, int Direction>
56class PartialReduxExpr : public internal::dense_xpr_base< PartialReduxExpr<MatrixType, MemberOp, Direction> >::type,
57 internal::no_assignment_operator
58{
59 public:
60
61 typedef typename internal::dense_xpr_base<PartialReduxExpr>::type Base;
62 EIGEN_DENSE_PUBLIC_INTERFACE(PartialReduxExpr)
63
64 EIGEN_DEVICE_FUNC
65 explicit PartialReduxExpr(const MatrixType& mat, const MemberOp& func = MemberOp())
66 : m_matrix(mat), m_functor(func) {}
67
68 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
69 Index rows() const EIGEN_NOEXCEPT { return (Direction==Vertical ? 1 : m_matrix.rows()); }
70 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
71 Index cols() const EIGEN_NOEXCEPT { return (Direction==Horizontal ? 1 : m_matrix.cols()); }
72
73 EIGEN_DEVICE_FUNC
74 typename MatrixType::Nested nestedExpression() const { return m_matrix; }
75
76 EIGEN_DEVICE_FUNC
77 const MemberOp& functor() const { return m_functor; }
78
79 protected:
80 typename MatrixType::Nested m_matrix;
81 const MemberOp m_functor;
82};
83
84template<typename A,typename B> struct partial_redux_dummy_func;
85
86#define EIGEN_MAKE_PARTIAL_REDUX_FUNCTOR(MEMBER,COST,VECTORIZABLE,BINARYOP) \
87 template <typename ResultType,typename Scalar> \
88 struct member_##MEMBER { \
89 EIGEN_EMPTY_STRUCT_CTOR(member_##MEMBER) \
90 typedef ResultType result_type; \
91 typedef BINARYOP<Scalar,Scalar> BinaryOp; \
92 template<int Size> struct Cost { enum { value = COST }; }; \
93 enum { Vectorizable = VECTORIZABLE }; \
94 template<typename XprType> \
95 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE \
96 ResultType operator()(const XprType& mat) const \
97 { return mat.MEMBER(); } \
98 BinaryOp binaryFunc() const { return BinaryOp(); } \
99 }
100
101#define EIGEN_MEMBER_FUNCTOR(MEMBER,COST) \
102 EIGEN_MAKE_PARTIAL_REDUX_FUNCTOR(MEMBER,COST,0,partial_redux_dummy_func)
103
104namespace internal {
105
106EIGEN_MEMBER_FUNCTOR(norm, (Size+5) * NumTraits<Scalar>::MulCost + (Size-1)*NumTraits<Scalar>::AddCost);
107EIGEN_MEMBER_FUNCTOR(stableNorm, (Size+5) * NumTraits<Scalar>::MulCost + (Size-1)*NumTraits<Scalar>::AddCost);
108EIGEN_MEMBER_FUNCTOR(blueNorm, (Size+5) * NumTraits<Scalar>::MulCost + (Size-1)*NumTraits<Scalar>::AddCost);
109EIGEN_MEMBER_FUNCTOR(hypotNorm, (Size-1) * functor_traits<scalar_hypot_op<Scalar> >::Cost );
110EIGEN_MEMBER_FUNCTOR(all, (Size-1)*NumTraits<Scalar>::AddCost);
111EIGEN_MEMBER_FUNCTOR(any, (Size-1)*NumTraits<Scalar>::AddCost);
112EIGEN_MEMBER_FUNCTOR(count, (Size-1)*NumTraits<Scalar>::AddCost);
113
114EIGEN_MAKE_PARTIAL_REDUX_FUNCTOR(sum, (Size-1)*NumTraits<Scalar>::AddCost, 1, internal::scalar_sum_op);
115EIGEN_MAKE_PARTIAL_REDUX_FUNCTOR(minCoeff, (Size-1)*NumTraits<Scalar>::AddCost, 1, internal::scalar_min_op);
116EIGEN_MAKE_PARTIAL_REDUX_FUNCTOR(maxCoeff, (Size-1)*NumTraits<Scalar>::AddCost, 1, internal::scalar_max_op);
117EIGEN_MAKE_PARTIAL_REDUX_FUNCTOR(prod, (Size-1)*NumTraits<Scalar>::MulCost, 1, internal::scalar_product_op);
118
119template <int p, typename ResultType,typename Scalar>
120struct member_lpnorm {
121 typedef ResultType result_type;
122 enum { Vectorizable = 0 };
123 template<int Size> struct Cost
124 { enum { value = (Size+5) * NumTraits<Scalar>::MulCost + (Size-1)*NumTraits<Scalar>::AddCost }; };
125 EIGEN_DEVICE_FUNC member_lpnorm() {}
126 template<typename XprType>
127 EIGEN_DEVICE_FUNC inline ResultType operator()(const XprType& mat) const
128 { return mat.template lpNorm<p>(); }
129};
130
131template <typename BinaryOpT, typename Scalar>
132struct member_redux {
133 typedef BinaryOpT BinaryOp;
134 typedef typename result_of<
135 BinaryOp(const Scalar&,const Scalar&)
136 >::type result_type;
137
138 enum { Vectorizable = functor_traits<BinaryOp>::PacketAccess };
139 template<int Size> struct Cost { enum { value = (Size-1) * functor_traits<BinaryOp>::Cost }; };
140 EIGEN_DEVICE_FUNC explicit member_redux(const BinaryOp func) : m_functor(func) {}
141 template<typename Derived>
142 EIGEN_DEVICE_FUNC inline result_type operator()(const DenseBase<Derived>& mat) const
143 { return mat.redux(m_functor); }
144 const BinaryOp& binaryFunc() const { return m_functor; }
145 const BinaryOp m_functor;
146};
147}
148
186template<typename ExpressionType, int Direction> class VectorwiseOp
187{
188 public:
189
190 typedef typename ExpressionType::Scalar Scalar;
191 typedef typename ExpressionType::RealScalar RealScalar;
193 typedef typename internal::ref_selector<ExpressionType>::non_const_type ExpressionTypeNested;
194 typedef typename internal::remove_all<ExpressionTypeNested>::type ExpressionTypeNestedCleaned;
195
196 template<template<typename OutScalar,typename InputScalar> class Functor,
197 typename ReturnScalar=Scalar> struct ReturnType
198 {
199 typedef PartialReduxExpr<ExpressionType,
200 Functor<ReturnScalar,Scalar>,
201 Direction
202 > Type;
203 };
204
205 template<typename BinaryOp> struct ReduxReturnType
206 {
207 typedef PartialReduxExpr<ExpressionType,
208 internal::member_redux<BinaryOp,Scalar>,
209 Direction
210 > Type;
211 };
212
213 enum {
214 isVertical = (Direction==Vertical) ? 1 : 0,
215 isHorizontal = (Direction==Horizontal) ? 1 : 0
216 };
217
218 protected:
219
220 template<typename OtherDerived> struct ExtendedType {
221 typedef Replicate<OtherDerived,
222 isVertical ? 1 : ExpressionType::RowsAtCompileTime,
223 isHorizontal ? 1 : ExpressionType::ColsAtCompileTime> Type;
224 };
225
228 template<typename OtherDerived>
229 EIGEN_DEVICE_FUNC
230 typename ExtendedType<OtherDerived>::Type
231 extendedTo(const DenseBase<OtherDerived>& other) const
232 {
233 EIGEN_STATIC_ASSERT(EIGEN_IMPLIES(isVertical, OtherDerived::MaxColsAtCompileTime==1),
234 YOU_PASSED_A_ROW_VECTOR_BUT_A_COLUMN_VECTOR_WAS_EXPECTED)
235 EIGEN_STATIC_ASSERT(EIGEN_IMPLIES(isHorizontal, OtherDerived::MaxRowsAtCompileTime==1),
236 YOU_PASSED_A_COLUMN_VECTOR_BUT_A_ROW_VECTOR_WAS_EXPECTED)
237 return typename ExtendedType<OtherDerived>::Type
238 (other.derived(),
239 isVertical ? 1 : m_matrix.rows(),
240 isHorizontal ? 1 : m_matrix.cols());
241 }
242
243 template<typename OtherDerived> struct OppositeExtendedType {
244 typedef Replicate<OtherDerived,
245 isHorizontal ? 1 : ExpressionType::RowsAtCompileTime,
246 isVertical ? 1 : ExpressionType::ColsAtCompileTime> Type;
247 };
248
251 template<typename OtherDerived>
252 EIGEN_DEVICE_FUNC
253 typename OppositeExtendedType<OtherDerived>::Type
254 extendedToOpposite(const DenseBase<OtherDerived>& other) const
255 {
256 EIGEN_STATIC_ASSERT(EIGEN_IMPLIES(isHorizontal, OtherDerived::MaxColsAtCompileTime==1),
257 YOU_PASSED_A_ROW_VECTOR_BUT_A_COLUMN_VECTOR_WAS_EXPECTED)
258 EIGEN_STATIC_ASSERT(EIGEN_IMPLIES(isVertical, OtherDerived::MaxRowsAtCompileTime==1),
259 YOU_PASSED_A_COLUMN_VECTOR_BUT_A_ROW_VECTOR_WAS_EXPECTED)
260 return typename OppositeExtendedType<OtherDerived>::Type
261 (other.derived(),
262 isHorizontal ? 1 : m_matrix.rows(),
263 isVertical ? 1 : m_matrix.cols());
264 }
265
266 public:
267 EIGEN_DEVICE_FUNC
268 explicit inline VectorwiseOp(ExpressionType& matrix) : m_matrix(matrix) {}
269
271 EIGEN_DEVICE_FUNC
272 inline const ExpressionType& _expression() const { return m_matrix; }
273
274 #ifdef EIGEN_PARSED_BY_DOXYGEN
278 random_access_iterator_type iterator;
280 random_access_iterator_type const_iterator;
281 #else
282 typedef internal::subvector_stl_iterator<ExpressionType, DirectionType(Direction)> iterator;
283 typedef internal::subvector_stl_iterator<const ExpressionType, DirectionType(Direction)> const_iterator;
284 typedef internal::subvector_stl_reverse_iterator<ExpressionType, DirectionType(Direction)> reverse_iterator;
285 typedef internal::subvector_stl_reverse_iterator<const ExpressionType, DirectionType(Direction)> const_reverse_iterator;
286 #endif
287
291 iterator begin() { return iterator (m_matrix, 0); }
293 const_iterator begin() const { return const_iterator(m_matrix, 0); }
295 const_iterator cbegin() const { return const_iterator(m_matrix, 0); }
296
300 reverse_iterator rbegin() { return reverse_iterator (m_matrix, m_matrix.template subVectors<DirectionType(Direction)>()-1); }
302 const_reverse_iterator rbegin() const { return const_reverse_iterator (m_matrix, m_matrix.template subVectors<DirectionType(Direction)>()-1); }
304 const_reverse_iterator crbegin() const { return const_reverse_iterator (m_matrix, m_matrix.template subVectors<DirectionType(Direction)>()-1); }
305
309 iterator end() { return iterator (m_matrix, m_matrix.template subVectors<DirectionType(Direction)>()); }
311 const_iterator end() const { return const_iterator(m_matrix, m_matrix.template subVectors<DirectionType(Direction)>()); }
313 const_iterator cend() const { return const_iterator(m_matrix, m_matrix.template subVectors<DirectionType(Direction)>()); }
314
318 reverse_iterator rend() { return reverse_iterator (m_matrix, -1); }
320 const_reverse_iterator rend() const { return const_reverse_iterator (m_matrix, -1); }
322 const_reverse_iterator crend() const { return const_reverse_iterator (m_matrix, -1); }
323
334 template<typename BinaryOp>
335 EIGEN_DEVICE_FUNC
336 const typename ReduxReturnType<BinaryOp>::Type
337 redux(const BinaryOp& func = BinaryOp()) const
338 {
339 eigen_assert(redux_length()>0 && "you are using an empty matrix");
340 return typename ReduxReturnType<BinaryOp>::Type(_expression(), internal::member_redux<BinaryOp,Scalar>(func));
341 }
342
343 typedef typename ReturnType<internal::member_minCoeff>::Type MinCoeffReturnType;
344 typedef typename ReturnType<internal::member_maxCoeff>::Type MaxCoeffReturnType;
345 typedef PartialReduxExpr<const CwiseUnaryOp<internal::scalar_abs2_op<Scalar>, const ExpressionTypeNestedCleaned>,internal::member_sum<RealScalar,RealScalar>,Direction> SquaredNormReturnType;
346 typedef CwiseUnaryOp<internal::scalar_sqrt_op<RealScalar>, const SquaredNormReturnType> NormReturnType;
347 typedef typename ReturnType<internal::member_blueNorm,RealScalar>::Type BlueNormReturnType;
348 typedef typename ReturnType<internal::member_stableNorm,RealScalar>::Type StableNormReturnType;
349 typedef typename ReturnType<internal::member_hypotNorm,RealScalar>::Type HypotNormReturnType;
350 typedef typename ReturnType<internal::member_sum>::Type SumReturnType;
351 typedef EIGEN_EXPR_BINARYOP_SCALAR_RETURN_TYPE(SumReturnType,Scalar,quotient) MeanReturnType;
352 typedef typename ReturnType<internal::member_all>::Type AllReturnType;
353 typedef typename ReturnType<internal::member_any>::Type AnyReturnType;
355 typedef typename ReturnType<internal::member_prod>::Type ProdReturnType;
356 typedef Reverse<const ExpressionType, Direction> ConstReverseReturnType;
357 typedef Reverse<ExpressionType, Direction> ReverseReturnType;
358
359 template<int p> struct LpNormReturnType {
361 };
362
375 EIGEN_DEVICE_FUNC
377 {
378 eigen_assert(redux_length()>0 && "you are using an empty matrix");
379 return MinCoeffReturnType(_expression());
380 }
381
394 EIGEN_DEVICE_FUNC
396 {
397 eigen_assert(redux_length()>0 && "you are using an empty matrix");
398 return MaxCoeffReturnType(_expression());
399 }
400
409 EIGEN_DEVICE_FUNC
411 { return SquaredNormReturnType(m_matrix.cwiseAbs2()); }
412
421 EIGEN_DEVICE_FUNC
422 const NormReturnType norm() const
423 { return NormReturnType(squaredNorm()); }
424
433 template<int p>
434 EIGEN_DEVICE_FUNC
435 const typename LpNormReturnType<p>::Type lpNorm() const
436 { return typename LpNormReturnType<p>::Type(_expression()); }
437
438
445 EIGEN_DEVICE_FUNC
447 { return BlueNormReturnType(_expression()); }
448
449
456 EIGEN_DEVICE_FUNC
458 { return StableNormReturnType(_expression()); }
459
460
467 EIGEN_DEVICE_FUNC
469 { return HypotNormReturnType(_expression()); }
470
478 EIGEN_DEVICE_FUNC
479 const SumReturnType sum() const
480 { return SumReturnType(_expression()); }
481
486 EIGEN_DEVICE_FUNC
487 const MeanReturnType mean() const
488 { return sum() / Scalar(Direction==Vertical?m_matrix.rows():m_matrix.cols()); }
489
495 EIGEN_DEVICE_FUNC
496 const AllReturnType all() const
497 { return AllReturnType(_expression()); }
498
504 EIGEN_DEVICE_FUNC
505 const AnyReturnType any() const
506 { return AnyReturnType(_expression()); }
507
517 EIGEN_DEVICE_FUNC
518 const CountReturnType count() const
519 { return CountReturnType(_expression()); }
520
528 EIGEN_DEVICE_FUNC
529 const ProdReturnType prod() const
530 { return ProdReturnType(_expression()); }
531
532
540 EIGEN_DEVICE_FUNC
542 { return ConstReverseReturnType( _expression() ); }
543
548 EIGEN_DEVICE_FUNC
550 { return ReverseReturnType( _expression() ); }
551
552 typedef Replicate<ExpressionType,(isVertical?Dynamic:1),(isHorizontal?Dynamic:1)> ReplicateReturnType;
553 EIGEN_DEVICE_FUNC
554 const ReplicateReturnType replicate(Index factor) const;
555
564 // NOTE implemented here because of sunstudio's compilation errors
565 // isVertical*Factor+isHorizontal instead of (isVertical?Factor:1) to handle CUDA bug with ternary operator
567 EIGEN_DEVICE_FUNC
568 replicate(Index factor = Factor) const
569 {
571 (_expression(),isVertical?factor:1,isHorizontal?factor:1);
572 }
573
575
577 template<typename OtherDerived>
578 EIGEN_DEVICE_FUNC
579 ExpressionType& operator=(const DenseBase<OtherDerived>& other)
580 {
581 EIGEN_STATIC_ASSERT_VECTOR_ONLY(OtherDerived)
582 EIGEN_STATIC_ASSERT_SAME_XPR_KIND(ExpressionType, OtherDerived)
583 //eigen_assert((m_matrix.isNull()) == (other.isNull())); FIXME
584 return m_matrix = extendedTo(other.derived());
585 }
586
588 template<typename OtherDerived>
589 EIGEN_DEVICE_FUNC
590 ExpressionType& operator+=(const DenseBase<OtherDerived>& other)
591 {
592 EIGEN_STATIC_ASSERT_VECTOR_ONLY(OtherDerived)
593 EIGEN_STATIC_ASSERT_SAME_XPR_KIND(ExpressionType, OtherDerived)
594 return m_matrix += extendedTo(other.derived());
595 }
596
598 template<typename OtherDerived>
599 EIGEN_DEVICE_FUNC
600 ExpressionType& operator-=(const DenseBase<OtherDerived>& other)
601 {
602 EIGEN_STATIC_ASSERT_VECTOR_ONLY(OtherDerived)
603 EIGEN_STATIC_ASSERT_SAME_XPR_KIND(ExpressionType, OtherDerived)
604 return m_matrix -= extendedTo(other.derived());
605 }
606
608 template<typename OtherDerived>
609 EIGEN_DEVICE_FUNC
610 ExpressionType& operator*=(const DenseBase<OtherDerived>& other)
611 {
612 EIGEN_STATIC_ASSERT_VECTOR_ONLY(OtherDerived)
613 EIGEN_STATIC_ASSERT_ARRAYXPR(ExpressionType)
614 EIGEN_STATIC_ASSERT_SAME_XPR_KIND(ExpressionType, OtherDerived)
615 m_matrix *= extendedTo(other.derived());
616 return m_matrix;
617 }
618
620 template<typename OtherDerived>
621 EIGEN_DEVICE_FUNC
622 ExpressionType& operator/=(const DenseBase<OtherDerived>& other)
623 {
624 EIGEN_STATIC_ASSERT_VECTOR_ONLY(OtherDerived)
625 EIGEN_STATIC_ASSERT_ARRAYXPR(ExpressionType)
626 EIGEN_STATIC_ASSERT_SAME_XPR_KIND(ExpressionType, OtherDerived)
627 m_matrix /= extendedTo(other.derived());
628 return m_matrix;
629 }
630
632 template<typename OtherDerived> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC
634 const ExpressionTypeNestedCleaned,
635 const typename ExtendedType<OtherDerived>::Type>
637 {
638 EIGEN_STATIC_ASSERT_VECTOR_ONLY(OtherDerived)
639 EIGEN_STATIC_ASSERT_SAME_XPR_KIND(ExpressionType, OtherDerived)
640 return m_matrix + extendedTo(other.derived());
641 }
642
644 template<typename OtherDerived>
645 EIGEN_DEVICE_FUNC
647 const ExpressionTypeNestedCleaned,
648 const typename ExtendedType<OtherDerived>::Type>
650 {
651 EIGEN_STATIC_ASSERT_VECTOR_ONLY(OtherDerived)
652 EIGEN_STATIC_ASSERT_SAME_XPR_KIND(ExpressionType, OtherDerived)
653 return m_matrix - extendedTo(other.derived());
654 }
655
658 template<typename OtherDerived> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC
660 const ExpressionTypeNestedCleaned,
661 const typename ExtendedType<OtherDerived>::Type>
662 EIGEN_DEVICE_FUNC
664 {
665 EIGEN_STATIC_ASSERT_VECTOR_ONLY(OtherDerived)
666 EIGEN_STATIC_ASSERT_ARRAYXPR(ExpressionType)
667 EIGEN_STATIC_ASSERT_SAME_XPR_KIND(ExpressionType, OtherDerived)
668 return m_matrix * extendedTo(other.derived());
669 }
670
673 template<typename OtherDerived>
674 EIGEN_DEVICE_FUNC
676 const ExpressionTypeNestedCleaned,
677 const typename ExtendedType<OtherDerived>::Type>
679 {
680 EIGEN_STATIC_ASSERT_VECTOR_ONLY(OtherDerived)
681 EIGEN_STATIC_ASSERT_ARRAYXPR(ExpressionType)
682 EIGEN_STATIC_ASSERT_SAME_XPR_KIND(ExpressionType, OtherDerived)
683 return m_matrix / extendedTo(other.derived());
684 }
685
690 EIGEN_DEVICE_FUNC
692 const ExpressionTypeNestedCleaned,
693 const typename OppositeExtendedType<NormReturnType>::Type>
694 normalized() const { return m_matrix.cwiseQuotient(extendedToOpposite(this->norm())); }
695
696
700 EIGEN_DEVICE_FUNC void normalize() {
701 m_matrix = this->normalized();
702 }
703
704 EIGEN_DEVICE_FUNC inline void reverseInPlace();
705
707
708 typedef Homogeneous<ExpressionType,Direction> HomogeneousReturnType;
709 EIGEN_DEVICE_FUNC
710 HomogeneousReturnType homogeneous() const;
711
712 typedef typename ExpressionType::PlainObject CrossReturnType;
713 template<typename OtherDerived>
714 EIGEN_DEVICE_FUNC
715 const CrossReturnType cross(const MatrixBase<OtherDerived>& other) const;
716
717 enum {
718 HNormalized_Size = Direction==Vertical ? internal::traits<ExpressionType>::RowsAtCompileTime
719 : internal::traits<ExpressionType>::ColsAtCompileTime,
720 HNormalized_SizeMinusOne = HNormalized_Size==Dynamic ? Dynamic : HNormalized_Size-1
721 };
722 typedef Block<const ExpressionType,
723 Direction==Vertical ? int(HNormalized_SizeMinusOne)
724 : int(internal::traits<ExpressionType>::RowsAtCompileTime),
725 Direction==Horizontal ? int(HNormalized_SizeMinusOne)
726 : int(internal::traits<ExpressionType>::ColsAtCompileTime)>
727 HNormalized_Block;
728 typedef Block<const ExpressionType,
729 Direction==Vertical ? 1 : int(internal::traits<ExpressionType>::RowsAtCompileTime),
730 Direction==Horizontal ? 1 : int(internal::traits<ExpressionType>::ColsAtCompileTime)>
731 HNormalized_Factors;
732 typedef CwiseBinaryOp<internal::scalar_quotient_op<typename internal::traits<ExpressionType>::Scalar>,
733 const HNormalized_Block,
734 const Replicate<HNormalized_Factors,
735 Direction==Vertical ? HNormalized_SizeMinusOne : 1,
736 Direction==Horizontal ? HNormalized_SizeMinusOne : 1> >
737 HNormalizedReturnType;
738
739 EIGEN_DEVICE_FUNC
740 const HNormalizedReturnType hnormalized() const;
741
742# ifdef EIGEN_VECTORWISEOP_PLUGIN
743# include EIGEN_VECTORWISEOP_PLUGIN
744# endif
745
746 protected:
747 Index redux_length() const
748 {
749 return Direction==Vertical ? m_matrix.rows() : m_matrix.cols();
750 }
751 ExpressionTypeNested m_matrix;
752};
753
754//const colwise moved to DenseBase.h due to CUDA compiler bug
755
756
761template<typename Derived>
762EIGEN_DEVICE_FUNC inline typename DenseBase<Derived>::ColwiseReturnType
764{
765 return ColwiseReturnType(derived());
766}
767
768//const rowwise moved to DenseBase.h due to CUDA compiler bug
769
770
775template<typename Derived>
776EIGEN_DEVICE_FUNC inline typename DenseBase<Derived>::RowwiseReturnType
778{
779 return RowwiseReturnType(derived());
780}
781
782} // end namespace Eigen
783
784#endif // EIGEN_PARTIAL_REDUX_H
Generic expression where a coefficient-wise binary operator is applied to two expressions.
Definition: CwiseBinaryOp.h:84
Generic expression where a coefficient-wise unary operator is applied to an expression.
Definition: CwiseUnaryOp.h:56
Base class for all dense matrices, vectors, and arrays.
Definition: DenseBase.h:47
Derived & derived()
Definition: EigenBase.h:46
Expression of one (or a set of) homogeneous vector(s)
Definition: Homogeneous.h:62
Base class for all dense matrices, vectors, and expressions.
Definition: MatrixBase.h:50
Generic expression of a partially reduxed matrix.
Definition: VectorwiseOp.h:58
Expression of the multiple replication of a matrix or vector.
Definition: Replicate.h:63
Expression of the reverse of a vector or matrix.
Definition: Reverse.h:65
Pseudo expression providing broadcasting and partial reduction operations.
Definition: VectorwiseOp.h:187
const HypotNormReturnType hypotNorm() const
Definition: VectorwiseOp.h:468
const SquaredNormReturnType squaredNorm() const
Definition: VectorwiseOp.h:410
const Replicate< ExpressionType, isVertical *Factor+isHorizontal, isHorizontal *Factor+isVertical > replicate(Index factor=Factor) const
Definition: VectorwiseOp.h:568
const ProdReturnType prod() const
Definition: VectorwiseOp.h:529
CwiseBinaryOp< internal::scalar_product_op< Scalar >, const ExpressionTypeNestedCleaned, const typename ExtendedType< OtherDerived >::Type > operator*(const DenseBase< OtherDerived > &other) const
Definition: VectorwiseOp.h:663
CwiseBinaryOp< internal::scalar_difference_op< Scalar, typename OtherDerived::Scalar >, const ExpressionTypeNestedCleaned, const typename ExtendedType< OtherDerived >::Type > operator-(const DenseBase< OtherDerived > &other) const
Definition: VectorwiseOp.h:649
const_reverse_iterator crend() const
Definition: VectorwiseOp.h:322
ExpressionType & operator-=(const DenseBase< OtherDerived > &other)
Definition: VectorwiseOp.h:600
CwiseBinaryOp< internal::scalar_sum_op< Scalar, typename OtherDerived::Scalar >, const ExpressionTypeNestedCleaned, const typename ExtendedType< OtherDerived >::Type > operator+(const DenseBase< OtherDerived > &other) const
Definition: VectorwiseOp.h:636
Eigen::Index Index
Definition: VectorwiseOp.h:192
const BlueNormReturnType blueNorm() const
Definition: VectorwiseOp.h:446
CwiseBinaryOp< internal::scalar_quotient_op< Scalar >, const ExpressionTypeNestedCleaned, const typename OppositeExtendedType< NormReturnType >::Type > normalized() const
Definition: VectorwiseOp.h:694
ExpressionType & operator/=(const DenseBase< OtherDerived > &other)
Definition: VectorwiseOp.h:622
const_iterator cend() const
Definition: VectorwiseOp.h:313
random_access_iterator_type const_iterator
Definition: VectorwiseOp.h:280
const CountReturnType count() const
Definition: VectorwiseOp.h:518
const ReplicateReturnType replicate(Index factor) const
Definition: Replicate.h:134
const MaxCoeffReturnType maxCoeff() const
Definition: VectorwiseOp.h:395
const SumReturnType sum() const
Definition: VectorwiseOp.h:479
ExpressionType & operator*=(const DenseBase< OtherDerived > &other)
Definition: VectorwiseOp.h:610
const_iterator cbegin() const
Definition: VectorwiseOp.h:295
const MinCoeffReturnType minCoeff() const
Definition: VectorwiseOp.h:376
reverse_iterator rbegin()
Definition: VectorwiseOp.h:300
const AnyReturnType any() const
Definition: VectorwiseOp.h:505
CwiseBinaryOp< internal::scalar_quotient_op< Scalar >, const ExpressionTypeNestedCleaned, const typename ExtendedType< OtherDerived >::Type > operator/(const DenseBase< OtherDerived > &other) const
Definition: VectorwiseOp.h:678
const AllReturnType all() const
Definition: VectorwiseOp.h:496
iterator end()
Definition: VectorwiseOp.h:309
random_access_iterator_type iterator
Definition: VectorwiseOp.h:278
const_reverse_iterator rbegin() const
Definition: VectorwiseOp.h:302
reverse_iterator rend()
Definition: VectorwiseOp.h:318
const MeanReturnType mean() const
Definition: VectorwiseOp.h:487
const_iterator begin() const
Definition: VectorwiseOp.h:293
ExpressionType & operator=(const DenseBase< OtherDerived > &other)
Definition: VectorwiseOp.h:579
ReverseReturnType reverse()
Definition: VectorwiseOp.h:549
void reverseInPlace()
Definition: Reverse.h:210
iterator begin()
Definition: VectorwiseOp.h:291
ExpressionType & operator+=(const DenseBase< OtherDerived > &other)
Definition: VectorwiseOp.h:590
const_iterator end() const
Definition: VectorwiseOp.h:311
void normalize()
Definition: VectorwiseOp.h:700
const StableNormReturnType stableNorm() const
Definition: VectorwiseOp.h:457
const LpNormReturnType< p >::Type lpNorm() const
Definition: VectorwiseOp.h:435
const ConstReverseReturnType reverse() const
Definition: VectorwiseOp.h:541
const_reverse_iterator crbegin() const
Definition: VectorwiseOp.h:304
const NormReturnType norm() const
Definition: VectorwiseOp.h:422
const ReduxReturnType< BinaryOp >::Type redux(const BinaryOp &func=BinaryOp()) const
Definition: VectorwiseOp.h:337
const_reverse_iterator rend() const
Definition: VectorwiseOp.h:320
static const Eigen::internal::all_t all
Definition: IndexedViewHelper.h:171
const HNormalizedReturnType hnormalized() const
column or row-wise homogeneous normalization
Definition: Homogeneous.h:198
const CrossReturnType cross(const MatrixBase< OtherDerived > &other) const
Definition: OrthoMethods.h:111
HomogeneousReturnType homogeneous() const
Definition: Homogeneous.h:150
DirectionType
Definition: Constants.h:261
@ Horizontal
Definition: Constants.h:267
@ Vertical
Definition: Constants.h:264
const unsigned int RowMajorBit
Definition: Constants.h:66
Namespace containing all symbols from the Eigen library.
Definition: Core:141
EIGEN_DEFAULT_DENSE_INDEX_TYPE Index
The Index type as used for the API.
Definition: Meta.h:74
const int Dynamic
Definition: Constants.h:22
Holds information about the various numeric (i.e. scalar) types allowed by Eigen.
Definition: NumTraits.h:233