newlib/libm/math/s_frexp.c

   1
   2 /* @(#)s_frexp.c 5.1 93/09/24 */
   3 /*
   4  * ====================================================
   5  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
   6  *
   7  * Developed at SunPro, a Sun Microsystems, Inc. business.
   8  * Permission to use, copy, modify, and distribute this
   9  * software is freely granted, provided that this notice
  10  * is preserved.
  11  * ====================================================
  12  */
  13
  14 /*
  15 FUNCTION
  16        <<frexp>>, <<frexpf>>---split floating-point number
  17 INDEX
  18         frexp
  19 INDEX
  20         frexpf
  21
  22 ANSI_SYNOPSIS
  23         #include <math.h>
  24         double frexp(double <[val]>, int *<[exp]>);
  25         float frexpf(float <[val]>, int *<[exp]>);
  26
  27 TRAD_SYNOPSIS
  28         #include <math.h>
  29         double frexp(<[val]>, <[exp]>)
  30         double <[val]>;
  31         int *<[exp]>;
  32
  33         float frexpf(<[val]>, <[exp]>)
  34         float <[val]>;
  35         int *<[exp]>;
  36
  37
  38 DESCRIPTION
  39         All non zero, normal numbers can be described as <[m]> * 2**<[p]>.
  40         <<frexp>> represents the double <[val]> as a mantissa <[m]>
  41         and a power of two <[p]>. The resulting mantissa will always
  42         be greater than or equal to <<0.5>>, and less than <<1.0>> (as
  43         long as <[val]> is nonzero). The power of two will be stored
  44         in <<*>><[exp]>.
  45
  46 @ifinfo
  47 <[m]> and <[p]> are calculated so that
  48 <[val]> is <[m]> times <<2>> to the power <[p]>.
  49 @end ifinfo
  50 @tex
  51 <[m]> and <[p]> are calculated so that
  52 $ val = m \times 2^p $.
  53 @end tex
  54
  55 <<frexpf>> is identical, other than taking and returning
  56 floats rather than doubles.
  57
  58 RETURNS
  59 <<frexp>> returns the mantissa <[m]>. If <[val]> is <<0>>, infinity,
  60 or Nan, <<frexp>> will set <<*>><[exp]> to <<0>> and return <[val]>.
  61
  62 PORTABILITY
  63 <<frexp>> is ANSI.
  64 <<frexpf>> is an extension.
  65
  66
  67 */
  68
  69 /*
  70  * for non-zero x
  71  *      x = frexp(arg,&exp);
  72  * return a double fp quantity x such that 0.5 <= |x| <1.0
  73  * and the corresponding binary exponent "exp". That is
  74  *      arg = x*2^exp.
  75  * If arg is inf, 0.0, or NaN, then frexp(arg,&exp) returns arg
  76  * with *exp=0.
  77  */
  78
  79 #include "fdlibm.h"
  80
  81 #ifndef _DOUBLE_IS_32BITS
  82
  83 #ifdef __STDC__
  84 static const double
  85 #else
  86 static double
  87 #endif
  88 two54 =  1.80143985094819840000e+16; /* 0x43500000, 0x00000000 */
  89
  90 #ifdef __STDC__
  91         double frexp(double x, int *eptr)
  92 #else
  93         double frexp(x, eptr)
  94         double x; int *eptr;
  95 #endif
  96 {
  97         __int32_t hx, ix, lx;
  98         EXTRACT_WORDS(hx,lx,x);
  99         ix = 0x7fffffff&hx;
 100         *eptr = 0;
 101         if(ix>=0x7ff00000||((ix|lx)==0)) return x;      /* 0,inf,nan */
 102         if (ix<0x00100000) {            /* subnormal */
 103             x *= two54;
 104             GET_HIGH_WORD(hx,x);
 105             ix = hx&0x7fffffff;
 106             *eptr = -54;
 107         }
 108         *eptr += (ix>>20)-1022;
 109         hx = (hx&0x800fffff)|0x3fe00000;
 110         SET_HIGH_WORD(x,hx);
 111         return x;
 112 }
 113
 114 #endif /* _DOUBLE_IS_32BITS */