Hi, there:<br><br>I was stucked again in the subroutine regterg in the file regterg.f90 when I read the pwscf code. There are two things bothering me. <br><br>One is in the subroutine, the variables psi and hpsi are declared as complex(dp), however, they are used as parameters to call the  Lapack subroutine DGEMM, which requires real numbers. For example, in the line 136, <br>
<br> CALL DGEMM( 'T', 'N', nbase, nbase, ndim2, 2.D0 , &<br>               psi, ndmx2, hpsi, ndmx2, 0.D0, hr, nvecx )<br><br>The other one is about the diagonalization of the reduced hamiltonian hr, which is declared as a real symmetric matrix in the subroutine. But in general, the wavefunction psi is complex, how can we guarentee all the elements in the matrix is real? Can the matrix hr be a Hermitian one?<br>
<br><br>Best,<br><br>Wei<br>