* LDA (input) INTEGER
* The leading dimension of the array A. LDA >= max(1,N).
*
-* WORK REAL array, dimension (2*N)
+* WORK (workspace) REAL array, dimension (2*N)
* Work space.
*
* INFO (output) INTEGER
* will be used. The algorithm terminates at the (K-1)st step
* if the pivot <= TOL.
*
-* WORK REAL array, dimension (2*N)
+* WORK (workspace) REAL array, dimension (2*N)
* Work space.
*
* INFO (output) INTEGER
* even and is N when is odd.
* See the Note below for more details. Unchanged on exit.
*
-* B (input/ouptut) COMPLEX array, dimension (LDB,N)
+* B (input/output) COMPLEX array, dimension (LDB,N)
* Before entry, the leading m by n part of the array B must
* contain the right-hand side matrix B, and on exit is
* overwritten by the solution matrix X.
* The number of right hand sides, i.e., the number of columns
* of the matrix B. NRHS >= 0.
*
-* A (input or input/ouptut) DOUBLE PRECISION array,
+* A (input/output) DOUBLE PRECISION array,
* dimension (LDA,N)
* On entry, the N-by-N coefficient matrix A.
* On exit, if iterative refinement has been successfully used
* The number of right hand sides, i.e., the number of columns
* of the matrix B. NRHS >= 0.
*
-* A (input or input/ouptut) DOUBLE PRECISION array,
+* A (input/output) DOUBLE PRECISION array,
* dimension (LDA,N)
* On entry, the symmetric matrix A. If UPLO = 'U', the leading
* N-by-N upper triangular part of A contains the upper
* LDA (input) INTEGER
* The leading dimension of the array A. LDA >= max(1,N).
*
-* WORK REAL array, dimension (2*N)
+* WORK (workspace) REAL array, dimension (2*N)
* Work space.
*
* INFO (output) INTEGER
* will be used. The algorithm terminates at the (K-1)st step
* if the pivot <= TOL.
*
-* WORK REAL array, dimension (2*N)
+* WORK (workspace) REAL array, dimension (2*N)
* Work space.
*
* INFO (output) INTEGER
* even and is N when is odd.
* See the Note below for more details. Unchanged on exit.
*
-* B (input/ouptut) REAL array, DIMENSION (LDB,N)
+* B (input/output) REAL array, DIMENSION (LDB,N)
* Before entry, the leading m by n part of the array B must
* contain the right-hand side matrix B, and on exit is
* overwritten by the solution matrix X.
* The number of right hand sides, i.e., the number of columns
* of the matrix B. NRHS >= 0.
*
-* A (input or input/ouptut) COMPLEX*16 array,
+* A (input/output) COMPLEX*16 array,
* dimension (LDA,N)
* On entry, the N-by-N coefficient matrix A.
* On exit, if iterative refinement has been successfully used
* The number of right hand sides, i.e., the number of columns
* of the matrix B. NRHS >= 0.
*
-* A (input or input/ouptut) COMPLEX*16 array,
+* A (input/output) COMPLEX*16 array,
* dimension (LDA,N)
* On entry, the Hermitian matrix A. If UPLO = 'U', the leading
* N-by-N upper triangular part of A contains the upper
* LDA (input) INTEGER
* The leading dimension of the array A. LDA >= max(1,N).
*
-* WORK DOUBLE PRECISION array, dimension (2*N)
+* WORK (workspace) DOUBLE PRECISION array, dimension (2*N)
* Work space.
*
* INFO (output) INTEGER
* will be used. The algorithm terminates at the (K-1)st step
* if the pivot <= TOL.
*
-* WORK DOUBLE PRECISION array, dimension (2*N)
+* WORK (workspace) DOUBLE PRECISION array, dimension (2*N)
* Work space.
*
* INFO (output) INTEGER