2: static char help[] = "Tests catching of floating point exceptions.\n\n";
4: #include <petscsys.h>
6: int CreateError(PetscReal x)
7: {
8: x = 1.0 / x;
9: PetscPrintf(PETSC_COMM_SELF, "x = %g\n", (double)x);
10: return 0;
11: }
13: int main(int argc, char **argv)
14: {
16: PetscInitialize(&argc, &argv, (char *)0, help);
17: PetscPrintf(PETSC_COMM_SELF, "This is a contrived example to test floating pointing\n");
18: PetscPrintf(PETSC_COMM_SELF, "It is not a true error.\n");
19: PetscPrintf(PETSC_COMM_SELF, "Run with -fp_trap to catch the floating point error\n");
20: CreateError(0.0);
21: return 0;
22: }
24: /*
26: Because this example may produce different output on different machines we filter out everything.
27: This makes the test ineffective but currently we don't have a good way to know which machines should handle
28: the floating point exceptions properly.
30: */
31: /*TEST
33: test:
34: args: -fp_trap -error_output_stdout
35: filter: Error: true
37: TEST*/