std::riemann_zeta, std::riemann_zetaf, std::riemann_zetal

来自cppreference.com
 
 
实验性
技术规范
文件系统库 (文件系统 TS)
库基础 (库基础 TS)
库基础 2 (库基础 TS v2)
库基础 3 (库基础 TS v3)
并行扩展 (并行 TS)
并行扩展 2 (并行 TS v2)
并发扩展 (并发 TS)
并发扩展 2 (并发 TS v2)
概念 (概念 TS)
范围 (范围 TS)
反射 (反射 TS)
数学特殊函数 (特殊函数 TR)
实验性非 TS 功能特性
模式匹配
线性代数
std::execution
契约
2D 图形
 
 
double      riemann_zeta( double arg );

double      riemann_zeta( float arg );
double      riemann_zeta( long double arg );
float       riemann_zetaf( float arg );

long double riemann_zetal( long double arg );
(1)
double      riemann_zeta( IntegralType arg );
(2)
1) 计算 arg黎曼 zeta 函数的值。
2) 一组接受任何整数类型实参的重载或函数模板。等价于将实参转型为 double 后的 (1)

与所有特殊函数一样,仅当实现将 __STDCPP_MATH_SPEC_FUNCS__ 定义为至少 201003L 的值,且用户在包含任何标准库头文件前定义了 __STDCPP_WANT_MATH_SPEC_FUNCS__ 时,riemann_zeta 才保证在 <cmath> 中可用。

参数

n - 多项式的度数,无符号整数类型的值
arg - 浮点或整数类型的值

返回值

如果未发生错误,则返回 arg 的黎曼 zeta 函数的值,即 ζ(arg),在整个实轴上定义:

  • 对于 arg > 1Σ
    n=1
    n-arg
  • 对于 0 ≤ arg ≤ 1
    1
    1 - 21-arg
    Σ
    n=1
    (-1)n-1
    n-arg
  • 对于 arg < 02arg
    πarg-1
    sin(
    πarg
    2
    )Γ(1 − arg)ζ(1 − arg)

错误处理

根据 math_errhandling 的规定进行错误报告。

  • 如果实参为 NaN,则返回 NaN 但不报告定义域错误。

注解

不支持 TR 29124 但支持 TR 19768 的实现,在头文件 tr1/cmath 和命名空间 std::tr1 中提供此函数。

此函数的实现也在 boost.math 中可用

示例

(以 gcc 6.0 运行)

#define __STDCPP_WANT_MATH_SPEC_FUNCS__ 1
#include <cmath>
#include <iostream>
 
int main()
{
    // spot checks for well-known values
    std::cout << "ζ(-1) = " << std::riemann_zeta(-1) << '\n'
              << "ζ(0) = " << std::riemann_zeta(0) << '\n'
              << "ζ(1) = " << std::riemann_zeta(1) << '\n'
              << "ζ(0.5) = " << std::riemann_zeta(0.5) << '\n'
              << "ζ(2) = " << std::riemann_zeta(2) << ' '
              << "(π²/6 = " << std::pow(std::acos(-1), 2) / 6 << ")\n";
}

输出:

ζ(-1) = -0.0833333
ζ(0) = -0.5
ζ(1) = inf
ζ(0.5) = -1.46035
ζ(2) = 1.64493 (π²/6 = 1.64493)

外部链接

Weisstein, Eric W. "Riemann Zeta Function." From MathWorld--A Wolfram Web Resource.