diff options
Diffstat (limited to 'Modules/fpetestmodule.c')
| -rw-r--r-- | Modules/fpetestmodule.c | 228 |
1 files changed, 228 insertions, 0 deletions
diff --git a/Modules/fpetestmodule.c b/Modules/fpetestmodule.c new file mode 100644 index 0000000..2156afb --- /dev/null +++ b/Modules/fpetestmodule.c @@ -0,0 +1,228 @@ +/* + --------------------------------------------------------------------- + / Copyright (c) 1996. \ + | The Regents of the University of California. | + | All rights reserved. | + | | + | Permission to use, copy, modify, and distribute this software for | + | any purpose without fee is hereby granted, provided that this en- | + | tire notice is included in all copies of any software which is or | + | includes a copy or modification of this software and in all | + | copies of the supporting documentation for such software. | + | | + | This work was produced at the University of California, Lawrence | + | Livermore National Laboratory under contract no. W-7405-ENG-48 | + | between the U.S. Department of Energy and The Regents of the | + | University of California for the operation of UC LLNL. | + | | + | DISCLAIMER | + | | + | This software was prepared as an account of work sponsored by an | + | agency of the United States Government. Neither the United States | + | Government nor the University of California nor any of their em- | + | ployees, makes any warranty, express or implied, or assumes any | + | liability or responsibility for the accuracy, completeness, or | + | usefulness of any information, apparatus, product, or process | + | disclosed, or represents that its use would not infringe | + | privately-owned rights. Reference herein to any specific commer- | + | cial products, process, or service by trade name, trademark, | + | manufacturer, or otherwise, does not necessarily constitute or | + | imply its endorsement, recommendation, or favoring by the United | + | States Government or the University of California. The views and | + | opinions of authors expressed herein do not necessarily state or | + | reflect those of the United States Government or the University | + | of California, and shall not be used for advertising or product | + \ endorsement purposes. / + --------------------------------------------------------------------- +*/ + +/* + Floating point exception test module. + + */ + +#include "Python.h" + +static PyObject *fpe_error; +void initfpetest(void); +static PyObject *test(PyObject *self,PyObject *args); +static int db0(void); +static int overflow(void); +static int nest1(double); +static int nest2(double); +static double nest3(double); + +static PyMethodDef fpetest_methods[] = { + {"test", (PyCFunction) test, 1}, + {0,0} +}; + +static PyObject *test(PyObject *self,PyObject *args) +{ + int i = 0, r; + + fprintf(stderr,"Test trapping overflow\n"); + r = overflow(); + if(r){ + fprintf(stderr,"(Note: No exception was raised.)\n"); + PyErr_Clear(); + }else{ + fprintf(stderr,"Trapped:: "); + PyErr_Print(); + PyErr_Clear(); + } + i += r; + + fprintf(stderr,"Test trapping division by zero\n"); + r = db0(); + if(r){ + fprintf(stderr,"(Note: No exception was raised.)\n"); + PyErr_Clear(); + }else{ + fprintf(stderr,"Trapped:: "); + PyErr_Print(); + PyErr_Clear(); + } + i += r; + + fprintf(stderr,"Test nested protection zones, outer zone\n"); + r = nest1(0.0); + if(r){ + fprintf(stderr,"(Note: No exception was raised.)\n"); + PyErr_Clear(); + }else{ + fprintf(stderr,"Trapped:: "); + PyErr_Print(); + PyErr_Clear(); + } + i += r; + + fprintf(stderr,"Test nested protection zones, inner zone\n"); + fprintf(stderr,"(Note: Return will apparently come from outer zone.)\n"); + r = nest1(1.0); + if(r){ + fprintf(stderr,"(Note: No exception was raised.)\n"); + PyErr_Clear(); + }else{ + fprintf(stderr,"Trapped:: "); + PyErr_Print(); + PyErr_Clear(); + } + i += r; + + fprintf(stderr,"Test nested protection zones, trailing outer zone\n"); + r = nest1(2.0); + if(r){ + fprintf(stderr,"(Note: No exception was raised.)\n"); + PyErr_Clear(); + }else{ + fprintf(stderr,"Trapped:: "); + PyErr_Print(); + PyErr_Clear(); + } + i += r; + + fprintf(stderr,"Test nested function calls, prior error\n"); + r = nest2(0.0); + if(r){ + fprintf(stderr,"(Note: No exception was raised.)\n"); + PyErr_Clear(); + }else{ + fprintf(stderr,"Trapped:: "); + PyErr_Print(); + PyErr_Clear(); + } + i += r; + + fprintf(stderr,"Test nested function calls, interior error\n"); + r = nest2(1.0); + if(r){ + fprintf(stderr,"(Note: No exception was raised.)\n"); + PyErr_Clear(); + }else{ + fprintf(stderr,"Trapped:: "); + PyErr_Print(); + PyErr_Clear(); + } + i += r; + + fprintf(stderr,"Test nested function calls, trailing error\n"); + r = nest2(2.0); + if(r){ + fprintf(stderr,"(Note: No exception was raised.)\n"); + PyErr_Clear(); + }else{ + fprintf(stderr,"Trapped:: "); + PyErr_Print(); + PyErr_Clear(); + } + i += r; + + fprintf(stderr,"Number of tests failed: %d\n", i); + Py_INCREF (Py_None); + return Py_None; +} + +static int nest1(double x) +{ + double a = 1.0; + PyFPE_START_PROTECT("Division by zero, outer zone", return 0) + a = 1./x; + PyFPE_START_PROTECT("Division by zero, inner zone", return 0) + a = 1./(1. - x); + PyFPE_END_PROTECT + a = 1./(2. - x); + PyFPE_END_PROTECT + return(1); +} + +static int nest2(double x) +{ + double a = 1.0; + PyFPE_START_PROTECT("Division by zero, prior error", return 0) + a = 1./x; + a = nest3(x); + a = 1./(2. - x); + PyFPE_END_PROTECT + return(1); +} + +static double nest3(double x) +{ + double result; + PyFPE_START_PROTECT("Division by zero, nest3 error", return 0) + result = 1./(1. - x); + PyFPE_END_PROTECT + return result; +} + +static int db0(void) +{ + double a = 1.0; + PyFPE_START_PROTECT("Division by zero", return 0) + a = 1./(a - 1.); + PyFPE_END_PROTECT + return(1); +} + +static int overflow(void) +{ + double a, b = 1.e200; + PyFPE_START_PROTECT("Overflow", return 0) + a = b*b; + PyFPE_END_PROTECT + return(1); +} + +void initfpetest(void) +{ + PyObject *m, *d; + + m = Py_InitModule("fpetest", fpetest_methods); + d = PyModule_GetDict(m); + fpe_error = PyString_FromString("fpetest.error"); + PyDict_SetItemString(d, "error", fpe_error); + + if (PyErr_Occurred()) + Py_FatalError("Cannot initialize module fpetest"); +} |