devtools/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2026-04-10 11:34:33 +08:00
Code Issues Packages Projects Releases Wiki Activity

Add Neon psqrt<Packet2d> and pexp<Packet2d>

Browse Source
This commit is contained in:
Guoqiang QI
2020-09-08 09:04:03 +00:00
committed by David Tellenbach
parent 5272106826
commit 85428a3440
4 changed files with 59 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
Eigen/src/Core/arch/Default/GenericPacketMathFunctions.h
View File

@@ -39,6 +39,16 @@ pldexp_float(Packet a, Packet exponent)
return pmul(a, preinterpret<Packet>(plogical_shift_left<23>(ei)));
}
template<typename Packet> EIGEN_STRONG_INLINE Packet
pldexp_double(Packet a, Packet exponent)
{
typedef typename unpacket_traits<Packet>::integer_packet PacketI;
const Packet cst_1023 = pset1<Packet>(1023.0);
// return a * 2^exponent
PacketI ei = pcast<Packet,PacketI>(padd(exponent, cst_1023));
return pmul(a, preinterpret<Packet>(plogical_shift_left<52>(ei)));
}
// Natural logarithm
// Computes log(x) as log(2^e * m) = C*e + log(m), where the constant C =log(2)
// and m is in the range [sqrt(1/2),sqrt(2)). In this range, the logarithm can