For symmetric eigensolvers only upper or lower triangle of a matrix should be set

and  remaining  part of matrix can contain any values including NaNs but  these LAPACKE
symmetric eigensolver drivers checked the whole matrix for NaNs and they used
transposition  of the whole matrix. The new versions check and transpose upper or lower
triangle of matrix.

Author: sergey-v-kuznetsov

LAPACKE/src/lapacke_cheev_work.c

@@ -70,15 +70,15 @@ lapack_int LAPACKE_cheev_work( int matrix_layout, char jobz, char uplo,
             goto exit_level_0;
         }
         /* Transpose input matrices */
-        LAPACKE_cge_trans( matrix_layout, n, n, a, lda, a_t, lda_t );
+        LAPACKE_che_trans( matrix_layout, uplo, n, a, lda, a_t, lda_t );
         /* Call LAPACK function and adjust info */
         LAPACK_cheev( &jobz, &uplo, &n, a_t, &lda_t, w, work, &lwork, rwork,
                       &info );
         if( info < 0 ) {
             info = info - 1;
         }
         /* Transpose output matrices */
-        LAPACKE_cge_trans( LAPACK_COL_MAJOR, n, n, a_t, lda_t, a, lda );
+        LAPACKE_che_trans( LAPACK_COL_MAJOR, uplo, n, a_t, lda_t, a, lda );
         /* Release memory and exit */
         LAPACKE_free( a_t );
 exit_level_0: