fpetestmodule.c 5.5 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46
/*
     ---------------------------------------------------------------------  
    /                       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;
47
PyMODINIT_FUNC initfpetest(void);
48
static PyObject *test(PyObject *self,PyObject *args);
49 50 51 52
static double db0(double);
static double overflow(double);
static double nest1(int, double);
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 59 60 61 62
    {0,0}
};

static PyObject *test(PyObject *self,PyObject *args)
{
63
    double r;
64

65 66 67
    fprintf(stderr,"overflow");
    r = overflow(1.e160);
    printerr(r);
68

69 70 71
    fprintf(stderr,"\ndiv by 0");
    r = db0(0.0);
    printerr(r);
72

73 74 75
    fprintf(stderr,"\nnested outer");
    r = nest1(0, 0.0);
    printerr(r);
76

77 78 79
    fprintf(stderr,"\nnested inner");
    r = nest1(1, 1.0);
    printerr(r);
80

81 82 83
    fprintf(stderr,"\ntrailing outer");
    r = nest1(2, 2.0);
    printerr(r);
84

85 86 87
    fprintf(stderr,"\nnested prior");
    r = nest2(0, 0.0);
    printerr(r);
88

89 90 91 92 93 94 95
    fprintf(stderr,"\nnested interior");
    r = nest2(1, 1.0);
    printerr(r);

    fprintf(stderr,"\nnested trailing");
    r = nest2(2, 2.0);
    printerr(r);
96 97 98 99 100

    Py_INCREF (Py_None);
    return Py_None;
}

101 102 103 104 105 106 107 108 109 110 111 112
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)
113 114
{
  double a = 1.0;
115 116 117 118 119 120 121

  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)
122
    a = 1./(1. - x);
123
    PyFPE_END_PROTECT(a)
124 125 126
  }else if(i == 2){
    a = 1./(2. - x);
  }
127
  PyFPE_END_PROTECT(a)
128 129

  return a;
130 131
}

132
static double nest2(int i, double x)
133 134
{
  double a = 1.0;
135 136 137 138 139 140 141 142
  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);
  }
143
  PyFPE_END_PROTECT(a)
144
  return a;
145 146 147 148 149
}

static double nest3(double x)
{
  double result;
150 151
  /* This (following) message is never seen. */
  PyFPE_START_PROTECT("Division by zero, nest3 error", return 3.1416)
152
  result = 1./(1. - x);
153
  PyFPE_END_PROTECT(result)
154 155 156
  return result;
}

157
static double db0(double x)
158
{
159 160 161
  double a;
  PyFPE_START_PROTECT("Division by zero", return 3.1416)
  a = 1./x;
162
  PyFPE_END_PROTECT(a)
163
  return a;
164 165
}

166
static double overflow(double b)
167
{
168 169
  double a;
  PyFPE_START_PROTECT("Overflow", return 3.1416)
170
  a = b*b;
171
  PyFPE_END_PROTECT(a)
172
  return a;
173 174
}

175
PyMODINIT_FUNC initfpetest(void)
176 177 178 179 180
{
    PyObject *m, *d;

    m = Py_InitModule("fpetest", fpetest_methods);
    d = PyModule_GetDict(m);
181 182 183
    fpe_error = PyErr_NewException("fpetest.error", NULL, NULL);
    if (fpe_error != NULL)
	    PyDict_SetItemString(d, "error", fpe_error);
184
}