我想为两个变量函数PDE解决以下f(q,y)
d f(q,y) / dq + 1/2 (d^2f(q,y)/dy^2 + x(q)*(df(q,y)/dy)^2) = 0,
在区间
-\inf < y < \inf, 0<=q<=1中边界条件为
f(1,y) = g(y)，其中g(y)是已知函数。用数字方式求解该方程的最佳C / C++软件包是什么？
最佳答案
如果可以访问NAG库(价格昂贵)，请尝试使用它们。如果您的大学/公司/研究所拥有正确的许可证类型，NAG的人员将为您提供帮助。我遇到了NAG的一个人，他们非常重视技术支持。有些人向NAG图书馆宣誓；我几乎没有使用过它们，因为它们是如此昂贵，而且我始终担心我的下一份工作将无法使用它们，因此我不想依赖它们。
在免费库方面，Netlib中有一些PDE求解器，但我从未使用过。另一个选择是C语言的数字食谱，我实际上会建议不要这样做。我见过的最差的数字代码已将数值食谱例程用作黑匣子。 GSL是免费的，我已经非常成功地将其用于求解耦合ODE的系统，但是GSL中没有PDE求解器。