1c4762a1bSJed Brown static char help[] = "Tests catching of floating point exceptions.\n\n";
2c4762a1bSJed Brown
3c4762a1bSJed Brown #include <petscsys.h>
4c4762a1bSJed Brown
CreateError(PetscReal x)53ba16761SJacob Faibussowitsch PetscErrorCode CreateError(PetscReal x)
6d71ae5a4SJacob Faibussowitsch {
7c4762a1bSJed Brown PetscFunctionBegin;
8c4762a1bSJed Brown x = 1.0 / x;
99566063dSJacob Faibussowitsch PetscCall(PetscPrintf(PETSC_COMM_SELF, "x = %g\n", (double)x));
103ba16761SJacob Faibussowitsch PetscFunctionReturn(PETSC_SUCCESS);
11c4762a1bSJed Brown }
12c4762a1bSJed Brown
main(int argc,char ** argv)13d71ae5a4SJacob Faibussowitsch int main(int argc, char **argv)
14d71ae5a4SJacob Faibussowitsch {
15327415f7SBarry Smith PetscFunctionBeginUser;
16c8025a54SPierre Jolivet PetscCall(PetscInitialize(&argc, &argv, NULL, help));
179566063dSJacob Faibussowitsch PetscCall(PetscPrintf(PETSC_COMM_SELF, "This is a contrived example to test floating pointing\n"));
189566063dSJacob Faibussowitsch PetscCall(PetscPrintf(PETSC_COMM_SELF, "It is not a true error.\n"));
199566063dSJacob Faibussowitsch PetscCall(PetscPrintf(PETSC_COMM_SELF, "Run with -fp_trap to catch the floating point error\n"));
209566063dSJacob Faibussowitsch PetscCall(CreateError(0.0));
21c4762a1bSJed Brown return 0;
22c4762a1bSJed Brown }
23c4762a1bSJed Brown
24c4762a1bSJed Brown /*
25c4762a1bSJed Brown
26c4762a1bSJed Brown Because this example may produce different output on different machines we filter out everything.
27c4762a1bSJed Brown This makes the test ineffective but currently we don't have a good way to know which machines should handle
28c4762a1bSJed Brown the floating point exceptions properly.
29c4762a1bSJed Brown
30c4762a1bSJed Brown */
31c4762a1bSJed Brown /*TEST
32c4762a1bSJed Brown
33c4762a1bSJed Brown test:
34c4762a1bSJed Brown args: -fp_trap -error_output_stdout
35c4762a1bSJed Brown filter: Error: true
36*3886731fSPierre Jolivet output_file: output/empty.out
37c4762a1bSJed Brown
38c4762a1bSJed Brown TEST*/
39