diff --git a/reference/cmath/lgamma.md b/reference/cmath/lgamma.md index 60de4faba..03027e27d 100644 --- a/reference/cmath/lgamma.md +++ b/reference/cmath/lgamma.md @@ -60,7 +60,8 @@ namespace std { - `x = -∞` の場合、戻り値は `+∞` となる。 - `x = +∞` の場合、戻り値は `+∞` となる。 - C++23では、(1)、(2)、(3)が(4)に統合され、拡張浮動小数点数型を含む浮動小数点数型へのオーバーロードとして定義された - +- この関数はガンマ関数 ([`tgamma`](tgamma.md)) がオーバーフローするような場合に使う。 + 具体例については[ガンマ関数の備考](tgamma.md#remarks-lgamma)を参照のこと。 ## 例 ```cpp example diff --git a/reference/cmath/tgamma.md b/reference/cmath/tgamma.md index 4dcbd333a..df93b11bb 100644 --- a/reference/cmath/tgamma.md +++ b/reference/cmath/tgamma.md @@ -64,6 +64,32 @@ $$ \Gamma (x) = \int_0^\infty t^{x-1} e^{-t} dt $$ - `gamma` という関数は既にあったが処理系によって定義が違ったため、本当の (true) ガンマ関数 `tgamma` と名付けられた。 - C++23では、(1)、(2)、(3)が(4)に統合され、拡張浮動小数点数型を含む浮動小数点数型へのオーバーロードとして定義された +### lgamma との使い分け +ガンマ関数は急激に増加し容易にオーバーフローするので、代わりにガンマ関数の結果を自然対数で返す関数 [`lgamma`](lgamma.md) を用いた方が良いことが多くある。 +例えばガンマ関数の比を計算する場合には、 ガンマ関数の対数の差を取ってから指数関数 [`std::exp`](exp.md) を適用するのが懸命である。 + +$$ \frac{\Gamma(2026)}{\Gamma(2025)} = \exp[\ln\Gamma(2026) - \ln\Gamma(2025)] $$ + +```cpp example +#include +#include +int main() { + std::cout << std::tgamma(2026.0) / std::tgamma(2025.0) << std::endl; + std::cout << std::exp(std::lgamma(2026.0) - std::lgamma(2025.0)) << std::endl; +} +``` +* std::tgamma[color ff0000] +* std::lgamma[color 0000ff][link tgamma.md] + +出力例 +``` +-nan +2025 +``` + +上の結果では、直接ガンマ関数を計算した場合はオーバーフローによって inf / inf となり最終結果が -nan になっているが、`lgamma` を使った場合には正しい値が計算できている。 +ただし、`lgamma` は飽くまでガンマ関数の「絶対値」の対数であることに注意する。 +ガンマ関数の引数が負になる場合はガンマ関数が負の値を取りうるので符号は別に求める必要がある。 ## 例 ```cpp example