summaryrefslogtreecommitdiffstats
path: root/Modules/fpetestmodule.c
diff options
context:
space:
mode:
Diffstat (limited to 'Modules/fpetestmodule.c')
-rw-r--r--Modules/fpetestmodule.c228
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");
+}