1c4762a1bSJed Brown static char help[] = "Tests catching of floating point exceptions.\n\n"; 2c4762a1bSJed Brown 3c4762a1bSJed Brown #include <petscsys.h> 4c4762a1bSJed Brown 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 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