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 {
using numext::conj;
return 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> >
120 Cost = NumTraits<Scalar>::IsComplex ? NumTraits<Scalar>::AddCost : 0,
121 PacketAccess = packet_traits<Scalar>::HasConj
130 template<
typename Scalar>
struct scalar_arg_op {
131 EIGEN_EMPTY_STRUCT_CTOR(scalar_arg_op)
132 typedef typename NumTraits<Scalar>::Real result_type;
133 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const result_type operator() (
const Scalar& a)
const {
using numext::arg;
return arg(a); }
134 template<
typename Packet>
135 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Packet packetOp(
const Packet& a)
const 136 {
return internal::parg(a); }
138 template<
typename Scalar>
139 struct functor_traits<scalar_arg_op<Scalar> >
142 Cost = NumTraits<Scalar>::IsComplex ? 5 * NumTraits<Scalar>::MulCost : NumTraits<Scalar>::AddCost,
143 PacketAccess = packet_traits<Scalar>::HasArg
151 template<
typename Scalar,
typename NewType>
152 struct scalar_cast_op {
153 EIGEN_EMPTY_STRUCT_CTOR(scalar_cast_op)
154 typedef NewType result_type;
155 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const NewType operator() (
const Scalar& a)
const {
return cast<Scalar, NewType>(a); }
157 template<
typename Scalar,
typename NewType>
158 struct functor_traits<scalar_cast_op<Scalar,NewType> >
159 {
enum { Cost = is_same<Scalar, NewType>::value ? 0 : NumTraits<NewType>::AddCost, PacketAccess =
false }; };
166 template<
typename Scalar>
167 struct scalar_real_op {
168 EIGEN_EMPTY_STRUCT_CTOR(scalar_real_op)
169 typedef typename NumTraits<Scalar>::Real result_type;
171 EIGEN_STRONG_INLINE result_type operator() (
const Scalar& a)
const {
return numext::real(a); }
173 template<
typename Scalar>
174 struct functor_traits<scalar_real_op<Scalar> >
175 {
enum { Cost = 0, PacketAccess =
false }; };
182 template<
typename Scalar>
183 struct scalar_imag_op {
184 EIGEN_EMPTY_STRUCT_CTOR(scalar_imag_op)
185 typedef typename NumTraits<Scalar>::Real result_type;
187 EIGEN_STRONG_INLINE result_type operator() (
const Scalar& a)
const {
return numext::imag(a); }
189 template<
typename Scalar>
190 struct functor_traits<scalar_imag_op<Scalar> >
191 {
enum { Cost = 0, PacketAccess =
false }; };
198 template<
typename Scalar>
199 struct scalar_real_ref_op {
200 EIGEN_EMPTY_STRUCT_CTOR(scalar_real_ref_op)
201 typedef typename NumTraits<Scalar>::Real result_type;
203 EIGEN_STRONG_INLINE result_type& operator() (
const Scalar& a)
const {
return numext::real_ref(*const_cast<Scalar*>(&a)); }
205 template<
typename Scalar>
206 struct functor_traits<scalar_real_ref_op<Scalar> >
207 {
enum { Cost = 0, PacketAccess =
false }; };
214 template<
typename Scalar>
215 struct scalar_imag_ref_op {
216 EIGEN_EMPTY_STRUCT_CTOR(scalar_imag_ref_op)
217 typedef typename NumTraits<Scalar>::Real result_type;
219 EIGEN_STRONG_INLINE result_type& operator() (
const Scalar& a)
const {
return numext::imag_ref(*const_cast<Scalar*>(&a)); }
221 template<
typename Scalar>
222 struct functor_traits<scalar_imag_ref_op<Scalar> >
223 {
enum { Cost = 0, PacketAccess =
false }; };
231 template<
typename Scalar>
struct scalar_exp_op {
232 EIGEN_EMPTY_STRUCT_CTOR(scalar_exp_op)
233 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::exp(a); }
234 template <
typename Packet>
235 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pexp(a); }
237 template <
typename Scalar>
238 struct functor_traits<scalar_exp_op<Scalar> > {
240 PacketAccess = packet_traits<Scalar>::HasExp,
242 #ifdef EIGEN_VECTORIZE_FMA 247 ? (8 * NumTraits<Scalar>::AddCost + 6 * NumTraits<Scalar>::MulCost)
249 : (14 * NumTraits<Scalar>::AddCost +
250 6 * NumTraits<Scalar>::MulCost +
251 scalar_div_cost<Scalar,packet_traits<Scalar>::HasDiv>::value))
256 ? (21 * NumTraits<Scalar>::AddCost + 13 * NumTraits<Scalar>::MulCost)
258 : (23 * NumTraits<Scalar>::AddCost +
259 12 * NumTraits<Scalar>::MulCost +
260 scalar_div_cost<Scalar,packet_traits<Scalar>::HasDiv>::value))
271 template<
typename Scalar>
struct scalar_log_op {
272 EIGEN_EMPTY_STRUCT_CTOR(scalar_log_op)
273 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::log(a); }
274 template <
typename Packet>
275 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::plog(a); }
277 template <
typename Scalar>
278 struct functor_traits<scalar_log_op<Scalar> > {
280 PacketAccess = packet_traits<Scalar>::HasLog,
284 #ifdef EIGEN_VECTORIZE_FMA 286 ? (20 * NumTraits<Scalar>::AddCost + 7 * NumTraits<Scalar>::MulCost)
289 ? (36 * NumTraits<Scalar>::AddCost + 14 * NumTraits<Scalar>::MulCost)
292 :
sizeof(Scalar)==4 ? 40 : 85)
302 template<
typename Scalar>
struct scalar_log1p_op {
303 EIGEN_EMPTY_STRUCT_CTOR(scalar_log1p_op)
304 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::log1p(a); }
305 template <
typename Packet>
306 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::plog1p(a); }
308 template <
typename Scalar>
309 struct functor_traits<scalar_log1p_op<Scalar> > {
311 PacketAccess = packet_traits<Scalar>::HasLog1p,
312 Cost = functor_traits<scalar_log_op<Scalar> >::Cost
322 template<
typename Scalar>
struct scalar_log10_op {
323 EIGEN_EMPTY_STRUCT_CTOR(scalar_log10_op)
324 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const { EIGEN_USING_STD_MATH(
log10)
return log10(a); }
325 template <
typename Packet>
326 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::plog10(a); }
328 template<
typename Scalar>
329 struct functor_traits<scalar_log10_op<Scalar> >
330 {
enum { Cost = 5 * NumTraits<Scalar>::MulCost, PacketAccess = packet_traits<Scalar>::HasLog10 }; };
336 template<
typename Scalar>
struct scalar_sqrt_op {
337 EIGEN_EMPTY_STRUCT_CTOR(scalar_sqrt_op)
338 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::sqrt(a); }
339 template <
typename Packet>
340 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::psqrt(a); }
342 template <
typename Scalar>
343 struct functor_traits<scalar_sqrt_op<Scalar> > {
347 Cost = (
sizeof(Scalar) == 8 ? 28
349 : (3 * NumTraits<Scalar>::AddCost +
350 5 * NumTraits<Scalar>::MulCost)),
353 Cost = (
sizeof(Scalar) == 8 ? 28 : 14),
355 PacketAccess = packet_traits<Scalar>::HasSqrt
363 template<
typename Scalar>
struct scalar_rsqrt_op {
364 EIGEN_EMPTY_STRUCT_CTOR(scalar_rsqrt_op)
365 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return Scalar(1)/numext::sqrt(a); }
366 template <
typename Packet>
367 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::prsqrt(a); }
370 template<
typename Scalar>
371 struct functor_traits<scalar_rsqrt_op<Scalar> >
373 Cost = 5 * NumTraits<Scalar>::MulCost,
374 PacketAccess = packet_traits<Scalar>::HasRsqrt
382 template<
typename Scalar>
struct scalar_cos_op {
383 EIGEN_EMPTY_STRUCT_CTOR(scalar_cos_op)
384 EIGEN_DEVICE_FUNC
inline Scalar operator() (
const Scalar& a)
const {
return numext::cos(a); }
385 template <
typename Packet>
386 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pcos(a); }
388 template<
typename Scalar>
389 struct functor_traits<scalar_cos_op<Scalar> >
392 Cost = 5 * NumTraits<Scalar>::MulCost,
393 PacketAccess = packet_traits<Scalar>::HasCos
401 template<
typename Scalar>
struct scalar_sin_op {
402 EIGEN_EMPTY_STRUCT_CTOR(scalar_sin_op)
403 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::sin(a); }
404 template <
typename Packet>
405 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::psin(a); }
407 template<
typename Scalar>
408 struct functor_traits<scalar_sin_op<Scalar> >
411 Cost = 5 * NumTraits<Scalar>::MulCost,
412 PacketAccess = packet_traits<Scalar>::HasSin
421 template<
typename Scalar>
struct scalar_tan_op {
422 EIGEN_EMPTY_STRUCT_CTOR(scalar_tan_op)
423 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::tan(a); }
424 template <
typename Packet>
425 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::ptan(a); }
427 template<
typename Scalar>
428 struct functor_traits<scalar_tan_op<Scalar> >
431 Cost = 5 * NumTraits<Scalar>::MulCost,
432 PacketAccess = packet_traits<Scalar>::HasTan
440 template<
typename Scalar>
struct scalar_acos_op {
441 EIGEN_EMPTY_STRUCT_CTOR(scalar_acos_op)
442 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::acos(a); }
443 template <
typename Packet>
444 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pacos(a); }
446 template<
typename Scalar>
447 struct functor_traits<scalar_acos_op<Scalar> >
450 Cost = 5 * NumTraits<Scalar>::MulCost,
451 PacketAccess = packet_traits<Scalar>::HasACos
459 template<
typename Scalar>
struct scalar_asin_op {
460 EIGEN_EMPTY_STRUCT_CTOR(scalar_asin_op)
461 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::asin(a); }
462 template <
typename Packet>
463 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pasin(a); }
465 template<
typename Scalar>
466 struct functor_traits<scalar_asin_op<Scalar> >
469 Cost = 5 * NumTraits<Scalar>::MulCost,
470 PacketAccess = packet_traits<Scalar>::HasASin
479 template<
typename Scalar>
struct scalar_atan_op {
480 EIGEN_EMPTY_STRUCT_CTOR(scalar_atan_op)
481 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::atan(a); }
482 template <
typename Packet>
483 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::patan(a); }
485 template<
typename Scalar>
486 struct functor_traits<scalar_atan_op<Scalar> >
489 Cost = 5 * NumTraits<Scalar>::MulCost,
490 PacketAccess = packet_traits<Scalar>::HasATan
498 template <
typename Scalar>
499 struct scalar_tanh_op {
500 EIGEN_EMPTY_STRUCT_CTOR(scalar_tanh_op)
501 EIGEN_DEVICE_FUNC
inline const Scalar operator()(
const Scalar& a)
const {
return numext::tanh(a); }
502 template <
typename Packet>
503 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& x)
const {
return ptanh(x); }
506 template <
typename Scalar>
507 struct functor_traits<scalar_tanh_op<Scalar> > {
509 PacketAccess = packet_traits<Scalar>::HasTanh,
510 Cost = ( (EIGEN_FAST_MATH && is_same<Scalar,float>::value)
512 #ifdef EIGEN_VECTORIZE_FMA
515 ? (2 * NumTraits<Scalar>::AddCost +
516 6 * NumTraits<Scalar>::MulCost +
517 scalar_div_cost<Scalar,packet_traits<Scalar>::HasDiv>::value)
519 ? (11 * NumTraits<Scalar>::AddCost +
520 11 * NumTraits<Scalar>::MulCost +
521 scalar_div_cost<Scalar,packet_traits<Scalar>::HasDiv>::value)
524 : (6 * NumTraits<Scalar>::AddCost +
525 3 * NumTraits<Scalar>::MulCost +
526 2 * scalar_div_cost<Scalar,packet_traits<Scalar>::HasDiv>::value +
527 functor_traits<scalar_exp_op<Scalar> >::Cost))
535 template<
typename Scalar>
struct scalar_sinh_op {
536 EIGEN_EMPTY_STRUCT_CTOR(scalar_sinh_op)
537 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::sinh(a); }
538 template <
typename Packet>
539 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::psinh(a); }
541 template<
typename Scalar>
542 struct functor_traits<scalar_sinh_op<Scalar> >
545 Cost = 5 * NumTraits<Scalar>::MulCost,
546 PacketAccess = packet_traits<Scalar>::HasSinh
554 template<
typename Scalar>
struct scalar_cosh_op {
555 EIGEN_EMPTY_STRUCT_CTOR(scalar_cosh_op)
556 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const {
return numext::cosh(a); }
557 template <
typename Packet>
558 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pcosh(a); }
560 template<
typename Scalar>
561 struct functor_traits<scalar_cosh_op<Scalar> >
564 Cost = 5 * NumTraits<Scalar>::MulCost,
565 PacketAccess = packet_traits<Scalar>::HasCosh
573 template<
typename Scalar>
574 struct scalar_inverse_op {
575 EIGEN_EMPTY_STRUCT_CTOR(scalar_inverse_op)
576 EIGEN_DEVICE_FUNC
inline Scalar operator() (
const Scalar& a)
const {
return Scalar(1)/a; }
577 template<
typename Packet>
578 EIGEN_DEVICE_FUNC
inline const Packet packetOp(
const Packet& a)
const 579 {
return internal::pdiv(pset1<Packet>(Scalar(1)),a); }
581 template<
typename Scalar>
582 struct functor_traits<scalar_inverse_op<Scalar> >
583 {
enum { Cost = NumTraits<Scalar>::MulCost, PacketAccess = packet_traits<Scalar>::HasDiv }; };
589 template<
typename Scalar>
590 struct scalar_square_op {
591 EIGEN_EMPTY_STRUCT_CTOR(scalar_square_op)
592 EIGEN_DEVICE_FUNC
inline Scalar operator() (
const Scalar& a)
const {
return a*a; }
593 template<
typename Packet>
594 EIGEN_DEVICE_FUNC
inline const Packet packetOp(
const Packet& a)
const 595 {
return internal::pmul(a,a); }
597 template<
typename Scalar>
598 struct functor_traits<scalar_square_op<Scalar> >
599 {
enum { Cost = NumTraits<Scalar>::MulCost, PacketAccess = packet_traits<Scalar>::HasMul }; };
605 template<
typename Scalar>
606 struct scalar_cube_op {
607 EIGEN_EMPTY_STRUCT_CTOR(scalar_cube_op)
608 EIGEN_DEVICE_FUNC
inline Scalar operator() (
const Scalar& a)
const {
return a*a*a; }
609 template<
typename Packet>
610 EIGEN_DEVICE_FUNC
inline const Packet packetOp(
const Packet& a)
const 611 {
return internal::pmul(a,pmul(a,a)); }
613 template<
typename Scalar>
614 struct functor_traits<scalar_cube_op<Scalar> >
615 {
enum { Cost = 2*NumTraits<Scalar>::MulCost, PacketAccess = packet_traits<Scalar>::HasMul }; };
621 template<
typename Scalar>
struct scalar_round_op {
622 EIGEN_EMPTY_STRUCT_CTOR(scalar_round_op)
623 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Scalar operator() (
const Scalar& a)
const {
return numext::round(a); }
624 template <
typename Packet>
625 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pround(a); }
627 template<
typename Scalar>
628 struct functor_traits<scalar_round_op<Scalar> >
631 Cost = NumTraits<Scalar>::MulCost,
632 PacketAccess = packet_traits<Scalar>::HasRound
640 template<
typename Scalar>
struct scalar_floor_op {
641 EIGEN_EMPTY_STRUCT_CTOR(scalar_floor_op)
642 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Scalar operator() (
const Scalar& a)
const {
return numext::floor(a); }
643 template <
typename Packet>
644 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pfloor(a); }
646 template<
typename Scalar>
647 struct functor_traits<scalar_floor_op<Scalar> >
650 Cost = NumTraits<Scalar>::MulCost,
651 PacketAccess = packet_traits<Scalar>::HasFloor
659 template<
typename Scalar>
struct scalar_ceil_op {
660 EIGEN_EMPTY_STRUCT_CTOR(scalar_ceil_op)
661 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
const Scalar operator() (
const Scalar& a)
const {
return numext::ceil(a); }
662 template <
typename Packet>
663 EIGEN_DEVICE_FUNC
inline Packet packetOp(
const Packet& a)
const {
return internal::pceil(a); }
665 template<
typename Scalar>
666 struct functor_traits<scalar_ceil_op<Scalar> >
669 Cost = NumTraits<Scalar>::MulCost,
670 PacketAccess = packet_traits<Scalar>::HasCeil
678 template<
typename Scalar>
struct scalar_isnan_op {
679 EIGEN_EMPTY_STRUCT_CTOR(scalar_isnan_op)
680 typedef bool result_type;
681 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE result_type operator() (
const Scalar& a)
const {
return (numext::isnan)(a); }
683 template<
typename Scalar>
684 struct functor_traits<scalar_isnan_op<Scalar> >
687 Cost = NumTraits<Scalar>::MulCost,
696 template<
typename Scalar>
struct scalar_isinf_op {
697 EIGEN_EMPTY_STRUCT_CTOR(scalar_isinf_op)
698 typedef bool result_type;
699 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE result_type operator() (
const Scalar& a)
const {
return (numext::isinf)(a); }
701 template<
typename Scalar>
702 struct functor_traits<scalar_isinf_op<Scalar> >
705 Cost = NumTraits<Scalar>::MulCost,
714 template<
typename Scalar>
struct scalar_isfinite_op {
715 EIGEN_EMPTY_STRUCT_CTOR(scalar_isfinite_op)
716 typedef bool result_type;
717 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE result_type operator() (
const Scalar& a)
const {
return (numext::isfinite)(a); }
719 template<
typename Scalar>
720 struct functor_traits<scalar_isfinite_op<Scalar> >
723 Cost = NumTraits<Scalar>::MulCost,
733 template<
typename Scalar>
struct scalar_boolean_not_op {
734 EIGEN_EMPTY_STRUCT_CTOR(scalar_boolean_not_op)
735 EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
bool operator() (
const bool& a)
const {
return !a; }
737 template<
typename Scalar>
738 struct functor_traits<scalar_boolean_not_op<Scalar> > {
740 Cost = NumTraits<bool>::AddCost,
749 template<typename Scalar,bool iscpx=(NumTraits<Scalar>::IsComplex!=0) >
struct scalar_sign_op;
750 template<
typename Scalar>
751 struct scalar_sign_op<Scalar,false> {
752 EIGEN_EMPTY_STRUCT_CTOR(scalar_sign_op)
753 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const 755 return Scalar( (a>Scalar(0)) - (a<Scalar(0)) );
761 template<
typename Scalar>
762 struct scalar_sign_op<Scalar,true> {
763 EIGEN_EMPTY_STRUCT_CTOR(scalar_sign_op)
764 EIGEN_DEVICE_FUNC
inline const Scalar operator() (
const Scalar& a)
const 766 typedef typename NumTraits<Scalar>::Real real_type;
767 real_type aa = numext::abs(a);
768 if (aa==real_type(0))
770 aa = real_type(1)/aa;
771 return Scalar(
real(a)*aa,
imag(a)*aa );
777 template<
typename Scalar>
778 struct functor_traits<scalar_sign_op<Scalar> >
781 NumTraits<Scalar>::IsComplex
782 ? ( 8*NumTraits<Scalar>::MulCost )
783 : ( 3*NumTraits<Scalar>::AddCost),
784 PacketAccess = packet_traits<Scalar>::HasSign
792 #endif // EIGEN_FUNCTORS_H const Eigen::CwiseUnaryOp< Eigen::internal::scalar_conjugate_op< typename Derived::Scalar >, const Derived > conj(const Eigen::ArrayBase< Derived > &x)
const Eigen::CwiseUnaryOp< Eigen::internal::scalar_arg_op< typename Derived::Scalar >, const Derived > arg(const Eigen::ArrayBase< Derived > &x)
Namespace containing all symbols from the Eigen library.
Definition: Core:287
const Eigen::CwiseUnaryOp< Eigen::internal::scalar_imag_op< typename Derived::Scalar >, const Derived > imag(const Eigen::ArrayBase< Derived > &x)
const Eigen::CwiseUnaryOp< Eigen::internal::scalar_real_op< typename Derived::Scalar >, const Derived > real(const Eigen::ArrayBase< Derived > &x)
Definition: Eigen_Colamd.h:50
const Eigen::CwiseUnaryOp< Eigen::internal::scalar_log10_op< typename Derived::Scalar >, const Derived > log10(const Eigen::ArrayBase< Derived > &x)