GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-09-13 01:43:13 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fix generic ceil for SSE2.

Browse Source
This commit is contained in:
Antonio Sánchez 2024-09-14 01:31:21 +00:00 committed by Charles Schlosser
parent 84282c42fc
commit 132f281f50
1 changed files with 3 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -2598,11 +2598,14 @@ template <typename Packet>
EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE Packet generic_ceil(const Packet& a) { EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE Packet generic_ceil(const Packet& a) {
using Scalar = typename unpacket_traits<Packet>::type; using Scalar = typename unpacket_traits<Packet>::type;
const Packet cst_1 = pset1<Packet>(Scalar(1)); const Packet cst_1 = pset1<Packet>(Scalar(1));
const Packet sign_mask = pset1<Packet>(static_cast<Scalar>(-0.0));
Packet rint_a = generic_rint(a); Packet rint_a = generic_rint(a);
// if rint(a) < a, then rint(a) == floor(a) // if rint(a) < a, then rint(a) == floor(a)
Packet mask = pcmp_lt(rint_a, a); Packet mask = pcmp_lt(rint_a, a);
Packet offset = pand(cst_1, mask); Packet offset = pand(cst_1, mask);
Packet result = padd(rint_a, offset); Packet result = padd(rint_a, offset);
// Signed zero must remain signed (e.g. ceil(-0.02) == -0).
result = por(result, pand(sign_mask, a));
return result; return result;
} }