*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is COMPLEX array dimension (min(M,N))
+*> TAUQ is COMPLEX array, dimension (min(M,N))
*> The scalar factors of the elementary reflectors which
*> represent the unitary matrix Q. See Further Details.
*> \endverbatim
*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is COMPLEX array dimension (min(M,N))
+*> TAUQ is COMPLEX array, dimension (min(M,N))
*> The scalar factors of the elementary reflectors which
*> represent the unitary matrix Q. See Further Details.
*> \endverbatim
*>
*> \param[in,out] CWORK
*> \verbatim
-*> CWORK is COMPLEX array, dimension max(1,LWORK).
+*> CWORK is COMPLEX array, dimension (max(1,LWORK))
*> Used as workspace.
*> If on entry LWORK .EQ. -1, then a workspace query is assumed and
*> no computation is done; CWORK(1) is set to the minial (and optimal)
*>
*> \param[in,out] RWORK
*> \verbatim
-*> RWORK is REAL array, dimension max(6,LRWORK).
+*> RWORK is REAL array, dimension (max(6,LRWORK))
*> On entry,
*> If JOBU .EQ. 'C' :
*> RWORK(1) = CTOL, where CTOL defines the threshold for convergence.
*>
*> \param[in] A
*> \verbatim
-*> A is COMPLEX array of DIMENSION ( LDA, n ).
+*> A is COMPLEX array, dimension ( LDA, n ).
*> Before entry, the leading m by n part of the array A must
*> contain the matrix of coefficients.
*> Unchanged on exit.
*>
*> \param[in] A
*> \verbatim
-*> A is COMPLEX array of DIMENSION ( LDA, n ).
+*> A is COMPLEX array, dimension ( LDA, n ).
*> Before entry, the leading m by n part of the array A must
*> contain the matrix of coefficients.
*> Unchanged on exit.
*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is COMPLEX array dimension (NB)
+*> TAUQ is COMPLEX array, dimension (NB)
*> The scalar factors of the elementary reflectors which
*> represent the unitary matrix Q. See Further Details.
*> \endverbatim
*>
*> \param[in] H
*> \verbatim
-*> H is COMPLEX array of dimension (LDH,N)
+*> H is COMPLEX array, dimension (LDH,N)
*> The 2-by-2 or 3-by-3 matrix H in (*).
*> \endverbatim
*>
*>
*> \param[out] V
*> \verbatim
-*> V is COMPLEX array of dimension N
+*> V is COMPLEX array, dimension (N)
*> A scalar multiple of the first column of the
*> matrix K in (*).
*> \endverbatim
*>
*> \param[in,out] S
*> \verbatim
-*> S is COMPLEX array of size (NSHFTS)
+*> S is COMPLEX array, dimension (NSHFTS)
*> S contains the shifts of origin that define the multi-
*> shift QR sweep. On output S may be reordered.
*> \endverbatim
*>
*> \param[in,out] H
*> \verbatim
-*> H is COMPLEX array of size (LDH,N)
+*> H is COMPLEX array, dimension (LDH,N)
*> On input H contains a Hessenberg matrix. On output a
*> multi-shift QR sweep with shifts SR(J)+i*SI(J) is applied
*> to the isolated diagonal block in rows and columns KTOP
*>
*> \param[in,out] Z
*> \verbatim
-*> Z is COMPLEX array of size (LDZ,IHIZ)
+*> Z is COMPLEX array, dimension (LDZ,IHIZ)
*> If WANTZ = .TRUE., then the QR Sweep unitary
*> similarity transformation is accumulated into
*> Z(ILOZ:IHIZ,ILOZ:IHIZ) from the right.
*>
*> \param[out] V
*> \verbatim
-*> V is COMPLEX array of size (LDV,NSHFTS/2)
+*> V is COMPLEX array, dimension (LDV,NSHFTS/2)
*> \endverbatim
*>
*> \param[in] LDV
*>
*> \param[out] U
*> \verbatim
-*> U is COMPLEX array of size
-*> (LDU,3*NSHFTS-3)
+*> U is COMPLEX array, dimension (LDU,3*NSHFTS-3)
*> \endverbatim
*>
*> \param[in] LDU
*>
*> \param[out] WH
*> \verbatim
-*> WH is COMPLEX array of size (LDWH,NH)
+*> WH is COMPLEX array, dimension (LDWH,NH)
*> \endverbatim
*>
*> \param[in] LDWH
*>
*> \param[out] WV
*> \verbatim
-*> WV is COMPLEX array of size
-*> (LDWV,3*NSHFTS-3)
+*> WV is COMPLEX array, dimension (LDWV,3*NSHFTS-3)
*> \endverbatim
*>
*> \param[in] LDWV
*>
*> \param[out] ISUPPZ
*> \verbatim
-*> ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*> ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
*> The support of the eigenvectors in Z, i.e., the indices
*> indicating the nonzero elements in Z. The i-th computed eigenvector
*> is nonzero only in elements ISUPPZ( 2*i-1 ) through
*>
*> \param[out] ISUPPZ
*> \verbatim
-*> ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*> ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
*> The support of the eigenvectors in Z, i.e., the indices
*> indicating the nonzero elements in Z. The i-th computed eigenvector
*> is nonzero only in elements ISUPPZ( 2*i-1 ) through
*>
*> \param[in,out] A
*> \verbatim
-*> A is COMPLEX arrays, dimensions (LDA,N)
+*> A is COMPLEX array, dimension (LDA,N)
*> On entry, the matrix A in the pair (A, B).
*> On exit, the updated matrix A.
*> \endverbatim
*>
*> \param[in,out] B
*> \verbatim
-*> B is COMPLEX arrays, dimensions (LDB,N)
+*> B is COMPLEX array, dimension (LDB,N)
*> On entry, the matrix B in the pair (A, B).
*> On exit, the updated matrix B.
*> \endverbatim
*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is DOUBLE PRECISION array dimension (min(M,N))
+*> TAUQ is DOUBLE PRECISION array, dimension (min(M,N))
*> The scalar factors of the elementary reflectors which
*> represent the orthogonal matrix Q. See Further Details.
*> \endverbatim
*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is DOUBLE PRECISION array dimension (min(M,N))
+*> TAUQ is DOUBLE PRECISION array, dimension (min(M,N))
*> The scalar factors of the elementary reflectors which
*> represent the orthogonal matrix Q. See Further Details.
*> \endverbatim
*>
*> \param[in,out] WORK
*> \verbatim
-*> WORK is DOUBLE PRECISION array, dimension MAX(6,M+N).
+*> WORK is DOUBLE PRECISION array, dimension (max(6,M+N))
*> On entry :
*> If JOBU .EQ. 'C' :
*> WORK(1) = CTOL, where CTOL defines the threshold for convergence.
*>
*> \param[in] AB
*> \verbatim
-*> AB is DOUBLE PRECISION array of DIMENSION ( LDAB, n )
+*> AB is DOUBLE PRECISION array, dimension ( LDAB, n )
*> Before entry, the leading m by n part of the array AB must
*> contain the matrix of coefficients.
*> Unchanged on exit.
*>
*> \param[in] A
*> \verbatim
-*> A is DOUBLE PRECISION array of DIMENSION ( LDA, n )
+*> A is DOUBLE PRECISION array, dimension ( LDA, n )
*> Before entry, the leading m by n part of the array A must
*> contain the matrix of coefficients.
*> Unchanged on exit.
*>
*> \param[in] A
*> \verbatim
-*> A is DOUBLE PRECISION array of DIMENSION ( LDA, n ).
+*> A is DOUBLE PRECISION array, dimension ( LDA, n ).
*> Before entry, the leading m by n part of the array A must
*> contain the matrix of coefficients.
*> Unchanged on exit.
*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is DOUBLE PRECISION array dimension (NB)
+*> TAUQ is DOUBLE PRECISION array, dimension (NB)
*> The scalar factors of the elementary reflectors which
*> represent the orthogonal matrix Q. See Further Details.
*> \endverbatim
*>
*> \param[in] H
*> \verbatim
-*> H is DOUBLE PRECISION array of dimension (LDH,N)
+*> H is DOUBLE PRECISION array, dimension (LDH,N)
*> The 2-by-2 or 3-by-3 matrix H in (*).
*> \endverbatim
*>
*>
*> \param[out] V
*> \verbatim
-*> V is DOUBLE PRECISION array of dimension N
+*> V is DOUBLE PRECISION array, dimension (N)
*> A scalar multiple of the first column of the
*> matrix K in (*).
*> \endverbatim
*>
*> \param[in,out] SR
*> \verbatim
-*> SR is DOUBLE PRECISION array of size (NSHFTS)
+*> SR is DOUBLE PRECISION array, dimension (NSHFTS)
*> \endverbatim
*>
*> \param[in,out] SI
*> \verbatim
-*> SI is DOUBLE PRECISION array of size (NSHFTS)
+*> SI is DOUBLE PRECISION array, dimension (NSHFTS)
*> SR contains the real parts and SI contains the imaginary
*> parts of the NSHFTS shifts of origin that define the
*> multi-shift QR sweep. On output SR and SI may be
*>
*> \param[in,out] H
*> \verbatim
-*> H is DOUBLE PRECISION array of size (LDH,N)
+*> H is DOUBLE PRECISION array, dimension (LDH,N)
*> On input H contains a Hessenberg matrix. On output a
*> multi-shift QR sweep with shifts SR(J)+i*SI(J) is applied
*> to the isolated diagonal block in rows and columns KTOP
*>
*> \param[in,out] Z
*> \verbatim
-*> Z is DOUBLE PRECISION array of size (LDZ,IHIZ)
+*> Z is DOUBLE PRECISION array, dimension (LDZ,IHIZ)
*> If WANTZ = .TRUE., then the QR Sweep orthogonal
*> similarity transformation is accumulated into
*> Z(ILOZ:IHIZ,ILOZ:IHIZ) from the right.
*>
*> \param[out] V
*> \verbatim
-*> V is DOUBLE PRECISION array of size (LDV,NSHFTS/2)
+*> V is DOUBLE PRECISION array, dimension (LDV,NSHFTS/2)
*> \endverbatim
*>
*> \param[in] LDV
*>
*> \param[out] U
*> \verbatim
-*> U is DOUBLE PRECISION array of size
-*> (LDU,3*NSHFTS-3)
+*> U is DOUBLE PRECISION array, dimension (LDU,3*NSHFTS-3)
*> \endverbatim
*>
*> \param[in] LDU
*>
*> \param[out] WH
*> \verbatim
-*> WH is DOUBLE PRECISION array of size (LDWH,NH)
+*> WH is DOUBLE PRECISION array, dimension (LDWH,NH)
*> \endverbatim
*>
*> \param[in] LDWH
*>
*> \param[out] WV
*> \verbatim
-*> WV is DOUBLE PRECISION array of size
-*> (LDWV,3*NSHFTS-3)
+*> WV is DOUBLE PRECISION array, dimension (LDWV,3*NSHFTS-3)
*> \endverbatim
*>
*> \param[in] LDWV
*>
*> \param[out] IDXP
*> \verbatim
-*> IDXP is INTEGER array dimension(N)
+*> IDXP is INTEGER array, dimension(N)
*> This will contain the permutation used to place deflated
*> values of D at the end of the array. On output IDXP(2:K)
*> points to the nondeflated D-values and IDXP(K+1:N)
*>
*> \param[out] IDX
*> \verbatim
-*> IDX is INTEGER array dimension(N)
+*> IDX is INTEGER array, dimension(N)
*> This will contain the permutation used to sort the contents of
*> D into ascending order.
*> \endverbatim
*>
*> \param[out] IDXC
*> \verbatim
-*> IDXC is INTEGER array dimension(N)
+*> IDXC is INTEGER array, dimension(N)
*> This will contain the permutation used to arrange the columns
*> of the deflated U matrix into three groups: the first group
*> contains non-zero entries only at and above NL, the second
*>
*> \param[in,out] IDXQ
*> \verbatim
-*> IDXQ is INTEGER array dimension(N)
+*> IDXQ is INTEGER array, dimension(N)
*> This contains the permutation which separately sorts the two
*> sub-problems in D into ascending order. Note that entries in
*> the first hlaf of this permutation must first be moved one
*>
*> \param[out] COLTYP
*> \verbatim
-*> COLTYP is INTEGER array dimension(N)
+*> COLTYP is INTEGER array, dimension(N)
*> As workspace, this will contain a label which will indicate
*> which of the following types a column in the U2 matrix or a
*> row in the VT2 matrix is:
*>
*> \param[out] ISUPPZ
*> \verbatim
-*> ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*> ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
*> The support of the eigenvectors in Z, i.e., the indices
*> indicating the nonzero elements in Z. The i-th computed eigenvector
*> is nonzero only in elements ISUPPZ( 2*i-1 ) through
*>
*> \param[out] ISUPPZ
*> \verbatim
-*> ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*> ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
*> The support of the eigenvectors in Z, i.e., the indices
*> indicating the nonzero elements in Z. The i-th computed eigenvector
*> is nonzero only in elements ISUPPZ( 2*i-1 ) through
*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is REAL array dimension (min(M,N))
+*> TAUQ is REAL array, dimension (min(M,N))
*> The scalar factors of the elementary reflectors which
*> represent the orthogonal matrix Q. See Further Details.
*> \endverbatim
*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is REAL array dimension (min(M,N))
+*> TAUQ is REAL array, dimension (min(M,N))
*> The scalar factors of the elementary reflectors which
*> represent the orthogonal matrix Q. See Further Details.
*> \endverbatim
*>
*> \param[in,out] WORK
*> \verbatim
-*> WORK is REAL array, dimension MAX(6,M+N).
+*> WORK is REAL array, dimension (max(6,M+N))
*> On entry,
*> If JOBU .EQ. 'C' :
*> WORK(1) = CTOL, where CTOL defines the threshold for convergence.
*>
*> \param[in] AB
*> \verbatim
-*> AB is REAL array of DIMENSION ( LDAB, n )
+*> AB is REAL array, dimension ( LDAB, n )
*> Before entry, the leading m by n part of the array AB must
*> contain the matrix of coefficients.
*> Unchanged on exit.
*>
*> \param[in] A
*> \verbatim
-*> A is REAL array of DIMENSION ( LDA, n )
+*> A is REAL array, dimension ( LDA, n )
*> Before entry, the leading m by n part of the array A must
*> contain the matrix of coefficients.
*> Unchanged on exit.
*>
*> \param[in] A
*> \verbatim
-*> A is REAL array of DIMENSION ( LDA, n ).
+*> A is REAL array, dimension ( LDA, n ).
*> Before entry, the leading m by n part of the array A must
*> contain the matrix of coefficients.
*> Unchanged on exit.
*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is REAL array dimension (NB)
+*> TAUQ is REAL array, dimension (NB)
*> The scalar factors of the elementary reflectors which
*> represent the orthogonal matrix Q. See Further Details.
*> \endverbatim
*>
*> \param[in] H
*> \verbatim
-*> H is REAL array of dimension (LDH,N)
+*> H is REAL array, dimension (LDH,N)
*> The 2-by-2 or 3-by-3 matrix H in (*).
*> \endverbatim
*>
*>
*> \param[out] V
*> \verbatim
-*> V is REAL array of dimension N
+*> V is REAL array, dimension (N)
*> A scalar multiple of the first column of the
*> matrix K in (*).
*> \endverbatim
*>
*> \param[in,out] SR
*> \verbatim
-*> SR is REAL array of size (NSHFTS)
+*> SR is REAL array, dimension (NSHFTS)
*> \endverbatim
*>
*> \param[in,out] SI
*> \verbatim
-*> SI is REAL array of size (NSHFTS)
+*> SI is REAL array, dimension (NSHFTS)
*> SR contains the real parts and SI contains the imaginary
*> parts of the NSHFTS shifts of origin that define the
*> multi-shift QR sweep. On output SR and SI may be
*>
*> \param[in,out] H
*> \verbatim
-*> H is REAL array of size (LDH,N)
+*> H is REAL array, dimension (LDH,N)
*> On input H contains a Hessenberg matrix. On output a
*> multi-shift QR sweep with shifts SR(J)+i*SI(J) is applied
*> to the isolated diagonal block in rows and columns KTOP
*>
*> \param[in,out] Z
*> \verbatim
-*> Z is REAL array of size (LDZ,IHIZ)
+*> Z is REAL array, dimension (LDZ,IHIZ)
*> If WANTZ = .TRUE., then the QR Sweep orthogonal
*> similarity transformation is accumulated into
*> Z(ILOZ:IHIZ,ILOZ:IHIZ) from the right.
*>
*> \param[out] V
*> \verbatim
-*> V is REAL array of size (LDV,NSHFTS/2)
+*> V is REAL array, dimension (LDV,NSHFTS/2)
*> \endverbatim
*>
*> \param[in] LDV
*>
*> \param[out] U
*> \verbatim
-*> U is REAL array of size
-*> (LDU,3*NSHFTS-3)
+*> U is REAL array, dimension (LDU,3*NSHFTS-3)
*> \endverbatim
*>
*> \param[in] LDU
*>
*> \param[out] WH
*> \verbatim
-*> WH is REAL array of size (LDWH,NH)
+*> WH is REAL array, dimension (LDWH,NH)
*> \endverbatim
*>
*> \param[in] LDWH
*>
*> \param[out] WV
*> \verbatim
-*> WV is REAL array of size
-*> (LDWV,3*NSHFTS-3)
+*> WV is REAL array, dimension (LDWV,3*NSHFTS-3)
*> \endverbatim
*>
*> \param[in] LDWV
*>
*> \param[in] A
*> \verbatim
-*> A is REAL array of DIMENSION (LDA,ka)
+*> A is REAL array, dimension (LDA,ka)
*> where KA
*> is K when TRANS = 'N' or 'n', and is N otherwise. Before
*> entry with TRANS = 'N' or 'n', the leading N--by--K part of
*>
*> \param[out] ISUPPZ
*> \verbatim
-*> ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*> ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
*> The support of the eigenvectors in Z, i.e., the indices
*> indicating the nonzero elements in Z. The i-th computed eigenvector
*> is nonzero only in elements ISUPPZ( 2*i-1 ) through
*>
*> \param[out] ISUPPZ
*> \verbatim
-*> ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*> ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
*> The support of the eigenvectors in Z, i.e., the indices
*> indicating the nonzero elements in Z. The i-th computed eigenvector
*> is nonzero only in elements ISUPPZ( 2*i-1 ) through
*>
*> \param[in,out] B
*> \verbatim
-*> B is REAL array, DIMENSION (LDB,N)
+*> B is 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.
*>
*> \param[in,out] A
*> \verbatim
-*> A is REAL arrays, dimensions (LDA,N)
+*> A is REAL array, dimension (LDA,N)
*> On entry, the matrix A in the pair (A, B).
*> On exit, the updated matrix A.
*> \endverbatim
*>
*> \param[in,out] B
*> \verbatim
-*> B is REAL arrays, dimensions (LDB,N)
+*> B is REAL array, dimension (LDB,N)
*> On entry, the matrix B in the pair (A, B).
*> On exit, the updated matrix B.
*> \endverbatim
*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is COMPLEX*16 array dimension (min(M,N))
+*> TAUQ is COMPLEX*16 array, dimension (min(M,N))
*> The scalar factors of the elementary reflectors which
*> represent the unitary matrix Q. See Further Details.
*> \endverbatim
*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is COMPLEX*16 array dimension (min(M,N))
+*> TAUQ is COMPLEX*16 array, dimension (min(M,N))
*> The scalar factors of the elementary reflectors which
*> represent the unitary matrix Q. See Further Details.
*> \endverbatim
*>
*> \param[in,out] CWORK
*> \verbatim
-*> CWORK is COMPLEX*16 array, dimension max(1,LWORK).
+*> CWORK is COMPLEX*16 array, dimension (max(1,LWORK))
*> Used as workspace.
*> If on entry LWORK .EQ. -1, then a workspace query is assumed and
*> no computation is done; CWORK(1) is set to the minial (and optimal)
*>
*> \param[in,out] RWORK
*> \verbatim
-*> RWORK is DOUBLE PRECISION array, dimension max(6,LRWORK).
+*> RWORK is DOUBLE PRECISION array, dimension (max(6,LRWORK))
*> On entry,
*> If JOBU .EQ. 'C' :
*> RWORK(1) = CTOL, where CTOL defines the threshold for convergence.
*>
*> \param[in] A
*> \verbatim
-*> A is COMPLEX*16 array of DIMENSION (LDA,ka)
+*> A is COMPLEX*16 array, dimension (LDA,ka)
*> where KA
*> is K when TRANS = 'N' or 'n', and is N otherwise. Before
*> entry with TRANS = 'N' or 'n', the leading N--by--K part of
*>
*> \param[in] AB
*> \verbatim
-*> AB is COMPLEX*16 array of DIMENSION ( LDAB, n )
+*> AB is COMPLEX*16 array, dimension ( LDAB, n )
*> Before entry, the leading m by n part of the array AB must
*> contain the matrix of coefficients.
*> Unchanged on exit.
*>
*> \param[in] A
*> \verbatim
-*> A is COMPLEX*16 array of DIMENSION ( LDA, n )
+*> A is COMPLEX*16 array, dimension ( LDA, n )
*> Before entry, the leading m by n part of the array A must
*> contain the matrix of coefficients.
*> Unchanged on exit.
*>
*> \param[in] X
*> \verbatim
-*> X is COMPLEX*16 array of DIMENSION at least
+*> X is COMPLEX*16 array, dimension at least
*> ( 1 + ( n - 1 )*abs( INCX ) ) when TRANS = 'N' or 'n'
*> and at least
*> ( 1 + ( m - 1 )*abs( INCX ) ) otherwise.
*>
*> \param[in] A
*> \verbatim
-*> A is COMPLEX*16 array, DIMENSION ( LDA, n ).
+*> A is COMPLEX*16 array, dimension ( LDA, n ).
*> Before entry, the leading m by n part of the array A must
*> contain the matrix of coefficients.
*> Unchanged on exit.
*>
*> \param[in] X
*> \verbatim
-*> X is COMPLEX*16 array, DIMENSION at least
+*> X is COMPLEX*16 array, dimension at least
*> ( 1 + ( n - 1 )*abs( INCX ) )
*> Before entry, the incremented array X must contain the
*> vector x.
*>
*> \param[in] A
*> \verbatim
-*> A is COMPLEX*16 array, DIMENSION ( LDA, n ).
+*> A is COMPLEX*16 array, dimension ( LDA, n ).
*> Before entry, the leading m by n part of the array A must
*> contain the matrix of coefficients.
*> Unchanged on exit.
*>
*> \param[in] X
*> \verbatim
-*> X is COMPLEX*16 array, DIMENSION at least
+*> X is COMPLEX*16 array, dimension at least
*> ( 1 + ( n - 1 )*abs( INCX ) )
*> Before entry, the incremented array X must contain the
*> vector x.
*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is COMPLEX*16 array dimension (NB)
+*> TAUQ is COMPLEX*16 array, dimension (NB)
*> The scalar factors of the elementary reflectors which
*> represent the unitary matrix Q. See Further Details.
*> \endverbatim
*>
*> \param[in] H
*> \verbatim
-*> H is COMPLEX*16 array of dimension (LDH,N)
+*> H is COMPLEX*16 array, dimension (LDH,N)
*> The 2-by-2 or 3-by-3 matrix H in (*).
*> \endverbatim
*>
*>
*> \param[out] V
*> \verbatim
-*> V is COMPLEX*16 array of dimension N
+*> V is COMPLEX*16 array, dimension (N)
*> A scalar multiple of the first column of the
*> matrix K in (*).
*> \endverbatim
*>
*> \param[in,out] S
*> \verbatim
-*> S is COMPLEX*16 array of size (NSHFTS)
+*> S is COMPLEX*16 array, dimension (NSHFTS)
*> S contains the shifts of origin that define the multi-
*> shift QR sweep. On output S may be reordered.
*> \endverbatim
*>
*> \param[in,out] H
*> \verbatim
-*> H is COMPLEX*16 array of size (LDH,N)
+*> H is COMPLEX*16 array, dimension (LDH,N)
*> On input H contains a Hessenberg matrix. On output a
*> multi-shift QR sweep with shifts SR(J)+i*SI(J) is applied
*> to the isolated diagonal block in rows and columns KTOP
*>
*> \param[in,out] Z
*> \verbatim
-*> Z is COMPLEX*16 array of size (LDZ,IHIZ)
+*> Z is COMPLEX*16 array, dimension (LDZ,IHIZ)
*> If WANTZ = .TRUE., then the QR Sweep unitary
*> similarity transformation is accumulated into
*> Z(ILOZ:IHIZ,ILOZ:IHIZ) from the right.
*>
*> \param[out] V
*> \verbatim
-*> V is COMPLEX*16 array of size (LDV,NSHFTS/2)
+*> V is COMPLEX*16 array, dimension (LDV,NSHFTS/2)
*> \endverbatim
*>
*> \param[in] LDV
*>
*> \param[out] U
*> \verbatim
-*> U is COMPLEX*16 array of size
-*> (LDU,3*NSHFTS-3)
+*> U is COMPLEX*16 array, dimension (LDU,3*NSHFTS-3)
*> \endverbatim
*>
*> \param[in] LDU
*>
*> \param[out] WH
*> \verbatim
-*> WH is COMPLEX*16 array of size (LDWH,NH)
+*> WH is COMPLEX*16 array, dimension (LDWH,NH)
*> \endverbatim
*>
*> \param[in] LDWH
*>
*> \param[out] WV
*> \verbatim
-*> WV is COMPLEX*16 array of size
-*> (LDWV,3*NSHFTS-3)
+*> WV is COMPLEX*16 array, dimension (LDWV,3*NSHFTS-3)
*> \endverbatim
*>
*> \param[in] LDWV
*>
*> \param[out] ISUPPZ
*> \verbatim
-*> ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*> ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
*> The support of the eigenvectors in Z, i.e., the indices
*> indicating the nonzero elements in Z. The i-th computed eigenvector
*> is nonzero only in elements ISUPPZ( 2*i-1 ) through
*>
*> \param[out] ISUPPZ
*> \verbatim
-*> ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*> ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
*> The support of the eigenvectors in Z, i.e., the indices
*> indicating the nonzero elements in Z. The i-th computed eigenvector
*> is nonzero only in elements ISUPPZ( 2*i-1 ) through
*>
*> \param[in,out] A
*> \verbatim
-*> A is COMPLEX*16 arrays, dimensions (LDA,N)
+*> A is COMPLEX*16 array, dimensions (LDA,N)
*> On entry, the matrix A in the pair (A, B).
*> On exit, the updated matrix A.
*> \endverbatim
*>
*> \param[in,out] B
*> \verbatim
-*> B is COMPLEX*16 arrays, dimensions (LDB,N)
+*> B is COMPLEX*16 array, dimensions (LDB,N)
*> On entry, the matrix B in the pair (A, B).
*> On exit, the updated matrix B.
*> \endverbatim