.PRIVECのデータがないパターンの検知

kohei-noda-qcrg commented 8 months ago

.PRIVECを指定してもなぜかvector printのデータがないときがある
- .ANALYZEがないとき？
- 条件を調査して判明したら、DIRACのインプットを読めば.PRIVECのデータがないことは分かるので、エラーにしたりWarningにしたり、無視してeigenvaluesのデータだけ読むように変更できるため、要調査

kohei-noda-qcrg commented 7 months ago

H₂OとF₂で検証

結論
- *.inpレベルでは全く同じものを試しても、PRIVECのデータが出る場合と出ない場合がある
- 原因は現時点で不明
- C2vとC32hなので、それで異なる?

インプットの中身

PRIVECするための最低限の設定

.VECPRIの設定はなし(デフォルト値)

noda@relqc01 ~/test/dirac/privec/dirac19
$ cat H2O/H2O.inp
**DIRAC
.TITLE
H2O
.WAVE FUNCTION
.ANALYZE
**HAMILTONIAN
.DOSSSS
**WAVE FUNCTIONS
.SCF
**ANALYZE
.PRIVEC
**MOLECULE
*BASIS
.DEFAULT
6-31G
*END OF

インプットの差

タイトルだけ

noda@relqc01 ~/test/dirac/privec/dirac19
$ diff F2/F2.inp H2O/H2O.inp
3c3
< F2
---
> H2O

H₂OのPRIVEC付近

デフォルト値のoccupied electronic solutionsの範囲のVector printが出力される


noda@relqc01 ~/test/dirac/privec/dirac19
$ grep -A 15 "\*\* Vector print" H2O/H2O_H2O.out
****************************** Vector print ******************************
**************************************************************************

Coefficients from DFCOEF

                          Fermion ircop E1
                          ----------------

Electronic eigenvalue no. 1: -20.577173840434

F₂のPRIVEC付近

Vector print に何も出ない


noda@relqc01 ~/test/dirac/privec/dirac19
$ grep -A 15 "\*\* Vector print" F2/F2_F2.out
****************************** Vector print ******************************
**************************************************************************      

                                                                                                                                                                                                                                                                                        Coefficients from DFCOEF
------------------------

** E N D of D I R A C output **

Date and time (Linux) : Fri Mar 1 23:05:01 2024

kohei-noda-qcrg commented 7 months ago

現時点のDIRAC最新コミット https://gitlab.com/dirac/dirac/-/tree/76204c28f3622253d4444d8b16d991b1edb2092b で試してもHDF5 CheckpointファイルからVector printを出力しようとする動作以外は同じ傾向
等核二原子分子と原子だと出ない?
- COとArを試した
- COはVector printの中身が出力される
- Arは出力されない

kohei-noda-qcrg commented 7 months ago

DIRAC側のバグの可能性が出てきたので、DIRACのソースコードを調査予定

kohei-noda-qcrg commented 7 months ago

仮にDIRAC側のバグだとしても、インプットによってはエラーが発生する場合があるので、その対処は必要

例: .MULPOPを付けていた場合、MULPOPのeigenvalueの表示がPRIVECのそれに酷似しているが、微妙に異なるのでエラーが発生してしまう

kohei-noda-qcrg commented 7 months ago

CALL TITLER('Vector print','*',130) のようなsubroutineがDIRACにはあり、DIRACの各機能が始まるときにこのサブルーチンを呼んでいるようなので

PRIVECを読もうとしているのに ****** **** Vector print * ******

のようなセクション、つまり.PRIVECのデータを読むコードが始まった後に*しかない行が見つかったら .PRIVECのデータが空であることが分かる

kohei-noda-qcrg commented 7 months ago

https://github.com/RQC-HU/sum_dirac_dfcoef/issues/91#issuecomment-1973304319

について、実際にDIRACのバグであることが判明したので

issueを https://gitlab.com/dirac/dirac/-/issues/92 に、merge requestを https://gitlab.com/dirac/dirac/-/merge_requests/165 に作成した。マージされればDIRAC25でこのDIRACのバグは解消されるはず。

RQC-HU / sum_dirac_dfcoef