From 35219cea68f8ac9f5624072492cab1e370155b1f Mon Sep 17 00:00:00 2001
From: Christoph Hertzberg <chtz@informatik.uni-bremen.de>
Date: Sat, 11 Jan 2020 14:57:22 +0100
Subject: Bug #1790: Make `areApprox` check `numext::isnan` instead of bitwise
 equality (NaNs don't have to be bitwise equal).

---
 test/packetmath_test_shared.h | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

(limited to 'test/packetmath_test_shared.h')

diff --git a/test/packetmath_test_shared.h b/test/packetmath_test_shared.h
index 046fd8104..de1a05508 100644
--- a/test/packetmath_test_shared.h
+++ b/test/packetmath_test_shared.h
@@ -100,8 +100,12 @@ template<typename Scalar> bool areApprox(const Scalar* a, const Scalar* b, int s
 {
   for (int i=0; i<size; ++i)
   {
-    if ((!internal::biteq(a[i],b[i])) && a[i]!=b[i] && !internal::isApprox(a[i],b[i]))
+    if (a[i]!=b[i] && !internal::isApprox(a[i],b[i]))
     {
+      if((numext::isnan)(a[i]) && (numext::isnan)(b[i]))
+      {
+        continue;
+      }
       std::cout << "ref: [" << Map<const Matrix<Scalar,1,Dynamic> >(a,size) << "]" << " != vec: [" << Map<const Matrix<Scalar,1,Dynamic> >(b,size) << "]\n";
       return false;
     }
-- 
cgit v1.2.3