• Home
  • Line#
  • Scopes#
  • Navigate#
  • Raw
  • Download
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