1*0Sstevel@tonic-gate /*
2*0Sstevel@tonic-gate * CDDL HEADER START
3*0Sstevel@tonic-gate *
4*0Sstevel@tonic-gate * The contents of this file are subject to the terms of the
5*0Sstevel@tonic-gate * Common Development and Distribution License, Version 1.0 only
6*0Sstevel@tonic-gate * (the "License"). You may not use this file except in compliance
7*0Sstevel@tonic-gate * with the License.
8*0Sstevel@tonic-gate *
9*0Sstevel@tonic-gate * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
10*0Sstevel@tonic-gate * or http://www.opensolaris.org/os/licensing.
11*0Sstevel@tonic-gate * See the License for the specific language governing permissions
12*0Sstevel@tonic-gate * and limitations under the License.
13*0Sstevel@tonic-gate *
14*0Sstevel@tonic-gate * When distributing Covered Code, include this CDDL HEADER in each
15*0Sstevel@tonic-gate * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
16*0Sstevel@tonic-gate * If applicable, add the following below this CDDL HEADER, with the
17*0Sstevel@tonic-gate * fields enclosed by brackets "[]" replaced with your own identifying
18*0Sstevel@tonic-gate * information: Portions Copyright [yyyy] [name of copyright owner]
19*0Sstevel@tonic-gate *
20*0Sstevel@tonic-gate * CDDL HEADER END
21*0Sstevel@tonic-gate */
22*0Sstevel@tonic-gate /*
23*0Sstevel@tonic-gate * Copyright (c) 1994-1997, by Sun Microsystems, Inc.
24*0Sstevel@tonic-gate * All rights reserved.
25*0Sstevel@tonic-gate */
26*0Sstevel@tonic-gate
27*0Sstevel@tonic-gate #pragma ident "%Z%%M% %I% %E% SMI"
28*0Sstevel@tonic-gate
29*0Sstevel@tonic-gate #include "quad.h"
30*0Sstevel@tonic-gate
31*0Sstevel@tonic-gate /*
32*0Sstevel@tonic-gate * _Qp_qtoux(x) returns (unsigned long)*x.
33*0Sstevel@tonic-gate */
34*0Sstevel@tonic-gate unsigned long
_Qp_qtoux(const union longdouble * x)35*0Sstevel@tonic-gate _Qp_qtoux(const union longdouble *x)
36*0Sstevel@tonic-gate {
37*0Sstevel@tonic-gate union longdouble z;
38*0Sstevel@tonic-gate unsigned long i, round;
39*0Sstevel@tonic-gate unsigned int xm, fsr;
40*0Sstevel@tonic-gate
41*0Sstevel@tonic-gate xm = x->l.msw & 0x7fffffff;
42*0Sstevel@tonic-gate
43*0Sstevel@tonic-gate __quad_getfsrp(&fsr);
44*0Sstevel@tonic-gate
45*0Sstevel@tonic-gate /* handle nan, inf, and out-of-range cases */
46*0Sstevel@tonic-gate if (xm >= 0x403e0000) {
47*0Sstevel@tonic-gate if (x->l.msw < 0x403f0000) {
48*0Sstevel@tonic-gate i = 0x8000000000000000ul |
49*0Sstevel@tonic-gate ((long) (xm & 0xffff) << 47) |
50*0Sstevel@tonic-gate ((long) x->l.frac2 << 15) | (x->l.frac3 >> 17);
51*0Sstevel@tonic-gate if ((x->l.frac3 & 0x1ffff) | x->l.frac4) {
52*0Sstevel@tonic-gate /* signal inexact */
53*0Sstevel@tonic-gate if (fsr & FSR_NXM) {
54*0Sstevel@tonic-gate /* z = x - 2^63 */
55*0Sstevel@tonic-gate if (xm & 0xffff ||
56*0Sstevel@tonic-gate x->l.frac2 & 0xffff0000) {
57*0Sstevel@tonic-gate z.l.msw = xm & 0xffff;
58*0Sstevel@tonic-gate z.l.frac2 = x->l.frac2;
59*0Sstevel@tonic-gate z.l.frac3 = x->l.frac3;
60*0Sstevel@tonic-gate z.l.frac4 = x->l.frac4;
61*0Sstevel@tonic-gate } else if (x->l.frac2 & 0xffff ||
62*0Sstevel@tonic-gate x->l.frac3 & 0xffff0000) {
63*0Sstevel@tonic-gate z.l.msw = x->l.frac2;
64*0Sstevel@tonic-gate z.l.frac2 = x->l.frac3;
65*0Sstevel@tonic-gate z.l.frac3 = x->l.frac4;
66*0Sstevel@tonic-gate z.l.frac4 = 0;
67*0Sstevel@tonic-gate } else if (x->l.frac3 & 0xffff ||
68*0Sstevel@tonic-gate x->l.frac4 & 0xffff0000) {
69*0Sstevel@tonic-gate z.l.msw = x->l.frac3;
70*0Sstevel@tonic-gate z.l.frac2 = x->l.frac4;
71*0Sstevel@tonic-gate z.l.frac3 = z.l.frac4 = 0;
72*0Sstevel@tonic-gate } else {
73*0Sstevel@tonic-gate z.l.msw = x->l.frac4;
74*0Sstevel@tonic-gate z.l.frac2 = z.l.frac3 =
75*0Sstevel@tonic-gate z.l.frac4 = 0;
76*0Sstevel@tonic-gate }
77*0Sstevel@tonic-gate xm = 0x403e;
78*0Sstevel@tonic-gate while ((z.l.msw & 0x10000) == 0) {
79*0Sstevel@tonic-gate z.l.msw = (z.l.msw << 1) |
80*0Sstevel@tonic-gate (z.l.frac2 >> 31);
81*0Sstevel@tonic-gate z.l.frac2 = (z.l.frac2 << 1) |
82*0Sstevel@tonic-gate (z.l.frac3 >> 31);
83*0Sstevel@tonic-gate z.l.frac3 = (z.l.frac3 << 1) |
84*0Sstevel@tonic-gate (z.l.frac4 >> 31);
85*0Sstevel@tonic-gate z.l.frac4 <<= 1;
86*0Sstevel@tonic-gate xm--;
87*0Sstevel@tonic-gate }
88*0Sstevel@tonic-gate z.l.msw |= (xm << 16);
89*0Sstevel@tonic-gate __quad_fqtox(&z, (long *)&i);
90*0Sstevel@tonic-gate i |= 0x8000000000000000ul;
91*0Sstevel@tonic-gate } else {
92*0Sstevel@tonic-gate fsr = (fsr & ~FSR_CEXC) | FSR_NXA |
93*0Sstevel@tonic-gate FSR_NXC;
94*0Sstevel@tonic-gate __quad_setfsrp(&fsr);
95*0Sstevel@tonic-gate }
96*0Sstevel@tonic-gate }
97*0Sstevel@tonic-gate return (i);
98*0Sstevel@tonic-gate }
99*0Sstevel@tonic-gate if (x->l.msw == 0xc03e0000 && x->l.frac2 == 0 &&
100*0Sstevel@tonic-gate (x->l.frac3 & 0xfffe0000) == 0) {
101*0Sstevel@tonic-gate /* return largest negative 64 bit int */
102*0Sstevel@tonic-gate i = 0x8000000000000000ul;
103*0Sstevel@tonic-gate if ((x->l.frac3 & 0x1ffff) | x->l.frac4) {
104*0Sstevel@tonic-gate /* signal inexact */
105*0Sstevel@tonic-gate if (fsr & FSR_NXM) {
106*0Sstevel@tonic-gate __quad_fqtox(x, (long *)&i);
107*0Sstevel@tonic-gate } else {
108*0Sstevel@tonic-gate fsr = (fsr & ~FSR_CEXC) | FSR_NXA |
109*0Sstevel@tonic-gate FSR_NXC;
110*0Sstevel@tonic-gate __quad_setfsrp(&fsr);
111*0Sstevel@tonic-gate }
112*0Sstevel@tonic-gate }
113*0Sstevel@tonic-gate return (i);
114*0Sstevel@tonic-gate }
115*0Sstevel@tonic-gate i = 0x7fffffffffffffffl;
116*0Sstevel@tonic-gate if (fsr & FSR_NVM) {
117*0Sstevel@tonic-gate __quad_fqtox(x, (long *)&i);
118*0Sstevel@tonic-gate } else {
119*0Sstevel@tonic-gate fsr = (fsr & ~FSR_CEXC) | FSR_NVA | FSR_NVC;
120*0Sstevel@tonic-gate __quad_setfsrp(&fsr);
121*0Sstevel@tonic-gate }
122*0Sstevel@tonic-gate return (i);
123*0Sstevel@tonic-gate }
124*0Sstevel@tonic-gate if (xm < 0x3fff0000) {
125*0Sstevel@tonic-gate i = 0l;
126*0Sstevel@tonic-gate if (xm | x->l.frac2 | x->l.frac3 | x->l.frac4) {
127*0Sstevel@tonic-gate /* signal inexact */
128*0Sstevel@tonic-gate if (fsr & FSR_NXM) {
129*0Sstevel@tonic-gate __quad_fqtox(x, (long *)&i);
130*0Sstevel@tonic-gate } else {
131*0Sstevel@tonic-gate fsr = (fsr & ~FSR_CEXC) | FSR_NXA | FSR_NXC;
132*0Sstevel@tonic-gate __quad_setfsrp(&fsr);
133*0Sstevel@tonic-gate }
134*0Sstevel@tonic-gate }
135*0Sstevel@tonic-gate return (i);
136*0Sstevel@tonic-gate }
137*0Sstevel@tonic-gate
138*0Sstevel@tonic-gate /* now x is in the range of long */
139*0Sstevel@tonic-gate i = 0x4000000000000000l | ((long) (xm & 0xffff) << 46) |
140*0Sstevel@tonic-gate ((long) x->l.frac2 << 14) | (x->l.frac3 >> 18);
141*0Sstevel@tonic-gate round = i & ((1l << (0x403d - (xm >> 16))) - 1);
142*0Sstevel@tonic-gate i >>= (0x403d - (xm >> 16));
143*0Sstevel@tonic-gate if (x->l.msw & 0x80000000)
144*0Sstevel@tonic-gate i = -i;
145*0Sstevel@tonic-gate if (round | (x->l.frac3 & 0x3ffff) | x->l.frac4) {
146*0Sstevel@tonic-gate /* signal inexact */
147*0Sstevel@tonic-gate if (fsr & FSR_NXM) {
148*0Sstevel@tonic-gate __quad_fqtox(x, (long *)&i);
149*0Sstevel@tonic-gate } else {
150*0Sstevel@tonic-gate fsr = (fsr & ~FSR_CEXC) | FSR_NXA | FSR_NXC;
151*0Sstevel@tonic-gate __quad_setfsrp(&fsr);
152*0Sstevel@tonic-gate }
153*0Sstevel@tonic-gate }
154*0Sstevel@tonic-gate return (i);
155*0Sstevel@tonic-gate }
156