/* | |
--------------------------------------------------------------------- | |
/ 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; | |
PyMODINIT_FUNC initfpetest(void); | |
static PyObject *test(PyObject *self,PyObject *args); | |
static double db0(double); | |
static double overflow(double); | |
static double nest1(int, double); | |
static double nest2(int, double); | |
static double nest3(double); | |
static void printerr(double); | |
static PyMethodDef fpetest_methods[] = { | |
{"test", (PyCFunction) test, METH_VARARGS}, | |
{0,0} | |
}; | |
static PyObject *test(PyObject *self,PyObject *args) | |
{ | |
double r; | |
fprintf(stderr,"overflow"); | |
r = overflow(1.e160); | |
printerr(r); | |
fprintf(stderr,"\ndiv by 0"); | |
r = db0(0.0); | |
printerr(r); | |
fprintf(stderr,"\nnested outer"); | |
r = nest1(0, 0.0); | |
printerr(r); | |
fprintf(stderr,"\nnested inner"); | |
r = nest1(1, 1.0); | |
printerr(r); | |
fprintf(stderr,"\ntrailing outer"); | |
r = nest1(2, 2.0); | |
printerr(r); | |
fprintf(stderr,"\nnested prior"); | |
r = nest2(0, 0.0); | |
printerr(r); | |
fprintf(stderr,"\nnested interior"); | |
r = nest2(1, 1.0); | |
printerr(r); | |
fprintf(stderr,"\nnested trailing"); | |
r = nest2(2, 2.0); | |
printerr(r); | |
Py_INCREF (Py_None); | |
return Py_None; | |
} | |
static void printerr(double r) | |
{ | |
if(r == 3.1416){ | |
fprintf(stderr,"\tPASS\n"); | |
PyErr_Print(); | |
}else{ | |
fprintf(stderr,"\tFAIL\n"); | |
} | |
PyErr_Clear(); | |
} | |
static double nest1(int i, double x) | |
{ | |
double a = 1.0; | |
PyFPE_START_PROTECT("Division by zero, outer zone", return 3.1416) | |
if(i == 0){ | |
a = 1./x; | |
}else if(i == 1){ | |
/* This (following) message is never seen. */ | |
PyFPE_START_PROTECT("Division by zero, inner zone", return 3.1416) | |
a = 1./(1. - x); | |
PyFPE_END_PROTECT(a) | |
}else if(i == 2){ | |
a = 1./(2. - x); | |
} | |
PyFPE_END_PROTECT(a) | |
return a; | |
} | |
static double nest2(int i, double x) | |
{ | |
double a = 1.0; | |
PyFPE_START_PROTECT("Division by zero, prior error", return 3.1416) | |
if(i == 0){ | |
a = 1./x; | |
}else if(i == 1){ | |
a = nest3(x); | |
}else if(i == 2){ | |
a = 1./(2. - x); | |
} | |
PyFPE_END_PROTECT(a) | |
return a; | |
} | |
static double nest3(double x) | |
{ | |
double result; | |
/* This (following) message is never seen. */ | |
PyFPE_START_PROTECT("Division by zero, nest3 error", return 3.1416) | |
result = 1./(1. - x); | |
PyFPE_END_PROTECT(result) | |
return result; | |
} | |
static double db0(double x) | |
{ | |
double a; | |
PyFPE_START_PROTECT("Division by zero", return 3.1416) | |
a = 1./x; | |
PyFPE_END_PROTECT(a) | |
return a; | |
} | |
static double overflow(double b) | |
{ | |
double a; | |
PyFPE_START_PROTECT("Overflow", return 3.1416) | |
a = b*b; | |
PyFPE_END_PROTECT(a) | |
return a; | |
} | |
PyMODINIT_FUNC initfpetest(void) | |
{ | |
PyObject *m, *d; | |
m = Py_InitModule("fpetest", fpetest_methods); | |
if (m == NULL) | |
return; | |
d = PyModule_GetDict(m); | |
fpe_error = PyErr_NewException("fpetest.error", NULL, NULL); | |
if (fpe_error != NULL) | |
PyDict_SetItemString(d, "error", fpe_error); | |
} |