sysdeps/sparc/sparc64/soft-fp/s_frexpl.c

   1 /* Software floating-point emulation.
   2    frexpl(x, exp)
   3    Copyright (C) 1999 Free Software Foundation, Inc.
   4    This file is part of the GNU C Library.
   5    Contributed by Jakub Jelinek (jj@ultra.linux.cz).
   6
   7    The GNU C Library is free software; you can redistribute it and/or
   8    modify it under the terms of the GNU Lesser General Public
   9    License as published by the Free Software Foundation; either
  10    version 2.1 of the License, or (at your option) any later version.
  11
  12    The GNU C Library is distributed in the hope that it will be useful,
  13    but WITHOUT ANY WARRANTY; without even the implied warranty of
  14    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  15    Lesser General Public License for more details.
  16
  17    You should have received a copy of the GNU Lesser General Public
  18    License along with the GNU C Library; if not, write to the Free
  19    Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
  20    02111-1307 USA.  */
  21
  22 /*
  23  * for non-zero x
  24  *      x = frexpl(arg,&exp);
  25  * return a long double fp quantity x such that 0.5 <= |x| <1.0
  26  * and the corresponding binary exponent "exp". That is
  27  *      arg = x*2^exp.
  28  * If arg is inf, 0.0, or NaN, then frexpl(arg,&exp) returns arg
  29  * with *exp=0.
  30  */
  31
  32 #include "soft-fp.h"
  33 #include "quad.h"
  34
  35 long double __frexpl(long double arg, int *exp)
  36 {
  37   FP_DECL_EX;
  38   FP_DECL_Q(A);
  39   long double r;
  40
  41   *exp = 0;
  42   FP_UNPACK_Q(A, arg);
  43   if (A_c != FP_CLS_NORMAL)
  44     return arg;
  45   *exp = A_e + 1;
  46   A_e = -1;
  47   FP_PACK_Q(r, A);
  48
  49   return r;
  50 }
  51
  52 weak_alias (__frexpl, frexpl)