Modules/fpetestmodule.c

/*
     ---------------------------------------------------------------------  
    /                       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;
void initfpetest(void);
static PyObject *test(PyObject *self,PyObject *args);
static int db0(void);
static int overflow(void);
static int nest1(double);
static int nest2(double);
static double nest3(double);

static PyMethodDef fpetest_methods[] = {
    {"test",		 (PyCFunction) test,		 1},
    {0,0}
};

static PyObject *test(PyObject *self,PyObject *args)
{
    int i = 0, r;

    fprintf(stderr,"Test trapping overflow\n");
    r = overflow();
    if(r){
      fprintf(stderr,"(Note: No exception was raised.)\n");
      PyErr_Clear();
    }else{
      fprintf(stderr,"Trapped:: ");
      PyErr_Print();
      PyErr_Clear();
    }
    i += r;

    fprintf(stderr,"Test trapping division by zero\n");
    r = db0();
    if(r){
      fprintf(stderr,"(Note: No exception was raised.)\n");
      PyErr_Clear();
    }else{
      fprintf(stderr,"Trapped:: ");
      PyErr_Print();
      PyErr_Clear();
    }
    i += r;

    fprintf(stderr,"Test nested protection zones, outer zone\n");
    r = nest1(0.0);
    if(r){
      fprintf(stderr,"(Note: No exception was raised.)\n");
      PyErr_Clear();
    }else{
      fprintf(stderr,"Trapped:: ");
      PyErr_Print();
      PyErr_Clear();
    }
    i += r;

    fprintf(stderr,"Test nested protection zones, inner zone\n");
    fprintf(stderr,"(Note: Return will apparently come from outer zone.)\n");
    r = nest1(1.0);
    if(r){
      fprintf(stderr,"(Note: No exception was raised.)\n");
      PyErr_Clear();
    }else{
      fprintf(stderr,"Trapped:: ");
      PyErr_Print();
      PyErr_Clear();
    }
    i += r;

    fprintf(stderr,"Test nested protection zones, trailing outer zone\n");
    r = nest1(2.0);
    if(r){
      fprintf(stderr,"(Note: No exception was raised.)\n");
      PyErr_Clear();
    }else{
      fprintf(stderr,"Trapped:: ");
      PyErr_Print();
      PyErr_Clear();
    }
    i += r;

    fprintf(stderr,"Test nested function calls, prior error\n");
    r = nest2(0.0);
    if(r){
      fprintf(stderr,"(Note: No exception was raised.)\n");
      PyErr_Clear();
    }else{
      fprintf(stderr,"Trapped:: ");
      PyErr_Print();
      PyErr_Clear();
    }
    i += r;

    fprintf(stderr,"Test nested function calls, interior error\n");
    r = nest2(1.0);
    if(r){
      fprintf(stderr,"(Note: No exception was raised.)\n");
      PyErr_Clear();
    }else{
      fprintf(stderr,"Trapped:: ");
      PyErr_Print();
      PyErr_Clear();
    }
    i += r;

    fprintf(stderr,"Test nested function calls, trailing error\n");
    r = nest2(2.0);
    if(r){
      fprintf(stderr,"(Note: No exception was raised.)\n");
      PyErr_Clear();
    }else{
      fprintf(stderr,"Trapped:: ");
      PyErr_Print();
      PyErr_Clear();
    }
    i += r;

    fprintf(stderr,"Number of tests failed: %d\n", i);
    Py_INCREF (Py_None);
    return Py_None;
}

static int nest1(double x)
{
  double a = 1.0;
  PyFPE_START_PROTECT("Division by zero, outer zone", return 0)
  a = 1./x;
    PyFPE_START_PROTECT("Division by zero, inner zone", return 0)
    a = 1./(1. - x);
    PyFPE_END_PROTECT
  a = 1./(2. - x);
  PyFPE_END_PROTECT
  return(1);
}

static int nest2(double x)
{
  double a = 1.0;
  PyFPE_START_PROTECT("Division by zero, prior error", return 0)
  a = 1./x;
  a = nest3(x);
  a = 1./(2. - x);
  PyFPE_END_PROTECT
  return(1);
}

static double nest3(double x)
{
  double result;
  PyFPE_START_PROTECT("Division by zero, nest3 error", return 0)
  result = 1./(1. - x);
  PyFPE_END_PROTECT
  return result;
}

static int db0(void)
{
  double a = 1.0;
  PyFPE_START_PROTECT("Division by zero", return 0)
  a = 1./(a - 1.);
  PyFPE_END_PROTECT
  return(1);
}

static int overflow(void)
{
  double a, b = 1.e200;
  PyFPE_START_PROTECT("Overflow", return 0)
  a = b*b;
  PyFPE_END_PROTECT
  return(1);
}

void initfpetest(void)
{
    PyObject *m, *d;

    m = Py_InitModule("fpetest", fpetest_methods);
    d = PyModule_GetDict(m);
    fpe_error = PyString_FromString("fpetest.error");
    PyDict_SetItemString(d, "error", fpe_error);

    if (PyErr_Occurred())
	Py_FatalError("Cannot initialize module fpetest");
}