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

Определено в заголовочном файле `<cmath>`
double comp_ellint_1 ( double k ); float comp_ellint_1 ( float k ); long double comp_ellint_1 ( long double k );	(1)	(начиная с C++17) (до C++23)
/* floating-point-type / comp_ellint_1( / floating-point-type */ k );	(1)		(начиная с 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. Библиотека предоставляет перегрузки std::comp_ellint_1 для всех типов с плавающей точкой без квалификаторов cv в качестве типа параметра k.(начиная с C++23)

A) Предоставляются дополнительные перегрузки для всех целочисленных типов, которые рассматриваются как double.

Параметры

эллиптический модуль или эксцентриситет (значение с плавающей точкой или целочисленное значение)

Возвращаемое значение

Если ошибок не возникает, возвращается значение полного эллиптического интеграла первого рода для 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)), где K — std::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

См. также

ellint_1ellint_1fellint_1l

(C++17)(C++17)(C++17)

(неполный) эллиптический интеграл первого рода
(функция)

Внешние ссылки

Вайсстайн, Эрик У. «Полный эллиптический интеграл первого рода». Из MathWorld — веб-ресурса Wolfram.

© cppreference.com
Licensed under the Creative Commons Attribution-ShareAlike Unported License v3.0.
https://en.cppreference.com/cpp/numeric/special_functions/comp_ellint_1