From 9aba527405b40132a308f5f782dacadf6ef50acd Mon Sep 17 00:00:00 2001 From: Rasmus Munk Larsen Date: Tue, 27 Aug 2019 15:35:29 -0700 Subject: Revert changes to std_falback::log1p that broke handling of arguments less than -1. Fix packet op accordingly. --- Eigen/src/Core/arch/Default/GenericPacketMathFunctions.h | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) (limited to 'Eigen/src/Core/arch/Default/GenericPacketMathFunctions.h') diff --git a/Eigen/src/Core/arch/Default/GenericPacketMathFunctions.h b/Eigen/src/Core/arch/Default/GenericPacketMathFunctions.h index 640aae05a..505a0eec8 100644 --- a/Eigen/src/Core/arch/Default/GenericPacketMathFunctions.h +++ b/Eigen/src/Core/arch/Default/GenericPacketMathFunctions.h @@ -137,10 +137,8 @@ Packet generic_plog1p(const Packet& x) Packet xp1 = padd(x, one); Packet small_mask = pcmp_eq(xp1, one); Packet log1 = plog(xp1); - // Add a check to handle x == +inf. - Packet pos_inf_mask = pcmp_eq(x, log1); Packet log_large = pmul(x, pdiv(log1, psub(xp1, one))); - return pselect(por(small_mask, pos_inf_mask), x, log_large); + return pselect(small_mask, x, log_large); } /** \internal \returns exp(x)-1 computed using W. Kahan's formula. -- cgit v1.2.3