1*25c28e83SPiotr Jasiukajtis /*
2*25c28e83SPiotr Jasiukajtis * CDDL HEADER START
3*25c28e83SPiotr Jasiukajtis *
4*25c28e83SPiotr Jasiukajtis * The contents of this file are subject to the terms of the
5*25c28e83SPiotr Jasiukajtis * Common Development and Distribution License (the "License").
6*25c28e83SPiotr Jasiukajtis * You may not use this file except in compliance with the License.
7*25c28e83SPiotr Jasiukajtis *
8*25c28e83SPiotr Jasiukajtis * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
9*25c28e83SPiotr Jasiukajtis * or http://www.opensolaris.org/os/licensing.
10*25c28e83SPiotr Jasiukajtis * See the License for the specific language governing permissions
11*25c28e83SPiotr Jasiukajtis * and limitations under the License.
12*25c28e83SPiotr Jasiukajtis *
13*25c28e83SPiotr Jasiukajtis * When distributing Covered Code, include this CDDL HEADER in each
14*25c28e83SPiotr Jasiukajtis * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
15*25c28e83SPiotr Jasiukajtis * If applicable, add the following below this CDDL HEADER, with the
16*25c28e83SPiotr Jasiukajtis * fields enclosed by brackets "[]" replaced with your own identifying
17*25c28e83SPiotr Jasiukajtis * information: Portions Copyright [yyyy] [name of copyright owner]
18*25c28e83SPiotr Jasiukajtis *
19*25c28e83SPiotr Jasiukajtis * CDDL HEADER END
20*25c28e83SPiotr Jasiukajtis */
21*25c28e83SPiotr Jasiukajtis
22*25c28e83SPiotr Jasiukajtis /*
23*25c28e83SPiotr Jasiukajtis * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
24*25c28e83SPiotr Jasiukajtis */
25*25c28e83SPiotr Jasiukajtis /*
26*25c28e83SPiotr Jasiukajtis * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
27*25c28e83SPiotr Jasiukajtis * Use is subject to license terms.
28*25c28e83SPiotr Jasiukajtis */
29*25c28e83SPiotr Jasiukajtis
30*25c28e83SPiotr Jasiukajtis #pragma weak frexpl = __frexpl
31*25c28e83SPiotr Jasiukajtis
32*25c28e83SPiotr Jasiukajtis #include "libm.h"
33*25c28e83SPiotr Jasiukajtis
34*25c28e83SPiotr Jasiukajtis #if defined(__sparc)
35*25c28e83SPiotr Jasiukajtis
36*25c28e83SPiotr Jasiukajtis long double
__frexpl(long double x,int * exp)37*25c28e83SPiotr Jasiukajtis __frexpl(long double x, int *exp) {
38*25c28e83SPiotr Jasiukajtis union {
39*25c28e83SPiotr Jasiukajtis unsigned i[4];
40*25c28e83SPiotr Jasiukajtis long double q;
41*25c28e83SPiotr Jasiukajtis } xx;
42*25c28e83SPiotr Jasiukajtis unsigned hx;
43*25c28e83SPiotr Jasiukajtis int e, s;
44*25c28e83SPiotr Jasiukajtis
45*25c28e83SPiotr Jasiukajtis xx.q = x;
46*25c28e83SPiotr Jasiukajtis hx = xx.i[0] & ~0x80000000;
47*25c28e83SPiotr Jasiukajtis
48*25c28e83SPiotr Jasiukajtis if (hx >= 0x7fff0000) { /* x is infinite or NaN */
49*25c28e83SPiotr Jasiukajtis *exp = 0;
50*25c28e83SPiotr Jasiukajtis return (x);
51*25c28e83SPiotr Jasiukajtis }
52*25c28e83SPiotr Jasiukajtis
53*25c28e83SPiotr Jasiukajtis e = 0;
54*25c28e83SPiotr Jasiukajtis if (hx < 0x00010000) { /* x is subnormal or zero */
55*25c28e83SPiotr Jasiukajtis if ((hx | xx.i[1] | xx.i[2] | xx.i[3]) == 0) {
56*25c28e83SPiotr Jasiukajtis *exp = 0;
57*25c28e83SPiotr Jasiukajtis return (x);
58*25c28e83SPiotr Jasiukajtis }
59*25c28e83SPiotr Jasiukajtis
60*25c28e83SPiotr Jasiukajtis /* normalize x */
61*25c28e83SPiotr Jasiukajtis s = xx.i[0] & 0x80000000;
62*25c28e83SPiotr Jasiukajtis while ((hx | (xx.i[1] & 0xffff0000)) == 0) {
63*25c28e83SPiotr Jasiukajtis hx = xx.i[1];
64*25c28e83SPiotr Jasiukajtis xx.i[1] = xx.i[2];
65*25c28e83SPiotr Jasiukajtis xx.i[2] = xx.i[3];
66*25c28e83SPiotr Jasiukajtis xx.i[3] = 0;
67*25c28e83SPiotr Jasiukajtis e -= 32;
68*25c28e83SPiotr Jasiukajtis }
69*25c28e83SPiotr Jasiukajtis while (hx < 0x10000) {
70*25c28e83SPiotr Jasiukajtis hx = (hx << 1) | (xx.i[1] >> 31);
71*25c28e83SPiotr Jasiukajtis xx.i[1] = (xx.i[1] << 1) | (xx.i[2] >> 31);
72*25c28e83SPiotr Jasiukajtis xx.i[2] = (xx.i[2] << 1) | (xx.i[3] >> 31);
73*25c28e83SPiotr Jasiukajtis xx.i[3] <<= 1;
74*25c28e83SPiotr Jasiukajtis e--;
75*25c28e83SPiotr Jasiukajtis }
76*25c28e83SPiotr Jasiukajtis xx.i[0] = s | hx;
77*25c28e83SPiotr Jasiukajtis }
78*25c28e83SPiotr Jasiukajtis
79*25c28e83SPiotr Jasiukajtis /* now xx.q is normal */
80*25c28e83SPiotr Jasiukajtis xx.i[0] = (xx.i[0] & ~0x7fff0000) | 0x3ffe0000;
81*25c28e83SPiotr Jasiukajtis *exp = e + (hx >> 16) - 0x3ffe;
82*25c28e83SPiotr Jasiukajtis return (xx.q);
83*25c28e83SPiotr Jasiukajtis }
84*25c28e83SPiotr Jasiukajtis
85*25c28e83SPiotr Jasiukajtis #elif defined(__x86)
86*25c28e83SPiotr Jasiukajtis
87*25c28e83SPiotr Jasiukajtis long double
__frexpl(long double x,int * exp)88*25c28e83SPiotr Jasiukajtis __frexpl(long double x, int *exp) {
89*25c28e83SPiotr Jasiukajtis union {
90*25c28e83SPiotr Jasiukajtis unsigned i[3];
91*25c28e83SPiotr Jasiukajtis long double e;
92*25c28e83SPiotr Jasiukajtis } xx;
93*25c28e83SPiotr Jasiukajtis unsigned hx;
94*25c28e83SPiotr Jasiukajtis int e;
95*25c28e83SPiotr Jasiukajtis
96*25c28e83SPiotr Jasiukajtis xx.e = x;
97*25c28e83SPiotr Jasiukajtis hx = xx.i[2] & 0x7fff;
98*25c28e83SPiotr Jasiukajtis
99*25c28e83SPiotr Jasiukajtis if (hx >= 0x7fff) { /* x is infinite or NaN */
100*25c28e83SPiotr Jasiukajtis *exp = 0;
101*25c28e83SPiotr Jasiukajtis return (x);
102*25c28e83SPiotr Jasiukajtis }
103*25c28e83SPiotr Jasiukajtis
104*25c28e83SPiotr Jasiukajtis e = 0;
105*25c28e83SPiotr Jasiukajtis if (hx < 0x0001) { /* x is subnormal or zero */
106*25c28e83SPiotr Jasiukajtis if ((xx.i[0] | xx.i[1]) == 0) {
107*25c28e83SPiotr Jasiukajtis *exp = 0;
108*25c28e83SPiotr Jasiukajtis return (x);
109*25c28e83SPiotr Jasiukajtis }
110*25c28e83SPiotr Jasiukajtis
111*25c28e83SPiotr Jasiukajtis /* normalize x */
112*25c28e83SPiotr Jasiukajtis xx.e *= 18446744073709551616.0L; /* 2^64 */
113*25c28e83SPiotr Jasiukajtis hx = xx.i[2] & 0x7fff;
114*25c28e83SPiotr Jasiukajtis e = -64;
115*25c28e83SPiotr Jasiukajtis }
116*25c28e83SPiotr Jasiukajtis
117*25c28e83SPiotr Jasiukajtis /* now xx.e is normal */
118*25c28e83SPiotr Jasiukajtis xx.i[2] = (xx.i[2] & 0x8000) | 0x3ffe;
119*25c28e83SPiotr Jasiukajtis *exp = e + hx - 0x3ffe;
120*25c28e83SPiotr Jasiukajtis return (xx.e);
121*25c28e83SPiotr Jasiukajtis }
122*25c28e83SPiotr Jasiukajtis
123*25c28e83SPiotr Jasiukajtis #else
124*25c28e83SPiotr Jasiukajtis #error Unknown architecture
125*25c28e83SPiotr Jasiukajtis #endif
126