10 #ifndef EIGEN_UNARY_FUNCTORS_H 11 #define EIGEN_UNARY_FUNCTORS_H 22 template<
typename Scalar>
struct scalar_opposite_op {
23 EIGEN_EMPTY_STRUCT_CTOR(scalar_opposite_op)
24 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Scalar operator() (
const Scalar& a)
const {
return -a; }
25 template<
typename Packet>
26 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Packet packetOp(
const Packet& a)
const 27 {
return internal::pnegate(a); }
29 template<
typename Scalar>
30 struct functor_traits<scalar_opposite_op<Scalar> >
32 Cost = NumTraits<Scalar>::AddCost,
33 PacketAccess = packet_traits<Scalar>::HasNegate };
41 template<
typename Scalar>
struct scalar_abs_op {
42 EIGEN_EMPTY_STRUCT_CTOR(scalar_abs_op)
43 typedef typename NumTraits<Scalar>::Real result_type;
44 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const result_type operator() (
const Scalar& a)
const {
return numext::abs(a); }
45 template<
typename Packet>
46 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Packet packetOp(
const Packet& a)
const 47 {
return internal::pabs(a); }
49 template<
typename Scalar>
50 struct functor_traits<scalar_abs_op<Scalar> >
53 Cost = NumTraits<Scalar>::AddCost,
54 PacketAccess = packet_traits<Scalar>::HasAbs
63 template<
typename Scalar>
struct scalar_score_coeff_op : scalar_abs_op<Scalar>
65 typedef void Score_is_abs;
67 template<
typename Scalar>
68 struct functor_traits<scalar_score_coeff_op<Scalar> > : functor_traits<scalar_abs_op<Scalar> > {};
71 template<
typename Scalar,
typename=
void>
struct abs_knowing_score
73 EIGEN_EMPTY_STRUCT_CTOR(abs_knowing_score)
74 typedef typename NumTraits<Scalar>::Real result_type;
75 template<
typename Score>
76 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const result_type operator() (
const Scalar& a,
const Score&)
const {
return numext::abs(a); }
78 template<
typename Scalar>
struct abs_knowing_score<Scalar, typename scalar_score_coeff_op<Scalar>::Score_is_abs>
80 EIGEN_EMPTY_STRUCT_CTOR(abs_knowing_score)
81 typedef typename NumTraits<Scalar>::Real result_type;
82 template<
typename Scal>
83 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const result_type operator() (
const Scal&,
const result_type& a)
const {
return a; }
91 template<
typename Scalar>
struct scalar_abs2_op {
92 EIGEN_EMPTY_STRUCT_CTOR(scalar_abs2_op)
93 typedef typename NumTraits<Scalar>::Real result_type;
95 EIGEN_STRONG_INLINE
const result_type operator() (
const Scalar& a)
const {
return numext::abs2(a); }
96 template<
typename Packet>
97 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Packet packetOp(
const Packet& a)
const 98 {
return internal::pmul(a,a); }
100 template<
typename Scalar>
101 struct functor_traits<scalar_abs2_op<Scalar> >
102 {
enum { Cost = NumTraits<Scalar>::MulCost, PacketAccess = packet_traits<Scalar>::HasAbs2 }; };
109 template<
typename Scalar>
struct scalar_conjugate_op {
110 EIGEN_EMPTY_STRUCT_CTOR(scalar_conjugate_op)
112 EIGEN_STRONG_INLINE
const Scalar operator() (
const Scalar& a)
const {
return numext::conj(a); }
113 template<
typename Packet>
114 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Packet packetOp(
const Packet& a)
const {
return internal::pconj(a); }
116 template<
typename Scalar>
117 struct functor_traits<scalar_conjugate_op<Scalar> >
129 PacketAccess = packet_traits<Scalar>::HasConj
138 template<
typename Scalar>
struct scalar_arg_op {
139 EIGEN_EMPTY_STRUCT_CTOR(scalar_arg_op)
140 typedef typename NumTraits<Scalar>::Real result_type;
141 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const result_type operator() (
const Scalar& a)
const {
return numext::arg(a); }
142 template<
typename Packet>
143 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Packet packetOp(
const Packet& a)
const 144 {
return internal::parg(a); }
146 template<
typename Scalar>
147 struct functor_traits<scalar_arg_op<Scalar> >
150 Cost = NumTraits<Scalar>::IsComplex ? 5 * NumTraits<Scalar>::MulCost : NumTraits<Scalar>::AddCost,
151 PacketAccess = packet_traits<Scalar>::HasArg
159 template<
typename Scalar,
typename NewType>
160 struct scalar_cast_op {
161 EIGEN_EMPTY_STRUCT_CTOR(scalar_cast_op)
162 typedef NewType result_type;
163 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const NewType operator() (
const Scalar& a)
const {
return cast<Scalar, NewType>(a); }
165 template<
typename Scalar,
typename NewType>
166 struct functor_traits<scalar_cast_op<Scalar,NewType> >
167 {
enum { Cost = is_same<Scalar, NewType>::value ? 0 : NumTraits<NewType>::AddCost, PacketAccess =
false }; };
174 template<
typename Scalar,
int N>
175 struct scalar_shift_right_op {
176 EIGEN_EMPTY_STRUCT_CTOR(scalar_shift_right_op)
178 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Scalar operator() (
const Scalar& a)
const 180 template<
typename Packet>
181 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Packet packetOp(
const Packet& a)
const 182 {
return internal::parithmetic_shift_right<N>(a); }
184 template<
typename Scalar,
int N>
185 struct functor_traits<scalar_shift_right_op<Scalar,N> >
186 {
enum { Cost = NumTraits<Scalar>::AddCost, PacketAccess = packet_traits<Scalar>::HasShift }; };
193 template<
typename Scalar,
int N>
194 struct scalar_shift_left_op {
195 EIGEN_EMPTY_STRUCT_CTOR(scalar_shift_left_op)
197 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Scalar operator() (
const Scalar& a)
const 199 template<
typename Packet>
200 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Packet packetOp(
const Packet& a)
const 201 {
return internal::plogical_shift_left<N>(a); }
203 template<
typename Scalar,
int N>
204 struct functor_traits<scalar_shift_left_op<Scalar,N> >
205 {
enum { Cost = NumTraits<Scalar>::AddCost, PacketAccess = packet_traits<Scalar>::HasShift }; };
212 template<
typename Scalar>
213 struct scalar_real_op {
214 EIGEN_EMPTY_STRUCT_CTOR(scalar_real_op)
215 typedef typename NumTraits<Scalar>::Real result_type;
217 EIGEN_STRONG_INLINE result_type operator() (
const Scalar& a)
const {
return numext::real(a); }
219 template<
typename Scalar>
220 struct functor_traits<scalar_real_op<Scalar> >
221 {
enum { Cost = 0, PacketAccess =
false }; };
228 template<
typename Scalar>
229 struct scalar_imag_op {
230 EIGEN_EMPTY_STRUCT_CTOR(scalar_imag_op)
231 typedef typename NumTraits<Scalar>::Real result_type;
233 EIGEN_STRONG_INLINE result_type operator() (
const Scalar& a)
const {
return numext::imag(a); }
235 template<
typename Scalar>
236 struct functor_traits<scalar_imag_op<Scalar> >
237 {
enum { Cost = 0, PacketAccess =
false }; };
244 template<
typename Scalar>
245 struct scalar_real_ref_op {
246 EIGEN_EMPTY_STRUCT_CTOR(scalar_real_ref_op)
247 typedef typename NumTraits<Scalar>::Real result_type;
249 EIGEN_STRONG_INLINE result_type& operator() (
const Scalar& a)
const {
return numext::real_ref(*const_cast<Scalar*>(&a)); }
251 template<
typename Scalar>
252 struct functor_traits<scalar_real_ref_op<Scalar> >
253 {
enum { Cost = 0, PacketAccess =
false }; };
260 template<
typename Scalar>
261 struct scalar_imag_ref_op {
262 EIGEN_EMPTY_STRUCT_CTOR(scalar_imag_ref_op)
263 typedef typename NumTraits<Scalar>::Real result_type;
265 EIGEN_STRONG_INLINE result_type& operator() (
const Scalar& a)
const {
return numext::imag_ref(*const_cast<Scalar*>(&a)); }
267 template<
typename Scalar>
268 struct functor_traits<scalar_imag_ref_op<Scalar> >
269 {
enum { Cost = 0, PacketAccess =
false }; };
277 template<
typename Scalar>
struct scalar_exp_op {
278 EIGEN_EMPTY_STRUCT_CTOR(scalar_exp_op)
279 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::exp(a); }
280 template <
typename Packet>
281 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pexp(a); }
283 template <
typename Scalar>
284 struct functor_traits<scalar_exp_op<Scalar> > {
286 PacketAccess = packet_traits<Scalar>::HasExp,
288 #ifdef EIGEN_VECTORIZE_FMA 293 ? (8 * NumTraits<Scalar>::AddCost + 6 * NumTraits<Scalar>::MulCost)
295 : (14 * NumTraits<Scalar>::AddCost +
296 6 * NumTraits<Scalar>::MulCost +
297 scalar_div_cost<Scalar,packet_traits<Scalar>::HasDiv>::value))
302 ? (21 * NumTraits<Scalar>::AddCost + 13 * NumTraits<Scalar>::MulCost)
304 : (23 * NumTraits<Scalar>::AddCost +
305 12 * NumTraits<Scalar>::MulCost +
306 scalar_div_cost<Scalar,packet_traits<Scalar>::HasDiv>::value))
317 template<
typename Scalar>
struct scalar_expm1_op {
318 EIGEN_EMPTY_STRUCT_CTOR(scalar_expm1_op)
319 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::expm1(a); }
320 template <
typename Packet>
321 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pexpm1(a); }
323 template <
typename Scalar>
324 struct functor_traits<scalar_expm1_op<Scalar> > {
326 PacketAccess = packet_traits<Scalar>::HasExpm1,
327 Cost = functor_traits<scalar_exp_op<Scalar> >::Cost
337 template<
typename Scalar>
struct scalar_log_op {
338 EIGEN_EMPTY_STRUCT_CTOR(scalar_log_op)
339 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::log(a); }
340 template <
typename Packet>
341 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::plog(a); }
343 template <
typename Scalar>
344 struct functor_traits<scalar_log_op<Scalar> > {
346 PacketAccess = packet_traits<Scalar>::HasLog,
350 #ifdef EIGEN_VECTORIZE_FMA 352 ? (20 * NumTraits<Scalar>::AddCost + 7 * NumTraits<Scalar>::MulCost)
355 ? (36 * NumTraits<Scalar>::AddCost + 14 * NumTraits<Scalar>::MulCost)
358 :
sizeof(Scalar)==4 ? 40 : 85)
368 template<
typename Scalar>
struct scalar_log1p_op {
369 EIGEN_EMPTY_STRUCT_CTOR(scalar_log1p_op)
370 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::log1p(a); }
371 template <
typename Packet>
372 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::plog1p(a); }
374 template <
typename Scalar>
375 struct functor_traits<scalar_log1p_op<Scalar> > {
377 PacketAccess = packet_traits<Scalar>::HasLog1p,
378 Cost = functor_traits<scalar_log_op<Scalar> >::Cost
388 template<
typename Scalar>
struct scalar_log10_op {
389 EIGEN_EMPTY_STRUCT_CTOR(scalar_log10_op)
390 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const { EIGEN_USING_STD(
log10)
return log10(a); }
391 template <
typename Packet>
392 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::plog10(a); }
394 template<
typename Scalar>
395 struct functor_traits<scalar_log10_op<Scalar> >
396 {
enum { Cost = 5 * NumTraits<Scalar>::MulCost, PacketAccess = packet_traits<Scalar>::HasLog10 }; };
404 template<
typename Scalar>
struct scalar_log2_op {
405 EIGEN_EMPTY_STRUCT_CTOR(scalar_log2_op)
406 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return Scalar(EIGEN_LOG2E) * numext::log(a); }
407 template <
typename Packet>
408 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::plog2(a); }
410 template<
typename Scalar>
411 struct functor_traits<scalar_log2_op<Scalar> >
412 {
enum { Cost = 5 * NumTraits<Scalar>::MulCost, PacketAccess = packet_traits<Scalar>::HasLog }; };
418 template<
typename Scalar>
struct scalar_sqrt_op {
419 EIGEN_EMPTY_STRUCT_CTOR(scalar_sqrt_op)
420 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::sqrt(a); }
421 template <
typename Packet>
422 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::psqrt(a); }
424 template <
typename Scalar>
425 struct functor_traits<scalar_sqrt_op<Scalar> > {
429 Cost = (
sizeof(Scalar) == 8 ? 28
431 : (3 * NumTraits<Scalar>::AddCost +
432 5 * NumTraits<Scalar>::MulCost)),
435 Cost = (
sizeof(Scalar) == 8 ? 28 : 14),
437 PacketAccess = packet_traits<Scalar>::HasSqrt
442 template<>
struct scalar_sqrt_op<bool> {
443 EIGEN_EMPTY_STRUCT_CTOR(scalar_sqrt_op)
444 EIGEN_DEPRECATED EIGEN_DEVICE_FUNC
inline bool operator() (
const bool& a)
const {
return a; }
445 template <
typename Packet>
446 EIGEN_DEPRECATED EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return a; }
449 struct functor_traits<scalar_sqrt_op<bool> > {
450 enum { Cost = 1, PacketAccess = packet_traits<bool>::Vectorizable };
457 template<
typename Scalar>
struct scalar_rsqrt_op {
458 EIGEN_EMPTY_STRUCT_CTOR(scalar_rsqrt_op)
459 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::rsqrt(a); }
460 template <
typename Packet>
461 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::prsqrt(a); }
464 template<
typename Scalar>
465 struct functor_traits<scalar_rsqrt_op<Scalar> >
467 Cost = 5 * NumTraits<Scalar>::MulCost,
468 PacketAccess = packet_traits<Scalar>::HasRsqrt
476 template<
typename Scalar>
struct scalar_cos_op {
477 EIGEN_EMPTY_STRUCT_CTOR(scalar_cos_op)
478 EIGEN_DEVICE_FUNC
inline Scalar operator() (
const Scalar& a)
const {
return numext::cos(a); }
479 template <
typename Packet>
480 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pcos(a); }
482 template<
typename Scalar>
483 struct functor_traits<scalar_cos_op<Scalar> >
486 Cost = 5 * NumTraits<Scalar>::MulCost,
487 PacketAccess = packet_traits<Scalar>::HasCos
495 template<
typename Scalar>
struct scalar_sin_op {
496 EIGEN_EMPTY_STRUCT_CTOR(scalar_sin_op)
497 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::sin(a); }
498 template <
typename Packet>
499 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::psin(a); }
501 template<
typename Scalar>
502 struct functor_traits<scalar_sin_op<Scalar> >
505 Cost = 5 * NumTraits<Scalar>::MulCost,
506 PacketAccess = packet_traits<Scalar>::HasSin
515 template<
typename Scalar>
struct scalar_tan_op {
516 EIGEN_EMPTY_STRUCT_CTOR(scalar_tan_op)
517 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::tan(a); }
518 template <
typename Packet>
519 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::ptan(a); }
521 template<
typename Scalar>
522 struct functor_traits<scalar_tan_op<Scalar> >
525 Cost = 5 * NumTraits<Scalar>::MulCost,
526 PacketAccess = packet_traits<Scalar>::HasTan
534 template<
typename Scalar>
struct scalar_acos_op {
535 EIGEN_EMPTY_STRUCT_CTOR(scalar_acos_op)
536 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::acos(a); }
537 template <
typename Packet>
538 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pacos(a); }
540 template<
typename Scalar>
541 struct functor_traits<scalar_acos_op<Scalar> >
544 Cost = 5 * NumTraits<Scalar>::MulCost,
545 PacketAccess = packet_traits<Scalar>::HasACos
553 template<
typename Scalar>
struct scalar_asin_op {
554 EIGEN_EMPTY_STRUCT_CTOR(scalar_asin_op)
555 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::asin(a); }
556 template <
typename Packet>
557 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pasin(a); }
559 template<
typename Scalar>
560 struct functor_traits<scalar_asin_op<Scalar> >
563 Cost = 5 * NumTraits<Scalar>::MulCost,
564 PacketAccess = packet_traits<Scalar>::HasASin
573 template<
typename Scalar>
struct scalar_atan_op {
574 EIGEN_EMPTY_STRUCT_CTOR(scalar_atan_op)
575 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::atan(a); }
576 template <
typename Packet>
577 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::patan(a); }
579 template<
typename Scalar>
580 struct functor_traits<scalar_atan_op<Scalar> >
583 Cost = 5 * NumTraits<Scalar>::MulCost,
584 PacketAccess = packet_traits<Scalar>::HasATan
592 template <
typename Scalar>
593 struct scalar_tanh_op {
594 EIGEN_EMPTY_STRUCT_CTOR(scalar_tanh_op)
595 EIGEN_DEVICE_FUNC
inline const Scalar operator()(
const Scalar& a)
const {
return numext::tanh(a); }
596 template <
typename Packet>
597 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& x)
const {
return ptanh(x); }
600 template <
typename Scalar>
601 struct functor_traits<scalar_tanh_op<Scalar> > {
603 PacketAccess = packet_traits<Scalar>::HasTanh,
604 Cost = ( (EIGEN_FAST_MATH && is_same<Scalar,float>::value)
606 #ifdef EIGEN_VECTORIZE_FMA
609 ? (2 * NumTraits<Scalar>::AddCost +
610 6 * NumTraits<Scalar>::MulCost +
611 scalar_div_cost<Scalar,packet_traits<Scalar>::HasDiv>::value)
613 ? (11 * NumTraits<Scalar>::AddCost +
614 11 * NumTraits<Scalar>::MulCost +
615 scalar_div_cost<Scalar,packet_traits<Scalar>::HasDiv>::value)
618 : (6 * NumTraits<Scalar>::AddCost +
619 3 * NumTraits<Scalar>::MulCost +
620 2 * scalar_div_cost<Scalar,packet_traits<Scalar>::HasDiv>::value +
621 functor_traits<scalar_exp_op<Scalar> >::Cost))
625 #if EIGEN_HAS_CXX11_MATH 630 template <
typename Scalar>
631 struct scalar_atanh_op {
632 EIGEN_EMPTY_STRUCT_CTOR(scalar_atanh_op)
633 EIGEN_DEVICE_FUNC
inline const Scalar operator()(
const Scalar& a)
const {
return numext::atanh(a); }
636 template <
typename Scalar>
637 struct functor_traits<scalar_atanh_op<Scalar> > {
638 enum { Cost = 5 * NumTraits<Scalar>::MulCost, PacketAccess =
false };
646 template<
typename Scalar>
struct scalar_sinh_op {
647 EIGEN_EMPTY_STRUCT_CTOR(scalar_sinh_op)
648 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::sinh(a); }
649 template <
typename Packet>
650 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::psinh(a); }
652 template<
typename Scalar>
653 struct functor_traits<scalar_sinh_op<Scalar> >
656 Cost = 5 * NumTraits<Scalar>::MulCost,
657 PacketAccess = packet_traits<Scalar>::HasSinh
661 #if EIGEN_HAS_CXX11_MATH 666 template <
typename Scalar>
667 struct scalar_asinh_op {
668 EIGEN_EMPTY_STRUCT_CTOR(scalar_asinh_op)
669 EIGEN_DEVICE_FUNC
inline const Scalar operator()(
const Scalar& a)
const {
return numext::asinh(a); }
672 template <
typename Scalar>
673 struct functor_traits<scalar_asinh_op<Scalar> > {
674 enum { Cost = 5 * NumTraits<Scalar>::MulCost, PacketAccess =
false };
682 template<
typename Scalar>
struct scalar_cosh_op {
683 EIGEN_EMPTY_STRUCT_CTOR(scalar_cosh_op)
684 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::cosh(a); }
685 template <
typename Packet>
686 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pcosh(a); }
688 template<
typename Scalar>
689 struct functor_traits<scalar_cosh_op<Scalar> >
692 Cost = 5 * NumTraits<Scalar>::MulCost,
693 PacketAccess = packet_traits<Scalar>::HasCosh
697 #if EIGEN_HAS_CXX11_MATH 702 template <
typename Scalar>
703 struct scalar_acosh_op {
704 EIGEN_EMPTY_STRUCT_CTOR(scalar_acosh_op)
705 EIGEN_DEVICE_FUNC
inline const Scalar operator()(
const Scalar& a)
const {
return numext::acosh(a); }
708 template <
typename Scalar>
709 struct functor_traits<scalar_acosh_op<Scalar> > {
710 enum { Cost = 5 * NumTraits<Scalar>::MulCost, PacketAccess =
false };
718 template<
typename Scalar>
719 struct scalar_inverse_op {
720 EIGEN_EMPTY_STRUCT_CTOR(scalar_inverse_op)
721 EIGEN_DEVICE_FUNC
inline Scalar operator() (
const Scalar& a)
const {
return Scalar(1)/a; }
722 template<
typename Packet>
723 EIGEN_DEVICE_FUNC
inline const Packet packetOp(
const Packet& a)
const 724 {
return internal::pdiv(pset1<Packet>(Scalar(1)),a); }
726 template <
typename Scalar>
727 struct functor_traits<scalar_inverse_op<Scalar> > {
729 PacketAccess = packet_traits<Scalar>::HasDiv,
730 Cost = scalar_div_cost<Scalar, PacketAccess>::value
738 template<
typename Scalar>
739 struct scalar_square_op {
740 EIGEN_EMPTY_STRUCT_CTOR(scalar_square_op)
741 EIGEN_DEVICE_FUNC
inline Scalar operator() (
const Scalar& a)
const {
return a*a; }
742 template<
typename Packet>
743 EIGEN_DEVICE_FUNC
inline const Packet packetOp(
const Packet& a)
const 744 {
return internal::pmul(a,a); }
746 template<
typename Scalar>
747 struct functor_traits<scalar_square_op<Scalar> >
748 {
enum { Cost = NumTraits<Scalar>::MulCost, PacketAccess = packet_traits<Scalar>::HasMul }; };
752 struct scalar_square_op<bool> {
753 EIGEN_EMPTY_STRUCT_CTOR(scalar_square_op)
754 EIGEN_DEPRECATED EIGEN_DEVICE_FUNC
inline bool operator() (
const bool& a)
const {
return a; }
755 template<
typename Packet>
756 EIGEN_DEPRECATED EIGEN_DEVICE_FUNC
inline const Packet packetOp(
const Packet& a)
const 760 struct functor_traits<scalar_square_op<bool> >
761 {
enum { Cost = 0, PacketAccess = packet_traits<bool>::Vectorizable }; };
767 template<
typename Scalar>
768 struct scalar_cube_op {
769 EIGEN_EMPTY_STRUCT_CTOR(scalar_cube_op)
770 EIGEN_DEVICE_FUNC
inline Scalar operator() (
const Scalar& a)
const {
return a*a*a; }
771 template<
typename Packet>
772 EIGEN_DEVICE_FUNC
inline const Packet packetOp(
const Packet& a)
const 773 {
return internal::pmul(a,pmul(a,a)); }
775 template<
typename Scalar>
776 struct functor_traits<scalar_cube_op<Scalar> >
777 {
enum { Cost = 2*NumTraits<Scalar>::MulCost, PacketAccess = packet_traits<Scalar>::HasMul }; };
781 struct scalar_cube_op<bool> {
782 EIGEN_EMPTY_STRUCT_CTOR(scalar_cube_op)
783 EIGEN_DEPRECATED EIGEN_DEVICE_FUNC
inline bool operator() (
const bool& a)
const {
return a; }
784 template<
typename Packet>
785 EIGEN_DEPRECATED EIGEN_DEVICE_FUNC
inline const Packet packetOp(
const Packet& a)
const 789 struct functor_traits<scalar_cube_op<bool> >
790 {
enum { Cost = 0, PacketAccess = packet_traits<bool>::Vectorizable }; };
796 template<
typename Scalar>
struct scalar_round_op {
797 EIGEN_EMPTY_STRUCT_CTOR(scalar_round_op)
798 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Scalar operator() (
const Scalar& a)
const {
return numext::round(a); }
799 template <
typename Packet>
800 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pround(a); }
802 template<
typename Scalar>
803 struct functor_traits<scalar_round_op<Scalar> >
806 Cost = NumTraits<Scalar>::MulCost,
807 PacketAccess = packet_traits<Scalar>::HasRound
815 template<
typename Scalar>
struct scalar_floor_op {
816 EIGEN_EMPTY_STRUCT_CTOR(scalar_floor_op)
817 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Scalar operator() (
const Scalar& a)
const {
return numext::floor(a); }
818 template <
typename Packet>
819 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pfloor(a); }
821 template<
typename Scalar>
822 struct functor_traits<scalar_floor_op<Scalar> >
825 Cost = NumTraits<Scalar>::MulCost,
826 PacketAccess = packet_traits<Scalar>::HasFloor
834 template<
typename Scalar>
struct scalar_rint_op {
835 EIGEN_EMPTY_STRUCT_CTOR(scalar_rint_op)
836 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Scalar operator() (
const Scalar& a)
const {
return numext::rint(a); }
837 template <
typename Packet>
838 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::print(a); }
840 template<
typename Scalar>
841 struct functor_traits<scalar_rint_op<Scalar> >
844 Cost = NumTraits<Scalar>::MulCost,
845 PacketAccess = packet_traits<Scalar>::HasRint
853 template<
typename Scalar>
struct scalar_ceil_op {
854 EIGEN_EMPTY_STRUCT_CTOR(scalar_ceil_op)
855 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Scalar operator() (
const Scalar& a)
const {
return numext::ceil(a); }
856 template <
typename Packet>
857 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pceil(a); }
859 template<
typename Scalar>
860 struct functor_traits<scalar_ceil_op<Scalar> >
863 Cost = NumTraits<Scalar>::MulCost,
864 PacketAccess = packet_traits<Scalar>::HasCeil
872 template<
typename Scalar>
struct scalar_isnan_op {
873 EIGEN_EMPTY_STRUCT_CTOR(scalar_isnan_op)
874 typedef bool result_type;
875 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE result_type operator() (
const Scalar& a)
const {
876 #if defined(SYCL_DEVICE_ONLY) 877 return numext::isnan(a);
879 return (numext::isnan)(a);
883 template<
typename Scalar>
884 struct functor_traits<scalar_isnan_op<Scalar> >
887 Cost = NumTraits<Scalar>::MulCost,
896 template<
typename Scalar>
struct scalar_isinf_op {
897 EIGEN_EMPTY_STRUCT_CTOR(scalar_isinf_op)
898 typedef bool result_type;
899 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE result_type operator() (
const Scalar& a)
const {
900 #if defined(SYCL_DEVICE_ONLY) 901 return numext::isinf(a);
903 return (numext::isinf)(a);
907 template<
typename Scalar>
908 struct functor_traits<scalar_isinf_op<Scalar> >
911 Cost = NumTraits<Scalar>::MulCost,
920 template<
typename Scalar>
struct scalar_isfinite_op {
921 EIGEN_EMPTY_STRUCT_CTOR(scalar_isfinite_op)
922 typedef bool result_type;
923 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE result_type operator() (
const Scalar& a)
const {
924 #if defined(SYCL_DEVICE_ONLY) 925 return numext::isfinite(a);
927 return (numext::isfinite)(a);
931 template<
typename Scalar>
932 struct functor_traits<scalar_isfinite_op<Scalar> >
935 Cost = NumTraits<Scalar>::MulCost,
945 template<
typename Scalar>
struct scalar_boolean_not_op {
946 EIGEN_EMPTY_STRUCT_CTOR(scalar_boolean_not_op)
947 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
bool operator() (
const bool& a)
const {
return !a; }
949 template<
typename Scalar>
950 struct functor_traits<scalar_boolean_not_op<Scalar> > {
952 Cost = NumTraits<bool>::AddCost,
961 template<typename Scalar,bool is_complex=(NumTraits<Scalar>::IsComplex!=0),
bool is_integer=(NumTraits<Scalar>::IsInteger!=0) >
struct scalar_sign_op;
962 template<
typename Scalar>
963 struct scalar_sign_op<Scalar, false, true> {
964 EIGEN_EMPTY_STRUCT_CTOR(scalar_sign_op)
965 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const 967 return Scalar( (a>Scalar(0)) - (a<Scalar(0)) );
974 template<
typename Scalar>
975 struct scalar_sign_op<Scalar, false, false> {
976 EIGEN_EMPTY_STRUCT_CTOR(scalar_sign_op)
977 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const 979 return (numext::isnan)(a) ? a : Scalar( (a>Scalar(0)) - (a<Scalar(0)) );
986 template<
typename Scalar,
bool is_
integer>
987 struct scalar_sign_op<Scalar,true, is_integer> {
988 EIGEN_EMPTY_STRUCT_CTOR(scalar_sign_op)
989 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const 991 typedef typename NumTraits<Scalar>::Real real_type;
992 real_type aa = numext::abs(a);
993 if (aa==real_type(0))
995 aa = real_type(1)/aa;
996 return Scalar(a.real()*aa, a.imag()*aa );
1002 template<
typename Scalar>
1003 struct functor_traits<scalar_sign_op<Scalar> >
1006 NumTraits<Scalar>::IsComplex
1007 ? ( 8*NumTraits<Scalar>::MulCost )
1008 : ( 3*NumTraits<Scalar>::AddCost),
1009 PacketAccess = packet_traits<Scalar>::HasSign
1017 template <
typename T>
1018 struct scalar_logistic_op {
1019 EIGEN_EMPTY_STRUCT_CTOR(scalar_logistic_op)
1020 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE T operator()(
const T& x)
const {
1024 template <
typename Packet> EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
1025 Packet packetOp(
const Packet& x)
const {
1026 const Packet one = pset1<Packet>(T(1));
1027 return pdiv(one, padd(one, pexp(pnegate(x))));
1031 #ifndef EIGEN_GPU_COMPILE_PHASE 1044 struct scalar_logistic_op<float> {
1045 EIGEN_EMPTY_STRUCT_CTOR(scalar_logistic_op)
1046 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
float operator()(
const float& x)
const {
1050 template <
typename Packet> EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
1051 Packet packetOp(
const Packet& _x)
const {
1052 const Packet cutoff_lower = pset1<Packet>(-9.f);
1053 const Packet lt_mask = pcmp_lt<Packet>(_x, cutoff_lower);
1054 const bool any_small = predux_any(lt_mask);
1058 #ifdef EIGEN_VECTORIZE_FMA 1059 const Packet cutoff_upper = pset1<Packet>(15.7243833541870117f);
1061 const Packet cutoff_upper = pset1<Packet>(15.6437711715698242f);
1063 const Packet x = pmin(_x, cutoff_upper);
1066 const Packet alpha_1 = pset1<Packet>(2.48287947061529e-01f);
1067 const Packet alpha_3 = pset1<Packet>(8.51377133304701e-03f);
1068 const Packet alpha_5 = pset1<Packet>(6.08574864600143e-05f);
1069 const Packet alpha_7 = pset1<Packet>(1.15627324459942e-07f);
1070 const Packet alpha_9 = pset1<Packet>(4.37031012579801e-11f);
1073 const Packet beta_0 = pset1<Packet>(9.93151921023180e-01f);
1074 const Packet beta_2 = pset1<Packet>(1.16817656904453e-01f);
1075 const Packet beta_4 = pset1<Packet>(1.70198817374094e-03f);
1076 const Packet beta_6 = pset1<Packet>(6.29106785017040e-06f);
1077 const Packet beta_8 = pset1<Packet>(5.76102136993427e-09f);
1078 const Packet beta_10 = pset1<Packet>(6.10247389755681e-13f);
1081 const Packet x2 = pmul(x, x);
1084 Packet p = pmadd(x2, alpha_9, alpha_7);
1085 p = pmadd(x2, p, alpha_5);
1086 p = pmadd(x2, p, alpha_3);
1087 p = pmadd(x2, p, alpha_1);
1091 Packet q = pmadd(x2, beta_10, beta_8);
1092 q = pmadd(x2, q, beta_6);
1093 q = pmadd(x2, q, beta_4);
1094 q = pmadd(x2, q, beta_2);
1095 q = pmadd(x2, q, beta_0);
1097 const Packet
logistic = padd(pdiv(p, q), pset1<Packet>(0.5f));
1098 if (EIGEN_PREDICT_FALSE(any_small)) {
1099 const Packet exponential = pexp(_x);
1100 return pselect(lt_mask, exponential, logistic);
1106 #endif // #ifndef EIGEN_GPU_COMPILE_PHASE 1108 template <
typename T>
1109 struct functor_traits<scalar_logistic_op<T> > {
1113 Cost = scalar_div_cost<T, packet_traits<T>::HasDiv>::value +
1114 (internal::is_same<T, float>::value
1115 ? NumTraits<T>::AddCost * 15 + NumTraits<T>::MulCost * 11
1116 : NumTraits<T>::AddCost * 2 +
1117 functor_traits<scalar_exp_op<T> >::Cost),
1119 packet_traits<T>::HasAdd && packet_traits<T>::HasDiv &&
1120 (internal::is_same<T, float>::value
1121 ? packet_traits<T>::HasMul && packet_traits<T>::HasMax &&
1122 packet_traits<T>::HasMin
1123 : packet_traits<T>::HasNegate && packet_traits<T>::HasExp)
1131 #endif // EIGEN_FUNCTORS_H Namespace containing all symbols from the Eigen library.
Definition: Core:141
Definition: Eigen_Colamd.h:50
const Eigen::CwiseUnaryOp< Eigen::internal::scalar_log10_op< typename Derived::Scalar >, const Derived > log10(const Eigen::ArrayBase< Derived > &x)
const Eigen::CwiseUnaryOp< Eigen::internal::scalar_logistic_op< typename Derived::Scalar >, const Derived > logistic(const Eigen::ArrayBase< Derived > &x)