25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * CDDL HEADER START
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * The contents of this file are subject to the terms of the
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Common Development and Distribution License (the "License").
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * You may not use this file except in compliance with the License.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * or http://www.opensolaris.org/os/licensing.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * See the License for the specific language governing permissions
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * and limitations under the License.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * When distributing Covered Code, include this CDDL HEADER in each
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * If applicable, add the following below this CDDL HEADER, with the
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * fields enclosed by brackets "[]" replaced with your own identifying
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * information: Portions Copyright [yyyy] [name of copyright owner]
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * CDDL HEADER END
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * Use is subject to license terms.
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis/* INDENT OFF */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtisstatic const int
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtisstatic const long double zero = 0.0L, one = 1.0L;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis/* INDENT ON */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis * On entrance: *quo is initialized to 0, x finite and y non-zero & ordered
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtisstatic long double
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtisfmodquol(long double x, long double y, int *quo) {
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis long double a, b;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis unsigned x1, x2, x3, y1, y2, y3, z1, z2, z3;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis /* determine ix = ilogbl(x) */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis while (x0 == 0) {
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis /* determine iy = ilogbl(y) */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis while (y0 == 0) {
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis /* fix point fmod */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis /* 0: we are done */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis m <<= (1 + n);
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis m &= ~0x80000000;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis /* x = z << 1 */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis m &= ~0x80000000;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis /* convert back to floating value and restore the sign */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis /* normalize output */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis __H0(a) = sx | (x0 - iu) | ((iy + 16383) << 16);
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis } else { /* subnormal output */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis if (k <= 16) {
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis while (n >= 32) {
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtisremquol(long double x, long double y, int *quo) {
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis long double v;
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis /* purge off exception values */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis /* y=0, y is NaN, x is NaN or inf */
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis if (y == 0.0L || y != y || hx >= 0x7fff0000)
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return ((x * y) / (x * y));
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis if (x + x > y) {
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis if (x + x >= y) {
25c28e83beb90e7c80452a7c818c5e6f73a07dc8Piotr Jasiukajtis return (sx == 0 ? x : -x);