私は重力ポテンシャル式(から15)を得て、軌道を計算するためにローカル勾配である重力を評価する必要があります。 P230、P4、P6の導関数を数万回の単値で評価する。Pythonでルジャンドル多項式の派生を取得する
私はthis questionで式を使用して、それを計算することができますが、明示的に有限の差として微分を評価する私を必要としない誘導体についてのpythonを依頼する方法がある場合、私は思ったんだけど。
SciPyでこれを自動的に行うことができませんでした。 numpy.polynomial.legendre.Legendre
にはderiv()
メソッドがありますが、多項式クラスでの操作はありません。
低次のルジャンドル多項式の一次導関数を数値的に積分するために一度に一度評価するのが最も速い方法は何でしょうか?
私が間違っていないとSciPyの一部である 'sympy'モジュールを試すことができます。それはあなたが本当に簡単な差別化や統合などを行うことを可能にしますが、まず無限の和を扱う必要があります.SymPyはそれを処理できるだけでなく、限界やそのようなものを評価するのにもかなり良いです。 – ForceBru
@ForceBruありがとう私は 'sympy'がここで仕事のために正しいとは思わない。私は自分自身で象徴的な部分を行うことができます - 答えはリンクされた質問の中にあります。私が必要とするのは速い数値評価であり、 'sympy'はそれをしません。 – uhoh