math/w_lgammal_r.c

   1 /* w_lgammal_r.c -- long double version of w_lgamma_r.c.
   2  * Conversion to long double by Ulrich Drepper,
   3  * Cygnus Support, drepper@cygnus.com.
   4  */
   5
   6 /*
   7  * ====================================================
   8  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
   9  *
  10  * Developed at SunPro, a Sun Microsystems, Inc. business.
  11  * Permission to use, copy, modify, and distribute this
  12  * software is freely granted, provided that this notice
  13  * is preserved.
  14  * ====================================================
  15  */
  16
  17 #if defined(LIBM_SCCS) && !defined(lint)
  18 static char rcsid[] = "$NetBSD: $";
  19 #endif
  20
  21 /*
  22  * wrapper long double lgammal_r(long double x, int *signgamp)
  23  */
  24
  25 #include <math.h>
  26 #include "math_private.h"
  27
  28
  29 #ifdef __STDC__
  30         long double __lgammal_r(long double x, int *signgamp)
  31                 /* wrapper lgamma_r */
  32 #else
  33         long double __lgammal_r(x,signgamp)             /* wrapper lgamma_r */
  34         long double x; int *signgamp;
  35 #endif
  36 {
  37 #ifdef _IEEE_LIBM
  38         return __ieee754_lgammal_r(x,signgamp);
  39 #else
  40         long double y;
  41         y = __ieee754_lgammal_r(x,signgamp);
  42         if(_LIB_VERSION == _IEEE_) return y;
  43         if(!__finitel(y)&&__finitel(x)) {
  44             if(__floorl(x)==x&&x<=0.0)
  45                 return __kernel_standard(x,x,215); /* lgamma pole */
  46             else
  47                 return __kernel_standard(x,x,214); /* lgamma overflow */
  48         } else
  49             return y;
  50 #endif
  51 }
  52 weak_alias (__lgammal_r, lgammal_r)