奇数サイト三角格子スピン模型のMPI並列がおかしい about hphi HOT 2 CLOSED

@ultimatile
報告ありがとうございます。
本件、全対角化で確認してみたところ、基底状態と第一励起状態が縮退しておりました。
縮退している場合は線形結合を適切に取らないと、固有ベクトルが一致しなくなります。

i=    0 Energy= -5.250000 N=  9.000000 Sz=  0.500000 S2=  0.750000 Doublon=  0.000000

i=    1 Energy= -5.250000 N=  9.000000 Sz=  0.500000 S2=  0.750000 Doublon=  0.000000

乱数の種(初期状態)が一致している場合には一致するのですが、
HΦではMPIを使用する場合に各プロセス毎に乱数を発生させる仕様となっています。
そのため、実質的にシングルプロセスとMPIプロセスで異なった初期状態を取るため、この現象が起きてしまいます。
サイズが大きい場合に、縮退しているか判定するため、CG法の実装をHΦでは行っておりますので、よろしければご活用ください。

from hphi.

ご返事ありがとうございます.縮退の可能性を失念しておりました.お手数をおかけして申し訳ありませんでした.

from hphi.