2016-12-22 14 views
0

私はスプラインを自分のデータにフィットさせるためにpythonを使用しています。私はノットと係数のセットを通常の使用のためのCプログラムにエクスポートできる必要があります。私は、スプラインに合うようにscipy.interpolate.SmoothBivariateSplineを使い、メソッドget_coeffs()とget_knots()を使ってノットと係数のコピーを作成しています。Scipys SmoothBivariateSplineによって与えられるノットと係数のセットからスプライン値を計算する方法

最終的に、C環境でスプラインを使用できるように、評価メソッド(ev())を書き直す必要があります。このメソッドを書く方法について誰かが私にいくつかのリソースを案内できますか?コードを探し

答えて

0

dfitpack.bispevに二変量FITPACK基づくスプラインデリゲートの評価、 https://github.com/scipy/scipy/blob/master/scipy/interpolate/fitpack2.py#L844

後者の機能は、https://github.com/scipy/scipy/blob/master/scipy/interpolate/fitpack/fpbisp.f

を呼び出す同じ名前のFortranのルーチン、 https://github.com/scipy/scipy/blob/master/scipy/interpolate/fitpack/bispev.fのf2pyラッパー、 https://github.com/scipy/scipy/blob/master/scipy/interpolate/src/fitpack.pyf#L352、あります

一つの落とし穴は、エンドノットがget_knotsget_coefsによって処理されている正確にどのようにあると注意するhttps://github.com/scipy/scipy/blob/master/scipy/interpolate/fitpack2.py#L780

これはおそらくDierckxの本のすべての細部の詳細で議論されているでしょう。あるいは、計算ルーチンは、2Dのためのde Boorの計画のかなり簡単な一般化のように見えるかもしれません。しかし、悪魔は細部にいる。

関連する問題