lib/msun/src/s_modff.c

   1 /* s_modff.c -- float version of s_modf.c.
   2  * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
   3  *
   4  * $FreeBSD: src/lib/msun/src/s_modff.c,v 1.5 1999/08/28 00:06:53 peter Exp $
   5  * $DragonFly: src/lib/msun/src/Attic/s_modff.c,v 1.2 2003/06/17 04:26:53 dillon Exp $
   6  */
   7
   8 /*
   9  * ====================================================
  10  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
  11  *
  12  * Developed at SunPro, a Sun Microsystems, Inc. business.
  13  * Permission to use, copy, modify, and distribute this
  14  * software is freely granted, provided that this notice
  15  * is preserved.
  16  * ====================================================
  17  */
  18
  19 #include "math.h"
  20 #include "math_private.h"
  21
  22 #ifdef __STDC__
  23 static const float one = 1.0;
  24 #else
  25 static float one = 1.0;
  26 #endif
  27
  28 #ifdef __STDC__
  29         float modff(float x, float *iptr)
  30 #else
  31         float modff(x, iptr)
  32         float x,*iptr;
  33 #endif
  34 {
  35         int32_t i0,j0;
  36         u_int32_t i;
  37         GET_FLOAT_WORD(i0,x);
  38         j0 = ((i0>>23)&0xff)-0x7f;      /* exponent of x */
  39         if(j0<23) {                     /* integer part in x */
  40             if(j0<0) {                  /* |x|<1 */
  41                 SET_FLOAT_WORD(*iptr,i0&0x80000000);    /* *iptr = +-0 */
  42                 return x;
  43             } else {
  44                 i = (0x007fffff)>>j0;
  45                 if((i0&i)==0) {                 /* x is integral */
  46                     u_int32_t ix;
  47                     *iptr = x;
  48                     GET_FLOAT_WORD(ix,x);
  49                     SET_FLOAT_WORD(x,ix&0x80000000);    /* return +-0 */
  50                     return x;
  51                 } else {
  52                     SET_FLOAT_WORD(*iptr,i0&(~i));
  53                     return x - *iptr;
  54                 }
  55             }
  56         } else {                        /* no fraction part */
  57             u_int32_t ix;
  58             *iptr = x*one;
  59             GET_FLOAT_WORD(ix,x);
  60             SET_FLOAT_WORD(x,ix&0x80000000);    /* return +-0 */
  61             return x;
  62         }
  63 }