Possible bug in h_psi.f90
In the following line https://gitlab.com/QEF/q-e/blob/develop/PW/src/h_psi.f90#L143 it is checked for
IF ( real_space .and. nkb > 0 )
so the code inside will be executed assuming that we are in real espace
but then there is an ELSE
https://gitlab.com/QEF/q-e/blob/develop/PW/src/h_psi.f90#L166
and the calculation there is assuming that we are in reciprocal space, at least from what it
is inferred in the comment.
But could be the case that real_space
is TRUE but nkb = 0
, in which case the ELSE
will be executed but we will be in real space.
Is this scenario possible?