Make test more robust: Check tolerance rather than ==

dealii · Jun 15, 2022 · 7633e80 · 7633e80
1 parent 0ab370f
commit 7633e80
Showing 1 changed file with 15 additions and 9 deletions.
diff --git a/tests/fe/fe_values_view_25b_nonsymmetric.cc b/tests/fe/fe_values_view_25b_nonsymmetric.cc
@@ -165,19 +165,21 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
                 dealii::Tensor<2, dim>::unrolled_to_component_indices(k);
               if (ind_k == nonzero_ind)
                 {
-                  AssertThrow(val_q[ind_k] == scalar_values[i_node][q],
+                  AssertThrow(std::abs(val_q[ind_k] -
+                                       scalar_values[i_node][q]) < 1e-12,
                               ExcInternalError());
-                  AssertThrow((grad_q[ind_k[0]][ind_k[1]] ==
-                               scalar_gradients[i_node][q]),
+                  AssertThrow((grad_q[ind_k[0]][ind_k[1]] -
+                               scalar_gradients[i_node][q])
+                                  .norm() < 1e-11,
                               ExcInternalError());
                 }
               else
                 {
-                  AssertThrow(val_q[ind_k] == 0.,
+                  AssertThrow(std::abs(val_q[ind_k]) < 1e-12,
                               ExcMessage(std::to_string(k) + " " +
                                          std::to_string(ind_k[0]) + " " +
                                          std::to_string(ind_k[1])));
-                  AssertThrow((grad_q[ind_k[0]][ind_k[1]] == Tensor<1, dim>()),
+                  AssertThrow(grad_q[ind_k[0]][ind_k[1]].norm() < 1e-11,
                               ExcInternalError());
                 }
             }
@@ -190,17 +192,21 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
               for (unsigned int j = 0; j < dim; ++j)
                 div_tmp += grad_q[i][j][j];
 
-              AssertThrow(div_tmp == div_q[i], ExcInternalError());
+              AssertThrow(std::abs(div_tmp - div_q[i]) < 1e-11,
+                          ExcInternalError());
             }
         }
     }
 
   // compare function values/gradients/divergences:
   for (unsigned int q = 0; q < quadrature.size(); ++q)
     {
-      AssertThrow(values_manual[q] == values[q], ExcInternalError());
-      AssertThrow(gradients_manual[q] == gradients[q], ExcInternalError());
-      AssertThrow(divergences_manual[q] == divergences[q], ExcInternalError());
+      AssertThrow((values_manual[q] - values[q]).norm() < 1e-12,
+                  ExcInternalError());
+      AssertThrow((gradients_manual[q] - gradients[q]).norm() < 1e-11,
+                  ExcInternalError());
+      AssertThrow((divergences_manual[q] - divergences[q]).norm() < 1e-11,
+                  ExcInternalError());
     }
 }