使用 GSL 进行复杂数值积分

计算科学 C++ 一体化 复分析
2021-12-15 07:42:40

我想使用 GSL 库在 C++ 中编写一个集成例程,但要使用复杂的函数。

我应该如何拆分我的被积函数以在其上应用 gsl_integration_qag 函数。

我猜只是将实部和虚部自己整合是不正确的。

我有以下函数: 其中是一些实常数。我稍后打算集成的函数使用谐波和作为基础,但只是它的总和。

S(a,N)=S(a,)n=11n+N
S(a,)

我打算做的积分是梅林变换,它由以下公式给出:

f(x)=12πiCdnxnfM(n).

任何帮助都非常受欢迎。

编辑:关于复杂数值积分的 C++ 程序有什么想法吗?

1个回答

我认为要走的路是用实参数对积分进行参数化,并分别整合实部和虚部。实现这一点很微妙,实现取决于具体情况。