std::lgamma, std::lgammaf, std::lgammal
| ヘッダー <cmath> で定義 |
||
| (1) | ||
float lgamma ( float num ); double lgamma ( double num ); |
(C++23まで) | |
| /*浮動小数点数型*/ lgamma ( /*浮動小数点型*/ num ); |
(C++23から) (C++26 以降 constexpr) |
|
float lgammaf( float num ); |
(2) | (C++11以降) (C++26 以降 constexpr) |
long double lgammal( long double num ); |
(3) | (C++11以降) (C++26 以降 constexpr) |
| SIMDオーバーロード (C++26以降) |
||
| ヘッダー <simd> で定義 |
||
| template< /*math-floating-point*/ V > constexpr /*deduced-simd-t*/<V> |
(S) | (C++26以降) |
| 追加のオーバーロード (C++11以降) |
||
| ヘッダー <cmath> で定義 |
||
template< class Integer > double lgamma ( Integer num ); |
(A) | (C++26 以降 constexpr) |
std::lgamma のオーバーロードを提供します。(since C++23)|
S) SIMD オーバーロードは、v_num に対して要素ごとの
std::lgamma を実行します。
|
(C++26以降) |
|
A) すべての整数型に対する追加のオーバーロードが提供されます。これらは double として扱われます。
|
(C++11以降) |
目次 |
[編集] パラメータ
| num | - | 浮動小数点数または整数値 |
[編集] 戻り値
エラーが発生しなかった場合、num のガンマ関数の対数の値、すなわち loge|∫∞
0tnum-1
e-t dt|が返されます。
極エラーが発生した場合、+HUGE_VAL、+HUGE_VALF、または +HUGE_VALL が返されます。
オーバーフローによる範囲エラーが発生した場合、±HUGE_VAL、±HUGE_VALF、または ±HUGE_VALL のいずれかを返します。
[編集] エラー処理
エラーは math_errhandling で指定された通りに報告される。
num がゼロまたは負の整数である場合、極エラーが発生する可能性があります。
実装がIEEE浮動小数点算術 (IEC 60559) をサポートしている場合、
- If the argument is 1, +0 is returned.
- 引数が 2 の場合、+0 が返されます。
- 引数が ±0 の場合、+∞ が返され、FE_DIVBYZERO が発生します。
- 引数が負の整数の場合、+∞ が返され、FE_DIVBYZERO が発生します。
- 引数が ±∞ の場合、+∞ が返されます。
- 引数が NaN の場合、NaN が返されます。
[編集] 注
num が自然数である場合、std::lgamma(num) は num - 1 の階乗の対数です。
POSIX 版の lgamma はスレッドセーフではありません。関数の各実行は、ガンマ関数の num における符号を静的外部変数 signgam に格納します。一部の実装では、signgam 用のユーザー提供ストレージへのポインタを 2 番目のパラメータとして受け取る lgamma_r が提供されており、これはスレッドセーフです。
様々な実装に gamma という名前の非標準関数がありますが、その定義は一貫性がありません。例えば、glibc と 4.2BSD 版の gamma は lgamma を実行しますが、4.4BSD 版の gamma は tgamma を実行します。
追加のオーバーロードは、(A) とまったく同じように提供される必要はありません。整数型の引数 num に対して、std::lgamma(num) が std::lgamma(static_cast<double>(num)) と同じ効果を持つことを保証するだけで十分です。
[編集] 例
#include <cerrno> #include <cfenv> #include <cmath> #include <cstring> #include <iostream> // #pragma STDC FENV_ACCESS ON const double pi = std::acos(-1); // or std::numbers::pi since C++20 int main() { std::cout << "lgamma(10) = " << std::lgamma(10) << ", log(9!) = " << std::log(std::tgamma(10)) << ", exp(lgamma(10)) = " << std::exp(std::lgamma(10)) << '\n' << "lgamma(0.5) = " << std::lgamma(0.5) << ", log(sqrt(pi)) = " << std::log(std::sqrt(pi)) << '\n'; // special values std::cout << "lgamma(1) = " << std::lgamma(1) << '\n' << "lgamma(+Inf) = " << std::lgamma(INFINITY) << '\n'; // error handling errno = 0; std::feclearexcept(FE_ALL_EXCEPT); std::cout << "lgamma(0) = " << std::lgamma(0) << '\n'; if (errno == ERANGE) std::cout << " errno == ERANGE: " << std::strerror(errno) << '\n'; if (std::fetestexcept(FE_DIVBYZERO)) std::cout << " FE_DIVBYZERO raised\n"; }
出力
lgamma(10) = 12.8018, log(9!) = 12.8018, exp(lgamma(10)) = 362880
lgamma(0.5) = 0.572365, log(sqrt(pi)) = 0.572365
lgamma(1) = 0
lgamma(+Inf) = inf
lgamma(0) = inf
errno == ERANGE: Numerical result out of range
FE_DIVBYZERO raised[編集] 関連項目
| (C++11)(C++11)(C++11) |
ガンマ関数 (関数) |
| C ドキュメント (lgamma)
| |
[編集] 外部リンク
| Weisstein, Eric W. "Log Gamma Function." MathWorld — A Wolfram Web Resource より。 |