Actual source code: ex3.c


  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*/