Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.
7036 Discussions

Are eigenvectors computed by zhpevd ( ) orthogonal ?

Customer__Intel4
Beginner
310 Views
Are eigenvectors computed by the functionLAPACKE_zhpevd ( ) for a hermitian matrix guaranteed to be orthogonal ? If not, is there a way to orthogonalize them ?
0 Kudos
1 Reply
Victor_K_Intel1
Employee
310 Views
Yes, the vectors must be orthogonal. Here it is an excerpt from the manual:

If jobz = 'V', then this array is overwritten by the unitary matrix Z which contains the eigenvectors of A. If

jobz = 'N', then z is not referenced.
So, if in your case you don't observe this feature this may be a bug. But you also should keep in mind that Divide-and-Conquer method is not the best algorithm for matrices with multiple eivenvalues.

Thanks
Victor

0 Kudos
Reply