Hi Prof Julian,

Could you please let me know as to :

(1) Which subroutine of LAPACK package is GULP using in order to calculate the eigenvectors and eigenvalues ?

(2) As I am not familiar enough with fortran so for the same subroutine (having the same dynamical matrix ) can I use Python(scipy) in order to calculate exactly the same set of eigenvectors as we are getting from GULP ?

I believe I need to call some specific subroutine inorder to calculate the eigenvectors for my Hermitian dynamical matrix.

(3) If not python then could you suggest me other ways by which I could get the same set of eigenvectors?

Thank you in advance.