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
305b2ba9d3SPiotr Jasiukajtis /*
315b2ba9d3SPiotr Jasiukajtis * ceill(x) return the biggest integral value below x
325b2ba9d3SPiotr Jasiukajtis * floorl(x) return the least integral value above x
335b2ba9d3SPiotr Jasiukajtis *
345b2ba9d3SPiotr Jasiukajtis * NOTE: aintl(x), anintl(x), ceill(x), floorl(x), and rintl(x) return result
355b2ba9d3SPiotr Jasiukajtis * with the same sign as x's, including 0.0.
365b2ba9d3SPiotr Jasiukajtis */
375b2ba9d3SPiotr Jasiukajtis
38*a9d3dcd5SRichard Lowe #pragma weak __ceill = ceill
39*a9d3dcd5SRichard Lowe #pragma weak __floorl = floorl
405b2ba9d3SPiotr Jasiukajtis
415b2ba9d3SPiotr Jasiukajtis #include "libm.h"
425b2ba9d3SPiotr Jasiukajtis #include "longdouble.h"
435b2ba9d3SPiotr Jasiukajtis
445b2ba9d3SPiotr Jasiukajtis static const long double qone = 1.0L;
455b2ba9d3SPiotr Jasiukajtis
465b2ba9d3SPiotr Jasiukajtis long double
ceill(long double x)475b2ba9d3SPiotr Jasiukajtis ceill(long double x) {
485b2ba9d3SPiotr Jasiukajtis long double t;
495b2ba9d3SPiotr Jasiukajtis
505b2ba9d3SPiotr Jasiukajtis if (!finitel(x))
515b2ba9d3SPiotr Jasiukajtis return (x + x);
525b2ba9d3SPiotr Jasiukajtis t = rintl(x);
535b2ba9d3SPiotr Jasiukajtis if (t >= x) /* already ceil(x) */
545b2ba9d3SPiotr Jasiukajtis return (t);
555b2ba9d3SPiotr Jasiukajtis else /* t < x case: return t+1 */
565b2ba9d3SPiotr Jasiukajtis return (copysignl(t + qone, x));
575b2ba9d3SPiotr Jasiukajtis }
585b2ba9d3SPiotr Jasiukajtis
595b2ba9d3SPiotr Jasiukajtis long double
floorl(long double x)605b2ba9d3SPiotr Jasiukajtis floorl(long double x) {
615b2ba9d3SPiotr Jasiukajtis long double t;
625b2ba9d3SPiotr Jasiukajtis
635b2ba9d3SPiotr Jasiukajtis if (!finitel(x))
645b2ba9d3SPiotr Jasiukajtis return (x + x);
655b2ba9d3SPiotr Jasiukajtis t = rintl(x);
665b2ba9d3SPiotr Jasiukajtis if (t <= x)
675b2ba9d3SPiotr Jasiukajtis return (t); /* already floor(x) */
685b2ba9d3SPiotr Jasiukajtis else /* x < t case: return t-1 */
695b2ba9d3SPiotr Jasiukajtis return (copysignl(t - qone, x));
705b2ba9d3SPiotr Jasiukajtis }
71