名前空間
変種
操作

std::lgamma, std::lgammaf, std::lgammal

From cppreference.com
< cpp‎ | numeric‎ | math
 
 
 
 
ヘッダー <cmath> で定義
(1)
float       lgamma ( float num );

double      lgamma ( double num );

long double lgamma ( long 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> で定義
template< /*math-floating-point*/ V >

constexpr /*deduced-simd-t*/<V>

            lgamma ( const V& v_num );
(S) (C++26以降)
ヘッダー <cmath> で定義
template< class Integer >
double      lgamma ( Integer num );
(A) (C++26 以降 constexpr)
1-3) numガンマ関数 の絶対値の自然対数を計算します。ライブラリは、パラメータの型として、すべての cv 修飾されていない浮動小数点型に対する std::lgamma のオーバーロードを提供します。(since C++23)
S) SIMD オーバーロードは、v_num に対して要素ごとの std::lgamma を実行します。
(定義については、math-floating-point および deduced-simd-t を参照してください。)
(C++26以降)
A) すべての整数型に対する追加のオーバーロードが提供されます。これらは double として扱われます。
(C++11以降)

目次

[編集] パラメータ

num - 浮動小数点数または整数値

[編集] 戻り値

エラーが発生しなかった場合、num のガンマ関数の対数の値、すなわち loge|
0
tnum-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 版の gammalgamma を実行しますが、4.4BSD 版の gammatgamma を実行します。

追加のオーバーロードは、(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)
ガンマ関数
(関数) [編集]

[編集] 外部リンク

Weisstein, Eric W. "Log Gamma Function." MathWorld — A Wolfram Web Resource より。
English 日本語 中文(简体) 中文(繁體)