命名空间
变体
操作

lgamma、lgammaf、lgammal

来自 cppreference.com
< c‎ | numeric‎ | math
 
 
 
常用数学函数
类型
(C99)(C99)    

(C99)(C99)    

函数
基本操作
(C99)
(C99)
(C99)
(C99)(C99)(C99)(C23)
最大/最小操作
(C99)
(C23)    
指数函数
(C23)
(C99)
(C99)
(C23)
(C23)
(C99)
(C99)(C23)
(C23)
(C23)
幂函数
(C99)
(C23)
(C23)
(C99)
(C23)
(C23)
三角函数和双曲函数
(C23)
(C23)
(C23)
(C23)
(C99)
(C99)
(C99)
误差函数和伽马函数
(C99)
(C99)
lgamma
(C99)
(C99)
最近整数浮点运算
(C99)(C99)(C99)
(C99)
(C99)(C99)(C99)
(C23)(C23)(C23)(C23)
浮点操作函数
(C99)(C99)
(C99)(C23)
(C99)
缩窄操作
(C23)
(C23)
(C23)
(C23)
(C23)
(C23)
量子和量子指数函数
十进制重新编码函数
全序和有效载荷函数
分类
(C99)
(C99)
(C99)
(C23)
宏常量
特殊浮点值
(C99)(C23)
参数和返回值
(C99)(C99)(C99)(C99)(C99)    
错误处理
(C99)    

 
定义在头文件 <math.h>
float       lgammaf( float arg );
(1) (自 C99)
double      lgamma( double arg );
(2) (自 C99)
long double lgammal( long double arg );
(3) (自 C99)
定义在头文件 <tgmath.h>
#define lgamma( arg )
(4) (自 C99)
1-3) 计算 arg伽马函数 的绝对值的自然对数。
4) 类型泛型宏:如果 arg 的类型为 long double,则调用 lgammal。否则,如果 arg 的类型为整数类型或类型 double,则调用 lgamma。否则,调用 lgammaf

内容

[编辑] 参数

arg - 浮点值

[编辑] 返回值

如果未发生错误,则返回 arg 的伽马函数的对数,即 log
e
|
0
targ-1
e-t dt|
的值。

如果发生极点错误,则返回 +HUGE_VAL+HUGE_VALF+HUGE_VALL

如果发生由于溢出导致的范围错误,则返回 ±HUGE_VAL±HUGE_VALF±HUGE_VALL

[编辑] 错误处理

错误报告方式如 math_errhandling 中所述。

如果 arg 为零或小于零的整数,则可能发生极点错误。

如果实现支持 IEEE 浮点运算 (IEC 60559),则

  • 如果参数为 1,则返回 +0。
  • 如果参数为 2,则返回 +0。
  • 如果参数为 ±0,则返回 +∞ 且引发 FE_DIVBYZERO
  • 如果参数为负整数,则返回 +∞ 且引发 FE_DIVBYZERO
  • 如果参数为 ±∞,则返回 +∞。
  • 如果参数为 NaN,则返回 NaN。

[编辑] 备注

如果 arg 为自然数,则 lgamma(arg)arg - 1 的阶乘的对数。

POSIX 版 lgamma 不是线程安全的:该函数的每次执行都会将 arg 的伽马函数的符号存储在静态外部变量 signgam 中。一些实现提供了 lgamma_r,它接受指向用户提供的 singgam 存储的指针作为第二个参数,并且是线程安全的。

在各种实现中,存在名为 gamma 的非标准函数,但其定义不一致。例如,glibc 和 4.2BSD 版 gamma 会执行 lgamma,但 4.4BSD 版 gamma 会执行 tgamma

[编辑] 示例

#include <errno.h>
#include <fenv.h>
#include <float.h>
#include <math.h>
#include <stdio.h>
// #pragma STDC FENV_ACCESS ON
 
int main(void)
{
    printf("lgamma(10) = %f, log(9!) = %f\n", lgamma(10),
                                              log(2 * 3 * 4 * 5 * 6 * 7 * 8 * 9));
    const double pi = acos(-1);
    printf("lgamma(0.5) = %f, log(sqrt(pi)) = %f\n", log(sqrt(pi)), lgamma(0.5));
    // special values
    printf("lgamma(1) = %f\n", lgamma(1));
    printf("lgamma(+Inf) = %f\n", lgamma(INFINITY));
 
    // error handling
    errno = 0; feclearexcept(FE_ALL_EXCEPT);
    printf("lgamma(0) = %f\n", lgamma(0));
    if (errno == ERANGE)
        perror("    errno == ERANGE");
    if (fetestexcept(FE_DIVBYZERO))
        puts("    FE_DIVBYZERO raised");
}

可能输出

lgamma(10) = 12.801827, log(9!) = 12.801827
lgamma(0.5) = 0.572365, log(sqrt(pi)) = 0.572365
lgamma(1) = 0.000000
lgamma(+Inf) = inf
lgamma(0) = inf
    errno == ERANGE: Numerical result out of range
    FE_DIVBYZERO raised

[编辑] 参考资料

  • C23 标准 (ISO/IEC 9899:2024)
  • 7.12.8.3 lgamma 函数 (p: 待定)
  • 7.25 类型泛型数学 <tgmath.h> (p: 待定)
  • F.10.5.3 lgamma 函数 (p: 待定)
  • C17 标准 (ISO/IEC 9899:2018)
  • 7.12.8.3 lgamma 函数 (p: 182)
  • 7.25 类型泛型数学 <tgmath.h> (p: 272-273)
  • F.10.5.3 lgamma 函数 (p: 383)
  • C11 标准 (ISO/IEC 9899:2011)
  • 7.12.8.3 lgamma 函数 (p: 250)
  • 7.25 类型泛型数学 <tgmath.h> (p: 373-375)
  • F.10.5.3 lgamma 函数 (p: 525)
  • C99 标准 (ISO/IEC 9899:1999)
  • 7.12.8.3 lgamma 函数 (p: 231)
  • 7.22 类型泛型数学 <tgmath.h> (p: 335-337)
  • F.9.5.3 lgamma 函数 (p: 462)

[编辑] 另请参阅

(C99)(C99)(C99)
计算伽玛函数
(函数) [编辑]
C++ 文档 for lgamma

[编辑] 外部链接

Weisstein, Eric W. "对数伽玛函数." 来自 Wolfram Web Resource 的 MathWorld。