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