Just taking a quick look on the new release I noticed some minor issues there, which probably caused by a generation script:
- If you explore diff for files (c/z)gbrfsx, cgbsvx, cgbsvxx, cheequb, clanhf, zheequb you could find that number of datatype descriptions in documentation were changed to incorrect one (DOUBLE PRECISION instead of COMPLEX, or COMPLEX*16 instead of just COMPLEX).
+ cunbdb has CMPLX instead of COMPLEX for X12 parameter.
+ zggevx RWORK become REAL instead of correct DOUBLE PRECISION
- Number of files got following string "/ output)", which seems meaningless. Just grep for it.
*>
*> \param[in] AB
*> \verbatim
-*> AB is DOUBLE PRECISION array, dimension (LDAB,N)
+*> AB is COMPLEX array, dimension (LDAB,N)
*> The original band matrix A, stored in rows 1 to KL+KU+1.
*> The j-th column of A is stored in the j-th column of the
*> array AB as follows:
*>
*> \param[in] AFB
*> \verbatim
-*> AFB is DOUBLE PRECISION array, dimension (LDAFB,N)
+*> AFB is COMPLEX array, dimension (LDAFB,N)
*> Details of the LU factorization of the band matrix A, as
*> computed by DGBTRF. U is stored as an upper triangular band
*> matrix with KL+KU superdiagonals in rows 1 to KL+KU+1, and
*>
*> \param[in] B
*> \verbatim
-*> B is REAL array, dimension (LDB,NRHS)
+*> B is COMPLEX array, dimension (LDB,NRHS)
*> The right hand side matrix B.
*> \endverbatim
*>
*>
*> \param[in,out] X
*> \verbatim
-*> X is REAL array, dimension (LDX,NRHS)
+*> X is COMPLEX array, dimension (LDX,NRHS)
*> On entry, the solution matrix X, as computed by SGETRS.
*> On exit, the improved solution matrix X.
*> \endverbatim
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] AB
*> \verbatim
-*> AB is REAL array, dimension (LDAB,N)
+*> AB is COMPLEX array, dimension (LDAB,N)
*> On entry, the matrix A in band storage, in rows 1 to KL+KU+1.
*> The j-th column of A is stored in the j-th column of the
*> array AB as follows:
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is REAL array, dimension (LDAFB,N)
+*> AFB is COMPLEX array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains details of the LU factorization of the band matrix
*> A, as computed by CGBTRF. U is stored as an upper triangular
*>
*> \param[in,out] B
*> \verbatim
-*> B is REAL array, dimension (LDB,NRHS)
+*> B is COMPLEX array, dimension (LDB,NRHS)
*> On entry, the N-by-NRHS right hand side matrix B.
*> On exit,
*> if EQUED = 'N', B is not modified;
*>
*> \param[out] X
*> \verbatim
-*> X is REAL array, dimension (LDX,NRHS)
+*> X is COMPLEX array, dimension (LDX,NRHS)
*> If INFO = 0, the N-by-NRHS solution matrix X to the original
*> system of equations. Note that A and B are modified on exit
*> if EQUED .ne. 'N', and the solution to the equilibrated system is
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[out] WORK
*> \verbatim
-*> WORK is DOUBLE PRECISION array, dimension (3*N)
+*> WORK is COMPLEX array, dimension (3*N)
*> \endverbatim
*>
*> \param[out] INFO
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in] A
*> \verbatim
-*> A is COMPLEX*16 array, dimension ( N*(N+1)/2 );
+*> A is COMPLEX array, dimension ( N*(N+1)/2 );
*> On entry, the matrix A in RFP Format.
*> RFP Format is described by TRANSR, UPLO and N as follows:
*> If TRANSR='N' then RFP A is (0:N,0:K-1) when N is even;
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] X12
*> \verbatim
-*> X12 is CMPLX array, dimension (LDX12,M-Q)
+*> X12 is COMPLEX array, dimension (LDX12,M-Q)
*> On entry, the top-right block of the unitary matrix to
*> be reduced. On exit, the form depends on TRANS:
*> If TRANS = 'N', then
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension (NPARAMS)
+*> PARAMS is DOUBLE PRECISION array, dimension (NPARAMS)
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension (NPARAMS)
+*> PARAMS is DOUBLE PRECISION array, dimension (NPARAMS)
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension (NPARAMS)
+*> PARAMS is DOUBLE PRECISION array, dimension (NPARAMS)
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension (NPARAMS)
+*> PARAMS is DOUBLE PRECISION array, dimension (NPARAMS)
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension (NPARAMS)
+*> PARAMS is DOUBLE PRECISION array, dimension (NPARAMS)
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension NPARAMS
+*> PARAMS is DOUBLE PRECISION array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension (NPARAMS)
+*> PARAMS is DOUBLE PRECISION array, dimension (NPARAMS)
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension (NPARAMS)
+*> PARAMS is DOUBLE PRECISION array, dimension (NPARAMS)
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) REAL array, dimension NPARAMS
+*> PARAMS is REAL array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in] AB
*> \verbatim
-*> AB is DOUBLE PRECISION array, dimension (LDAB,N)
+*> AB is COMPLEX*16 array, dimension (LDAB,N)
*> The original band matrix A, stored in rows 1 to KL+KU+1.
*> The j-th column of A is stored in the j-th column of the
*> array AB as follows:
*>
*> \param[in] AFB
*> \verbatim
-*> AFB is DOUBLE PRECISION array, dimension (LDAFB,N)
+*> AFB is COMPLEX*16 array, dimension (LDAFB,N)
*> Details of the LU factorization of the band matrix A, as
*> computed by DGBTRF. U is stored as an upper triangular band
*> matrix with KL+KU superdiagonals in rows 1 to KL+KU+1, and
*>
*> \param[in,out] X
*> \verbatim
-*> X is DOUBLE PRECISION array, dimension (LDX,NRHS)
+*> X is COMPLEX*16 array, dimension (LDX,NRHS)
*> On entry, the solution matrix X, as computed by DGETRS.
*> On exit, the improved solution matrix X.
*> \endverbatim
*>
*> \param[out] BERR
*> \verbatim
-*> BERR is DOUBLE PRECISION array, dimension (NRHS)
+*> BERR is COMPLEX*16 array, dimension (NRHS)
*> Componentwise relative backward error. This is the
*> componentwise relative backward error of each solution vector X(j)
*> (i.e., the smallest relative change in any element of A or B that
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension NPARAMS
+*> PARAMS is DOUBLE PRECISION array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension NPARAMS
+*> PARAMS is DOUBLE PRECISION array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension NPARAMS
+*> PARAMS is DOUBLE PRECISION array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension NPARAMS
+*> PARAMS is DOUBLE PRECISION array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[out] RWORK
*> \verbatim
-*> RWORK is REAL array, dimension (lrwork)
+*> RWORK is DOUBLE PRECISION array, dimension (lrwork)
*> lrwork must be at least max(1,6*N) if BALANC = 'S' or 'B',
*> and at least max(1,2*N) otherwise.
*> Real workspace.
*>
*> \param[out] WORK
*> \verbatim
-*> WORK is DOUBLE PRECISION array, dimension (3*N)
+*> WORK is COMPLEX*16 array, dimension (3*N)
*> \endverbatim
*>
*> \param[out] INFO
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension NPARAMS
+*> PARAMS is DOUBLE PRECISION array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension NPARAMS
+*> PARAMS is DOUBLE PRECISION array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension NPARAMS
+*> PARAMS is DOUBLE PRECISION array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension NPARAMS
+*> PARAMS is DOUBLE PRECISION array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension NPARAMS
+*> PARAMS is DOUBLE PRECISION array, dimension NPARAMS
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults