Seg fault using UPF pseudopotentials
Posted: Tue Feb 05, 2013 8:11 pm
Hello,
I've recently attempted to reproduce some Quantum Espresso calculations using Abinit, but have had little success performing calculations using the UPF format of pseudopotential. Versions 7.0.5 and 6.12.3 both segfault inside the SCF loop if a UPF pseudopotential is used. This is easily reproducible (for me at least) by taking the bulk Silicon example from the tutorials and switching the 14si.pspnc file for the 14-Si.nlcc.UPF file, also found in the bundled tutorial material.
The culprit seems to be the section of vtowfk.F90 at line 864 which is included via a compiler #if/#else directive. When using UPF pseudopotentials a segfault is encountered when calling cg_zdotc. This can be avoided by executing the presently excluded branch.
I've not delved into the mechanics of this fix, so I can't suggest how the current code is failing or even if this is a sensible thing to do. The results obtained after this procedure seem correct.
Yours,
Tom
I've recently attempted to reproduce some Quantum Espresso calculations using Abinit, but have had little success performing calculations using the UPF format of pseudopotential. Versions 7.0.5 and 6.12.3 both segfault inside the SCF loop if a UPF pseudopotential is used. This is easily reproducible (for me at least) by taking the bulk Silicon example from the tutorials and switching the 14si.pspnc file for the 14-Si.nlcc.UPF file, also found in the bundled tutorial material.
The culprit seems to be the section of vtowfk.F90 at line 864 which is included via a compiler #if/#else directive. When using UPF pseudopotentials a segfault is encountered when calling cg_zdotc. This can be avoided by executing the presently excluded branch.
I've not delved into the mechanics of this fix, so I can't suggest how the current code is failing or even if this is a sensible thing to do. The results obtained after this procedure seem correct.
Yours,
Tom