800 {
804
806
808 LDA = *n+1;
810
811 for( i=0; i<*n; i++ )
812 for( j=0; j<*n; j++ ) {
813 A[ LDA*i+j ].
real=a[ (*lda)*j+i ].
real;
814 A[ LDA*i+j ].
imag=a[ (*lda)*j+i ].
imag;
815 }
816
818 for( i=0; i<*n; i++ )
819 for( j=0; j<*n; j++ ) {
820 a[ (*lda)*j+i ].
real=A[ LDA*i+j ].
real;
821 a[ (*lda)*j+i ].
imag=A[ LDA*i+j ].
imag;
822 }
823 free(A);
824 }
827 else
829}
void cblas_zher(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const CBLAS_INT N, const double alpha, const void *X, const CBLAS_INT incX, void *A, const CBLAS_INT lda)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)