WebSVN - shark - Blame - Rev 2 - /shark/branches/pj/oslib/libm/msun/src/k

Rev	Author	Line No.	Line
2	pj	1	/* @(#)k_standard.c 5.1 93/09/24 */
		2	/*
		3	* ====================================================
		4	* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
		5	*
		6	* Developed at SunPro, a Sun Microsystems, Inc. business.
		7	* Permission to use, copy, modify, and distribute this
		8	* software is freely granted, provided that this notice
		9	* is preserved.
		10	* ====================================================
		11	*/
		12
		13	#ifndef lint
		14	static char rcsid[] = "$\Id: k_standard.c,v 1.2 1995/05/30 05:49:13 rgrimes Exp $";
		15	#endif
		16
		17	#include "math.h"
		18	#include "math_private.h"
		19	#include <errno.h>
		20
		21	/* Undefine these if you're actually going to use the FreeBSD libc. */
		22	#define _USE_WRITE
		23	#define _NO_FFLUSH
		24
		25	#ifndef _USE_WRITE
		26	#include <stdio.h> /* fputs(), stderr */
		27	#define WRITE2(u,v) fputs(u, stderr)
		28	#else /* !defined(_USE_WRITE) */
		29	#include <unistd.h> /* write */
		30	#define WRITE2(u,v) write(2, u, v)
		31	#undef fflush
		32	#endif /* !defined(_USE_WRITE) */
		33
		34	#ifdef __STDC__
		35	static const double zero = 0.0; /* used as const */
		36	#else
		37	static double zero = 0.0; /* used as const */
		38	#endif
		39
		40	/*
		41	* Standard conformance (non-IEEE) on exception cases.
		42	* Mapping:
		43	* 1 -- acos(\|x\|>1)
		44	* 2 -- asin(\|x\|>1)
		45	* 3 -- atan2(+-0,+-0)
		46	* 4 -- hypot overflow
		47	* 5 -- cosh overflow
		48	* 6 -- exp overflow
		49	* 7 -- exp underflow
		50	* 8 -- y0(0)
		51	* 9 -- y0(-ve)
		52	* 10-- y1(0)
		53	* 11-- y1(-ve)
		54	* 12-- yn(0)
		55	* 13-- yn(-ve)
		56	* 14-- lgamma(finite) overflow
		57	* 15-- lgamma(-integer)
		58	* 16-- log(0)
		59	* 17-- log(x<0)
		60	* 18-- log10(0)
		61	* 19-- log10(x<0)
		62	* 20-- pow(0.0,0.0)
		63	* 21-- pow(x,y) overflow
		64	* 22-- pow(x,y) underflow
		65	* 23-- pow(0,negative)
		66	* 24-- pow(neg,non-integral)
		67	* 25-- sinh(finite) overflow
		68	* 26-- sqrt(negative)
		69	* 27-- fmod(x,0)
		70	* 28-- remainder(x,0)
		71	* 29-- acosh(x<1)
		72	* 30-- atanh(\|x\|>1)
		73	* 31-- atanh(\|x\|=1)
		74	* 32-- scalb overflow
		75	* 33-- scalb underflow
		76	* 34-- j0(\|x\|>X_TLOSS)
		77	* 35-- y0(x>X_TLOSS)
		78	* 36-- j1(\|x\|>X_TLOSS)
		79	* 37-- y1(x>X_TLOSS)
		80	* 38-- jn(\|x\|>X_TLOSS, n)
		81	* 39-- yn(x>X_TLOSS, n)
		82	* 40-- gamma(finite) overflow
		83	* 41-- gamma(-integer)
		84	* 42-- pow(NaN,0.0)
		85	*/
		86
		87
		88	#ifdef __STDC__
		89	double __kernel_standard(double x, double y, int type)
		90	#else
		91	double __kernel_standard(x,y,type)
		92	double x,y; int type;
		93	#endif
		94	{
		95	struct exception exc;
		96	#ifndef HUGE_VAL /* this is the only routine that uses HUGE_VAL */
		97	#define HUGE_VAL inf
		98	double inf = 0.0;
		99
		100	SET_HIGH_WORD(inf,0x7ff00000); /* set inf to infinite */
		101	#endif
		102
		103	#if defined(_USE_WRITE) && !defined(_NO_FFLUSH)
		104	(void) fflush(stdout);
		105	#endif
		106	exc.arg1 = x;
		107	exc.arg2 = y;
		108	switch(type) {
		109	case 1:
		110	case 101:
		111	/* acos(\|x\|>1) */
		112	exc.type = DOMAIN;
		113	exc.name = type < 100 ? "acos" : "acosf";
		114	exc.retval = zero;
		115	if (_LIB_VERSION == _POSIX_)
		116	errno = EDOM;
		117	else if (!matherr(&exc)) {
		118	if(_LIB_VERSION == _SVID_) {
		119	(void) WRITE2("acos: DOMAIN error\n", 19);
		120	}
		121	errno = EDOM;
		122	}
		123	break;
		124	case 2:
		125	case 102:
		126	/* asin(\|x\|>1) */
		127	exc.type = DOMAIN;
		128	exc.name = type < 100 ? "asin" : "asinf";
		129	exc.retval = zero;
		130	if(_LIB_VERSION == _POSIX_)
		131	errno = EDOM;
		132	else if (!matherr(&exc)) {
		133	if(_LIB_VERSION == _SVID_) {
		134	(void) WRITE2("asin: DOMAIN error\n", 19);
		135	}
		136	errno = EDOM;
		137	}
		138	break;
		139	case 3:
		140	case 103:
		141	/* atan2(+-0,+-0) */
		142	exc.arg1 = y;
		143	exc.arg2 = x;
		144	exc.type = DOMAIN;
		145	exc.name = type < 100 ? "atan2" : "atan2f";
		146	exc.retval = zero;
		147	if(_LIB_VERSION == _POSIX_)
		148	errno = EDOM;
		149	else if (!matherr(&exc)) {
		150	if(_LIB_VERSION == _SVID_) {
		151	(void) WRITE2("atan2: DOMAIN error\n", 20);
		152	}
		153	errno = EDOM;
		154	}
		155	break;
		156	case 4:
		157	case 104:
		158	/* hypot(finite,finite) overflow */
		159	exc.type = OVERFLOW;
		160	exc.name = type < 100 ? "hypot" : "hypotf";
		161	if (_LIB_VERSION == _SVID_)
		162	exc.retval = HUGE;
		163	else
		164	exc.retval = HUGE_VAL;
		165	if (_LIB_VERSION == _POSIX_)
		166	errno = ERANGE;
		167	else if (!matherr(&exc)) {
		168	errno = ERANGE;
		169	}
		170	break;
		171	case 5:
		172	case 105:
		173	/* cosh(finite) overflow */
		174	exc.type = OVERFLOW;
		175	exc.name = type < 100 ? "cosh" : "coshf";
		176	if (_LIB_VERSION == _SVID_)
		177	exc.retval = HUGE;
		178	else
		179	exc.retval = HUGE_VAL;
		180	if (_LIB_VERSION == _POSIX_)
		181	errno = ERANGE;
		182	else if (!matherr(&exc)) {
		183	errno = ERANGE;
		184	}
		185	break;
		186	case 6:
		187	case 106:
		188	/* exp(finite) overflow */
		189	exc.type = OVERFLOW;
		190	exc.name = type < 100 ? "exp" : "expf";
		191	if (_LIB_VERSION == _SVID_)
		192	exc.retval = HUGE;
		193	else
		194	exc.retval = HUGE_VAL;
		195	if (_LIB_VERSION == _POSIX_)
		196	errno = ERANGE;
		197	else if (!matherr(&exc)) {
		198	errno = ERANGE;
		199	}
		200	break;
		201	case 7:
		202	case 107:
		203	/* exp(finite) underflow */
		204	exc.type = UNDERFLOW;
		205	exc.name = type < 100 ? "exp" : "expf";
		206	exc.retval = zero;
		207	if (_LIB_VERSION == _POSIX_)
		208	errno = ERANGE;
		209	else if (!matherr(&exc)) {
		210	errno = ERANGE;
		211	}
		212	break;
		213	case 8:
		214	case 108:
		215	/* y0(0) = -inf */
		216	exc.type = DOMAIN; /* should be SING for IEEE */
		217	exc.name = type < 100 ? "y0" : "y0f";
		218	if (_LIB_VERSION == _SVID_)
		219	exc.retval = -HUGE;
		220	else
		221	exc.retval = -HUGE_VAL;
		222	if (_LIB_VERSION == _POSIX_)
		223	errno = EDOM;
		224	else if (!matherr(&exc)) {
		225	if (_LIB_VERSION == _SVID_) {
		226	(void) WRITE2("y0: DOMAIN error\n", 17);
		227	}
		228	errno = EDOM;
		229	}
		230	break;
		231	case 9:
		232	case 109:
		233	/* y0(x<0) = NaN */
		234	exc.type = DOMAIN;
		235	exc.name = type < 100 ? "y0" : "y0f";
		236	if (_LIB_VERSION == _SVID_)
		237	exc.retval = -HUGE;
		238	else
		239	exc.retval = -HUGE_VAL;
		240	if (_LIB_VERSION == _POSIX_)
		241	errno = EDOM;
		242	else if (!matherr(&exc)) {
		243	if (_LIB_VERSION == _SVID_) {
		244	(void) WRITE2("y0: DOMAIN error\n", 17);
		245	}
		246	errno = EDOM;
		247	}
		248	break;
		249	case 10:
		250	case 110:
		251	/* y1(0) = -inf */
		252	exc.type = DOMAIN; /* should be SING for IEEE */
		253	exc.name = type < 100 ? "y1" : "y1f";
		254	if (_LIB_VERSION == _SVID_)
		255	exc.retval = -HUGE;
		256	else
		257	exc.retval = -HUGE_VAL;
		258	if (_LIB_VERSION == _POSIX_)
		259	errno = EDOM;
		260	else if (!matherr(&exc)) {
		261	if (_LIB_VERSION == _SVID_) {
		262	(void) WRITE2("y1: DOMAIN error\n", 17);
		263	}
		264	errno = EDOM;
		265	}
		266	break;
		267	case 11:
		268	case 111:
		269	/* y1(x<0) = NaN */
		270	exc.type = DOMAIN;
		271	exc.name = type < 100 ? "y1" : "y1f";
		272	if (_LIB_VERSION == _SVID_)
		273	exc.retval = -HUGE;
		274	else
		275	exc.retval = -HUGE_VAL;
		276	if (_LIB_VERSION == _POSIX_)
		277	errno = EDOM;
		278	else if (!matherr(&exc)) {
		279	if (_LIB_VERSION == _SVID_) {
		280	(void) WRITE2("y1: DOMAIN error\n", 17);
		281	}
		282	errno = EDOM;
		283	}
		284	break;
		285	case 12:
		286	case 112:
		287	/* yn(n,0) = -inf */
		288	exc.type = DOMAIN; /* should be SING for IEEE */
		289	exc.name = type < 100 ? "yn" : "ynf";
		290	if (_LIB_VERSION == _SVID_)
		291	exc.retval = -HUGE;
		292	else
		293	exc.retval = -HUGE_VAL;
		294	if (_LIB_VERSION == _POSIX_)
		295	errno = EDOM;
		296	else if (!matherr(&exc)) {
		297	if (_LIB_VERSION == _SVID_) {
		298	(void) WRITE2("yn: DOMAIN error\n", 17);
		299	}
		300	errno = EDOM;
		301	}
		302	break;
		303	case 13:
		304	case 113:
		305	/* yn(x<0) = NaN */
		306	exc.type = DOMAIN;
		307	exc.name = type < 100 ? "yn" : "ynf";
		308	if (_LIB_VERSION == _SVID_)
		309	exc.retval = -HUGE;
		310	else
		311	exc.retval = -HUGE_VAL;
		312	if (_LIB_VERSION == _POSIX_)
		313	errno = EDOM;
		314	else if (!matherr(&exc)) {
		315	if (_LIB_VERSION == _SVID_) {
		316	(void) WRITE2("yn: DOMAIN error\n", 17);
		317	}
		318	errno = EDOM;
		319	}
		320	break;
		321	case 14:
		322	case 114:
		323	/* lgamma(finite) overflow */
		324	exc.type = OVERFLOW;
		325	exc.name = type < 100 ? "lgamma" : "lgammaf";
		326	if (_LIB_VERSION == _SVID_)
		327	exc.retval = HUGE;
		328	else
		329	exc.retval = HUGE_VAL;
		330	if (_LIB_VERSION == _POSIX_)
		331	errno = ERANGE;
		332	else if (!matherr(&exc)) {
		333	errno = ERANGE;
		334	}
		335	break;
		336	case 15:
		337	case 115:
		338	/* lgamma(-integer) or lgamma(0) */
		339	exc.type = SING;
		340	exc.name = type < 100 ? "lgamma" : "lgammaf";
		341	if (_LIB_VERSION == _SVID_)
		342	exc.retval = HUGE;
		343	else
		344	exc.retval = HUGE_VAL;
		345	if (_LIB_VERSION == _POSIX_)
		346	errno = EDOM;
		347	else if (!matherr(&exc)) {
		348	if (_LIB_VERSION == _SVID_) {
		349	(void) WRITE2("lgamma: SING error\n", 19);
		350	}
		351	errno = EDOM;
		352	}
		353	break;
		354	case 16:
		355	case 116:
		356	/* log(0) */
		357	exc.type = SING;
		358	exc.name = type < 100 ? "log" : "logf";
		359	if (_LIB_VERSION == _SVID_)
		360	exc.retval = -HUGE;
		361	else
		362	exc.retval = -HUGE_VAL;
		363	if (_LIB_VERSION == _POSIX_)
		364	errno = ERANGE;
		365	else if (!matherr(&exc)) {
		366	if (_LIB_VERSION == _SVID_) {
		367	(void) WRITE2("log: SING error\n", 16);
		368	}
		369	errno = EDOM;
		370	}
		371	break;
		372	case 17:
		373	case 117:
		374	/* log(x<0) */
		375	exc.type = DOMAIN;
		376	exc.name = type < 100 ? "log" : "logf";
		377	if (_LIB_VERSION == _SVID_)
		378	exc.retval = -HUGE;
		379	else
		380	exc.retval = -HUGE_VAL;
		381	if (_LIB_VERSION == _POSIX_)
		382	errno = EDOM;
		383	else if (!matherr(&exc)) {
		384	if (_LIB_VERSION == _SVID_) {
		385	(void) WRITE2("log: DOMAIN error\n", 18);
		386	}
		387	errno = EDOM;
		388	}
		389	break;
		390	case 18:
		391	case 118:
		392	/* log10(0) */
		393	exc.type = SING;
		394	exc.name = type < 100 ? "log10" : "log10f";
		395	if (_LIB_VERSION == _SVID_)
		396	exc.retval = -HUGE;
		397	else
		398	exc.retval = -HUGE_VAL;
		399	if (_LIB_VERSION == _POSIX_)
		400	errno = ERANGE;
		401	else if (!matherr(&exc)) {
		402	if (_LIB_VERSION == _SVID_) {
		403	(void) WRITE2("log10: SING error\n", 18);
		404	}
		405	errno = EDOM;
		406	}
		407	break;
		408	case 19:
		409	case 119:
		410	/* log10(x<0) */
		411	exc.type = DOMAIN;
		412	exc.name = type < 100 ? "log10" : "log10f";
		413	if (_LIB_VERSION == _SVID_)
		414	exc.retval = -HUGE;
		415	else
		416	exc.retval = -HUGE_VAL;
		417	if (_LIB_VERSION == _POSIX_)
		418	errno = EDOM;
		419	else if (!matherr(&exc)) {
		420	if (_LIB_VERSION == _SVID_) {
		421	(void) WRITE2("log10: DOMAIN error\n", 20);
		422	}
		423	errno = EDOM;
		424	}
		425	break;
		426	case 20:
		427	case 120:
		428	/* pow(0.0,0.0) */
		429	/* error only if _LIB_VERSION == _SVID_ */
		430	exc.type = DOMAIN;
		431	exc.name = type < 100 ? "pow" : "powf";
		432	exc.retval = zero;
		433	if (_LIB_VERSION != _SVID_) exc.retval = 1.0;
		434	else if (!matherr(&exc)) {
		435	(void) WRITE2("pow(0,0): DOMAIN error\n", 23);
		436	errno = EDOM;
		437	}
		438	break;
		439	case 21:
		440	case 121:
		441	/* pow(x,y) overflow */
		442	exc.type = OVERFLOW;
		443	exc.name = type < 100 ? "pow" : "powf";
		444	if (_LIB_VERSION == _SVID_) {
		445	exc.retval = HUGE;
		446	y *= 0.5;
		447	if(x<zero&&rint(y)!=y) exc.retval = -HUGE;
		448	} else {
		449	exc.retval = HUGE_VAL;
		450	y *= 0.5;
		451	if(x<zero&&rint(y)!=y) exc.retval = -HUGE_VAL;
		452	}
		453	if (_LIB_VERSION == _POSIX_)
		454	errno = ERANGE;
		455	else if (!matherr(&exc)) {
		456	errno = ERANGE;
		457	}
		458	break;
		459	case 22:
		460	case 122:
		461	/* pow(x,y) underflow */
		462	exc.type = UNDERFLOW;
		463	exc.name = type < 100 ? "pow" : "powf";
		464	exc.retval = zero;
		465	if (_LIB_VERSION == _POSIX_)
		466	errno = ERANGE;
		467	else if (!matherr(&exc)) {
		468	errno = ERANGE;
		469	}
		470	break;
		471	case 23:
		472	case 123:
		473	/* 0*neg /
		474	exc.type = DOMAIN;
		475	exc.name = type < 100 ? "pow" : "powf";
		476	if (_LIB_VERSION == _SVID_)
		477	exc.retval = zero;
		478	else
		479	exc.retval = -HUGE_VAL;
		480	if (_LIB_VERSION == _POSIX_)
		481	errno = EDOM;
		482	else if (!matherr(&exc)) {
		483	if (_LIB_VERSION == _SVID_) {
		484	(void) WRITE2("pow(0,neg): DOMAIN error\n", 25);
		485	}
		486	errno = EDOM;
		487	}
		488	break;
		489	case 24:
		490	case 124:
		491	/* neg*non-integral /
		492	exc.type = DOMAIN;
		493	exc.name = type < 100 ? "pow" : "powf";
		494	if (_LIB_VERSION == _SVID_)
		495	exc.retval = zero;
		496	else
		497	exc.retval = zero/zero; /* X/Open allow NaN */
		498	if (_LIB_VERSION == _POSIX_)
		499	errno = EDOM;
		500	else if (!matherr(&exc)) {
		501	if (_LIB_VERSION == _SVID_) {
		502	(void) WRITE2("neg**non-integral: DOMAIN error\n", 32);
		503	}
		504	errno = EDOM;
		505	}
		506	break;
		507	case 25:
		508	case 125:
		509	/* sinh(finite) overflow */
		510	exc.type = OVERFLOW;
		511	exc.name = type < 100 ? "sinh" : "sinhf";
		512	if (_LIB_VERSION == _SVID_)
		513	exc.retval = ( (x>zero) ? HUGE : -HUGE);
		514	else
		515	exc.retval = ( (x>zero) ? HUGE_VAL : -HUGE_VAL);
		516	if (_LIB_VERSION == _POSIX_)
		517	errno = ERANGE;
		518	else if (!matherr(&exc)) {
		519	errno = ERANGE;
		520	}
		521	break;
		522	case 26:
		523	case 126:
		524	/* sqrt(x<0) */
		525	exc.type = DOMAIN;
		526	exc.name = type < 100 ? "sqrt" : "sqrtf";
		527	if (_LIB_VERSION == _SVID_)
		528	exc.retval = zero;
		529	else
		530	exc.retval = zero/zero;
		531	if (_LIB_VERSION == _POSIX_)
		532	errno = EDOM;
		533	else if (!matherr(&exc)) {
		534	if (_LIB_VERSION == _SVID_) {
		535	(void) WRITE2("sqrt: DOMAIN error\n", 19);
		536	}
		537	errno = EDOM;
		538	}
		539	break;
		540	case 27:
		541	case 127:
		542	/* fmod(x,0) */
		543	exc.type = DOMAIN;
		544	exc.name = type < 100 ? "fmod" : "fmodf";
		545	if (_LIB_VERSION == _SVID_)
		546	exc.retval = x;
		547	else
		548	exc.retval = zero/zero;
		549	if (_LIB_VERSION == _POSIX_)
		550	errno = EDOM;
		551	else if (!matherr(&exc)) {
		552	if (_LIB_VERSION == _SVID_) {
		553	(void) WRITE2("fmod: DOMAIN error\n", 20);
		554	}
		555	errno = EDOM;
		556	}
		557	break;
		558	case 28:
		559	case 128:
		560	/* remainder(x,0) */
		561	exc.type = DOMAIN;
		562	exc.name = type < 100 ? "remainder" : "remainderf";
		563	exc.retval = zero/zero;
		564	if (_LIB_VERSION == _POSIX_)
		565	errno = EDOM;
		566	else if (!matherr(&exc)) {
		567	if (_LIB_VERSION == _SVID_) {
		568	(void) WRITE2("remainder: DOMAIN error\n", 24);
		569	}
		570	errno = EDOM;
		571	}
		572	break;
		573	case 29:
		574	case 129:
		575	/* acosh(x<1) */
		576	exc.type = DOMAIN;
		577	exc.name = type < 100 ? "acosh" : "acoshf";
		578	exc.retval = zero/zero;
		579	if (_LIB_VERSION == _POSIX_)
		580	errno = EDOM;
		581	else if (!matherr(&exc)) {
		582	if (_LIB_VERSION == _SVID_) {
		583	(void) WRITE2("acosh: DOMAIN error\n", 20);
		584	}
		585	errno = EDOM;
		586	}
		587	break;
		588	case 30:
		589	case 130:
		590	/* atanh(\|x\|>1) */
		591	exc.type = DOMAIN;
		592	exc.name = type < 100 ? "atanh" : "atanhf";
		593	exc.retval = zero/zero;
		594	if (_LIB_VERSION == _POSIX_)
		595	errno = EDOM;
		596	else if (!matherr(&exc)) {
		597	if (_LIB_VERSION == _SVID_) {
		598	(void) WRITE2("atanh: DOMAIN error\n", 20);
		599	}
		600	errno = EDOM;
		601	}
		602	break;
		603	case 31:
		604	case 131:
		605	/* atanh(\|x\|=1) */
		606	exc.type = SING;
		607	exc.name = type < 100 ? "atanh" : "atanhf";
		608	exc.retval = x/zero; /* sign(x)inf /
		609	if (_LIB_VERSION == _POSIX_)
		610	errno = EDOM;
		611	else if (!matherr(&exc)) {
		612	if (_LIB_VERSION == _SVID_) {
		613	(void) WRITE2("atanh: SING error\n", 18);
		614	}
		615	errno = EDOM;
		616	}
		617	break;
		618	case 32:
		619	case 132:
		620	/* scalb overflow; SVID also returns +-HUGE_VAL */
		621	exc.type = OVERFLOW;
		622	exc.name = type < 100 ? "scalb" : "scalbf";
		623	exc.retval = x > zero ? HUGE_VAL : -HUGE_VAL;
		624	if (_LIB_VERSION == _POSIX_)
		625	errno = ERANGE;
		626	else if (!matherr(&exc)) {
		627	errno = ERANGE;
		628	}
		629	break;
		630	case 33:
		631	case 133:
		632	/* scalb underflow */
		633	exc.type = UNDERFLOW;
		634	exc.name = type < 100 ? "scalb" : "scalbf";
		635	exc.retval = copysign(zero,x);
		636	if (_LIB_VERSION == _POSIX_)
		637	errno = ERANGE;
		638	else if (!matherr(&exc)) {
		639	errno = ERANGE;
		640	}
		641	break;
		642	case 34:
		643	case 134:
		644	/* j0(\|x\|>X_TLOSS) */
		645	exc.type = TLOSS;
		646	exc.name = type < 100 ? "j0" : "j0f";
		647	exc.retval = zero;
		648	if (_LIB_VERSION == _POSIX_)
		649	errno = ERANGE;
		650	else if (!matherr(&exc)) {
		651	if (_LIB_VERSION == _SVID_) {
		652	(void) WRITE2(exc.name, 2);
		653	(void) WRITE2(": TLOSS error\n", 14);
		654	}
		655	errno = ERANGE;
		656	}
		657	break;
		658	case 35:
		659	case 135:
		660	/* y0(x>X_TLOSS) */
		661	exc.type = TLOSS;
		662	exc.name = type < 100 ? "y0" : "y0f";
		663	exc.retval = zero;
		664	if (_LIB_VERSION == _POSIX_)
		665	errno = ERANGE;
		666	else if (!matherr(&exc)) {
		667	if (_LIB_VERSION == _SVID_) {
		668	(void) WRITE2(exc.name, 2);
		669	(void) WRITE2(": TLOSS error\n", 14);
		670	}
		671	errno = ERANGE;
		672	}
		673	break;
		674	case 36:
		675	case 136:
		676	/* j1(\|x\|>X_TLOSS) */
		677	exc.type = TLOSS;
		678	exc.name = type < 100 ? "j1" : "j1f";
		679	exc.retval = zero;
		680	if (_LIB_VERSION == _POSIX_)
		681	errno = ERANGE;
		682	else if (!matherr(&exc)) {
		683	if (_LIB_VERSION == _SVID_) {
		684	(void) WRITE2(exc.name, 2);
		685	(void) WRITE2(": TLOSS error\n", 14);
		686	}
		687	errno = ERANGE;
		688	}
		689	break;
		690	case 37:
		691	case 137:
		692	/* y1(x>X_TLOSS) */
		693	exc.type = TLOSS;
		694	exc.name = type < 100 ? "y1" : "y1f";
		695	exc.retval = zero;
		696	if (_LIB_VERSION == _POSIX_)
		697	errno = ERANGE;
		698	else if (!matherr(&exc)) {
		699	if (_LIB_VERSION == _SVID_) {
		700	(void) WRITE2(exc.name, 2);
		701	(void) WRITE2(": TLOSS error\n", 14);
		702	}
		703	errno = ERANGE;
		704	}
		705	break;
		706	case 38:
		707	case 138:
		708	/* jn(\|x\|>X_TLOSS) */
		709	exc.type = TLOSS;
		710	exc.name = type < 100 ? "jn" : "jnf";
		711	exc.retval = zero;
		712	if (_LIB_VERSION == _POSIX_)
		713	errno = ERANGE;
		714	else if (!matherr(&exc)) {
		715	if (_LIB_VERSION == _SVID_) {
		716	(void) WRITE2(exc.name, 2);
		717	(void) WRITE2(": TLOSS error\n", 14);
		718	}
		719	errno = ERANGE;
		720	}
		721	break;
		722	case 39:
		723	case 139:
		724	/* yn(x>X_TLOSS) */
		725	exc.type = TLOSS;
		726	exc.name = type < 100 ? "yn" : "ynf";
		727	exc.retval = zero;
		728	if (_LIB_VERSION == _POSIX_)
		729	errno = ERANGE;
		730	else if (!matherr(&exc)) {
		731	if (_LIB_VERSION == _SVID_) {
		732	(void) WRITE2(exc.name, 2);
		733	(void) WRITE2(": TLOSS error\n", 14);
		734	}
		735	errno = ERANGE;
		736	}
		737	break;
		738	case 40:
		739	case 140:
		740	/* gamma(finite) overflow */
		741	exc.type = OVERFLOW;
		742	exc.name = type < 100 ? "gamma" : "gammaf";
		743	if (_LIB_VERSION == _SVID_)
		744	exc.retval = HUGE;
		745	else
		746	exc.retval = HUGE_VAL;
		747	if (_LIB_VERSION == _POSIX_)
		748	errno = ERANGE;
		749	else if (!matherr(&exc)) {
		750	errno = ERANGE;
		751	}
		752	break;
		753	case 41:
		754	case 141:
		755	/* gamma(-integer) or gamma(0) */
		756	exc.type = SING;
		757	exc.name = type < 100 ? "gamma" : "gammaf";
		758	if (_LIB_VERSION == _SVID_)
		759	exc.retval = HUGE;
		760	else
		761	exc.retval = HUGE_VAL;
		762	if (_LIB_VERSION == _POSIX_)
		763	errno = EDOM;
		764	else if (!matherr(&exc)) {
		765	if (_LIB_VERSION == _SVID_) {
		766	(void) WRITE2("gamma: SING error\n", 18);
		767	}
		768	errno = EDOM;
		769	}
		770	break;
		771	case 42:
		772	case 142:
		773	/* pow(NaN,0.0) */
		774	/* error only if _LIB_VERSION == _SVID_ & _XOPEN_ */
		775	exc.type = DOMAIN;
		776	exc.name = type < 100 ? "pow" : "powf";
		777	exc.retval = x;
		778	if (_LIB_VERSION == _IEEE_ \|\|
		779	_LIB_VERSION == _POSIX_) exc.retval = 1.0;
		780	else if (!matherr(&exc)) {
		781	errno = EDOM;
		782	}
		783	break;
		784	}
		785	return exc.retval;
		786	}

Subversion Repositories shark

(root)/shark/branches/pj/oslib/libm/msun/src/k_standa.c - Rev 2