命名空间
变体
操作

lgamma, lgammaf, lgammal

来自 cppreference.cn
< c‎ | numeric‎ | math
 
 
 
常用数学函数
函数
基本操作
(C99)
(C99)
(C99)
(C99)(C99)(C99)(C23)
最大值/最小值操作
(C99)
(C99)
指数函数
(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)(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)
(C99)
lgamma
(C99)
(C99)
类型
宏常量
特殊浮点值
(C99)(C23)
参数和返回值
错误处理
快速操作指示符
 
定义于头文件 <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 的伽玛函数的对数值,即 loge|
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 The lgamma functions (p: TBD)
  • 7.25 Type-generic math <tgmath.h> (p: TBD)
  • F.10.5.3 The lgamma functions (p: TBD)
  • C17 标准 (ISO/IEC 9899:2018)
  • 7.12.8.3 The lgamma functions (p: 182)
  • 7.25 Type-generic math <tgmath.h> (p: 272-273)
  • F.10.5.3 The lgamma functions (p: 383)
  • C11 标准 (ISO/IEC 9899:2011)
  • 7.12.8.3 The lgamma functions (p: 250)
  • 7.25 Type-generic math <tgmath.h> (p: 373-375)
  • F.10.5.3 The lgamma functions (p: 525)
  • C99 标准 (ISO/IEC 9899:1999)
  • 7.12.8.3 The lgamma functions (p: 231)
  • 7.22 Type-generic math <tgmath.h> (p: 335-337)
  • F.9.5.3 The lgamma functions (p: 462)

[编辑] 参见

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

[编辑] 外部链接

Weisstein, Eric W. "Log Gamma Function." 来自 MathWorld — Wolfram Web 资源。