具有多对数函数的库

计算科学 C++ 正则 C 图书馆
2021-12-21 08:53:00

我正在寻找一个可靠的库,我可以在其中找到用于求解非线性方程的多对数函数和牛顿/割线法基本上我有这样的事情:

f(x) = x - A*PolyLog(3/2, B*Exp(-t*x))
f(x) = 0

f(x,y,z) = 0
g(x,y,z) = 0
h(x,y,z) = 0

涉及多对数函数。Polylogarithm 函数的实现需要类似于 Mathematica 或 Python 的实现(可以返回复数)并定义为非整数值(这里我们有 3/2)。因此,牛顿法也应该适用于复数。

该库应使用FortranC/C++编写。

2个回答

实现多对数函数的一种方法可以在“数值食谱:科学计算的艺术”一书中找到。本书包括评估特殊函数和寻找方程零点的一些理论和实践方面。如果你愿意,你也可以在他们的网页上购买代码。

这里提供了一个很好的算法,包括分析延续我很少建议自己编写这些东西,但在这种情况下它非常简单。