source/wham/src/xdrf.org/xdr_float.c

   1 /* @(#)xdr_float.c      2.1 88/07/29 4.0 RPCSRC */
   2 /*
   3  * Sun RPC is a product of Sun Microsystems, Inc. and is provided for
   4  * unrestricted use provided that this legend is included on all tape
   5  * media and as a part of the software program in whole or part.  Users
   6  * may copy or modify Sun RPC without charge, but are not authorized
   7  * to license or distribute it to anyone else except as part of a product or
   8  * program developed by the user.
   9  *
  10  * SUN RPC IS PROVIDED AS IS WITH NO WARRANTIES OF ANY KIND INCLUDING THE
  11  * WARRANTIES OF DESIGN, MERCHANTIBILITY AND FITNESS FOR A PARTICULAR
  12  * PURPOSE, OR ARISING FROM A COURSE OF DEALING, USAGE OR TRADE PRACTICE.
  13  *
  14  * Sun RPC is provided with no support and without any obligation on the
  15  * part of Sun Microsystems, Inc. to assist in its use, correction,
  16  * modification or enhancement.
  17  *
  18  * SUN MICROSYSTEMS, INC. SHALL HAVE NO LIABILITY WITH RESPECT TO THE
  19  * INFRINGEMENT OF COPYRIGHTS, TRADE SECRETS OR ANY PATENTS BY SUN RPC
  20  * OR ANY PART THEREOF.
  21  *
  22  * In no event will Sun Microsystems, Inc. be liable for any lost revenue
  23  * or profits or other special, indirect and consequential damages, even if
  24  * Sun has been advised of the possibility of such damages.
  25  *
  26  * Sun Microsystems, Inc.
  27  * 2550 Garcia Avenue
  28  * Mountain View, California  94043
  29  */
  30 #if !defined(lint) && defined(SCCSIDS)
  31 static char sccsid[] = "@(#)xdr_float.c 1.12 87/08/11 Copyr 1984 Sun Micro";
  32 #endif
  33
  34 /*
  35  * xdr_float.c, Generic XDR routines implementation.
  36  *
  37  * Copyright (C) 1984, Sun Microsystems, Inc.
  38  *
  39  * These are the "floating point" xdr routines used to (de)serialize
  40  * most common data items.  See xdr.h for more info on the interface to
  41  * xdr.
  42  */
  43
  44 #include <stdio.h>
  45 #include <endian.h>
  46
  47 #include "types.h"
  48 #include "xdr.h"
  49
  50 /*
  51  * NB: Not portable.
  52  * This routine works on Suns (Sky / 68000's) and Vaxen.
  53  */
  54
  55 #define LSW     (__FLOAT_WORD_ORDER == __BIG_ENDIAN)
  56
  57 #ifdef vax
  58
  59 /* What IEEE single precision floating point looks like on a Vax */
  60 struct  ieee_single {
  61         unsigned int    mantissa: 23;
  62         unsigned int    exp     : 8;
  63         unsigned int    sign    : 1;
  64 };
  65
  66 /* Vax single precision floating point */
  67 struct  vax_single {
  68         unsigned int    mantissa1 : 7;
  69         unsigned int    exp       : 8;
  70         unsigned int    sign      : 1;
  71         unsigned int    mantissa2 : 16;
  72 };
  73
  74 #define VAX_SNG_BIAS    0x81
  75 #define IEEE_SNG_BIAS   0x7f
  76
  77 static struct sgl_limits {
  78         struct vax_single s;
  79         struct ieee_single ieee;
  80 } sgl_limits[2] = {
  81         {{ 0x7f, 0xff, 0x0, 0xffff },   /* Max Vax */
  82         { 0x0, 0xff, 0x0 }},            /* Max IEEE */
  83         {{ 0x0, 0x0, 0x0, 0x0 },        /* Min Vax */
  84         { 0x0, 0x0, 0x0 }}              /* Min IEEE */
  85 };
  86 #endif /* vax */
  87
  88 bool_t
  89 xdr_float(xdrs, fp)
  90      XDR *xdrs;
  91      float *fp;
  92 {
  93 #ifdef vax
  94         struct ieee_single is;
  95         struct vax_single vs, *vsp;
  96         struct sgl_limits *lim;
  97         int i;
  98 #endif
  99         switch (xdrs->x_op) {
 100
 101         case XDR_ENCODE:
 102 #ifdef vax
 103                 vs = *((struct vax_single *)fp);
 104                 for (i = 0, lim = sgl_limits;
 105                         i < sizeof(sgl_limits)/sizeof(struct sgl_limits);
 106                         i++, lim++) {
 107                         if ((vs.mantissa2 == lim->s.mantissa2) &&
 108                                 (vs.exp == lim->s.exp) &&
 109                                 (vs.mantissa1 == lim->s.mantissa1)) {
 110                                 is = lim->ieee;
 111                                 goto shipit;
 112                         }
 113                 }
 114                 is.exp = vs.exp - VAX_SNG_BIAS + IEEE_SNG_BIAS;
 115                 is.mantissa = (vs.mantissa1 << 16) | vs.mantissa2;
 116         shipit:
 117                 is.sign = vs.sign;
 118                 return (XDR_PUTLONG(xdrs, (long *)&is));
 119 #else
 120                 if (sizeof(float) == sizeof(long))
 121                         return (XDR_PUTLONG(xdrs, (long *)fp));
 122                 else if (sizeof(float) == sizeof(int)) {
 123                         long tmp = *(int *)fp;
 124                         return (XDR_PUTLONG(xdrs, &tmp));
 125                 }
 126                 break;
 127 #endif
 128
 129         case XDR_DECODE:
 130 #ifdef vax
 131                 vsp = (struct vax_single *)fp;
 132                 if (!XDR_GETLONG(xdrs, (long *)&is))
 133                         return (FALSE);
 134                 for (i = 0, lim = sgl_limits;
 135                         i < sizeof(sgl_limits)/sizeof(struct sgl_limits);
 136                         i++, lim++) {
 137                         if ((is.exp == lim->ieee.exp) &&
 138                                 (is.mantissa == lim->ieee.mantissa)) {
 139                                 *vsp = lim->s;
 140                                 goto doneit;
 141                         }
 142                 }
 143                 vsp->exp = is.exp - IEEE_SNG_BIAS + VAX_SNG_BIAS;
 144                 vsp->mantissa2 = is.mantissa;
 145                 vsp->mantissa1 = (is.mantissa >> 16);
 146         doneit:
 147                 vsp->sign = is.sign;
 148                 return (TRUE);
 149 #else
 150                 if (sizeof(float) == sizeof(long))
 151                         return (XDR_GETLONG(xdrs, (long *)fp));
 152                 else if (sizeof(float) == sizeof(int)) {
 153                         long tmp;
 154                         if (XDR_GETLONG(xdrs, &tmp)) {
 155                                 *(int *)fp = tmp;
 156                                 return (TRUE);
 157                         }
 158                 }
 159                 break;
 160 #endif
 161
 162         case XDR_FREE:
 163                 return (TRUE);
 164         }
 165         return (FALSE);
 166 }
 167
 168 /*
 169  * This routine works on Suns (Sky / 68000's) and Vaxen.
 170  */
 171
 172 #ifdef vax
 173 /* What IEEE double precision floating point looks like on a Vax */
 174 struct  ieee_double {
 175         unsigned int    mantissa1 : 20;
 176         unsigned int    exp       : 11;
 177         unsigned int    sign      : 1;
 178         unsigned int    mantissa2 : 32;
 179 };
 180
 181 /* Vax double precision floating point */
 182 struct  vax_double {
 183         unsigned int    mantissa1 : 7;
 184         unsigned int    exp       : 8;
 185         unsigned int    sign      : 1;
 186         unsigned int    mantissa2 : 16;
 187         unsigned int    mantissa3 : 16;
 188         unsigned int    mantissa4 : 16;
 189 };
 190
 191 #define VAX_DBL_BIAS    0x81
 192 #define IEEE_DBL_BIAS   0x3ff
 193 #define MASK(nbits)     ((1 << nbits) - 1)
 194
 195 static struct dbl_limits {
 196         struct  vax_double d;
 197         struct  ieee_double ieee;
 198 } dbl_limits[2] = {
 199         {{ 0x7f, 0xff, 0x0, 0xffff, 0xffff, 0xffff },   /* Max Vax */
 200         { 0x0, 0x7ff, 0x0, 0x0 }},                      /* Max IEEE */
 201         {{ 0x0, 0x0, 0x0, 0x0, 0x0, 0x0},               /* Min Vax */
 202         { 0x0, 0x0, 0x0, 0x0 }}                         /* Min IEEE */
 203 };
 204
 205 #endif /* vax */
 206
 207
 208 bool_t
 209 xdr_double(xdrs, dp)
 210      XDR *xdrs;
 211      double *dp;
 212 {
 213 #ifdef vax
 214         struct  ieee_double id;
 215         struct  vax_double vd;
 216         register struct dbl_limits *lim;
 217         int i;
 218 #endif
 219
 220         switch (xdrs->x_op) {
 221
 222         case XDR_ENCODE:
 223 #ifdef vax
 224                 vd = *((struct vax_double *)dp);
 225                 for (i = 0, lim = dbl_limits;
 226                         i < sizeof(dbl_limits)/sizeof(struct dbl_limits);
 227                         i++, lim++) {
 228                         if ((vd.mantissa4 == lim->d.mantissa4) &&
 229                                 (vd.mantissa3 == lim->d.mantissa3) &&
 230                                 (vd.mantissa2 == lim->d.mantissa2) &&
 231                                 (vd.mantissa1 == lim->d.mantissa1) &&
 232                                 (vd.exp == lim->d.exp)) {
 233                                 id = lim->ieee;
 234                                 goto shipit;
 235                         }
 236                 }
 237                 id.exp = vd.exp - VAX_DBL_BIAS + IEEE_DBL_BIAS;
 238                 id.mantissa1 = (vd.mantissa1 << 13) | (vd.mantissa2 >> 3);
 239                 id.mantissa2 = ((vd.mantissa2 & MASK(3)) << 29) |
 240                                 (vd.mantissa3 << 13) |
 241                                 ((vd.mantissa4 >> 3) & MASK(13));
 242         shipit:
 243                 id.sign = vd.sign;
 244                 dp = (double *)&id;
 245 #endif
 246                 if (2*sizeof(long) == sizeof(double)) {
 247                         long *lp = (long *)dp;
 248                         return (XDR_PUTLONG(xdrs, lp+!LSW) &&
 249                                 XDR_PUTLONG(xdrs, lp+LSW));
 250                 } else if (2*sizeof(int) == sizeof(double)) {
 251                         int *ip = (int *)dp;
 252                         long tmp[2];
 253                         tmp[0] = ip[!LSW];
 254                         tmp[1] = ip[LSW];
 255                         return (XDR_PUTLONG(xdrs, tmp) &&
 256                                 XDR_PUTLONG(xdrs, tmp+1));
 257                 }
 258                 break;
 259
 260         case XDR_DECODE:
 261 #ifdef vax
 262                 lp = (long *)&id;
 263                 if (!XDR_GETLONG(xdrs, lp++) || !XDR_GETLONG(xdrs, lp))
 264                         return (FALSE);
 265                 for (i = 0, lim = dbl_limits;
 266                         i < sizeof(dbl_limits)/sizeof(struct dbl_limits);
 267                         i++, lim++) {
 268                         if ((id.mantissa2 == lim->ieee.mantissa2) &&
 269                                 (id.mantissa1 == lim->ieee.mantissa1) &&
 270                                 (id.exp == lim->ieee.exp)) {
 271                                 vd = lim->d;
 272                                 goto doneit;
 273                         }
 274                 }
 275                 vd.exp = id.exp - IEEE_DBL_BIAS + VAX_DBL_BIAS;
 276                 vd.mantissa1 = (id.mantissa1 >> 13);
 277                 vd.mantissa2 = ((id.mantissa1 & MASK(13)) << 3) |
 278                                 (id.mantissa2 >> 29);
 279                 vd.mantissa3 = (id.mantissa2 >> 13);
 280                 vd.mantissa4 = (id.mantissa2 << 3);
 281         doneit:
 282                 vd.sign = id.sign;
 283                 *dp = *((double *)&vd);
 284                 return (TRUE);
 285 #else
 286                 if (2*sizeof(long) == sizeof(double)) {
 287                         long *lp = (long *)dp;
 288                         return (XDR_GETLONG(xdrs, lp+!LSW) &&
 289                                 XDR_GETLONG(xdrs, lp+LSW));
 290                 } else if (2*sizeof(int) == sizeof(double)) {
 291                         int *ip = (int *)dp;
 292                         long tmp[2];
 293                         if (XDR_GETLONG(xdrs, tmp+!LSW) &&
 294                             XDR_GETLONG(xdrs, tmp+LSW)) {
 295                                 ip[0] = tmp[0];
 296                                 ip[1] = tmp[1];
 297                                 return (TRUE);
 298                         }
 299                 }
 300                 break;
 301 #endif
 302
 303         case XDR_FREE:
 304                 return (TRUE);
 305         }
 306         return (FALSE);
 307 }