c++ - C/C++ で非線形反放物型偏微分方程式を解くためのパッケージ

Question

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++ パッケージは何ですか?

Answer 1

アクセスできる場合は、NAG ライブラリを試してください (高価です)。あなたの大学/会社/機関が適切な種類のライセンスを持っている場合、NAG のスタッフがお手伝いします。NAG の人に会いましたが、彼らは技術サポートに非常に真剣です。一部の人々は、NAG ライブラリを信頼しています。値段が高いのでほとんど使っていませんし、次の仕事で使えなくなるのではないかといつも心配しているので、頼りたくありません。

フリー ライブラリに関しては、Netlib に PDE ソルバーがいくつかありますが、私はそれらを使用したことがありません。もう 1 つのオプションは、C の Numerical Recipes です。これは実際にはお勧めしません。私が見た最悪の数値コードは、Numerical Recipes ルーチンをブラック ボックスとして使用していました。GSL は無料で、連立 ODE 系を解くために使用して非常に成功していますが、GSL には PDE ソルバーはありません。