4 * The contents of this file are subject to the terms of the
5 * Common Development and Distribution License (the "License").
6 * You may not use this file except in compliance with the License.
8 * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
9 * or http://www.opensolaris.org/os/licensing.
10 * See the License for the specific language governing permissions
11 * and limitations under the License.
13 * When distributing Covered Code, include this CDDL HEADER in each
14 * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
15 * If applicable, add the following below this CDDL HEADER, with the
16 * fields enclosed by brackets "[]" replaced with your own identifying
17 * information: Portions Copyright [yyyy] [name of copyright owner]
23 * Copyright 2008 Sun Microsystems, Inc. All rights reserved.
24 * Use is subject to license terms.
28 #include <sys/types.h>
29 #include "base_conversion.h"
30 #include <sys/isa_defs.h>
33 * Miscellaneous support routines used in base conversion
41 { 0x00000000u
, 0x00100000u
},
42 { 0x00000001u
, 0x7ff00000u
}
44 { 0x00100000u
, 0x00000000u
},
45 { 0x7ff00000u
, 0x00000001u
}
49 #define minnormal C[0].d
50 #define signalingnan C[1].d
52 /* raise the floating point exceptions indicated by ef */
54 __base_conversion_set_exception(fp_exception_field_type ef
)
57 volatile double tstored __unused
;
59 if (ef
== (1 << fp_inexact
)) {
60 t
= 9.999999962747097015E-1;
62 * 28 sig bits so product isn't inexact in extended
63 * accumulator, causing two inexact traps.
65 } else if ((ef
& (1 << fp_invalid
)) != 0) {
67 } else if ((ef
& (1 << fp_overflow
)) != 0) {
68 t
= 4.149515553422842866E+180;
70 * 28 sig bits so product isn't inexact in extended
71 * accumulator, causing inexact trap prior to overflow trap
74 } else if ((ef
& (1 << fp_underflow
)) != 0) {
79 /* Storage forces exception */
87 * The following routine is no longer used in libc, but we have
88 * to leave it for now because it's still used by Sun's old Fortran
89 * runtime libraries. Today this is a bug; in the days of SunOS 4.x,
90 * when the relevant design decisions were made, it was a feature.
93 __class_quadruple(quadruple
*x
)
95 quadruple_equivalence kluge
;
98 if (kluge
.f
.msw
.exponent
== 0) { /* 0 or sub */
99 if ((kluge
.f
.msw
.significand
== 0) &&
100 (kluge
.f
.significand2
== 0) &&
101 (kluge
.f
.significand3
== 0) &&
102 (kluge
.f
.significand4
== 0))
105 return (fp_subnormal
);
106 } else if (kluge
.f
.msw
.exponent
== 0x7fff) { /* inf or nan */
107 if ((kluge
.f
.msw
.significand
== 0) &&
108 (kluge
.f
.significand2
== 0) &&
109 (kluge
.f
.significand3
== 0) &&
110 (kluge
.f
.significand4
== 0))
111 return (fp_infinity
);
112 else if ((kluge
.f
.msw
.significand
& 0xffff) >=
113 (unsigned int)0x8000)
116 return (fp_signaling
);