?
This document uses PHP Chinese website manual Release
在頭文件<math.h>中定義 | ||
---|---|---|
float lgammaf( float arg ); | (1) | (since C99) |
double lgamma( double arg ); | (2) | (since C99) |
long double lgammal( long double arg ); | (3) | (since C99) |
Defined in header <tgmath.h> | ||
#define lgamma( arg ) | (4) | (since C99) |
1-3)計(jì)算的絕對(duì)值的自然對(duì)數(shù)伽馬函數(shù)的arg
。
4)類型 - 通用宏:如果arg
有類型long double
,lgammal
被調(diào)用。否則,如果arg
有整數(shù)類型或類型double
,lgamma
則調(diào)用。否則,lgammaf
被調(diào)用。
ARG | - | 浮點(diǎn)值 |
---|
如果沒(méi)有發(fā)生錯(cuò)誤,伽瑪函數(shù)的對(duì)數(shù)值arg
即log
e|∫∞
0_t_arg-1
e -td t |返回。
如果發(fā)生極錯(cuò)誤+HUGE_VAL
,+HUGE_VALF
或+HUGE_VALL
返回。
如果范圍誤差由于發(fā)生溢出,±HUGE_VAL
,±HUGE_VALF
,或±HUGE_VALL
返回。
按照math_errhandling中的指定報(bào)告錯(cuò)誤。
如果arg
是零或者是小于零的整數(shù),則可能會(huì)發(fā)生極點(diǎn)錯(cuò)誤。
如果實(shí)現(xiàn)支持IEEE浮點(diǎn)運(yùn)算(IEC 60559),
如果參數(shù)是1,則返回+0
如果參數(shù)是2,則返回+0
如果參數(shù)為±0,則返回+∞并將FE_DIVBYZERO
其升高
如果參數(shù)是負(fù)整數(shù),+∞被返回,并FE_DIVBYZERO
提高
如果參數(shù)是±∞,則返回+∞。
如果參數(shù)是NaN,則返回NaN
如果arg
是自然數(shù),lgamma(arg)
則是階乘的對(duì)數(shù)arg-1
。
lgamma的POSIX版本不是線程安全的:函數(shù)的每次執(zhí)行都會(huì)將gamma函數(shù)的符號(hào)存儲(chǔ)arg
在靜態(tài)外部變量中signgam
。一些實(shí)現(xiàn)提供lgamma_r
了將用戶提供的存儲(chǔ)指針作為第二個(gè)參數(shù)并且是線程安全的。
gamma
在各種實(shí)現(xiàn)中有一個(gè)非標(biāo)準(zhǔn)函數(shù),但其定義不一致。例如,glibc和4.2BSD版本gamma
執(zhí)行lgamma
,但4.4BSD版本gamma
執(zhí)行tgamma
。
#include <stdio.h>#include <math.h>#include <float.h>#include <errno.h>#include <fenv.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)); 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.801827lgamma(0.5) = 0.572365, log(sqrt(pi)) = 0.572365lgamma(1) = 0.000000lgamma(+Inf) = inflgamma(0) = inf errno == ERANGE: Numerical result out of range FE_DIVBYZERO raised
C11標(biāo)準(zhǔn)(ISO / IEC 9899:2011):
7.12.8.3 lgamma函數(shù)(p:250)
7.25類型通用數(shù)學(xué)<tgmath.h>(p:373-375)
F.10.5.3 lgamma函數(shù)(p:525)
C99標(biāo)準(zhǔn)(ISO / IEC 9899:1999):
7.12.8.3 lgamma函數(shù)(p:231)
7.22類型通用數(shù)學(xué)<tgmath.h>(p:335-337)
F.9.5.3 lgamma函數(shù)(p:462)