15b2ba9d3SPiotr Jasiukajtis /*
25b2ba9d3SPiotr Jasiukajtis * CDDL HEADER START
35b2ba9d3SPiotr Jasiukajtis *
45b2ba9d3SPiotr Jasiukajtis * The contents of this file are subject to the terms of the
55b2ba9d3SPiotr Jasiukajtis * Common Development and Distribution License (the "License").
65b2ba9d3SPiotr Jasiukajtis * You may not use this file except in compliance with the License.
75b2ba9d3SPiotr Jasiukajtis *
85b2ba9d3SPiotr Jasiukajtis * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
95b2ba9d3SPiotr Jasiukajtis * or http://www.opensolaris.org/os/licensing.
105b2ba9d3SPiotr Jasiukajtis * See the License for the specific language governing permissions
115b2ba9d3SPiotr Jasiukajtis * and limitations under the License.
125b2ba9d3SPiotr Jasiukajtis *
135b2ba9d3SPiotr Jasiukajtis * When distributing Covered Code, include this CDDL HEADER in each
145b2ba9d3SPiotr Jasiukajtis * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
155b2ba9d3SPiotr Jasiukajtis * If applicable, add the following below this CDDL HEADER, with the
165b2ba9d3SPiotr Jasiukajtis * fields enclosed by brackets "[]" replaced with your own identifying
175b2ba9d3SPiotr Jasiukajtis * information: Portions Copyright [yyyy] [name of copyright owner]
185b2ba9d3SPiotr Jasiukajtis *
195b2ba9d3SPiotr Jasiukajtis * CDDL HEADER END
205b2ba9d3SPiotr Jasiukajtis */
215b2ba9d3SPiotr Jasiukajtis
225b2ba9d3SPiotr Jasiukajtis /*
235b2ba9d3SPiotr Jasiukajtis * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
245b2ba9d3SPiotr Jasiukajtis */
255b2ba9d3SPiotr Jasiukajtis /*
265b2ba9d3SPiotr Jasiukajtis * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
275b2ba9d3SPiotr Jasiukajtis * Use is subject to license terms.
285b2ba9d3SPiotr Jasiukajtis */
295b2ba9d3SPiotr Jasiukajtis
30*a9d3dcd5SRichard Lowe #pragma weak __feclearexcept = feclearexcept
31*a9d3dcd5SRichard Lowe #pragma weak __feraiseexcept = feraiseexcept
32*a9d3dcd5SRichard Lowe #pragma weak __fetestexcept = fetestexcept
33*a9d3dcd5SRichard Lowe #pragma weak __fegetexceptflag = fegetexceptflag
34*a9d3dcd5SRichard Lowe #pragma weak __fesetexceptflag = fesetexceptflag
355b2ba9d3SPiotr Jasiukajtis
36*a9d3dcd5SRichard Lowe #pragma weak feclearexcept96 = feclearexcept
37*a9d3dcd5SRichard Lowe #pragma weak feraiseexcept96 = feraiseexcept
38*a9d3dcd5SRichard Lowe #pragma weak fetestexcept96 = fetestexcept
39*a9d3dcd5SRichard Lowe #pragma weak fegetexceptflag96 = fegetexceptflag
40*a9d3dcd5SRichard Lowe #pragma weak fesetexceptflag96 = fesetexceptflag
415b2ba9d3SPiotr Jasiukajtis
425b2ba9d3SPiotr Jasiukajtis #include <fenv.h>
435b2ba9d3SPiotr Jasiukajtis #include <sys/ieeefp.h>
445b2ba9d3SPiotr Jasiukajtis #include <ucontext.h>
455b2ba9d3SPiotr Jasiukajtis #include <thread.h>
465b2ba9d3SPiotr Jasiukajtis #include "fex_handler.h"
475b2ba9d3SPiotr Jasiukajtis #include "fenv_inlines.h"
485b2ba9d3SPiotr Jasiukajtis
495b2ba9d3SPiotr Jasiukajtis
feclearexcept(int e)505b2ba9d3SPiotr Jasiukajtis int feclearexcept(int e)
515b2ba9d3SPiotr Jasiukajtis {
525b2ba9d3SPiotr Jasiukajtis unsigned long fsr;
535b2ba9d3SPiotr Jasiukajtis
545b2ba9d3SPiotr Jasiukajtis __fenv_getfsr(&fsr);
555b2ba9d3SPiotr Jasiukajtis __fenv_set_ex(fsr, __fenv_get_ex(fsr) & ~e);
565b2ba9d3SPiotr Jasiukajtis __fenv_setfsr(&fsr);
575b2ba9d3SPiotr Jasiukajtis if (fex_get_log())
585b2ba9d3SPiotr Jasiukajtis __fex_update_te();
595b2ba9d3SPiotr Jasiukajtis return 0;
605b2ba9d3SPiotr Jasiukajtis }
615b2ba9d3SPiotr Jasiukajtis
625b2ba9d3SPiotr Jasiukajtis /*
635b2ba9d3SPiotr Jasiukajtis * note - __fex_hdlr depends on fetestexcept following feraiseexcept
645b2ba9d3SPiotr Jasiukajtis */
feraiseexcept(int e)655b2ba9d3SPiotr Jasiukajtis int feraiseexcept(int e)
665b2ba9d3SPiotr Jasiukajtis {
675b2ba9d3SPiotr Jasiukajtis volatile double t;
685b2ba9d3SPiotr Jasiukajtis unsigned long fsr;
695b2ba9d3SPiotr Jasiukajtis
705b2ba9d3SPiotr Jasiukajtis if (e & FE_INVALID) {
715b2ba9d3SPiotr Jasiukajtis t = 0.0;
725b2ba9d3SPiotr Jasiukajtis t /= 0.0;
735b2ba9d3SPiotr Jasiukajtis }
745b2ba9d3SPiotr Jasiukajtis if (e & FE_DIVBYZERO) {
755b2ba9d3SPiotr Jasiukajtis t = 1.0e300;
765b2ba9d3SPiotr Jasiukajtis t /= 0.0;
775b2ba9d3SPiotr Jasiukajtis }
785b2ba9d3SPiotr Jasiukajtis if (e & FE_OVERFLOW) {
795b2ba9d3SPiotr Jasiukajtis /* if overflow is not trapped, avoid raising inexact */
805b2ba9d3SPiotr Jasiukajtis __fenv_getfsr(&fsr);
815b2ba9d3SPiotr Jasiukajtis if (!(__fenv_get_te(fsr) & (1 << fp_trap_overflow))) {
825b2ba9d3SPiotr Jasiukajtis __fenv_set_ex(fsr, __fenv_get_ex(fsr) | FE_OVERFLOW);
835b2ba9d3SPiotr Jasiukajtis __fenv_setfsr(&fsr);
845b2ba9d3SPiotr Jasiukajtis }
855b2ba9d3SPiotr Jasiukajtis else {
865b2ba9d3SPiotr Jasiukajtis t = 1.0e300;
875b2ba9d3SPiotr Jasiukajtis t *= 1.0e300;
885b2ba9d3SPiotr Jasiukajtis }
895b2ba9d3SPiotr Jasiukajtis }
905b2ba9d3SPiotr Jasiukajtis if (e & FE_UNDERFLOW) {
915b2ba9d3SPiotr Jasiukajtis /* if underflow is not trapped, avoid raising inexact */
925b2ba9d3SPiotr Jasiukajtis __fenv_getfsr(&fsr);
935b2ba9d3SPiotr Jasiukajtis if (!(__fenv_get_te(fsr) & (1 << fp_trap_underflow))) {
945b2ba9d3SPiotr Jasiukajtis __fenv_set_ex(fsr, __fenv_get_ex(fsr) | FE_UNDERFLOW);
955b2ba9d3SPiotr Jasiukajtis __fenv_setfsr(&fsr);
965b2ba9d3SPiotr Jasiukajtis }
975b2ba9d3SPiotr Jasiukajtis else {
985b2ba9d3SPiotr Jasiukajtis t = 1.0e-307;
995b2ba9d3SPiotr Jasiukajtis t -= 1.001e-307;
1005b2ba9d3SPiotr Jasiukajtis }
1015b2ba9d3SPiotr Jasiukajtis }
1025b2ba9d3SPiotr Jasiukajtis if (e & FE_INEXACT) {
1035b2ba9d3SPiotr Jasiukajtis t = 1.0e300;
1045b2ba9d3SPiotr Jasiukajtis t += 1.0e-307;
1055b2ba9d3SPiotr Jasiukajtis }
1065b2ba9d3SPiotr Jasiukajtis return 0;
1075b2ba9d3SPiotr Jasiukajtis }
1085b2ba9d3SPiotr Jasiukajtis
fetestexcept(int e)1095b2ba9d3SPiotr Jasiukajtis int fetestexcept(int e)
1105b2ba9d3SPiotr Jasiukajtis {
1115b2ba9d3SPiotr Jasiukajtis unsigned long fsr;
1125b2ba9d3SPiotr Jasiukajtis
1135b2ba9d3SPiotr Jasiukajtis __fenv_getfsr(&fsr);
1145b2ba9d3SPiotr Jasiukajtis return (int)__fenv_get_ex(fsr) & e;
1155b2ba9d3SPiotr Jasiukajtis }
1165b2ba9d3SPiotr Jasiukajtis
fegetexceptflag(fexcept_t * p,int e)1175b2ba9d3SPiotr Jasiukajtis int fegetexceptflag(fexcept_t *p, int e)
1185b2ba9d3SPiotr Jasiukajtis {
1195b2ba9d3SPiotr Jasiukajtis unsigned long fsr;
1205b2ba9d3SPiotr Jasiukajtis
1215b2ba9d3SPiotr Jasiukajtis __fenv_getfsr(&fsr);
1225b2ba9d3SPiotr Jasiukajtis *p = (int)__fenv_get_ex(fsr) & e;
1235b2ba9d3SPiotr Jasiukajtis return 0;
1245b2ba9d3SPiotr Jasiukajtis }
1255b2ba9d3SPiotr Jasiukajtis
fesetexceptflag(const fexcept_t * p,int e)1265b2ba9d3SPiotr Jasiukajtis int fesetexceptflag(const fexcept_t *p, int e)
1275b2ba9d3SPiotr Jasiukajtis {
1285b2ba9d3SPiotr Jasiukajtis unsigned long fsr;
1295b2ba9d3SPiotr Jasiukajtis
1305b2ba9d3SPiotr Jasiukajtis __fenv_getfsr(&fsr);
1315b2ba9d3SPiotr Jasiukajtis __fenv_set_ex(fsr, (((int)__fenv_get_ex(fsr) & ~e) | (*p & e)) &
1325b2ba9d3SPiotr Jasiukajtis FE_ALL_EXCEPT);
1335b2ba9d3SPiotr Jasiukajtis __fenv_setfsr(&fsr);
1345b2ba9d3SPiotr Jasiukajtis if (fex_get_log())
1355b2ba9d3SPiotr Jasiukajtis __fex_update_te();
1365b2ba9d3SPiotr Jasiukajtis return 0;
1375b2ba9d3SPiotr Jasiukajtis }
138