std::tgamma, std::tgammaf, std::tgammal
定义于头文件 <cmath> |
||
(1) | ||
float tgamma ( float num ); double tgamma ( double num ); |
(直到 C++23) | |
/*floating-point-type*/ tgamma ( /*floating-point-type*/ num ); |
(自 C++23 起) (constexpr 自 C++26 起) |
|
float tgammaf( float num ); |
(2) | (自 C++11 起) (constexpr 自 C++26 起) |
long double tgammal( long double num ); |
(3) | (自 C++11 起) (constexpr 自 C++26 起) |
SIMD 重载 (自 C++26 起) |
||
定义于头文件 <simd> |
||
template< /*math-floating-point*/ V > constexpr /*deduced-simd-t*/<V> |
(S) | (自 C++26 起) |
附加重载 (自 C++11 起) |
||
定义于头文件 <cmath> |
||
template< class Integer > double tgamma ( Integer num ); |
(A) | (constexpr 自 C++26 起) |
S) SIMD 重载对 v_num 执行元素级
std::tgamma 。
|
(自 C++26 起) |
A) 为所有整数类型提供了附加重载,它们被视为 double。
|
(自 C++11 起) |
目录 |
[编辑] 参数
num | - | 浮点值或整数值 |
[编辑] 返回值
如果未发生错误,则返回 num 的伽玛函数值,即 ∫∞
0tnum-1
e-t dt。
如果发生域错误,则返回实现定义的数值(在支持的情况下为 NaN)。
如果发生极点错误,则返回 ±HUGE_VAL、±HUGE_VALF
或 ±HUGE_VALL
。
如果由于溢出而发生范围错误,则返回 ±HUGE_VAL、±HUGE_VALF
或 ±HUGE_VALL
。
如果由于下溢而发生范围错误,则返回正确的值(舍入后)。
[编辑] 错误处理
错误报告按 math_errhandling 中的规定进行。
如果 num 为零或小于零的整数,则可能发生极点错误或域错误。
如果实现支持 IEEE 浮点算术 (IEC 60559),
- 如果参数为 ±0,则返回 ±∞ 并引发 FE_DIVBYZERO。
- 如果参数为负整数,则返回 NaN 并引发 FE_INVALID。
- 如果参数为 -∞,则返回 NaN 并引发 FE_INVALID。
- 如果参数为 +∞,则返回 +∞。
- 如果参数为 NaN,则返回 NaN。
[编辑] 注意
如果 num 是自然数,则 std::tgamma(num) 是 num - 1 的阶乘。如果参数是足够小的整数,许多实现会计算精确的整数域阶乘。
对于 IEEE 兼容类型 double,如果 0 < num && num < 1 / DBL_MAX 或如果 num > 171.7,则会发生溢出。
POSIX 要求,如果参数为零,则会发生极点错误,但当参数为负整数时,会发生域错误。它还规定,将来,对于负整数参数,域错误可能会被极点错误取代(在这种情况下,这些情况下的返回值将从 NaN 更改为 ±∞)。
在各种实现中都有一个名为 gamma
的非标准函数,但其定义不一致。例如,glibc 和 4.2BSD 版本的 gamma
执行 lgamma
,但 4.4BSD 版本的 gamma
执行 tgamma
。
不需要完全按照 (A) 的形式提供附加重载。它们只需要足以确保对于整数类型的参数 num,std::tgamma(num) 与 std::tgamma(static_cast<double>(num)) 具有相同的效果。
[编辑] 示例
#include <cerrno> #include <cfenv> #include <cmath> #include <cstring> #include <iostream> // #pragma STDC FENV_ACCESS ON int main() { std::cout << "tgamma(10) = " << std::tgamma(10) << ", 9! = " << 2 * 3 * 4 * 5 * 6 * 7 * 8 * 9 << '\n' << "tgamma(0.5) = " << std::tgamma(0.5) << ", sqrt(pi) = " << std::sqrt(std::acos(-1)) << '\n'; // special values std::cout << "tgamma(1) = " << std::tgamma(1) << '\n' << "tgamma(+Inf) = " << std::tgamma(INFINITY) << '\n'; // error handling errno = 0; std::feclearexcept(FE_ALL_EXCEPT); std::cout << "tgamma(-1) = " << std::tgamma(-1) << '\n'; if (errno == EDOM) std::cout << " errno == EDOM: " << std::strerror(errno) << '\n'; if (std::fetestexcept(FE_INVALID)) std::cout << " FE_INVALID raised\n"; }
可能的输出
tgamma(10) = 362880, 9! = 362880 tgamma(0.5) = 1.77245, sqrt(pi) = 1.77245 tgamma(1) = 1 tgamma(+Inf) = inf tgamma(-1) = nan errno == EDOM: Numerical argument out of domain FE_INVALID raised
[编辑] 参见
(C++11)(C++11)(C++11) |
伽玛函数的自然对数 (函数) |
(C++17)(C++17)(C++17) |
贝塔函数 (函数) |
C 文档 关于 tgamma
|
[编辑] 外部链接
Weisstein, Eric W. "Gamma Function." 来自 MathWorld — Wolfram Web 资源。 |