User | Revision | Line number | New contents of line |
ykuroda |
0:13a5d365ba16
|
1
|
// This file is part of Eigen, a lightweight C++ template library
|
ykuroda |
0:13a5d365ba16
|
2
|
// for linear algebra.
|
ykuroda |
0:13a5d365ba16
|
3
|
//
|
ykuroda |
0:13a5d365ba16
|
4
|
// Copyright (C) 2010-2012 Gael Guennebaud <gael.guennebaud@inria.fr>
|
ykuroda |
0:13a5d365ba16
|
5
|
// Copyright (C) 2010 Benoit Jacob <jacob.benoit.1@gmail.com>
|
ykuroda |
0:13a5d365ba16
|
6
|
//
|
ykuroda |
0:13a5d365ba16
|
7
|
// This Source Code Form is subject to the terms of the Mozilla
|
ykuroda |
0:13a5d365ba16
|
8
|
// Public License v. 2.0. If a copy of the MPL was not distributed
|
ykuroda |
0:13a5d365ba16
|
9
|
// with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
|
ykuroda |
0:13a5d365ba16
|
10
|
|
ykuroda |
0:13a5d365ba16
|
11
|
#ifndef EIGEN_GLOBAL_FUNCTIONS_H
|
ykuroda |
0:13a5d365ba16
|
12
|
#define EIGEN_GLOBAL_FUNCTIONS_H
|
ykuroda |
0:13a5d365ba16
|
13
|
|
ykuroda |
0:13a5d365ba16
|
14
|
#define EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(NAME,FUNCTOR) \
|
ykuroda |
0:13a5d365ba16
|
15
|
template<typename Derived> \
|
ykuroda |
0:13a5d365ba16
|
16
|
inline const Eigen::CwiseUnaryOp<Eigen::internal::FUNCTOR<typename Derived::Scalar>, const Derived> \
|
ykuroda |
0:13a5d365ba16
|
17
|
NAME(const Eigen::ArrayBase<Derived>& x) { \
|
ykuroda |
0:13a5d365ba16
|
18
|
return x.derived(); \
|
ykuroda |
0:13a5d365ba16
|
19
|
}
|
ykuroda |
0:13a5d365ba16
|
20
|
|
ykuroda |
0:13a5d365ba16
|
21
|
#define EIGEN_ARRAY_DECLARE_GLOBAL_EIGEN_UNARY(NAME,FUNCTOR) \
|
ykuroda |
0:13a5d365ba16
|
22
|
\
|
ykuroda |
0:13a5d365ba16
|
23
|
template<typename Derived> \
|
ykuroda |
0:13a5d365ba16
|
24
|
struct NAME##_retval<ArrayBase<Derived> > \
|
ykuroda |
0:13a5d365ba16
|
25
|
{ \
|
ykuroda |
0:13a5d365ba16
|
26
|
typedef const Eigen::CwiseUnaryOp<Eigen::internal::FUNCTOR<typename Derived::Scalar>, const Derived> type; \
|
ykuroda |
0:13a5d365ba16
|
27
|
}; \
|
ykuroda |
0:13a5d365ba16
|
28
|
template<typename Derived> \
|
ykuroda |
0:13a5d365ba16
|
29
|
struct NAME##_impl<ArrayBase<Derived> > \
|
ykuroda |
0:13a5d365ba16
|
30
|
{ \
|
ykuroda |
0:13a5d365ba16
|
31
|
static inline typename NAME##_retval<ArrayBase<Derived> >::type run(const Eigen::ArrayBase<Derived>& x) \
|
ykuroda |
0:13a5d365ba16
|
32
|
{ \
|
ykuroda |
0:13a5d365ba16
|
33
|
return x.derived(); \
|
ykuroda |
0:13a5d365ba16
|
34
|
} \
|
ykuroda |
0:13a5d365ba16
|
35
|
};
|
ykuroda |
0:13a5d365ba16
|
36
|
|
ykuroda |
0:13a5d365ba16
|
37
|
|
ykuroda |
0:13a5d365ba16
|
38
|
namespace Eigen
|
ykuroda |
0:13a5d365ba16
|
39
|
{
|
ykuroda |
0:13a5d365ba16
|
40
|
EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(real,scalar_real_op)
|
ykuroda |
0:13a5d365ba16
|
41
|
EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(imag,scalar_imag_op)
|
ykuroda |
0:13a5d365ba16
|
42
|
EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(conj,scalar_conjugate_op)
|
ykuroda |
0:13a5d365ba16
|
43
|
EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(sin,scalar_sin_op)
|
ykuroda |
0:13a5d365ba16
|
44
|
EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(cos,scalar_cos_op)
|
ykuroda |
0:13a5d365ba16
|
45
|
EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(asin,scalar_asin_op)
|
ykuroda |
0:13a5d365ba16
|
46
|
EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(acos,scalar_acos_op)
|
ykuroda |
0:13a5d365ba16
|
47
|
EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(tan,scalar_tan_op)
|
ykuroda |
0:13a5d365ba16
|
48
|
EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(exp,scalar_exp_op)
|
ykuroda |
0:13a5d365ba16
|
49
|
EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(log,scalar_log_op)
|
ykuroda |
0:13a5d365ba16
|
50
|
EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(abs,scalar_abs_op)
|
ykuroda |
0:13a5d365ba16
|
51
|
EIGEN_ARRAY_DECLARE_GLOBAL_UNARY(sqrt,scalar_sqrt_op)
|
ykuroda |
0:13a5d365ba16
|
52
|
|
ykuroda |
0:13a5d365ba16
|
53
|
template<typename Derived>
|
ykuroda |
0:13a5d365ba16
|
54
|
inline const Eigen::CwiseUnaryOp<Eigen::internal::scalar_pow_op<typename Derived::Scalar>, const Derived>
|
ykuroda |
0:13a5d365ba16
|
55
|
pow(const Eigen::ArrayBase<Derived>& x, const typename Derived::Scalar& exponent) {
|
ykuroda |
0:13a5d365ba16
|
56
|
return x.derived().pow(exponent);
|
ykuroda |
0:13a5d365ba16
|
57
|
}
|
ykuroda |
0:13a5d365ba16
|
58
|
|
ykuroda |
0:13a5d365ba16
|
59
|
template<typename Derived>
|
ykuroda |
0:13a5d365ba16
|
60
|
inline const Eigen::CwiseBinaryOp<Eigen::internal::scalar_binary_pow_op<typename Derived::Scalar, typename Derived::Scalar>, const Derived, const Derived>
|
ykuroda |
0:13a5d365ba16
|
61
|
pow(const Eigen::ArrayBase<Derived>& x, const Eigen::ArrayBase<Derived>& exponents)
|
ykuroda |
0:13a5d365ba16
|
62
|
{
|
ykuroda |
0:13a5d365ba16
|
63
|
return Eigen::CwiseBinaryOp<Eigen::internal::scalar_binary_pow_op<typename Derived::Scalar, typename Derived::Scalar>, const Derived, const Derived>(
|
ykuroda |
0:13a5d365ba16
|
64
|
x.derived(),
|
ykuroda |
0:13a5d365ba16
|
65
|
exponents.derived()
|
ykuroda |
0:13a5d365ba16
|
66
|
);
|
ykuroda |
0:13a5d365ba16
|
67
|
}
|
ykuroda |
0:13a5d365ba16
|
68
|
|
ykuroda |
0:13a5d365ba16
|
69
|
/**
|
ykuroda |
0:13a5d365ba16
|
70
|
* \brief Component-wise division of a scalar by array elements.
|
ykuroda |
0:13a5d365ba16
|
71
|
**/
|
ykuroda |
0:13a5d365ba16
|
72
|
template <typename Derived>
|
ykuroda |
0:13a5d365ba16
|
73
|
inline const Eigen::CwiseUnaryOp<Eigen::internal::scalar_inverse_mult_op<typename Derived::Scalar>, const Derived>
|
ykuroda |
0:13a5d365ba16
|
74
|
operator/(const typename Derived::Scalar& s, const Eigen::ArrayBase<Derived>& a)
|
ykuroda |
0:13a5d365ba16
|
75
|
{
|
ykuroda |
0:13a5d365ba16
|
76
|
return Eigen::CwiseUnaryOp<Eigen::internal::scalar_inverse_mult_op<typename Derived::Scalar>, const Derived>(
|
ykuroda |
0:13a5d365ba16
|
77
|
a.derived(),
|
ykuroda |
0:13a5d365ba16
|
78
|
Eigen::internal::scalar_inverse_mult_op<typename Derived::Scalar>(s)
|
ykuroda |
0:13a5d365ba16
|
79
|
);
|
ykuroda |
0:13a5d365ba16
|
80
|
}
|
ykuroda |
0:13a5d365ba16
|
81
|
|
ykuroda |
0:13a5d365ba16
|
82
|
namespace internal
|
ykuroda |
0:13a5d365ba16
|
83
|
{
|
ykuroda |
0:13a5d365ba16
|
84
|
EIGEN_ARRAY_DECLARE_GLOBAL_EIGEN_UNARY(real,scalar_real_op)
|
ykuroda |
0:13a5d365ba16
|
85
|
EIGEN_ARRAY_DECLARE_GLOBAL_EIGEN_UNARY(imag,scalar_imag_op)
|
ykuroda |
0:13a5d365ba16
|
86
|
EIGEN_ARRAY_DECLARE_GLOBAL_EIGEN_UNARY(abs2,scalar_abs2_op)
|
ykuroda |
0:13a5d365ba16
|
87
|
}
|
ykuroda |
0:13a5d365ba16
|
88
|
}
|
ykuroda |
0:13a5d365ba16
|
89
|
|
ykuroda |
0:13a5d365ba16
|
90
|
// TODO: cleanly disable those functions that are not supported on Array (numext::real_ref, internal::random, internal::isApprox...)
|
ykuroda |
0:13a5d365ba16
|
91
|
|
ykuroda |
0:13a5d365ba16
|
92
|
#endif // EIGEN_GLOBAL_FUNCTIONS_H |