From: Timo Heister <timo.heister@gmail.com>
Date: Thu, 7 Dec 2023 14:29:43 +0000 (-0500)
Subject: make test more robust
X-Git-Tag: relicensing~157^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F16327%2Fhead;p=dealii.git

make test more robust

By avoiding a floating point equality comparison. This test fails for me
in release mode with march=native with a residual of 1e-14.
---

diff --git a/tests/mpi/affine_constraints_get_view_01.cc b/tests/mpi/affine_constraints_get_view_01.cc
index def3035761..763a8718f4 100644
--- a/tests/mpi/affine_constraints_get_view_01.cc
+++ b/tests/mpi/affine_constraints_get_view_01.cc
@@ -192,8 +192,15 @@ test()
     pressure_constraints.distribute(check_vector_2.block(1));
   }
 
-  // Now the assertion part: These vectors are the same:
-  Assert(check_vector_1 == check_vector_2, ExcInternalError());
+  // Finally check that these two vectors are the same. Note that
+  // when we compile with native optimizations we might have a slight
+  // difference in results:
+  {
+    TrilinosWrappers::MPI::BlockVector result = check_vector_1;
+    result -= check_vector_2;
+    Assert(result.l2_norm() / check_vector_1.l2_norm() < 1e-8,
+           ExcInternalError());
+  }
 
   deallog << "OK" << std::endl;
 }