スピン系で全対角化がセグメンテーション違反する

[yomichi]

スピン系の計算で Exchange をつかって非対角項を記述して全対角化計算を行うとセグメンテーション違反で落ちる。

• InterAll ですべて記述した場合は問題なく計算が実行される。
• Lanczos 法および TPQ 法では、pure InterAll による計算結果と InterALL + Exchange とによる計算結果は一致する（二体グリーン関数）

入力ファイルをtar でかためたものをdropbox においておきます。
https://dl.dropboxusercontent.com/u/61544927/HPhi-example-XXX-chain.tar.gz

$ HPhi -e namelist.def                                                         (git)-[master] -?
Start: Read File 'namelist.def'.
End: Read File 'namelist.def'.
Read File 'calcmod.def' for CalcMod.
Read File 'modpara.def' for ModPara.
Read File 'zlocspn.def' for LocSpin.
Read File 'zTrans.def' for Trans.
Read File 'zExchange.def' for Exchange.
Read File 'zInterAll.def' for InterAll.
Read File 'greenone.def' for OneBodyG.
Read File 'greentwo.def' for TwoBodyG.
Definition files are correct.
Read File 'zlocspn.def'.
Nsite= 8.
Read File 'zTrans.def'.
Read File 'zExchange.def'.
Read File 'zInterAll.def'.
Read File 'greenone.def'.
Read File 'greentwo.def'.
comb_sum= 70
MAX DIMENSION idim_max=70
APPROXIMATE REQUIRED MEMORY  max_mem=0.000007 GB
sdim=16 =2^4
Indices and Parameters of Definition files(*.def) are complete.
LARGE ALLOCATE FINISH !
EDTrans EDNTransfer=0
EDTrans EDNChemi=0
Start: Calc Sz.
idim=256.000000 irght=15 ilft=240 ihfbit=16
num_threads==4
Check.sdim=16, ihfbit=16
ib=0, num_up_rgh=4, jb=1
ib=1, num_up_rgh=3, jb=5
ib=2, num_up_rgh=3, jb=9
ib=3, num_up_rgh=2, jb=15
ib=4, num_up_rgh=3, jb=19
ib=5, num_up_rgh=2, jb=25
ib=6, num_up_rgh=2, jb=31
ib=7, num_up_rgh=1, jb=35
ib=8, num_up_rgh=3, jb=39
ib=9, num_up_rgh=2, jb=45
ib=10, num_up_rgh=2, jb=51
ib=11, num_up_rgh=1, jb=55
ib=12, num_up_rgh=2, jb=61
ib=13, num_up_rgh=1, jb=65
ib=14, num_up_rgh=1, jb=69
ib=15, num_up_rgh=0, jb=70
Xicnt=70
End  : Calc Sz.
Start: calc diagaonal components of Hamiltonian.
End  : calc diagaonal components of Hamiltonian.
Start Setting Hamiltonian.
zsh: segmentation fault  HPhi -e namelist.def

[yomichi]

環境はMac OS X Yosemite で、コンパイラはGCC 5.2.0 です（homebrew）

[k-yoshimi]

全対角化で非対角化成分を取得する際の処理に不備がありました。
修正の上、現在確認作業中です。
bugfixできましたら、修正版を早期リリースするようにします。

[k-yoshimi]

ver. 0.1.1で対応いたしましたので、ご確認ください。