summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorlangou <julien.langou@ucdenver.edu>2017-06-03 19:51:37 +0000
committerGitHub <noreply@github.com>2017-06-03 19:51:37 +0000
commit0b0dcee3f8c2d2cea5aba274f56018bbdc979e72 (patch)
tree3225276d7cffc79b2983c72c4bd389514c788bce
parenta3fb05f665cae999625f82dc91eccad8388a97c5 (diff)
parentbd1af2e5f3ea836aaaa638d9356a1ca91bc1c889 (diff)
downloadlapack-0b0dcee3f8c2d2cea5aba274f56018bbdc979e72.tar.gz
lapack-0b0dcee3f8c2d2cea5aba274f56018bbdc979e72.tar.bz2
lapack-0b0dcee3f8c2d2cea5aba274f56018bbdc979e72.zip
Merge pull request #160 from thielema/cleanup-comments
more cleaned up comments
-rw-r--r--SRC/cgebd2.f2
-rw-r--r--SRC/cgebrd.f2
-rw-r--r--SRC/cgesvj.f6
-rw-r--r--SRC/cla_heamv.f2
-rw-r--r--SRC/cla_syamv.f2
-rw-r--r--SRC/clabrd.f2
-rw-r--r--SRC/claqr1.f4
-rw-r--r--SRC/claqr5.f20
-rw-r--r--SRC/clarrv.f2
-rw-r--r--SRC/cstegr.f2
-rw-r--r--SRC/cstemr.f2
-rw-r--r--SRC/ctgex2.f4
-rw-r--r--SRC/cuncsd.f8
-rw-r--r--SRC/cunm22.f4
-rw-r--r--SRC/dgebal.f4
-rw-r--r--SRC/dgebd2.f2
-rw-r--r--SRC/dgebrd.f2
-rw-r--r--SRC/dgesvj.f4
-rw-r--r--SRC/dla_gbamv.f2
-rw-r--r--SRC/dla_geamv.f2
-rw-r--r--SRC/dla_syamv.f2
-rw-r--r--SRC/dlabrd.f2
-rw-r--r--SRC/dlaqr1.f4
-rw-r--r--SRC/dlaqr5.f22
-rw-r--r--SRC/dlasd0.f6
-rw-r--r--SRC/dlasd2.f10
-rw-r--r--SRC/dlasq5.f2
-rw-r--r--SRC/dorcsd.f8
-rw-r--r--SRC/dorm22.f4
-rw-r--r--SRC/dstegr.f2
-rw-r--r--SRC/dstemr.f2
-rw-r--r--SRC/dtgsen.f2
-rw-r--r--SRC/ilaclr.f2
-rw-r--r--SRC/sgebd2.f2
-rw-r--r--SRC/sgebrd.f2
-rw-r--r--SRC/sgees.f2
-rw-r--r--SRC/sgesvj.f4
-rw-r--r--SRC/sla_gbamv.f2
-rw-r--r--SRC/sla_geamv.f2
-rw-r--r--SRC/sla_syamv.f2
-rw-r--r--SRC/slabrd.f2
-rw-r--r--SRC/slaqr1.f4
-rw-r--r--SRC/slaqr5.f22
-rw-r--r--SRC/sorcsd.f8
-rw-r--r--SRC/sorm22.f4
-rw-r--r--SRC/ssfrk.f2
-rw-r--r--SRC/sstegr.f2
-rw-r--r--SRC/sstemr.f2
-rw-r--r--SRC/stfsm.f2
-rw-r--r--SRC/stgex2.f4
-rw-r--r--SRC/zgebal.f2
-rw-r--r--SRC/zgebd2.f2
-rw-r--r--SRC/zgebrd.f2
-rw-r--r--SRC/zgesvj.f6
-rw-r--r--SRC/zhfrk.f2
-rw-r--r--SRC/zla_gbamv.f2
-rw-r--r--SRC/zla_geamv.f4
-rw-r--r--SRC/zla_heamv.f4
-rw-r--r--SRC/zla_syamv.f4
-rw-r--r--SRC/zlabrd.f2
-rw-r--r--SRC/zlaqr1.f4
-rw-r--r--SRC/zlaqr5.f20
-rw-r--r--SRC/zstegr.f2
-rw-r--r--SRC/zstemr.f2
-rw-r--r--SRC/ztgex2.f4
-rw-r--r--SRC/zuncsd.f8
-rw-r--r--SRC/zunm22.f4
67 files changed, 140 insertions, 148 deletions
diff --git a/SRC/cgebd2.f b/SRC/cgebd2.f
index e8741202..7ce92a91 100644
--- a/SRC/cgebd2.f
+++ b/SRC/cgebd2.f
@@ -100,7 +100,7 @@
*>
*> \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
diff --git a/SRC/cgebrd.f b/SRC/cgebrd.f
index d01e228a..a15a2aa0 100644
--- a/SRC/cgebrd.f
+++ b/SRC/cgebrd.f
@@ -102,7 +102,7 @@
*>
*> \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
diff --git a/SRC/cgesvj.f b/SRC/cgesvj.f
index 0168c5c3..8cce26b4 100644
--- a/SRC/cgesvj.f
+++ b/SRC/cgesvj.f
@@ -52,7 +52,7 @@
*
*> \param[in] JOBA
*> \verbatim
-*> JOBA is CHARACTER* 1
+*> JOBA is CHARACTER*1
*> Specifies the structure of A.
*> = 'L': The input matrix A is lower triangular;
*> = 'U': The input matrix A is upper triangular;
@@ -206,7 +206,7 @@
*>
*> \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)
@@ -221,7 +221,7 @@
*>
*> \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.
diff --git a/SRC/cla_heamv.f b/SRC/cla_heamv.f
index 6a3eef1b..2b630631 100644
--- a/SRC/cla_heamv.f
+++ b/SRC/cla_heamv.f
@@ -89,7 +89,7 @@
*>
*> \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.
diff --git a/SRC/cla_syamv.f b/SRC/cla_syamv.f
index 362d4559..dc8bd5dc 100644
--- a/SRC/cla_syamv.f
+++ b/SRC/cla_syamv.f
@@ -90,7 +90,7 @@
*>
*> \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.
diff --git a/SRC/clabrd.f b/SRC/clabrd.f
index 87bcb1bc..0972010e 100644
--- a/SRC/clabrd.f
+++ b/SRC/clabrd.f
@@ -111,7 +111,7 @@
*>
*> \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
diff --git a/SRC/claqr1.f b/SRC/claqr1.f
index b129bbab..b448ddb2 100644
--- a/SRC/claqr1.f
+++ b/SRC/claqr1.f
@@ -56,7 +56,7 @@
*>
*> \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
*>
@@ -81,7 +81,7 @@
*>
*> \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
diff --git a/SRC/claqr5.f b/SRC/claqr5.f
index d5c532e4..5b88d7c4 100644
--- a/SRC/claqr5.f
+++ b/SRC/claqr5.f
@@ -47,14 +47,14 @@
*
*> \param[in] WANTT
*> \verbatim
-*> WANTT is logical scalar
+*> WANTT is LOGICAL
*> WANTT = .true. if the triangular Schur factor
*> is being computed. WANTT is set to .false. otherwise.
*> \endverbatim
*>
*> \param[in] WANTZ
*> \verbatim
-*> WANTZ is logical scalar
+*> WANTZ is LOGICAL
*> WANTZ = .true. if the unitary Schur factor is being
*> computed. WANTZ is set to .false. otherwise.
*> \endverbatim
@@ -107,14 +107,14 @@
*>
*> \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
@@ -142,7 +142,7 @@
*>
*> \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.
@@ -158,7 +158,7 @@
*>
*> \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
@@ -170,8 +170,7 @@
*>
*> \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
@@ -190,7 +189,7 @@
*>
*> \param[out] WH
*> \verbatim
-*> WH is COMPLEX array of size (LDWH,NH)
+*> WH is COMPLEX array, dimension (LDWH,NH)
*> \endverbatim
*>
*> \param[in] LDWH
@@ -209,8 +208,7 @@
*>
*> \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
diff --git a/SRC/clarrv.f b/SRC/clarrv.f
index d6503013..72ef0459 100644
--- a/SRC/clarrv.f
+++ b/SRC/clarrv.f
@@ -199,7 +199,7 @@
*>
*> \param[out] Z
*> \verbatim
-*> Z is array, dimension (LDZ, max(1,M) )
+*> Z is COMPLEX array, dimension (LDZ, max(1,M) )
*> If INFO = 0, the first M columns of Z contain the
*> orthonormal eigenvectors of the matrix T
*> corresponding to the input eigenvalues, with the i-th
diff --git a/SRC/cstegr.f b/SRC/cstegr.f
index 3209f27e..77a2c210 100644
--- a/SRC/cstegr.f
+++ b/SRC/cstegr.f
@@ -184,7 +184,7 @@
*>
*> \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
diff --git a/SRC/cstemr.f b/SRC/cstemr.f
index f7e0abbe..1e9df5c3 100644
--- a/SRC/cstemr.f
+++ b/SRC/cstemr.f
@@ -239,7 +239,7 @@
*>
*> \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
diff --git a/SRC/ctgex2.f b/SRC/ctgex2.f
index c487e8f1..63ca4781 100644
--- a/SRC/ctgex2.f
+++ b/SRC/ctgex2.f
@@ -76,7 +76,7 @@
*>
*> \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
@@ -89,7 +89,7 @@
*>
*> \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
diff --git a/SRC/cuncsd.f b/SRC/cuncsd.f
index 829c4bcf..a8637e35 100644
--- a/SRC/cuncsd.f
+++ b/SRC/cuncsd.f
@@ -188,7 +188,7 @@
*>
*> \param[out] U1
*> \verbatim
-*> U1 is COMPLEX array, dimension (P)
+*> U1 is COMPLEX array, dimension (P,P)
*> If JOBU1 = 'Y', U1 contains the P-by-P unitary matrix U1.
*> \endverbatim
*>
@@ -201,7 +201,7 @@
*>
*> \param[out] U2
*> \verbatim
-*> U2 is COMPLEX array, dimension (M-P)
+*> U2 is COMPLEX array, dimension (M-P,M-P)
*> If JOBU2 = 'Y', U2 contains the (M-P)-by-(M-P) unitary
*> matrix U2.
*> \endverbatim
@@ -215,7 +215,7 @@
*>
*> \param[out] V1T
*> \verbatim
-*> V1T is COMPLEX array, dimension (Q)
+*> V1T is COMPLEX array, dimension (Q,Q)
*> If JOBV1T = 'Y', V1T contains the Q-by-Q matrix unitary
*> matrix V1**H.
*> \endverbatim
@@ -229,7 +229,7 @@
*>
*> \param[out] V2T
*> \verbatim
-*> V2T is COMPLEX array, dimension (M-Q)
+*> V2T is COMPLEX array, dimension (M-Q,M-Q)
*> If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) unitary
*> matrix V2**H.
*> \endverbatim
diff --git a/SRC/cunm22.f b/SRC/cunm22.f
index 85c2269a..f62939bc 100644
--- a/SRC/cunm22.f
+++ b/SRC/cunm22.f
@@ -52,8 +52,8 @@
*> N2-by-N2 upper triangular matrix.
*> \endverbatim
*
-* Arguments
-* =========
+* Arguments:
+* ==========
*
*> \param[in] SIDE
*> \verbatim
diff --git a/SRC/dgebal.f b/SRC/dgebal.f
index 93efd289..0527220a 100644
--- a/SRC/dgebal.f
+++ b/SRC/dgebal.f
@@ -67,7 +67,7 @@
*>
*> \param[in,out] A
*> \verbatim
-*> A is DOUBLE array, dimension (LDA,N)
+*> A is DOUBLE PRECISION array, dimension (LDA,N)
*> On entry, the input matrix A.
*> On exit, A is overwritten by the balanced matrix.
*> If JOB = 'N', A is not referenced.
@@ -94,7 +94,7 @@
*>
*> \param[out] SCALE
*> \verbatim
-*> SCALE is DOUBLE array, dimension (N)
+*> SCALE is DOUBLE PRECISION array, dimension (N)
*> Details of the permutations and scaling factors applied to
*> A. If P(j) is the index of the row and column interchanged
*> with row and column j and D(j) is the scaling factor
diff --git a/SRC/dgebd2.f b/SRC/dgebd2.f
index bb4035db..c7a7d787 100644
--- a/SRC/dgebd2.f
+++ b/SRC/dgebd2.f
@@ -100,7 +100,7 @@
*>
*> \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
diff --git a/SRC/dgebrd.f b/SRC/dgebrd.f
index 885ad9bb..c8780be3 100644
--- a/SRC/dgebrd.f
+++ b/SRC/dgebrd.f
@@ -101,7 +101,7 @@
*>
*> \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
diff --git a/SRC/dgesvj.f b/SRC/dgesvj.f
index 2b259942..09025c23 100644
--- a/SRC/dgesvj.f
+++ b/SRC/dgesvj.f
@@ -54,7 +54,7 @@
*
*> \param[in] JOBA
*> \verbatim
-*> JOBA is CHARACTER* 1
+*> JOBA is CHARACTER*1
*> Specifies the structure of A.
*> = 'L': The input matrix A is lower triangular;
*> = 'U': The input matrix A is upper triangular;
@@ -208,7 +208,7 @@
*>
*> \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.
diff --git a/SRC/dla_gbamv.f b/SRC/dla_gbamv.f
index 577866cf..707177e9 100644
--- a/SRC/dla_gbamv.f
+++ b/SRC/dla_gbamv.f
@@ -106,7 +106,7 @@
*>
*> \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.
diff --git a/SRC/dla_geamv.f b/SRC/dla_geamv.f
index 9a91f6ff..5b6690b8 100644
--- a/SRC/dla_geamv.f
+++ b/SRC/dla_geamv.f
@@ -94,7 +94,7 @@
*>
*> \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.
diff --git a/SRC/dla_syamv.f b/SRC/dla_syamv.f
index 1f948a2d..127371ba 100644
--- a/SRC/dla_syamv.f
+++ b/SRC/dla_syamv.f
@@ -88,7 +88,7 @@
*>
*> \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.
diff --git a/SRC/dlabrd.f b/SRC/dlabrd.f
index 36c2e85b..0cad5089 100644
--- a/SRC/dlabrd.f
+++ b/SRC/dlabrd.f
@@ -110,7 +110,7 @@
*>
*> \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
diff --git a/SRC/dlaqr1.f b/SRC/dlaqr1.f
index 25021d5a..c3e6e8ab 100644
--- a/SRC/dlaqr1.f
+++ b/SRC/dlaqr1.f
@@ -61,7 +61,7 @@
*>
*> \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
*>
@@ -95,7 +95,7 @@
*>
*> \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
diff --git a/SRC/dlaqr5.f b/SRC/dlaqr5.f
index 9b8763f0..a5a49610 100644
--- a/SRC/dlaqr5.f
+++ b/SRC/dlaqr5.f
@@ -48,14 +48,14 @@
*
*> \param[in] WANTT
*> \verbatim
-*> WANTT is logical scalar
+*> WANTT is LOGICAL
*> WANTT = .true. if the quasi-triangular Schur factor
*> is being computed. WANTT is set to .false. otherwise.
*> \endverbatim
*>
*> \param[in] WANTZ
*> \verbatim
-*> WANTZ is logical scalar
+*> WANTZ is LOGICAL
*> WANTZ = .true. if the orthogonal Schur factor is being
*> computed. WANTZ is set to .false. otherwise.
*> \endverbatim
@@ -108,12 +108,12 @@
*>
*> \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
@@ -122,7 +122,7 @@
*>
*> \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
@@ -150,7 +150,7 @@
*>
*> \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.
@@ -166,7 +166,7 @@
*>
*> \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
@@ -178,8 +178,7 @@
*>
*> \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
@@ -198,7 +197,7 @@
*>
*> \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
@@ -217,8 +216,7 @@
*>
*> \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
diff --git a/SRC/dlasd0.f b/SRC/dlasd0.f
index ca0b3b98..832481a4 100644
--- a/SRC/dlasd0.f
+++ b/SRC/dlasd0.f
@@ -112,14 +112,12 @@
*>
*> \param[out] IWORK
*> \verbatim
-*> IWORK is INTEGER work array.
-*> Dimension must be at least (8 * N)
+*> IWORK is INTEGER array, dimension (8*N)
*> \endverbatim
*>
*> \param[out] WORK
*> \verbatim
-*> WORK is DOUBLE PRECISION work array.
-*> Dimension must be at least (3 * M**2 + 2 * M)
+*> WORK is DOUBLE PRECISION array, dimension (3*M**2+2*M)
*> \endverbatim
*>
*> \param[out] INFO
diff --git a/SRC/dlasd2.f b/SRC/dlasd2.f
index a7ced418..cbc10c42 100644
--- a/SRC/dlasd2.f
+++ b/SRC/dlasd2.f
@@ -190,7 +190,7 @@
*>
*> \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)
@@ -199,14 +199,14 @@
*>
*> \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
@@ -216,7 +216,7 @@
*>
*> \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
@@ -226,7 +226,7 @@
*>
*> \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:
diff --git a/SRC/dlasq5.f b/SRC/dlasq5.f
index 99d4f678..0eb1fd30 100644
--- a/SRC/dlasq5.f
+++ b/SRC/dlasq5.f
@@ -121,7 +121,7 @@
*> IEEE is LOGICAL
*> Flag for IEEE or non IEEE arithmetic.
*> \endverbatim
-*
+*>
*> \param[in] EPS
*> \verbatim
*> EPS is DOUBLE PRECISION
diff --git a/SRC/dorcsd.f b/SRC/dorcsd.f
index 340e16a5..34c1dbba 100644
--- a/SRC/dorcsd.f
+++ b/SRC/dorcsd.f
@@ -186,7 +186,7 @@
*>
*> \param[out] U1
*> \verbatim
-*> U1 is DOUBLE PRECISION array, dimension (P)
+*> U1 is DOUBLE PRECISION array, dimension (P,P)
*> If JOBU1 = 'Y', U1 contains the P-by-P orthogonal matrix U1.
*> \endverbatim
*>
@@ -199,7 +199,7 @@
*>
*> \param[out] U2
*> \verbatim
-*> U2 is DOUBLE PRECISION array, dimension (M-P)
+*> U2 is DOUBLE PRECISION array, dimension (M-P,M-P)
*> If JOBU2 = 'Y', U2 contains the (M-P)-by-(M-P) orthogonal
*> matrix U2.
*> \endverbatim
@@ -213,7 +213,7 @@
*>
*> \param[out] V1T
*> \verbatim
-*> V1T is DOUBLE PRECISION array, dimension (Q)
+*> V1T is DOUBLE PRECISION array, dimension (Q,Q)
*> If JOBV1T = 'Y', V1T contains the Q-by-Q matrix orthogonal
*> matrix V1**T.
*> \endverbatim
@@ -227,7 +227,7 @@
*>
*> \param[out] V2T
*> \verbatim
-*> V2T is DOUBLE PRECISION array, dimension (M-Q)
+*> V2T is DOUBLE PRECISION array, dimension (M-Q,M-Q)
*> If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) orthogonal
*> matrix V2**T.
*> \endverbatim
diff --git a/SRC/dorm22.f b/SRC/dorm22.f
index ac79e1e7..3e17d493 100644
--- a/SRC/dorm22.f
+++ b/SRC/dorm22.f
@@ -53,8 +53,8 @@
*> N2-by-N2 upper triangular matrix.
*> \endverbatim
*
-* Arguments
-* =========
+* Arguments:
+* ==========
*
*> \param[in] SIDE
*> \verbatim
diff --git a/SRC/dstegr.f b/SRC/dstegr.f
index f3286032..84c08be1 100644
--- a/SRC/dstegr.f
+++ b/SRC/dstegr.f
@@ -184,7 +184,7 @@
*>
*> \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
diff --git a/SRC/dstemr.f b/SRC/dstemr.f
index 924d738d..c50183b9 100644
--- a/SRC/dstemr.f
+++ b/SRC/dstemr.f
@@ -222,7 +222,7 @@
*>
*> \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
diff --git a/SRC/dtgsen.f b/SRC/dtgsen.f
index 9f492393..c61b613c 100644
--- a/SRC/dtgsen.f
+++ b/SRC/dtgsen.f
@@ -222,7 +222,7 @@
*> \verbatim
*> PL is DOUBLE PRECISION
*> \endverbatim
-
+*>
*> \param[out] PR
*> \verbatim
*> PR is DOUBLE PRECISION
diff --git a/SRC/ilaclr.f b/SRC/ilaclr.f
index c2e0584b..dcb4fe2d 100644
--- a/SRC/ilaclr.f
+++ b/SRC/ilaclr.f
@@ -53,7 +53,7 @@
*>
*> \param[in] A
*> \verbatim
-*> A is array, dimension (LDA,N)
+*> A is COMPLEX array, dimension (LDA,N)
*> The m by n matrix A.
*> \endverbatim
*>
diff --git a/SRC/sgebd2.f b/SRC/sgebd2.f
index eb125fba..d00d189d 100644
--- a/SRC/sgebd2.f
+++ b/SRC/sgebd2.f
@@ -100,7 +100,7 @@
*>
*> \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
diff --git a/SRC/sgebrd.f b/SRC/sgebrd.f
index 0f38156d..2666c233 100644
--- a/SRC/sgebrd.f
+++ b/SRC/sgebrd.f
@@ -101,7 +101,7 @@
*>
*> \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
diff --git a/SRC/sgees.f b/SRC/sgees.f
index 06319bf7..ad65a0c2 100644
--- a/SRC/sgees.f
+++ b/SRC/sgees.f
@@ -80,7 +80,7 @@
*>
*> \param[in] SELECT
*> \verbatim
-*> SELECT is LOGICAL FUNCTION of two REAL arguments
+*> SELECT is a LOGICAL FUNCTION of two REAL arguments
*> SELECT must be declared EXTERNAL in the calling subroutine.
*> If SORT = 'S', SELECT is used to select eigenvalues to sort
*> to the top left of the Schur form.
diff --git a/SRC/sgesvj.f b/SRC/sgesvj.f
index 5e53cea1..adce3451 100644
--- a/SRC/sgesvj.f
+++ b/SRC/sgesvj.f
@@ -54,7 +54,7 @@
*
*> \param[in] JOBA
*> \verbatim
-*> JOBA is CHARACTER* 1
+*> JOBA is CHARACTER*1
*> Specifies the structure of A.
*> = 'L': The input matrix A is lower triangular;
*> = 'U': The input matrix A is upper triangular;
@@ -208,7 +208,7 @@
*>
*> \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.
diff --git a/SRC/sla_gbamv.f b/SRC/sla_gbamv.f
index 0798bacd..073f5610 100644
--- a/SRC/sla_gbamv.f
+++ b/SRC/sla_gbamv.f
@@ -106,7 +106,7 @@
*>
*> \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.
diff --git a/SRC/sla_geamv.f b/SRC/sla_geamv.f
index 35ce8b80..083a823b 100644
--- a/SRC/sla_geamv.f
+++ b/SRC/sla_geamv.f
@@ -94,7 +94,7 @@
*>
*> \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.
diff --git a/SRC/sla_syamv.f b/SRC/sla_syamv.f
index 962e17ac..aad2fbf2 100644
--- a/SRC/sla_syamv.f
+++ b/SRC/sla_syamv.f
@@ -88,7 +88,7 @@
*>
*> \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.
diff --git a/SRC/slabrd.f b/SRC/slabrd.f
index 8073d003..6321c191 100644
--- a/SRC/slabrd.f
+++ b/SRC/slabrd.f
@@ -110,7 +110,7 @@
*>
*> \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
diff --git a/SRC/slaqr1.f b/SRC/slaqr1.f
index ca8e28f4..318b976e 100644
--- a/SRC/slaqr1.f
+++ b/SRC/slaqr1.f
@@ -61,7 +61,7 @@
*>
*> \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
*>
@@ -95,7 +95,7 @@
*>
*> \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
diff --git a/SRC/slaqr5.f b/SRC/slaqr5.f
index 5198c1fb..59692bd0 100644
--- a/SRC/slaqr5.f
+++ b/SRC/slaqr5.f
@@ -48,14 +48,14 @@
*
*> \param[in] WANTT
*> \verbatim
-*> WANTT is logical scalar
+*> WANTT is LOGICAL
*> WANTT = .true. if the quasi-triangular Schur factor
*> is being computed. WANTT is set to .false. otherwise.
*> \endverbatim
*>
*> \param[in] WANTZ
*> \verbatim
-*> WANTZ is logical scalar
+*> WANTZ is LOGICAL
*> WANTZ = .true. if the orthogonal Schur factor is being
*> computed. WANTZ is set to .false. otherwise.
*> \endverbatim
@@ -108,12 +108,12 @@
*>
*> \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
@@ -122,7 +122,7 @@
*>
*> \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
@@ -150,7 +150,7 @@
*>
*> \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.
@@ -166,7 +166,7 @@
*>
*> \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
@@ -178,8 +178,7 @@
*>
*> \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
@@ -198,7 +197,7 @@
*>
*> \param[out] WH
*> \verbatim
-*> WH is REAL array of size (LDWH,NH)
+*> WH is REAL array, dimension (LDWH,NH)
*> \endverbatim
*>
*> \param[in] LDWH
@@ -217,8 +216,7 @@
*>
*> \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
diff --git a/SRC/sorcsd.f b/SRC/sorcsd.f
index 5621d58e..9694ef54 100644
--- a/SRC/sorcsd.f
+++ b/SRC/sorcsd.f
@@ -186,7 +186,7 @@
*>
*> \param[out] U1
*> \verbatim
-*> U1 is REAL array, dimension (P)
+*> U1 is REAL array, dimension (P,P)
*> If JOBU1 = 'Y', U1 contains the P-by-P orthogonal matrix U1.
*> \endverbatim
*>
@@ -199,7 +199,7 @@
*>
*> \param[out] U2
*> \verbatim
-*> U2 is REAL array, dimension (M-P)
+*> U2 is REAL array, dimension (M-P,M-P)
*> If JOBU2 = 'Y', U2 contains the (M-P)-by-(M-P) orthogonal
*> matrix U2.
*> \endverbatim
@@ -213,7 +213,7 @@
*>
*> \param[out] V1T
*> \verbatim
-*> V1T is REAL array, dimension (Q)
+*> V1T is REAL array, dimension (Q,Q)
*> If JOBV1T = 'Y', V1T contains the Q-by-Q matrix orthogonal
*> matrix V1**T.
*> \endverbatim
@@ -227,7 +227,7 @@
*>
*> \param[out] V2T
*> \verbatim
-*> V2T is REAL array, dimension (M-Q)
+*> V2T is REAL array, dimension (M-Q,M-Q)
*> If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) orthogonal
*> matrix V2**T.
*> \endverbatim
diff --git a/SRC/sorm22.f b/SRC/sorm22.f
index fdb5cd8b..06f18d68 100644
--- a/SRC/sorm22.f
+++ b/SRC/sorm22.f
@@ -53,8 +53,8 @@
*> N2-by-N2 upper triangular matrix.
*> \endverbatim
*
-* Arguments
-* =========
+* Arguments:
+* ==========
*
*> \param[in] SIDE
*> \verbatim
diff --git a/SRC/ssfrk.f b/SRC/ssfrk.f
index b2cc4fca..a0c0e9ce 100644
--- a/SRC/ssfrk.f
+++ b/SRC/ssfrk.f
@@ -117,7 +117,7 @@
*>
*> \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
diff --git a/SRC/sstegr.f b/SRC/sstegr.f
index 0a28c584..99a4be62 100644
--- a/SRC/sstegr.f
+++ b/SRC/sstegr.f
@@ -184,7 +184,7 @@
*>
*> \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
diff --git a/SRC/sstemr.f b/SRC/sstemr.f
index cff89ef9..917e896c 100644
--- a/SRC/sstemr.f
+++ b/SRC/sstemr.f
@@ -222,7 +222,7 @@
*>
*> \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
diff --git a/SRC/stfsm.f b/SRC/stfsm.f
index b8b81a87..4ac3c3f5 100644
--- a/SRC/stfsm.f
+++ b/SRC/stfsm.f
@@ -159,7 +159,7 @@
*>
*> \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.
diff --git a/SRC/stgex2.f b/SRC/stgex2.f
index fc5f9330..d7795b14 100644
--- a/SRC/stgex2.f
+++ b/SRC/stgex2.f
@@ -77,7 +77,7 @@
*>
*> \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
@@ -90,7 +90,7 @@
*>
*> \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
diff --git a/SRC/zgebal.f b/SRC/zgebal.f
index 601d5431..94b8aadf 100644
--- a/SRC/zgebal.f
+++ b/SRC/zgebal.f
@@ -83,10 +83,12 @@
*>
*> \param[out] ILO
*> \verbatim
+*> ILO is INTEGER
*> \endverbatim
*>
*> \param[out] IHI
*> \verbatim
+*> IHI is INTEGER
*> ILO and IHI are set to INTEGER such that on exit
*> A(i,j) = 0 if i > j and j = 1,...,ILO-1 or I = IHI+1,...,N.
*> If JOB = 'N' or 'S', ILO = 1 and IHI = N.
diff --git a/SRC/zgebd2.f b/SRC/zgebd2.f
index d5752e04..77779d6c 100644
--- a/SRC/zgebd2.f
+++ b/SRC/zgebd2.f
@@ -100,7 +100,7 @@
*>
*> \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
diff --git a/SRC/zgebrd.f b/SRC/zgebrd.f
index 26879a75..47db32f4 100644
--- a/SRC/zgebrd.f
+++ b/SRC/zgebrd.f
@@ -101,7 +101,7 @@
*>
*> \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
diff --git a/SRC/zgesvj.f b/SRC/zgesvj.f
index 3670ec66..1e1be50c 100644
--- a/SRC/zgesvj.f
+++ b/SRC/zgesvj.f
@@ -52,7 +52,7 @@
*
*> \param[in] JOBA
*> \verbatim
-*> JOBA is CHARACTER* 1
+*> JOBA is CHARACTER*1
*> Specifies the structure of A.
*> = 'L': The input matrix A is lower triangular;
*> = 'U': The input matrix A is upper triangular;
@@ -206,7 +206,7 @@
*>
*> \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)
@@ -221,7 +221,7 @@
*>
*> \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.
diff --git a/SRC/zhfrk.f b/SRC/zhfrk.f
index cfc3e111..a4d55078 100644
--- a/SRC/zhfrk.f
+++ b/SRC/zhfrk.f
@@ -117,7 +117,7 @@
*>
*> \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
diff --git a/SRC/zla_gbamv.f b/SRC/zla_gbamv.f
index 9d5291b8..673535fa 100644
--- a/SRC/zla_gbamv.f
+++ b/SRC/zla_gbamv.f
@@ -107,7 +107,7 @@
*>
*> \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.
diff --git a/SRC/zla_geamv.f b/SRC/zla_geamv.f
index 8d221691..ffcf6890 100644
--- a/SRC/zla_geamv.f
+++ b/SRC/zla_geamv.f
@@ -96,7 +96,7 @@
*>
*> \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.
@@ -113,7 +113,7 @@
*>
*> \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.
diff --git a/SRC/zla_heamv.f b/SRC/zla_heamv.f
index 1fa0e0a8..cc6407d6 100644
--- a/SRC/zla_heamv.f
+++ b/SRC/zla_heamv.f
@@ -89,7 +89,7 @@
*>
*> \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.
@@ -106,7 +106,7 @@
*>
*> \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.
diff --git a/SRC/zla_syamv.f b/SRC/zla_syamv.f
index 01854a9b..0f1f4e05 100644
--- a/SRC/zla_syamv.f
+++ b/SRC/zla_syamv.f
@@ -90,7 +90,7 @@
*>
*> \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.
@@ -107,7 +107,7 @@
*>
*> \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.
diff --git a/SRC/zlabrd.f b/SRC/zlabrd.f
index 07b5e9fc..812a973b 100644
--- a/SRC/zlabrd.f
+++ b/SRC/zlabrd.f
@@ -111,7 +111,7 @@
*>
*> \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
diff --git a/SRC/zlaqr1.f b/SRC/zlaqr1.f
index 50b3a682..cc7200eb 100644
--- a/SRC/zlaqr1.f
+++ b/SRC/zlaqr1.f
@@ -56,7 +56,7 @@
*>
*> \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
*>
@@ -81,7 +81,7 @@
*>
*> \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
diff --git a/SRC/zlaqr5.f b/SRC/zlaqr5.f
index 3e463b4a..422d7b21 100644
--- a/SRC/zlaqr5.f
+++ b/SRC/zlaqr5.f
@@ -47,14 +47,14 @@
*
*> \param[in] WANTT
*> \verbatim
-*> WANTT is logical scalar
+*> WANTT is LOGICAL
*> WANTT = .true. if the triangular Schur factor
*> is being computed. WANTT is set to .false. otherwise.
*> \endverbatim
*>
*> \param[in] WANTZ
*> \verbatim
-*> WANTZ is logical scalar
+*> WANTZ is LOGICAL
*> WANTZ = .true. if the unitary Schur factor is being
*> computed. WANTZ is set to .false. otherwise.
*> \endverbatim
@@ -107,14 +107,14 @@
*>
*> \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
@@ -142,7 +142,7 @@
*>
*> \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.
@@ -158,7 +158,7 @@
*>
*> \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
@@ -170,8 +170,7 @@
*>
*> \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
@@ -190,7 +189,7 @@
*>
*> \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
@@ -209,8 +208,7 @@
*>
*> \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
diff --git a/SRC/zstegr.f b/SRC/zstegr.f
index 71f304c1..b8d20729 100644
--- a/SRC/zstegr.f
+++ b/SRC/zstegr.f
@@ -184,7 +184,7 @@
*>
*> \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
diff --git a/SRC/zstemr.f b/SRC/zstemr.f
index 681c87ad..7fa8f3cb 100644
--- a/SRC/zstemr.f
+++ b/SRC/zstemr.f
@@ -239,7 +239,7 @@
*>
*> \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
diff --git a/SRC/ztgex2.f b/SRC/ztgex2.f
index 26b83dec..3f484285 100644
--- a/SRC/ztgex2.f
+++ b/SRC/ztgex2.f
@@ -76,7 +76,7 @@
*>
*> \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
@@ -89,7 +89,7 @@
*>
*> \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
diff --git a/SRC/zuncsd.f b/SRC/zuncsd.f
index 77a83c09..18dce987 100644
--- a/SRC/zuncsd.f
+++ b/SRC/zuncsd.f
@@ -188,7 +188,7 @@
*>
*> \param[out] U1
*> \verbatim
-*> U1 is COMPLEX*16 array, dimension (P)
+*> U1 is COMPLEX*16 array, dimension (P,P)
*> If JOBU1 = 'Y', U1 contains the P-by-P unitary matrix U1.
*> \endverbatim
*>
@@ -201,7 +201,7 @@
*>
*> \param[out] U2
*> \verbatim
-*> U2 is COMPLEX*16 array, dimension (M-P)
+*> U2 is COMPLEX*16 array, dimension (M-P,M-P)
*> If JOBU2 = 'Y', U2 contains the (M-P)-by-(M-P) unitary
*> matrix U2.
*> \endverbatim
@@ -215,7 +215,7 @@
*>
*> \param[out] V1T
*> \verbatim
-*> V1T is COMPLEX*16 array, dimension (Q)
+*> V1T is COMPLEX*16 array, dimension (Q,Q)
*> If JOBV1T = 'Y', V1T contains the Q-by-Q matrix unitary
*> matrix V1**H.
*> \endverbatim
@@ -229,7 +229,7 @@
*>
*> \param[out] V2T
*> \verbatim
-*> V2T is COMPLEX*16 array, dimension (M-Q)
+*> V2T is COMPLEX*16 array, dimension (M-Q,M-Q)
*> If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) unitary
*> matrix V2**H.
*> \endverbatim
diff --git a/SRC/zunm22.f b/SRC/zunm22.f
index 468d7d8c..01f8b73c 100644
--- a/SRC/zunm22.f
+++ b/SRC/zunm22.f
@@ -52,8 +52,8 @@
*> N2-by-N2 upper triangular matrix.
*> \endverbatim
*
-* Arguments
-* =========
+* Arguments:
+* ==========
*
*> \param[in] SIDE
*> \verbatim