I’m trying to wrap gdma calculations into a python script, but I’m stuck with a systematic error/crash: when I launch the gdma calculation from within a python interpreter or from a python script through the command line (with this command: python gdma.py), it results in a systematic crash at the gdma step.
I found some temporary files after the crashes, notably a file fort.52 which contains the following, unique line: “Could not open psi4 output for writing.unit = 51”.
So I wonder what is the default psi4 output directory where psi4.gdma is trying to write into and how I could change it.
To be noted, the same calculation using the psi4 executable (psi4 -i gdma.dat) succeded (I don’t know if it produced any temporary files, though) with no error…
I am using psi4 1.2.1 with python 3.6. The dev version didn’t solve this issue.
No other message error printed at all.
All I got was a .dma file (looks like a regular gdma input file to me) and the aforementioned fort.52 file. The python script produced the same outputs as the psi4 binary for the energy step executed right before the gdma step in my script.
What really suprised me was that the code worked using the psi4 executable but not the python wrapping. To me, the GDMA binary used in background is the same in both cases, or am I wrong?
The binary is the same, but something about how information gets passed from Psi to GDMA must depend on the input method. It shouldn’t, so this is definitely a bug.
For future reference, I can reproduce the bug with the following input file: