問題タブ [mpmath]

https://code.google.com/p/mpmath/issues/detail?id=239

Python バックエンドで mpmath データをピクルし、gmp バックエンドでそれらを取得しようとする問題を経験しました。これにより、エラーが発生します。

これは入っlibmpf.pyていて、これが呼び出されたときですMPZ(man, 16)

リンクに従って、それを代入MPZしlongてもエラーは発生しなくなります。誰かがこれについての洞察を持っていますか。これは私が躊躇しているハックだったようです。

私の問題に対するベストプラクティスの解決策は何ですか?

ありがとうございました

mpmath.findroot の多次元ニュートン法を使用して、連立方程式を数値的に高精度に解こうとしています。システムの例を次に示します。

残念ながら、解決策は得られませんが、エラーがスローされます

ZeroDivisionError: 行列は数値的に特異です

これは、findroot がヤコビアンを計算しようとしているために発生していますか? これは、システムが過小決定されていることを意味しますか?

scipy の fmin_cg を使用して出発点を見つけましたが、解をより高い精度に磨きたいと考えています。fmin_cg の関数として、関数 f の 27 エントリの二乗和を最小化しました。

mpmath.findroot の問題を回避できない場合、このシステムを高精度で解決するためのより良い方法はありますか?

ベネットのアルゴリズムを使用してフィボナッチ数を計算しようとしている以前の質問を続けようとしています。任意の精度で作業するには、mpmath. ただし、実装は特定の値を超えると失敗するようです。たとえば、99 番目の値は次のようになります。

218922995834555891712

これは ( ref )である必要があります。

218922995834555169026

これが私のコードです：

これは、Python で JIT を使用する最初の試みであり、スピードアップしたいユース ケースです。numba について少し読んだところ、十分に単純に思えましたが、次のコードでは速度が向上しませんでした。私が犯しているかもしれない明らかな間違いを許してください。

また、cython の基本的なチュートリアルが示唆することを実行しようとしましたが、時間に違いはありませんでした。 http://docs.cython.org/src/tutorial/cython_tutorial.html

変数を宣言するようなことをしなければならないと思いますか？他のライブラリを使用しますか? すべてに排他的に for ループを使用しますか? 参考になるガイダンスや例をいただければ幸いです。

たとえば、以前の質問から、numpy と比較して mpmath の Elementwise operations が遅いことと、 mpmath の代わりに gmpy を使用した方がはるかに高速であるという解決策を知っています。

これらのmpmathの間に、anacondaがインストールされ、Pythonの多くのライブラリがインストールされています。powerlaw パッケージを実行しようとすると、次のエラーが表示されます。

さらに、単純にやろうとすると import mpmath、同じエラーが発生します。

私はFedoraを使用しています。すでに試しyum remove python-mpmathてからyum install python-mpmath. しかし、私は同じエラーが発生しています。

他の質問で、これは複数のパスがあり、それらの両方を sys.path=['' など] で追加する必要があることが原因である可能性があることを見てきました。

"sys.path=['', etc]" とはどういう意味ですか?

私がする時：

私は得る：