1 /*
2 ---------------------------------------------------------------------
3 / Copyright (c) 1996. \
4 | The Regents of the University of California. |
5 | All rights reserved. |
6 | |
7 | Permission to use, copy, modify, and distribute this software for |
8 | any purpose without fee is hereby granted, provided that this en- |
9 | tire notice is included in all copies of any software which is or |
10 | includes a copy or modification of this software and in all |
11 | copies of the supporting documentation for such software. |
12 | |
13 | This work was produced at the University of California, Lawrence |
14 | Livermore National Laboratory under contract no. W-7405-ENG-48 |
15 | between the U.S. Department of Energy and The Regents of the |
16 | University of California for the operation of UC LLNL. |
17 | |
18 | DISCLAIMER |
19 | |
20 | This software was prepared as an account of work sponsored by an |
21 | agency of the United States Government. Neither the United States |
22 | Government nor the University of California nor any of their em- |
23 | ployees, makes any warranty, express or implied, or assumes any |
24 | liability or responsibility for the accuracy, completeness, or |
25 | usefulness of any information, apparatus, product, or process |
26 | disclosed, or represents that its use would not infringe |
27 | privately-owned rights. Reference herein to any specific commer- |
28 | cial products, process, or service by trade name, trademark, |
29 | manufacturer, or otherwise, does not necessarily constitute or |
30 | imply its endorsement, recommendation, or favoring by the United |
31 | States Government or the University of California. The views and |
32 | opinions of authors expressed herein do not necessarily state or |
33 | reflect those of the United States Government or the University |
34 | of California, and shall not be used for advertising or product |
35 \ endorsement purposes. /
36 ---------------------------------------------------------------------
37 */
38
39 /*
40 Floating point exception test module.
41
42 */
43
44 #include "Python.h"
45
46 static PyObject *fpe_error;
47 PyMODINIT_FUNC initfpetest(void);
48 static PyObject *test(PyObject *self,PyObject *args);
49 static double db0(double);
50 static double overflow(double);
51 static double nest1(int, double);
52 static double nest2(int, double);
53 static double nest3(double);
54 static void printerr(double);
55
56 static PyMethodDef fpetest_methods[] = {
57 {"test", (PyCFunction) test, METH_VARARGS},
58 {0,0}
59 };
60
test(PyObject * self,PyObject * args)61 static PyObject *test(PyObject *self,PyObject *args)
62 {
63 double r;
64
65 fprintf(stderr,"overflow");
66 r = overflow(1.e160);
67 printerr(r);
68
69 fprintf(stderr,"\ndiv by 0");
70 r = db0(0.0);
71 printerr(r);
72
73 fprintf(stderr,"\nnested outer");
74 r = nest1(0, 0.0);
75 printerr(r);
76
77 fprintf(stderr,"\nnested inner");
78 r = nest1(1, 1.0);
79 printerr(r);
80
81 fprintf(stderr,"\ntrailing outer");
82 r = nest1(2, 2.0);
83 printerr(r);
84
85 fprintf(stderr,"\nnested prior");
86 r = nest2(0, 0.0);
87 printerr(r);
88
89 fprintf(stderr,"\nnested interior");
90 r = nest2(1, 1.0);
91 printerr(r);
92
93 fprintf(stderr,"\nnested trailing");
94 r = nest2(2, 2.0);
95 printerr(r);
96
97 Py_INCREF (Py_None);
98 return Py_None;
99 }
100
printerr(double r)101 static void printerr(double r)
102 {
103 if(r == 3.1416){
104 fprintf(stderr,"\tPASS\n");
105 PyErr_Print();
106 }else{
107 fprintf(stderr,"\tFAIL\n");
108 }
109 PyErr_Clear();
110 }
111
nest1(int i,double x)112 static double nest1(int i, double x)
113 {
114 double a = 1.0;
115
116 PyFPE_START_PROTECT("Division by zero, outer zone", return 3.1416)
117 if(i == 0){
118 a = 1./x;
119 }else if(i == 1){
120 /* This (following) message is never seen. */
121 PyFPE_START_PROTECT("Division by zero, inner zone", return 3.1416)
122 a = 1./(1. - x);
123 PyFPE_END_PROTECT(a)
124 }else if(i == 2){
125 a = 1./(2. - x);
126 }
127 PyFPE_END_PROTECT(a)
128
129 return a;
130 }
131
nest2(int i,double x)132 static double nest2(int i, double x)
133 {
134 double a = 1.0;
135 PyFPE_START_PROTECT("Division by zero, prior error", return 3.1416)
136 if(i == 0){
137 a = 1./x;
138 }else if(i == 1){
139 a = nest3(x);
140 }else if(i == 2){
141 a = 1./(2. - x);
142 }
143 PyFPE_END_PROTECT(a)
144 return a;
145 }
146
nest3(double x)147 static double nest3(double x)
148 {
149 double result;
150 /* This (following) message is never seen. */
151 PyFPE_START_PROTECT("Division by zero, nest3 error", return 3.1416)
152 result = 1./(1. - x);
153 PyFPE_END_PROTECT(result)
154 return result;
155 }
156
db0(double x)157 static double db0(double x)
158 {
159 double a;
160 PyFPE_START_PROTECT("Division by zero", return 3.1416)
161 a = 1./x;
162 PyFPE_END_PROTECT(a)
163 return a;
164 }
165
overflow(double b)166 static double overflow(double b)
167 {
168 double a;
169 PyFPE_START_PROTECT("Overflow", return 3.1416)
170 a = b*b;
171 PyFPE_END_PROTECT(a)
172 return a;
173 }
174
initfpetest(void)175 PyMODINIT_FUNC initfpetest(void)
176 {
177 PyObject *m, *d;
178
179 m = Py_InitModule("fpetest", fpetest_methods);
180 if (m == NULL)
181 return;
182 d = PyModule_GetDict(m);
183 fpe_error = PyErr_NewException("fpetest.error", NULL, NULL);
184 if (fpe_error != NULL)
185 PyDict_SetItemString(d, "error", fpe_error);
186 }
187