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

   1 /* Software floating-point emulation.
   2    scalblnl(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 Library General Public License as
   9    published by the Free Software Foundation; either version 2 of the
  10    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    Library General Public License for more details.
  16
  17    You should have received a copy of the GNU Library General Public
  18    License along with the GNU C Library; see the file COPYING.LIB.  If
  19    not, write to the Free Software Foundation, Inc.,
  20    59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.  */
  21
  22 /*
  23  * scalblnl (long double x, long int n)
  24  * scalblnl(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 __scalblnl(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 (__scalblnl, scalblnl)