MATHERR 3M "630 MTG" ""

NAME
matherr - error-handling function

SYNOPSIS

#include <ccs/math.h>

int matherr (x)
struct exception *x;

DESCRIPTION
Matherr is invoked by functions in the Math Library when errors are detected. Users may define their own procedures for handling errors, by including a function named matherr in their programs. Matherr must be of the form described above. When an error occurs, a pointer to the exception structure x will be passed to the user-supplied matherr function. This structure, which is defined in the <ccs/math.h> header file, is as follows:

    struct exception {
    	int type;
    	char *name;
    	double arg1, arg2, retval;
    };

The element type is an integer describing the type of error that has occurred, from the following list of constants (defined in the header file):

    DOMAIN	argument domain error
    SING	argument singularity
    OVERFLOW	overflow range error
    UNDERFLOW 	underflow range error
    TLOSS	total loss of significance
    PLOSS	partial loss of significance

The element name points to a string containing the name of the function that incurred the error. The variables arg1 and arg2 are the arguments with which the function was invoked. Retval is set to the default value that will be returned by the function unless the user's matherr sets it to a different value.

If the user's matherr function returns non-zero, no error message will be printed, and errno will not be set.

If matherr is not supplied by the user, the default error-handling procedures, described with the math functions involved, will be invoked upon error. These procedures are also summarized in the table below. In every case, errno is set to EDOM or ERANGE and the program continues.

EXAMPLE
#include <ccs/math.h>
int
matherr(x)
register struct exception *x;
{
	switch (x->type) {
	case DOMAIN:
		/* change sqrt to return sqrt(-arg1), not 0 */
		if (!strcmp(x->name, "sqrt")) {
			x->retval = sqrt(-x->arg1);
			return (0); /* print message and set errno */
		}
	case SING:
		/* all other domain or sing errors, print message
		   and abort */
		lprintf("domain error in %s\n", x->name);
		for(;;) sleep(300);
	case PLOSS:
		/* print detailed error message */
		lprintf("loss of significance in %s(%g) = %g\n",
		    x->name, x->arg1, x->retval);
		return (1); /* take no other action */
	}
	/* all other errors, execute default procedure */
	return (0);
}

DEFAULT ERROR HANDLING PROCEDURES
 | Types of Errors
type | DOMAIN | SING | OVERFLOW | UNDERFLOW | TLOSS | PLOSS
errno | EDOM | EDOM | ERANGE | ERANGE | ERANGE | ERANGE
BESSEL: | - | - | - | - | M, 0 | *
y0, y1, yn (arg ≤ 0) | M, -H | - | - | - | - | -
EXP: | - | - | H | 0 | - | -
LOG, LOG10:
(arg < 0) | M, -H | - | - | - | - | -
(arg = 0) | - | M, -H | - | - | - | -
POW: | - | - | ±H | 0 | - | -
neg ** non-int | M, 0 | - | - | - | - | -
0 ** non-pos
SQRT: | M, 0 | - | - | - | - | -
GAMMA: | - | M, H | H | - | - | -
HYPOT: | - | - | H | - | - | -
SINH: | - | - | ±H | - | - | -
COSH: | - | - | H | - | - | -
SIN, COS, TAN: - | - | - | - | M, 0 | *
ASIN, ACOS, ATAN2: M, 0 | - | - | - | - | -

ABBREVIATIONS
* | As much as possible of the value is returned.
M | Message is printed (EDOM error).
H | HUGE is returned.
-H | -HUGE is returned.
±H | HUGE or -HUGE is returned.
0 | 0 is returned.