問題タブ [qutip]

Qutip の関数 bloch() を使用して、ブロッホ球にライブ データをプロットしようとしています。

これまでのところ、そこに b.show() があると、コードは常に中断します。

同様の問題に対する多くの解決策をオンラインで見つけましたが、それらのほとんどは、bloch クラスでは機能しないように見える matplotlib.draw() のような直接の matplotlib コマンドを使用しています。次に、たとえば Tk や GTKagg を利用する他のソリューションがあります (たとえば、 https://stackoverflow.com/a/15742183/3276735またはmatplotlib を使用した while ループでのリアルタイム プロット) 。

ブロッホクラスで同じ問題に対処する方法を教えてください。

編集:これは最小限の例です:

基本的に、できればループで、一度に 1 つのポイントでプロットを更新したいと考えています。私の目標は、ファイルから読み取る必要があるプロットにライブ データを表示することです。

QuTiP でいくつかのコードを実行しようとしていますが、関数を並行して実行するとparforエラーが発生します。

エラー：

これが私のコードです：

ここで何が悪いのか誰か知っていますか？

量子計算の数値シミュレーションを行っています。大きなエルミート行列 (~2^14 行/列) の固有ベクトルを見つけたいと考えています。

24 コア/48 スレッドの XEON マシンで実行しています。コードはもともと Qutip ライブラリの助けを借りて作成されました。含まれているeigenstates()関数はマシン上で単一のスレッドしか使用しないことがわかったので、それを行うためのより高速な方法を見つけようとしています。

scipy.linalg eig()andeigh()関数とscipy.sparse.linalg eig()andを使用してみましeigh()たが、どちらも Qutip に組み込まれている関数よりも遅いようです。

slepc4py を使用すると速度が向上する可能性があるという提案をいくつか見ましたが、パッケージのドキュメントが非常に不足しているようです。numpy の複雑な配列を SLEPC 行列に変換する方法がわかりません。