grid: Check floating point exceptions in grid_collocate_replay

cp2k · Sep 19, 2020 · 4ce1476 · 4ce1476
1 parent 06ab5f4
commit 4ce1476
Show file tree

Hide file tree

Showing 3 changed files with 916 additions and 1 deletion.
diff --git a/src/grid/grid_collocate_replay.c b/src/grid/grid_collocate_replay.c
@@ -8,6 +8,7 @@
 #define _XOPEN_SOURCE 700 /* Enable POSIX 2008/13 */
 
 #include <assert.h>
+#include <fenv.h>
 #include <float.h>
 #include <limits.h>
 #include <math.h>
@@ -453,6 +454,20 @@ double grid_collocate_replay(const char *filename, const int cycles,
   free(grid_ref);
   free(grid_test);
 
+  // Check floating point exceptions.
+  if (fetestexcept(FE_INVALID) != 0) {
+    fprintf(stderr, "Error: Floating point exception FE_INVALID.\n");
+    exit(1);
+  }
+  if (fetestexcept(FE_DIVBYZERO) != 0) {
+    fprintf(stderr, "Error: Floating point exception FE_DIVBYZERO.\n");
+    exit(1);
+  }
+  if (fetestexcept(FE_OVERFLOW) != 0) {
+    fprintf(stderr, "Error: Floating point exception FE_OVERFLOW.\n");
+    exit(1);
+  }
+
   return max_diff;
 }
 

diff --git a/src/grid/grid_collocate_unittest.c b/src/grid/grid_collocate_unittest.c
@@ -65,12 +65,13 @@ int main(int argc, char *argv[]) {
   errors += run_test(argv[1], "collocate_ortho_density_l2200.task");
   errors += run_test(argv[1], "collocate_ortho_density_l3300.task");
   errors += run_test(argv[1], "collocate_ortho_density_l3333.task");
+  errors += run_test(argv[1], "collocate_ortho_non_periodic.task");
   errors += run_test(argv[1], "collocate_ortho_tau.task");
   errors += run_test(argv[1], "collocate_general_density.task");
   errors += run_test(argv[1], "collocate_general_tau.task");
   errors += run_test(argv[1], "collocate_general_subpatch0.task");
   errors += run_test(argv[1], "collocate_general_subpatch16.task");
-  errors += run_test(argv[1], "collocate_ortho_non_periodic.task");
+  errors += run_test(argv[1], "collocate_general_overflow.task");
 
   grid_library_print_stats(&mpi_sum_func, 0, &print_func, 0);
   grid_library_finalize();