名前空間
変種
操作

std::comp_ellint_1, std::comp_ellint_1f, std::comp_ellint_1l

From cppreference.com
 
 
 
 
ヘッダー <cmath> で定義
(1)
double      comp_ellint_1 ( double k );

float       comp_ellint_1 ( float k );

long double comp_ellint_1 ( long double k );
(C++17以降)
(C++23まで)
/* floating-point-type */ comp_ellint_1( /* floating-point-type */ k );
(C++23から)
float       comp_ellint_1f( float k );
(2) (C++17以降)
long double comp_ellint_1l( long double k );
(3) (C++17以降)
ヘッダー <cmath> で定義
template< class Integer >
double      comp_ellint_1 ( Integer k );
(A) (C++17以降)
1-3) k第一種完全楕円積分を計算します。ライブラリは、パラメータ k の型として、すべての cv-修飾されていない浮動小数点型に対する std::comp_ellint_1 のオーバーロードを提供します。(since C++23)
A) すべての整数型に対する追加のオーバーロードが提供されます。これらは double として扱われます。

目次

[編集] パラメータ

k - 楕円モジュラスまたは離心率 (浮動小数点数または整数値)

[編集] 戻り値

エラーが発生しなかった場合、k の第一種完全楕円積分の値、すなわち std::ellint_1(k, π/2) が返されます。

[編集] エラー処理

math_errhandling で指定されたとおりにエラーが報告される場合があります。

  • 引数がNaNの場合、NaNが返され、ドメインエラーは報告されません。
  • |k|>1 の場合、ドメインエラーが発生する可能性があります。

[編集] 注記

C++17をサポートしないが、ISO 29124:2010をサポートする実装では、実装によって__STDCPP_MATH_SPEC_FUNCS__が少なくとも201003L以上の値に定義され、ユーザーが標準ライブラリヘッダをインクルードする前に__STDCPP_WANT_MATH_SPEC_FUNCS__を定義した場合、この関数が提供されます。

ISO 29124:2010をサポートしないが、TR 19768:2007 (TR1) をサポートする実装では、tr1/cmathヘッダおよびstd::tr1名前空間でこの関数が提供されます。

この関数の実装は boost.math でも利用可能です。

追加のオーバーロードは (A) とまったく同じように提供される必要はありません。整数型の引数 num に対して、std::comp_ellint_1(num)std::comp_ellint_1(static_cast<double>(num)) と同じ効果を持つことを保証するだけで十分です。

[編集]

長さ l、重力加速度 g、初期角 θ の単振り子の周期は、4⋅l/g⋅K(sin(θ/2)) です。ここで Kstd::comp_ellint_1 です。

#include <cmath>
#include <iostream>
#include <numbers>
 
int main()
{
    constexpr double π{std::numbers::pi};
 
    std::cout << "K(0) ≈ " << std::comp_ellint_1(0) << '\n'
              << "π/2 ≈ " << π / 2 << '\n'
              << "K(0.5) ≈ " << std::comp_ellint_1(0.5) << '\n'
              << "F(0.5, π/2) ≈ " << std::ellint_1(0.5, π / 2) << '\n'
              << "The period of a pendulum length 1m at 10° initial angle ≈ "
              << 4 * std::sqrt(1 / 9.80665) * std::comp_ellint_1(std::sin(π / 18 / 2))
              << "s,\n" "whereas the linear approximation gives ≈ "
              << 2 * π * std::sqrt(1 / 9.80665) << '\n';
}

出力

K(0) ≈ 1.5708
π/2 ≈ 1.5708
K(0.5) ≈ 1.68575
F(0.5, π/2) ≈ 1.68575
The period of a pendulum length 1 m at 10° initial angle ≈ 2.01024s,
whereas the linear approximation gives ≈ 2.00641

[編集] 関連項目

(C++17)(C++17)(C++17)
(不完全)第一種楕円積分
(関数) [編集]

[編集] 外部リンク

Weisstein, Eric W. "Complete Elliptic Integral of the First Kind." From MathWorld — A Wolfram Web Resource.
English 日本語 中文(简体) 中文(繁體)