xdr_float.c revision 2
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz/*
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * CDDL HEADER START
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz *
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * The contents of this file are subject to the terms of the
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * Common Development and Distribution License (the "License").
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed * You may not use this file except in compliance with the License.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz *
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * or http://www.opensolaris.org/os/licensing.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * See the License for the specific language governing permissions
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * and limitations under the License.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz *
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * When distributing Covered Code, include this CDDL HEADER in each
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * If applicable, add the following below this CDDL HEADER, with the
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * fields enclosed by brackets "[]" replaced with your own identifying
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * information: Portions Copyright [yyyy] [name of copyright owner]
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz *
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * CDDL HEADER END
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz/*
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * Copyright (c) 1989, 2011, Oracle and/or its affiliates. All rights reserved.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz/* Copyright (c) 1983, 1984, 1985, 1986, 1987, 1988, 1989 AT&T */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz/* All Rights Reserved */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz/*
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * Portions of this source code were derived from Berkeley
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * 4.3 BSD under license from the Regents of the University of
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * California.
381a2a9a387f449fab7d0c7e97c4184c26963abfdr */
381a2a9a387f449fab7d0c7e97c4184c26963abfdr
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz/*
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * Generic XDR routines impelmentation.
193974072f41a843678abf5f61979c748687e66bSherry Moore *
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * These are the "floating point" xdr routines used to (de)serialize
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * most common data items. See xdr.h for more info on the interface to
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * xdr.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#include "mt.h"
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#include <sys/types.h>
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#include <stdio.h>
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#include <rpc/types.h>
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#include <rpc/xdr.h>
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz/*
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * This routine works on Suns, 3b2, 68000s, 386 and Vaxen in a manner
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * which is very efficient as bit twiddling is all that is needed. All
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * other machines can use this code but the code is inefficient as
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * various mathematical operations are used to generate the ieee format.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * In addition rounding errors may occur due to the calculations involved.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * To be most efficient, new machines should have their own ifdefs.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * The encoding routines will fail if the machines try to encode a
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * float/double whose value can not be represented by the ieee format,
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * e.g. the exponent is too big/small.
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * ieee largest float = (2 ^ 128) * 0x1.fffff
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * ieee smallest float = (2 ^ -127) * 0x1.00000
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * ieee largest double = (2 ^ 1024) * 0x1.fffff
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * ieee smallest double = (2 ^ -1023) * 0x1.00000
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * The decoding routines assumes that the receiving machine can handle
193974072f41a843678abf5f61979c748687e66bSherry Moore * floats/doubles as large/small as the values stated above. If you
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * use a machine which can not represent these values, you will need
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * to put ifdefs in the decode sections to identify areas of failure.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#if defined(vax)
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed/*
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * What IEEE single precision floating point looks like this on a
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * vax.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yzstruct ieee_single {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int mantissa: 23;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int exp : 8;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int sign : 1;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz};
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#define IEEE_SNG_BIAS 0x7f
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#define VAX_SNG_BIAS 0x81
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz/* Vax single precision floating point */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yzstruct vax_single {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int mantissa1 : 7;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int exp : 8;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int sign : 1;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int mantissa2 : 16;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz};
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#define VAX_SNG_BIAS 0x81
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yzstatic struct sgl_limits {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz struct vax_single s;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz struct ieee_single ieee;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz} sgl_limits[2] = {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz {{ 0x7f, 0xff, 0x0, 0xffff }, /* Max Vax */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz { 0x0, 0xff, 0x0 }}, /* Max IEEE */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz {{ 0x0, 0x0, 0x0, 0x0 }, /* Min Vax */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz { 0x0, 0x0, 0x0 }} /* Min IEEE */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz};
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#endif /* vax */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yzbool_t
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yzxdr_float(XDR *xdrs, float *fp)
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz{
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#if defined(vax)
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz struct ieee_single is;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz struct vax_single vs, *vsp;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz struct sgl_limits *lim;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz size_t i;
193974072f41a843678abf5f61979c748687e66bSherry Moore#endif
193974072f41a843678abf5f61979c748687e66bSherry Moore
193974072f41a843678abf5f61979c748687e66bSherry Moore switch (xdrs->x_op) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz case XDR_ENCODE:
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed#if defined(mc68000) || defined(sparc) || defined(u3b2) || \
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed defined(u3b15) || defined(i386) || defined(__amd64)
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed return (XDR_PUTINT32(xdrs, (int *)fp));
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed#else
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#if defined(vax)
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz vs = *((struct vax_single *)fp);
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz if ((vs.exp == 1) || (vs.exp == 2)) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz /* map these to subnormals */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz is.exp = 0;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz is.mantissa = (vs.mantissa1 << 16) | vs.mantissa2;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz /* lose some precision */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz is.mantissa >>= 3 - vs.exp;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz is.mantissa += (1 << (20 + vs.exp));
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz goto shipit;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz for (i = 0, lim = sgl_limits;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz i < (int)(sizeof (sgl_limits) /
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz sizeof (struct sgl_limits));
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz i++, lim++) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz if ((vs.mantissa2 == lim->s.mantissa2) &&
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz (vs.exp == lim->s.exp) &&
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz (vs.mantissa1 == lim->s.mantissa1)) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz is = lim->ieee;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz goto shipit;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz is.exp = vs.exp - VAX_SNG_BIAS + IEEE_SNG_BIAS;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz is.mantissa = (vs.mantissa1 << 16) | vs.mantissa2;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz shipit:
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz is.sign = vs.sign;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz return (XDR_PUTINT32(xdrs, (int32_t *)&is));
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#else
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz /*
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * Every machine can do this, its just not very efficient.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * In addtion, some rounding errors may occur do to the
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * calculations involved.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz float f;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz int neg = 0;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz int exp = 0;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz int32_t val;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz f = *fp;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz if (f == 0) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz val = 0;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz return (XDR_PUTINT32(xdrs, &val));
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz if (f < 0) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz f = 0 - f;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz neg = 1;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz while (f < 1) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz f = f * 2;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz --exp;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz while (f >= 2) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz f = f/2;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz ++exp;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz if ((exp > 128) || (exp < -127)) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz /* over or under flowing ieee exponent */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz return (FALSE);
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz val = neg;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz val = val << 8; /* for the exponent */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz val += 127 + exp; /* 127 is the bias */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz val = val << 23; /* for the mantissa */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz val += (int32_t)((f - 1) * 8388608); /* 2 ^ 23 */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh return (XDR_PUTINT32(xdrs, &val));
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#endif
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#endif
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz case XDR_DECODE:
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#if defined(mc68000) || defined(sparc) || defined(u3b2) || \
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz defined(u3b15) || defined(i386) || defined(__amd64)
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz return (XDR_GETINT32(xdrs, (int *)fp));
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#else
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#if defined(vax)
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz vsp = (struct vax_single *)fp;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz if (!XDR_GETINT32(xdrs, (int32_t *)&is))
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz return (FALSE);
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz for (i = 0, lim = sgl_limits;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz i < (int)(sizeof (sgl_limits) /
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz sizeof (struct sgl_limits));
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz i++, lim++) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz if ((is.exp == lim->ieee.exp) &&
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz (is.mantissa == lim->ieee.mantissa)) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz *vsp = lim->s;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz goto doneit;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz } else if ((is.exp == 0) && (lim->ieee.exp == 0)) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz /* Special Case */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned tmp = is.mantissa >> 20;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz if (tmp >= 4) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz vsp->exp = 2;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz } else if (tmp >= 2) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz vsp->exp = 1;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz } else {
381a2a9a387f449fab7d0c7e97c4184c26963abfdr *vsp = min.s;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz break;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz } /* else */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz tmp = is.mantissa - (1 << (20 + vsp->exp));
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed tmp <<= 3 - vsp->exp;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz vsp->mantissa2 = tmp;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed vsp->mantissa1 = (tmp >> 16);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed goto doneit;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed }
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed vsp->exp = is.exp - IEEE_SNG_BIAS + VAX_SNG_BIAS;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed vsp->mantissa2 = is.mantissa;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed vsp->mantissa1 = (is.mantissa >> 16);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed doneit:
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed vsp->sign = is.sign;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed return (TRUE);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed#else
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz {
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed /*
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed * Every machine can do this, its just not very
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed * efficient. It assumes that the decoding machine's
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed * float can represent any value in the range of
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed * ieee largest float = (2 ^ 128) * 0x1.fffff
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed * to
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed * ieee smallest float = (2 ^ -127) * 0x1.00000
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed * In addtion, some rounding errors may occur do to the
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed * calculations involved.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz float f;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz int neg = 0;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed int exp = 0;
193974072f41a843678abf5f61979c748687e66bSherry Moore int32_t val;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz if (!XDR_GETINT32(xdrs, (int32_t *)&val))
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz return (FALSE);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed neg = val & 0x80000000;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz exp = (val & 0x7f800000) >> 23;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed exp -= 127; /* subtract exponent base */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz f = (val & 0x007fffff) * 0.00000011920928955078125;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz /* 2 ^ -23 */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh f++;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz while (exp != 0) {
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh if (exp < 0) {
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed f = f/2.0;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh ++exp;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz } else {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz f = f * 2.0;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz --exp;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz }
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz if (neg)
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz f = 0 - f;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz *fp = f;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz return (TRUE);
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#endif
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh#endif
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh case XDR_FREE:
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz return (TRUE);
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz return (FALSE);
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz}
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz/*
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz * This routine works on Suns (Sky / 68000's) and Vaxen.
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#if defined(vax)
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz/* What IEEE double precision floating point looks like on a Vax */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yzstruct ieee_double {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int mantissa1 : 20;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int exp : 11;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int sign : 1;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int mantissa2 : 32;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz};
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz/* Vax double precision floating point */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yzstruct vax_double {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int mantissa1 : 7;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int exp : 8;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int sign : 1;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int mantissa2 : 16;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int mantissa3 : 16;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz unsigned int mantissa4 : 16;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz};
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#define VAX_DBL_BIAS 0x81
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#define IEEE_DBL_BIAS 0x3ff
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#define MASK(nbits) ((1 << nbits) - 1)
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yzstatic struct dbl_limits {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz struct vax_double d;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz struct ieee_double ieee;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz} dbl_limits[2] = {
193974072f41a843678abf5f61979c748687e66bSherry Moore {{ 0x7f, 0xff, 0x0, 0xffff, 0xffff, 0xffff }, /* Max Vax */
193974072f41a843678abf5f61979c748687e66bSherry Moore { 0x0, 0x7ff, 0x0, 0x0 }}, /* Max IEEE */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz {{ 0x0, 0x0, 0x0, 0x0, 0x0, 0x0}, /* Min Vax */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz { 0x0, 0x0, 0x0, 0x0 }} /* Min IEEE */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz};
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#endif /* vax */
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reedbool_t
193974072f41a843678abf5f61979c748687e66bSherry Moorexdr_double(XDR *xdrs, double *dp)
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz{
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz int *lp;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#if defined(vax)
193974072f41a843678abf5f61979c748687e66bSherry Moore struct ieee_double id;
193974072f41a843678abf5f61979c748687e66bSherry Moore struct vax_double vd;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz struct dbl_limits *lim;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz size_t i;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#endif
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz switch (xdrs->x_op) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz case XDR_ENCODE:
193974072f41a843678abf5f61979c748687e66bSherry Moore#if defined(mc68000) || defined(u3b2) || defined(u3b15) || \
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz defined(_LONG_LONG_HTOL)
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz lp = (int *)dp;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz return (XDR_PUTINT32(xdrs, lp++) && XDR_PUTINT32(xdrs, lp));
193974072f41a843678abf5f61979c748687e66bSherry Moore#else
193974072f41a843678abf5f61979c748687e66bSherry Moore#if defined(_LONG_LONG_LTOH)
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz lp = (int *)dp;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz lp++;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz return (XDR_PUTINT32(xdrs, lp--) && XDR_PUTINT32(xdrs, lp));
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#else
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz#if defined(vax)
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz vd = *((struct vax_double *)dp);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed for (i = 0, lim = dbl_limits;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz i < (int)(sizeof (dbl_limits) /
193974072f41a843678abf5f61979c748687e66bSherry Moore sizeof (struct dbl_limits));
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz i++, lim++) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz if ((vd.mantissa4 == lim->d.mantissa4) &&
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz (vd.mantissa3 == lim->d.mantissa3) &&
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz (vd.mantissa2 == lim->d.mantissa2) &&
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz (vd.mantissa1 == lim->d.mantissa1) &&
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz (vd.exp == lim->d.exp)) {
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz id = lim->ieee;
193974072f41a843678abf5f61979c748687e66bSherry Moore goto shipit;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed }
193974072f41a843678abf5f61979c748687e66bSherry Moore }
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz id.exp = vd.exp - VAX_DBL_BIAS + IEEE_DBL_BIAS;
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz id.mantissa1 = (vd.mantissa1 << 13) | (vd.mantissa2 >> 3);
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz id.mantissa2 = ((vd.mantissa2 & MASK(3)) << 29) |
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz (vd.mantissa3 << 13) |
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz ((vd.mantissa4 >> 3) & MASK(13));
ab25eeb551a4be927a4b6ae2cf8aff7ed17decb4yz shipit:
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh id.sign = vd.sign;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh lp = (int32_t *)&id;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh#else
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh {
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed /*
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * Every machine can do this, its just not very efficient.
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * In addtion, some rounding errors may occur do to the
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * calculations involved.
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed */
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed double d;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh int neg = 0;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh int exp = 0;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed int32_t val[2];
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh d = *dp;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh if (d == 0) {
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh val[0] = 0;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh val[1] = 0;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh lp = val;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh return (XDR_PUTINT32(xdrs, lp++) &&
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh XDR_PUTINT32(xdrs, lp));
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh }
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh if (d < 0) {
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh d = 0 - d;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh neg = 1;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh }
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh while (d < 1) {
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh d = d * 2;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh --exp;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh }
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh while (d >= 2) {
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed d = d/2;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed ++exp;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed }
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed if ((exp > 1024) || (exp < -1023)) {
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed /* over or under flowing ieee exponent */
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed return (FALSE);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed }
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh val[0] = neg;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh val[0] = val[0] << 11; /* for the exponent */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh val[0] += 1023 + exp; /* 1023 is the bias */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh val[0] = val[0] << 20; /* for the mantissa */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh val[0] += (int32_t)((d - 1) * 1048576); /* 2 ^ 20 */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh val[1] += (int32_t)((((d - 1) * 1048576) - val[0])
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * 4294967296);
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh /* 2 ^ 32 */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh lp = val;
23f4867fdff96a11dd674de6259a5a0d0a13251cnordmark }
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed#endif
23f4867fdff96a11dd674de6259a5a0d0a13251cnordmark return (XDR_PUTINT32(xdrs, lp++) && XDR_PUTINT32(xdrs, lp));
23f4867fdff96a11dd674de6259a5a0d0a13251cnordmark#endif
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed#endif
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed case XDR_DECODE:
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed#if defined(mc68000) || defined(u3b2) || defined(u3b15) || \
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed defined(_LONG_LONG_HTOL)
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed lp = (int *)dp;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed return (XDR_GETINT32(xdrs, lp++) && XDR_GETINT32(xdrs, lp));
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed#else
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh#if defined(_LONG_LONG_LTOH)
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh lp = (int *)dp;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh lp++;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed return (XDR_GETINT32(xdrs, lp--) && XDR_GETINT32(xdrs, lp));
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed#else
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed#if defined(vax)
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed lp = (int32_t *)&id;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed if (!XDR_GETINT32(xdrs, lp++) || !XDR_GETINT32(xdrs, lp))
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed return (FALSE);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed for (i = 0, lim = dbl_limits;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed i < sizeof (dbl_limits) /
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed sizeof (struct dbl_limits);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed i++, lim++) {
193974072f41a843678abf5f61979c748687e66bSherry Moore if ((id.mantissa2 == lim->ieee.mantissa2) &&
193974072f41a843678abf5f61979c748687e66bSherry Moore (id.mantissa1 == lim->ieee.mantissa1) &&
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed (id.exp == lim->ieee.exp)) {
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed vd = lim->d;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed goto doneit;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed }
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed }
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed vd.exp = id.exp - IEEE_DBL_BIAS + VAX_DBL_BIAS;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed vd.mantissa1 = (id.mantissa1 >> 13);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed vd.mantissa2 = ((id.mantissa1 & MASK(13)) << 3) |
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed (id.mantissa2 >> 29);
193974072f41a843678abf5f61979c748687e66bSherry Moore vd.mantissa3 = (id.mantissa2 >> 13);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed vd.mantissa4 = (id.mantissa2 << 3);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed doneit:
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed vd.sign = id.sign;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh *dp = *((double *)&vd);
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh return (TRUE);
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh#else
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh {
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh /*
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * Every machine can do this, its just not very
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * efficient. It assumes that the decoding machine's
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * double can represent any value in the range of
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * ieee largest double = (2 ^ 1024) * 0x1.fffffffffffff
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * to
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * ieee smallest double = (2 ^ -1023) * 0x1.0000000000000
193974072f41a843678abf5f61979c748687e66bSherry Moore * In addtion, some rounding errors may occur do to the
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * calculations involved.
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh double d;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh int neg = 0;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh int exp = 0;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh int32_t val[2];
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh lp = val;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh if (!XDR_GETINT32(xdrs, lp++) || !XDR_GETINT32(xdrs, lp))
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh return (FALSE);
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh neg = val[0] & 0x80000000;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh exp = (val[0] & 0x7ff00000) >> 20;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh exp -= 1023; /* subtract exponent base */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh d = (val[0] & 0x000fffff) * 0.00000095367431640625;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh /* 2 ^ -20 */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh d += (val[1] * 0.0000000000000002220446049250313);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed /* 2 ^ -52 */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh d++;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh while (exp != 0) {
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed if (exp < 0) {
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed d = d/2.0;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed ++exp;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed } else {
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed d = d * 2.0;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed --exp;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed }
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed }
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh if (neg)
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh d = 0 - d;
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed *dp = d;
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh }
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh#endif
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh#endif
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh#endif
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed case XDR_FREE:
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh return (TRUE);
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh }
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed return (FALSE);
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh}
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed/* ARGSUSED */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dhbool_t
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dhxdr_quadruple(XDR *xdrs, long double *fp)
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh{
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh/*
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh * The Sparc uses IEEE FP encoding, so just do a byte copy
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh */
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh#if !defined(sparc)
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh return (FALSE);
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh#else
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh switch (xdrs->x_op) {
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh case XDR_ENCODE:
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed return (XDR_PUTBYTES(xdrs, (char *)fp, sizeof (long double)));
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed case XDR_DECODE:
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh return (XDR_GETBYTES(xdrs, (char *)fp, sizeof (long double)));
f4b3ec61df05330d25f55a36b975b4d7519fdeb1dh case XDR_FREE:
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed return (TRUE);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed }
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed return (FALSE);
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed#endif
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed}
7ddc9b1afd18f260b9fb78ec7732facd91769131Darren Reed