std::comp_ellint_1, std::comp_ellint_1f, std::comp_ellint_1l
来自cppreference.com
< cpp | numeric | special functions
在标头 <cmath> 定义 | ||
(1) | ||
double comp_ellint_1 (double k ); float comp_ellint_1 (float k ); | (C++17 起) (C++23 前) | |
/* 浮点数类型 */ comp_ellint_1(/* 浮点数类型 */ k ); | (C++23 起) | |
float comp_ellint_1f(float k ); | (2) | (C++17 起) |
longdouble comp_ellint_1l(longdouble k ); | (3) | (C++17 起) |
在标头 <cmath> 定义 | ||
template<class Integer > double comp_ellint_1 ( Integer k ); | (A) | (C++17 起) |
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)) 的效果相同即可。
[编辑]示例
给定重力加速度 g,及初始角 θ,那么长度是 l 的钟摆周期等于 4⋅√l/g⋅K(sin(θ/2)),其中 K 是 std::comp_ellint_1
。
运行此代码
#include <cmath>#include <iostream>#include <numbers> int main(){constexprdouble π{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'<<"在 10° 初始角的长度 1m 的钟摆周期 ≈ "<<4*std::sqrt(1/9.80665)* std::comp_ellint_1(std::sin(π /18/2))<<"s,\n""作为对比,线性估计 ≈ "<<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 在 10° 初始角的长度 1m 的钟摆周期 ≈ 2.01024s, 作为对比,线性估计 ≈ 2.00641
[编辑]参阅
(C++17)(C++17)(C++17) | 第一类(不完全)椭圆积分 (函数) |
[编辑]外部链接
Weisstein, Eric W. “第一类完全椭圆积分”来自 MathWorld--A Wolfram Web Resource。 |