numpyとの正確な合計

Question

数値のxとyデータをnumpyの配列として取り、それをxの領域に（数値的に）積分するプログラムを書いています。私はこれをして、今は精度について疑問に思っています。

私は、math.fsum()がsum()でない特定のエラーを修正していることに気がつきました。numpy.sumもこれらの修正を行っていますか？

sum = numpy.sum(array)などと書いても、結果（合計）はできるだけ正確ですか？私はそれ自身がいくつかの操作の積である何千ものデータ点を合計しているので、小さな誤差が蓄積する可能性があります。私はPythonで十分精通していないので、小数点以下の桁までの精度が重要なときに使用する必要がある、より正確な数学モジュールがあるかどうかを知ることができます。

Answer 1

いいえ、numpy.sumはありませんは修正math.fsumを行います。それはpairwise summationを使用しています。これは、テクニックmath.fsumが使用するほど良くはありませんが、それほど悪くはありません。

Answer 2

scipyのダウンロードで、いくつかの数値積分の実装、numpyのLIBSは、あなたがポイントの配列に適用することができますが、より良い方法をしたい場合

一部（ペア）があります。 numpy.trapz

か、補間を使用することができ、小さな加数がたくさんある場合scipy.integrate