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