[QE-developers] Bug in Modules/read_pseudo.f90 (Quantum ESPRESSO v.6.8)
Lothar Brendel
lothar.brendel at uni-due.de
Fri Sep 24 10:16:12 CEST 2021
Dear espressoers,
I'm pretty sure to have found an inconsistency regarding the treatment
of the case of zero wavefunctions in a PP file (nwfc = 0):
In upflib/read_upf_v1.f90:subroutine read_pseudo_pswfc(), the allocation
of the chi array is protected against the zero case:
ALLOCATE( upf%chi( upf%mesh, MAX( upf%nwfc, 1 ) ) )
But for the MPI slave ranks in Modules/read_pseudo.f90:SUBROUTINE
upf_bcast(), it is *not*:
ALLOCATE( upf%chi(upf%mesh,upf%nwfc) )
Hence, when subsequently broadcasting the array, the slave ranks are
receiving the data into a zero sized array.
Kind regards
Lothar Brendel
More information about the developers
mailing list