xYYzzz

YY=

BD 	bidiagonal
DI 	diagonal
GB 	general band
GE 	general (i.e., unsymmetric, in some cases rectangular)
GG 	general matrices, generalized problem (i.e., a pair of general matrices)
GT 	general tridiagonal
HB 	(complex) Hermitian band
HE 	(complex) Hermitian
HG 	upper Hessenberg matrix, generalized problem (i.e a Hessenberg and a triangular matrix)
HP 	(complex) Hermitian, packed storage
HS 	upper Hessenberg
OP 	(real) orthogonal, packed storage
OR 	(real) orthogonal
PB 	symmetric or Hermitian positive definite band
PO 	symmetric or Hermitian positive definite
PP 	symmetric or Hermitian positive definite, packed storage
PT 	symmetric or Hermitian positive definite tridiagonal
SB 	(real) symmetric band
SP 	symmetric, packed storage
ST 	(real) symmetric tridiagonal
SY 	symmetric
TB 	triangular band
TG 	triangular matrices, generalized problem (i.e., a pair of triangular matrices)
TP 	triangular, packed storage
TR 	triangular (or in some cases quasi-triangular)
TZ 	trapezoidal
UN 	(complex) unitary
UP 	(complex) unitary, packed storage

_GEMM (             TRANSA, TRANSB,      M, N, K, ALPHA, A, LDA, B, LDB, BETA, C, LDC ) S, D, C, Z
_SYMM ( SIDE, UPLO,                      M, N,    ALPHA, A, LDA, B, LDB, BETA, C, LDC ) S, D, C, Z
_HEMM ( SIDE, UPLO,                      M, N,    ALPHA, A, LDA, B, LDB, BETA, C, LDC ) C, Z
_SYRK (       UPLO, TRANS,                  N, K, ALPHA, A, LDA,         BETA, C, LDC ) S, D, C, Z
_HERK (       UPLO, TRANS,                  N, K, ALPHA, A, LDA,         BETA, C, LDC ) C, Z
_SYR2K(       UPLO, TRANS,                  N, K, ALPHA, A, LDA, B, LDB, BETA, C, LDC ) S, D, C, Z
_HER2K(       UPLO, TRANS,                  N, K, ALPHA, A, LDA, B, LDB, BETA, C, LDC ) C, Z
_TRMM ( SIDE, UPLO, TRANSA,        DIAG, M, N,    ALPHA, A, LDA, B, LDB ) S, D, C, Z
_TRSM ( SIDE, UPLO, TRANSA,        DIAG, M, N,    ALPHA, A, LDA, B, LDB ) S, D, C, Z

_GEMM :  C = a op(A) op(B) + b C ; op(X) = X, X^T, X^H ; C~(m*n)
_SYMM :  C = a A B + b C
	 C = a B A + b C ; C~(m*n) ; A = A^T
_HEMM :  C = a A B + b C
	 C = a B A + b C ; C~(m*n) ; A = A^H
_SYRK :  C = a A A^T + b C 
	 C = a A^T A + b C ; C~(n*n)
_HERK :  C = a A A^H + b C 
	 C = a A^H A + b C ; C~(n*n)
_SYR2K : C = a A B^T + a B A^T + b C 
	 C = a A^T B + a B^T A + b C ; C~(n*n)
_HER2K : C = a A B^H + a* B A^H + b C 
	 C = a A^H B + a* B^H A + b C ; C~(n*n)
_TRMM :  B = a op(A) B
	 B = a B op(A) ; op(A) = A, A^T, A^H ; B~(m*n)
_TRSM :  B = a op(A^-1) B
	 B = a B op(A^-1) ; op(A) = A, A^T, A^H ; B~(m*n)

TRANS	= 'No transpose', 'Transpose', 'Conjugate transpose' ( X, X T, XC)
UPLO 	= 'Upper triangular', 'Lower triangular'
DIAG 	= 'Non-unit triangular', 'Unit triangular'
SIDE 	= 'Left', 'Right' (A or op(A) on the left, or A or op(A) on the right)

Symmetric Eigenproblems: Ax=lx, A=A^T -> A=VLV^T
	simple driver -EV computes all the eigenvalues and vectors
	expert driver -EVX computes subset the eigenvalues and optionally  vectors
	divide-and-conquer driver -EVD is faster than simple driver for large matrices
	relatively robust representation driver -EVR is the fastest


