Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 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; |
Guido van Rossum | 5de5420 | 1998-12-10 16:49:28 +0000 | [diff] [blame^] | 47 | DL_EXPORT(void) initfpetest(void); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 48 | static PyObject *test(PyObject *self,PyObject *args); |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 49 | static double db0(double); |
| 50 | static double overflow(double); |
| 51 | static double nest1(int, double); |
| 52 | static double nest2(int, double); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 53 | static double nest3(double); |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 54 | static void printerr(double); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 55 | |
| 56 | static PyMethodDef fpetest_methods[] = { |
| 57 | {"test", (PyCFunction) test, 1}, |
| 58 | {0,0} |
| 59 | }; |
| 60 | |
| 61 | static PyObject *test(PyObject *self,PyObject *args) |
| 62 | { |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 63 | double r; |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 64 | |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 65 | fprintf(stderr,"overflow"); |
| 66 | r = overflow(1.e160); |
| 67 | printerr(r); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 68 | |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 69 | fprintf(stderr,"\ndiv by 0"); |
| 70 | r = db0(0.0); |
| 71 | printerr(r); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 72 | |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 73 | fprintf(stderr,"\nnested outer"); |
| 74 | r = nest1(0, 0.0); |
| 75 | printerr(r); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 76 | |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 77 | fprintf(stderr,"\nnested inner"); |
| 78 | r = nest1(1, 1.0); |
| 79 | printerr(r); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 80 | |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 81 | fprintf(stderr,"\ntrailing outer"); |
| 82 | r = nest1(2, 2.0); |
| 83 | printerr(r); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 84 | |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 85 | fprintf(stderr,"\nnested prior"); |
| 86 | r = nest2(0, 0.0); |
| 87 | printerr(r); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 88 | |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 89 | fprintf(stderr,"\nnested interior"); |
| 90 | r = nest2(1, 1.0); |
| 91 | printerr(r); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 92 | |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 93 | fprintf(stderr,"\nnested trailing"); |
| 94 | r = nest2(2, 2.0); |
| 95 | printerr(r); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 96 | |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 97 | Py_INCREF (Py_None); |
| 98 | return Py_None; |
| 99 | } |
| 100 | |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 101 | static void printerr(double r) |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 102 | { |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 103 | if(r == 3.1416){ |
| 104 | fprintf(stderr,"\tPASS\n"); |
| 105 | PyErr_Print(); |
| 106 | }else{ |
| 107 | fprintf(stderr,"\tFAIL\n"); |
| 108 | } |
| 109 | PyErr_Clear(); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 110 | } |
| 111 | |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 112 | static double nest1(int i, double x) |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 113 | { |
| 114 | double a = 1.0; |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 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); |
Guido van Rossum | 45b8391 | 1997-03-14 04:32:50 +0000 | [diff] [blame] | 123 | PyFPE_END_PROTECT(a) |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 124 | }else if(i == 2){ |
| 125 | a = 1./(2. - x); |
| 126 | } |
Guido van Rossum | 45b8391 | 1997-03-14 04:32:50 +0000 | [diff] [blame] | 127 | PyFPE_END_PROTECT(a) |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 128 | |
| 129 | return a; |
| 130 | } |
| 131 | |
| 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 | } |
Guido van Rossum | 45b8391 | 1997-03-14 04:32:50 +0000 | [diff] [blame] | 143 | PyFPE_END_PROTECT(a) |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 144 | return a; |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 145 | } |
| 146 | |
| 147 | static double nest3(double x) |
| 148 | { |
| 149 | double result; |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 150 | /* This (following) message is never seen. */ |
| 151 | PyFPE_START_PROTECT("Division by zero, nest3 error", return 3.1416) |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 152 | result = 1./(1. - x); |
Guido van Rossum | 45b8391 | 1997-03-14 04:32:50 +0000 | [diff] [blame] | 153 | PyFPE_END_PROTECT(result) |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 154 | return result; |
| 155 | } |
| 156 | |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 157 | static double db0(double x) |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 158 | { |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 159 | double a; |
| 160 | PyFPE_START_PROTECT("Division by zero", return 3.1416) |
| 161 | a = 1./x; |
Guido van Rossum | 45b8391 | 1997-03-14 04:32:50 +0000 | [diff] [blame] | 162 | PyFPE_END_PROTECT(a) |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 163 | return a; |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 164 | } |
| 165 | |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 166 | static double overflow(double b) |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 167 | { |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 168 | double a; |
| 169 | PyFPE_START_PROTECT("Overflow", return 3.1416) |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 170 | a = b*b; |
Guido van Rossum | 45b8391 | 1997-03-14 04:32:50 +0000 | [diff] [blame] | 171 | PyFPE_END_PROTECT(a) |
Guido van Rossum | 05bc4af | 1997-02-21 01:22:47 +0000 | [diff] [blame] | 172 | return a; |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 173 | } |
| 174 | |
Guido van Rossum | 3886bb6 | 1998-12-04 18:50:17 +0000 | [diff] [blame] | 175 | DL_EXPORT(void) initfpetest(void) |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 176 | { |
| 177 | PyObject *m, *d; |
| 178 | |
| 179 | m = Py_InitModule("fpetest", fpetest_methods); |
| 180 | d = PyModule_GetDict(m); |
Guido van Rossum | 0cb96de | 1997-10-01 04:29:29 +0000 | [diff] [blame] | 181 | fpe_error = PyErr_NewException("fpetest.error", NULL, NULL); |
| 182 | if (fpe_error != NULL) |
| 183 | PyDict_SetItemString(d, "error", fpe_error); |
Guido van Rossum | 52fa3a6 | 1997-02-14 22:59:58 +0000 | [diff] [blame] | 184 | } |