
! Copyright (C) 2014 J. K. Dewhurst, S. Sharma and E. K. U. Gross.
! This file is distributed under the terms of the GNU General Public License.
! See the file COPYING for license details.

subroutine dynevs(ev,dq,wq)
use modmain
use modphonon
implicit none
! arguments
complex(8), intent(in) :: ev(nbph,nbph)
complex(8), intent(inout) :: dq(nbph,nbph)
real(8), intent(out) :: wq(nbph)
! local variables
integer i,j,k
real(8) t1,t2
complex(8) z1
! automatic arrays
real(8) wt(nbph)
! find the eigenvalues and eigenvectors of the matrix dq
call zheevdi(nbph,nbph,dq,wq)
! reorder eigenvalues so that the eigenvectors maximally overlap with ev
wt(:)=wq(:)
do i=1,nbph
  j=1
  t1=0.d0
  do k=1,nbph
    z1=dot_product(ev(:,i),dq(:,k))
    t2=z1%re**2+z1%im**2
    if (t2 > t1) then
      j=k
      t1=t2
    end if
  end do
  wq(i)=wt(j)
end do
end subroutine

