5 * Common Development and Distribution License (the "License").
6 * You may not use this file except in compliance with the License.
7 *
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.
12 *
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]
18 *
19 * CDDL HEADER END
20 */
21
22 /*
23 * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
24 */
25 /*
26 * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
27 * Use is subject to license terms.
28 */
29
30 #pragma weak __nextafterf = nextafterf
31
32 #include "libm.h"
33
34 float
35 nextafterf(float x, float y) {
36 float w;
37 int *pw = (int *) &w;
38 int *px = (int *) &x;
39 int *py = (int *) &y;
40 int ix, iy, iz;
41
42 ix = px[0];
43 iy = py[0];
44 if ((ix & ~0x80000000) > 0x7f800000)
45 return (x * y); /* + -> * for Cheetah */
46 if ((iy & ~0x80000000) > 0x7f800000)
47 return (y * x); /* + -> * for Cheetah */
48 if (ix == iy || (ix | iy) == 0x80000000)
49 return (y); /* C99 requirement */
50 if ((ix & ~0x80000000) == 0)
51 iz = 1 | (iy & 0x80000000);
52 else if (ix > 0) {
53 if (ix > iy)
54 iz = ix - 1;
55 else
56 iz = ix + 1;
57 } else {
58 if (iy < 0 && ix < iy)
59 iz = ix + 1;
60 else
61 iz = ix - 1;
62 }
63 pw[0] = iz;
64 ix = iz & 0x7f800000;
65 if (ix == 0x7f800000) {
66 /* raise overflow */
67 volatile float t;
68
69 *(int *) &t = 0x7f7fffff;
70 t *= t;
71 } else if (ix == 0) {
72 /* raise underflow */
73 volatile float t;
74
75 *(int *) &t = 0x00800000;
76 t *= t;
77 }
78 return (w);
79 }
|
5 * Common Development and Distribution License (the "License").
6 * You may not use this file except in compliance with the License.
7 *
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.
12 *
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]
18 *
19 * CDDL HEADER END
20 */
21
22 /*
23 * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
24 */
25
26 /*
27 * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
28 * Use is subject to license terms.
29 */
30
31 #pragma weak __nextafterf = nextafterf
32
33 #include "libm.h"
34
35 float
36 nextafterf(float x, float y)
37 {
38 float w;
39 int *pw = (int *)&w;
40 int *px = (int *)&x;
41 int *py = (int *)&y;
42 int ix, iy, iz;
43
44 ix = px[0];
45 iy = py[0];
46
47 if ((ix & ~0x80000000) > 0x7f800000)
48 return (x * y); /* + -> * for Cheetah */
49
50 if ((iy & ~0x80000000) > 0x7f800000)
51 return (y * x); /* + -> * for Cheetah */
52
53 if (ix == iy || (ix | iy) == 0x80000000)
54 return (y); /* C99 requirement */
55
56 if ((ix & ~0x80000000) == 0) {
57 iz = 1 | (iy & 0x80000000);
58 } else if (ix > 0) {
59 if (ix > iy)
60 iz = ix - 1;
61 else
62 iz = ix + 1;
63 } else {
64 if (iy < 0 && ix < iy)
65 iz = ix + 1;
66 else
67 iz = ix - 1;
68 }
69
70 pw[0] = iz;
71 ix = iz & 0x7f800000;
72
73 if (ix == 0x7f800000) {
74 /* raise overflow */
75 volatile float t;
76
77 *(int *)&t = 0x7f7fffff;
78 t *= t;
79 } else if (ix == 0) {
80 /* raise underflow */
81 volatile float t;
82
83 *(int *)&t = 0x00800000;
84 t *= t;
85 }
86
87 return (w);
88 }
|