1 

     2 /*

     3  * Copyright (c) 1998, 2001, Oracle and/or its affiliates. All rights reserved.

     4  * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER.

     5  *

     6  * This code is free software; you can redistribute it and/or modify it

     7  * under the terms of the GNU General Public License version 2 only, as

     8  * published by the Free Software Foundation.  Oracle designates this

     9  * particular file as subject to the "Classpath" exception as provided

    10  * by Oracle in the LICENSE file that accompanied this code.

    11  *

    12  * This code is distributed in the hope that it will be useful, but WITHOUT

    13  * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or

    14  * FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License

    15  * version 2 for more details (a copy is included in the LICENSE file that

    16  * accompanied this code).

    17  *

    18  * You should have received a copy of the GNU General Public License version

    19  * 2 along with this work; if not, write to the Free Software Foundation,

    20  * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA.

    21  *

    22  * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA

    23  * or visit www.oracle.com if you need additional information or have any

    24  * questions.

    25  */

    26 

    27 /* double lgamma(double x)

    28  * Return the logarithm of the Gamma function of x.

    29  *

    30  * Method: call __ieee754_lgamma_r

    31  */

    32 

    33 #include "fdlibm.h"

    34 

    35 extern int signgam;

    36 

    37 #ifdef __STDC__

    38         double lgamma(double x)

    39 #else

    40         double lgamma(x)

    41         double x;

    42 #endif

    43 {

    44 #ifdef _IEEE_LIBM

    45         return __ieee754_lgamma_r(x,&signgam);

    46 #else

    47         double y;

    48         y = __ieee754_lgamma_r(x,&signgam);

    49         if(_LIB_VERSION == _IEEE_) return y;

    50         if(!finite(y)&&finite(x)) {

    51             if(floor(x)==x&&x<=0.0)

    52                 return __kernel_standard(x,x,15); /* lgamma pole */

    53             else

    54                 return __kernel_standard(x,x,14); /* lgamma overflow */

    55         } else

    56             return y;

    57 #endif

    58 }