mirror/eigen
2
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-03-07 18:27:40 +08:00
Code Issues Packages Projects Releases Wiki Activity

make the ei_p* math functions overloads instead of template

Browse Source 
specializations
This commit is contained in:
Gael Guennebaud 2009-04-22 21:35:50 +00:00
parent e0904a70ce
commit c7bb7436f9
1 changed files with 7 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
Eigen/src/Core/arch/SSE/MathFunctions.h
View File

@ -30,7 +30,7 @@
#ifndef EIGEN_MATH_FUNCTIONS_SSE_H
#define EIGEN_MATH_FUNCTIONS_SSE_H
template<> EIGEN_DONT_INLINE Packet4f ei_plog(Packet4f x)
static EIGEN_DONT_INLINE Packet4f ei_plog(Packet4f x)
{
_EIGEN_DECLARE_CONST_Packet4f(1 , 1.0f);
_EIGEN_DECLARE_CONST_Packet4f(half, 0.5f);
@ -108,7 +108,7 @@ template<> EIGEN_DONT_INLINE Packet4f ei_plog(Packet4f x)
return _mm_or_ps(x, invalid_mask); // negative arg will be NAN
}
template<> EIGEN_DONT_INLINE Packet4f ei_pexp(Packet4f x)
static EIGEN_DONT_INLINE Packet4f ei_pexp(Packet4f x)
{
_EIGEN_DECLARE_CONST_Packet4f(1 , 1.0f);
_EIGEN_DECLARE_CONST_Packet4f(half, 0.5f);
@ -181,7 +181,7 @@ template<> EIGEN_DONT_INLINE Packet4f ei_pexp(Packet4f x)
surprising but correct result.
*/
template<> EIGEN_DONT_INLINE Packet4f ei_psin(Packet4f x)
static EIGEN_DONT_INLINE Packet4f ei_psin(Packet4f x)
{
_EIGEN_DECLARE_CONST_Packet4f(1 , 1.0f);
_EIGEN_DECLARE_CONST_Packet4f(half, 0.5f);
@ -211,6 +211,8 @@ template<> EIGEN_DONT_INLINE Packet4f ei_psin(Packet4f x)
/* take the absolute value */
x = ei_pabs(x);
/* take the modulo */
/* extract the sign bit (upper one) */
sign_bit = _mm_and_ps(sign_bit, ei_p4f_sign_mask);
@ -278,7 +280,7 @@ template<> EIGEN_DONT_INLINE Packet4f ei_psin(Packet4f x)
}
/* almost the same as ei_psin */
template<> Packet4f ei_pcos(Packet4f x)
static Packet4f ei_pcos(Packet4f x)
{
_EIGEN_DECLARE_CONST_Packet4f(1 , 1.0f);
_EIGEN_DECLARE_CONST_Packet4f(half, 0.5f);
@ -363,7 +365,7 @@ template<> Packet4f ei_pcos(Packet4f x)
return _mm_xor_ps(y, sign_bit);
}
template<> Packet4f ei_psqrt(Packet4f _x)
static Packet4f ei_psqrt(Packet4f _x)
{
Packet4f half = ei_pmul(_x, ei_pset1(.5f));
Packet4f x = _mm_rsqrt_ps(_x);