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 __ceilf = ceilf
31*a9d3dcd5SRichard Lowe #pragma weak __floorf = floorf
325b2ba9d3SPiotr Jasiukajtis
335b2ba9d3SPiotr Jasiukajtis /* INDENT OFF */
345b2ba9d3SPiotr Jasiukajtis /*
355b2ba9d3SPiotr Jasiukajtis * ceilf(x) return the biggest integral value (in float) below x
365b2ba9d3SPiotr Jasiukajtis * floorf(x) return the least integral value (in float) above x
375b2ba9d3SPiotr Jasiukajtis *
385b2ba9d3SPiotr Jasiukajtis * NOTE: ceilf(x) and floorf(x) return result
395b2ba9d3SPiotr Jasiukajtis * with the same sign as x's, including 0.0F.
405b2ba9d3SPiotr Jasiukajtis */
415b2ba9d3SPiotr Jasiukajtis
425b2ba9d3SPiotr Jasiukajtis #include "libm.h"
435b2ba9d3SPiotr Jasiukajtis
445b2ba9d3SPiotr Jasiukajtis static const float xf[] = {
455b2ba9d3SPiotr Jasiukajtis /* ZEROF */ 0.0f,
465b2ba9d3SPiotr Jasiukajtis /* ONEF */ 1.0f,
475b2ba9d3SPiotr Jasiukajtis /* MONEF */ -1.0f,
485b2ba9d3SPiotr Jasiukajtis /* HUGEF */ 1.0e30f,
495b2ba9d3SPiotr Jasiukajtis };
505b2ba9d3SPiotr Jasiukajtis
515b2ba9d3SPiotr Jasiukajtis #define ZEROF xf[0]
525b2ba9d3SPiotr Jasiukajtis #define ONEF xf[1]
535b2ba9d3SPiotr Jasiukajtis #define MONEF xf[2]
545b2ba9d3SPiotr Jasiukajtis #define HUGEF xf[3]
555b2ba9d3SPiotr Jasiukajtis /* INDENT ON */
565b2ba9d3SPiotr Jasiukajtis
575b2ba9d3SPiotr Jasiukajtis float
ceilf(float x)585b2ba9d3SPiotr Jasiukajtis ceilf(float x) {
595b2ba9d3SPiotr Jasiukajtis volatile float dummy;
605b2ba9d3SPiotr Jasiukajtis int hx, k, j, ix;
615b2ba9d3SPiotr Jasiukajtis
625b2ba9d3SPiotr Jasiukajtis hx = *(int *) &x;
635b2ba9d3SPiotr Jasiukajtis ix = hx & ~0x80000000;
645b2ba9d3SPiotr Jasiukajtis k = ix >> 23;
655b2ba9d3SPiotr Jasiukajtis if (((k - 127) ^ (k - 150)) < 0) {
665b2ba9d3SPiotr Jasiukajtis k = (1 << (150 - k)) - 1;
675b2ba9d3SPiotr Jasiukajtis if ((k & hx) != 0)
685b2ba9d3SPiotr Jasiukajtis dummy = HUGEF + x; /* raise inexact */
695b2ba9d3SPiotr Jasiukajtis j = k & (~(hx >> 31));
705b2ba9d3SPiotr Jasiukajtis *(int *) &x = (hx + j) & ~k;
715b2ba9d3SPiotr Jasiukajtis return (x);
725b2ba9d3SPiotr Jasiukajtis } else if (k <= 126) {
735b2ba9d3SPiotr Jasiukajtis dummy = HUGEF + x;
745b2ba9d3SPiotr Jasiukajtis if (hx > 0)
755b2ba9d3SPiotr Jasiukajtis return (ONEF);
765b2ba9d3SPiotr Jasiukajtis else if (ix == 0)
775b2ba9d3SPiotr Jasiukajtis return (x);
785b2ba9d3SPiotr Jasiukajtis else
795b2ba9d3SPiotr Jasiukajtis return (-ZEROF);
805b2ba9d3SPiotr Jasiukajtis } else
815b2ba9d3SPiotr Jasiukajtis /* signal invalid if x is a SNaN */
825b2ba9d3SPiotr Jasiukajtis return (x * ONEF); /* +0 -> *1 for Cheetah */
835b2ba9d3SPiotr Jasiukajtis }
845b2ba9d3SPiotr Jasiukajtis
855b2ba9d3SPiotr Jasiukajtis float
floorf(float x)865b2ba9d3SPiotr Jasiukajtis floorf(float x) {
875b2ba9d3SPiotr Jasiukajtis volatile float dummy;
885b2ba9d3SPiotr Jasiukajtis int hx, k, j, ix;
895b2ba9d3SPiotr Jasiukajtis
905b2ba9d3SPiotr Jasiukajtis hx = *(int *) &x;
915b2ba9d3SPiotr Jasiukajtis ix = hx & ~0x80000000;
925b2ba9d3SPiotr Jasiukajtis k = ix >> 23;
935b2ba9d3SPiotr Jasiukajtis if (((k - 127) ^ (k - 150)) < 0) {
945b2ba9d3SPiotr Jasiukajtis k = (1 << (150 - k)) - 1;
955b2ba9d3SPiotr Jasiukajtis if ((k & hx) != 0)
965b2ba9d3SPiotr Jasiukajtis dummy = HUGEF + x; /* raise inexact */
975b2ba9d3SPiotr Jasiukajtis j = k & (hx >> 31);
985b2ba9d3SPiotr Jasiukajtis *(int *) &x = (hx + j) & ~k;
995b2ba9d3SPiotr Jasiukajtis return (x);
1005b2ba9d3SPiotr Jasiukajtis } else if (k <= 126) {
1015b2ba9d3SPiotr Jasiukajtis dummy = HUGEF + x;
1025b2ba9d3SPiotr Jasiukajtis if (hx > 0)
1035b2ba9d3SPiotr Jasiukajtis return (ZEROF);
1045b2ba9d3SPiotr Jasiukajtis else if (ix == 0)
1055b2ba9d3SPiotr Jasiukajtis return (x);
1065b2ba9d3SPiotr Jasiukajtis else
1075b2ba9d3SPiotr Jasiukajtis return (MONEF);
1085b2ba9d3SPiotr Jasiukajtis } else
1095b2ba9d3SPiotr Jasiukajtis /* signal invalid if x is a SNaN */
1105b2ba9d3SPiotr Jasiukajtis return (x * ONEF); /* +0 -> *1 for Cheetah */
1115b2ba9d3SPiotr Jasiukajtis }
112