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

   1 /* Software floating-point emulation.
   2    scalbnl(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  * scalbnl (long double x, int n)
  24  * scalbnl(x,n) returns x* 2**n  computed by  exponent
  25  * manipulation rather than by actually performing an
  26  * exponentiation or a multiplication.
  27  */
  28
  29 #include "soft-fp.h"
  30 #include "quad.h"
  31
  32 long double __scalbnl(long double arg, int exp)
  33 {
  34   FP_DECL_EX;
  35   FP_DECL_Q(A);
  36   long double r;
  37
  38   FP_UNPACK_Q(A, arg);
  39   switch (A_c)
  40     {
  41     case FP_CLS_ZERO:
  42       return arg;
  43     case FP_CLS_NAN:
  44     case FP_CLS_INF:
  45       FP_HANDLE_EXCEPTIONS;
  46       return arg;
  47     }
  48   A_e += exp;
  49   FP_PACK_Q(r, A);
  50   FP_HANDLE_EXCEPTIONS;
  51
  52   return r;
  53 }
  54
  55 weak_alias (__scalbnl, scalbnl)