Merge pull request #162 from thielema/array-dimensions

more consistency of array dimensions in Fortran declarations and Doxygen comments

26 files changed, 40 insertions, 43 deletions

diff --git a/SRC/csytri2x.f b/SRC/csytri2x.f
index 046d6134..4201a4f7 100644
--- a/SRC/csytri2x.f
+++ b/SRC/csytri2x.f
@@ -87,7 +87,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is COMPLEX array, dimension (N+NNB+1,NNB+3)
+*>          WORK is COMPLEX array, dimension (N+NB+1,NB+3)
 *> \endverbatim
 *>
 *> \param[in] NB
diff --git a/SRC/ctgex2.f b/SRC/ctgex2.f
index 63ca4781..328b14a0 100644
--- a/SRC/ctgex2.f
+++ b/SRC/ctgex2.f
@@ -102,7 +102,7 @@
 *>
 *> \param[in,out] Q
 *> \verbatim
-*>          Q is COMPLEX array, dimension (LDZ,N)
+*>          Q is COMPLEX array, dimension (LDQ,N)
 *>          If WANTQ = .TRUE, on entry, the unitary matrix Q. On exit,
 *>          the updated matrix Q.
 *>          Not referenced if WANTQ = .FALSE..
diff --git a/SRC/ctgexc.f b/SRC/ctgexc.f
index f6ccdcb9..8f6fe5c8 100644
--- a/SRC/ctgexc.f
+++ b/SRC/ctgexc.f
@@ -102,7 +102,7 @@
 *>
 *> \param[in,out] Q
 *> \verbatim
-*>          Q is COMPLEX array, dimension (LDZ,N)
+*>          Q is COMPLEX array, dimension (LDQ,N)
 *>          On entry, if WANTQ = .TRUE., the unitary matrix Q.
 *>          On exit, the updated matrix Q.
 *>          If WANTQ = .FALSE., Q is not referenced.
diff --git a/SRC/cuncsd.f b/SRC/cuncsd.f
index a8637e35..b01036dd 100644
--- a/SRC/cuncsd.f
+++ b/SRC/cuncsd.f
@@ -188,7 +188,7 @@
 *>
 *> \param[out] U1
 *> \verbatim
-*>          U1 is COMPLEX array, dimension (P,P)
+*>          U1 is COMPLEX array, dimension (LDU1,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,M-P)
+*>          U2 is COMPLEX array, dimension (LDU2,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,Q)
+*>          V1T is COMPLEX array, dimension (LDV1T,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,M-Q)
+*>          V2T is COMPLEX array, dimension (LDV2T,M-Q)
 *>          If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) unitary
 *>          matrix V2**H.
 *> \endverbatim
diff --git a/SRC/dgejsv.f b/SRC/dgejsv.f
index 865e7c46..c230dfb1 100644
--- a/SRC/dgejsv.f
+++ b/SRC/dgejsv.f
@@ -271,7 +271,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is DOUBLE PRECISION array, dimension at least LWORK.
+*>          WORK is DOUBLE PRECISION array, dimension (LWORK)
 *>          On exit, if N.GT.0 .AND. M.GT.0 (else not referenced),
 *>          WORK(1) = SCALE = WORK(2) / WORK(1) is the scaling factor such
 *>                    that SCALE*SVA(1:N) are the computed singular values
diff --git a/SRC/dgesvj.f b/SRC/dgesvj.f
index 09025c23..1e8933b6 100644
--- a/SRC/dgesvj.f
+++ b/SRC/dgesvj.f
@@ -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 (LWORK)
 *>          On entry :
 *>          If JOBU .EQ. 'C' :
 *>          WORK(1) = CTOL, where CTOL defines the threshold for convergence.
diff --git a/SRC/dlasd0.f b/SRC/dlasd0.f
index 832481a4..f639f81c 100644
--- a/SRC/dlasd0.f
+++ b/SRC/dlasd0.f
@@ -81,7 +81,7 @@
 *>
 *> \param[out] U
 *> \verbatim
-*>          U is DOUBLE PRECISION array, dimension at least (LDQ, N)
+*>          U is DOUBLE PRECISION array, dimension (LDU, N)
 *>         On exit, U contains the left singular vectors.
 *> \endverbatim
 *>
@@ -93,7 +93,7 @@
 *>
 *> \param[out] VT
 *> \verbatim
-*>          VT is DOUBLE PRECISION array, dimension at least (LDVT, M)
+*>          VT is DOUBLE PRECISION array, dimension (LDVT, M)
 *>         On exit, VT**T contains the right singular vectors.
 *> \endverbatim
 *>
diff --git a/SRC/dlasd3.f b/SRC/dlasd3.f
index 57d0abd4..11be8872 100644
--- a/SRC/dlasd3.f
+++ b/SRC/dlasd3.f
@@ -94,8 +94,7 @@
 *>
 *> \param[out] Q
 *> \verbatim
-*>          Q is DOUBLE PRECISION array,
-*>                     dimension at least (LDQ,K).
+*>          Q is DOUBLE PRECISION array, dimension (LDQ,K)
 *> \endverbatim
 *>
 *> \param[in] LDQ
diff --git a/SRC/dlasd8.f b/SRC/dlasd8.f
index 245e814a..5aaab5ad 100644
--- a/SRC/dlasd8.f
+++ b/SRC/dlasd8.f
@@ -133,7 +133,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is DOUBLE PRECISION array, dimension at least 3 * K
+*>          WORK is DOUBLE PRECISION array, dimension (3*K)
 *> \endverbatim
 *>
 *> \param[out] INFO
diff --git a/SRC/dorcsd.f b/SRC/dorcsd.f
index 34c1dbba..61b70ffd 100644
--- a/SRC/dorcsd.f
+++ b/SRC/dorcsd.f
@@ -186,7 +186,7 @@
 *>
 *> \param[out] U1
 *> \verbatim
-*>          U1 is DOUBLE PRECISION array, dimension (P,P)
+*>          U1 is DOUBLE PRECISION array, dimension (LDU1,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,M-P)
+*>          U2 is DOUBLE PRECISION array, dimension (LDU2,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,Q)
+*>          V1T is DOUBLE PRECISION array, dimension (LDV1T,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,M-Q)
+*>          V2T is DOUBLE PRECISION array, dimension (LDV2T,M-Q)
 *>          If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) orthogonal
 *>          matrix V2**T.
 *> \endverbatim
diff --git a/SRC/dsytri2x.f b/SRC/dsytri2x.f
index bcd5c942..77fda2de 100644
--- a/SRC/dsytri2x.f
+++ b/SRC/dsytri2x.f
@@ -87,7 +87,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is DOUBLE PRECISION array, dimension (N+NNB+1,NNB+3)
+*>          WORK is DOUBLE PRECISION array, dimension (N+NB+1,NB+3)
 *> \endverbatim
 *>
 *> \param[in] NB
diff --git a/SRC/sgejsv.f b/SRC/sgejsv.f
index b302844e..98b1dcd0 100644
--- a/SRC/sgejsv.f
+++ b/SRC/sgejsv.f
@@ -271,7 +271,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is REAL array, dimension at least LWORK.
+*>          WORK is REAL array, dimension (LWORK)
 *>          On exit,
 *>          WORK(1) = SCALE = WORK(2) / WORK(1) is the scaling factor such
 *>                    that SCALE*SVA(1:N) are the computed singular values
diff --git a/SRC/sgesvj.f b/SRC/sgesvj.f
index adce3451..8763ab4b 100644
--- a/SRC/sgesvj.f
+++ b/SRC/sgesvj.f
@@ -208,7 +208,7 @@
 *>
 *> \param[in,out] WORK
 *> \verbatim
-*>          WORK is REAL array, dimension (max(6,M+N))
+*>          WORK is REAL array, dimension (LWORK)
 *>          On entry,
 *>          If JOBU .EQ. 'C' :
 *>          WORK(1) = CTOL, where CTOL defines the threshold for convergence.
diff --git a/SRC/slasd0.f b/SRC/slasd0.f
index b3eb0735..60b05eb4 100644
--- a/SRC/slasd0.f
+++ b/SRC/slasd0.f
@@ -81,7 +81,7 @@
 *>
 *> \param[out] U
 *> \verbatim
-*>          U is REAL array, dimension at least (LDQ, N)
+*>          U is REAL array, dimension (LDU, N)
 *>         On exit, U contains the left singular vectors.
 *> \endverbatim
 *>
@@ -93,7 +93,7 @@
 *>
 *> \param[out] VT
 *> \verbatim
-*>          VT is REAL array, dimension at least (LDVT, M)
+*>          VT is REAL array, dimension (LDVT, M)
 *>         On exit, VT**T contains the right singular vectors.
 *> \endverbatim
 *>
diff --git a/SRC/slasd3.f b/SRC/slasd3.f
index 6b1d0f00..f9b0c01d 100644
--- a/SRC/slasd3.f
+++ b/SRC/slasd3.f
@@ -94,8 +94,7 @@
 *>
 *> \param[out] Q
 *> \verbatim
-*>          Q is REAL array,
-*>                     dimension at least (LDQ,K).
+*>          Q is REAL array, dimension (LDQ,K)
 *> \endverbatim
 *>
 *> \param[in] LDQ
diff --git a/SRC/slasd8.f b/SRC/slasd8.f
index 81a8625a..578d569c 100644
--- a/SRC/slasd8.f
+++ b/SRC/slasd8.f
@@ -133,7 +133,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is REAL array, dimension at least 3 * K
+*>          WORK is REAL array, dimension (3*K)
 *> \endverbatim
 *>
 *> \param[out] INFO
diff --git a/SRC/sorcsd.f b/SRC/sorcsd.f
index 9694ef54..bba032f8 100644
--- a/SRC/sorcsd.f
+++ b/SRC/sorcsd.f
@@ -186,7 +186,7 @@
 *>
 *> \param[out] U1
 *> \verbatim
-*>          U1 is REAL array, dimension (P,P)
+*>          U1 is REAL array, dimension (LDU1,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,M-P)
+*>          U2 is REAL array, dimension (LDU2,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,Q)
+*>          V1T is REAL array, dimension (LDV1T,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,M-Q)
+*>          V2T is REAL array, dimension (LDV2T,M-Q)
 *>          If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) orthogonal
 *>          matrix V2**T.
 *> \endverbatim
diff --git a/SRC/ssygvx.f b/SRC/ssygvx.f
index 2982e172..0e8f6515 100644
--- a/SRC/ssygvx.f
+++ b/SRC/ssygvx.f
@@ -110,7 +110,7 @@
 *>
 *> \param[in,out] B
 *> \verbatim
-*>          B is REAL array, dimension (LDA, N)
+*>          B is REAL array, dimension (LDB, N)
 *>          On entry, the symmetric matrix B.  If UPLO = 'U', the
 *>          leading N-by-N upper triangular part of B contains the
 *>          upper triangular part of the matrix B.  If UPLO = 'L',
diff --git a/SRC/ssytri2x.f b/SRC/ssytri2x.f
index c218383a..6d15fad1 100644
--- a/SRC/ssytri2x.f
+++ b/SRC/ssytri2x.f
@@ -87,7 +87,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is REAL array, dimension (N+NNB+1,NNB+3)
+*>          WORK is REAL array, dimension (N+NB+1,NB+3)
 *> \endverbatim
 *>
 *> \param[in] NB
diff --git a/SRC/stgex2.f b/SRC/stgex2.f
index d7795b14..7fb67fcc 100644
--- a/SRC/stgex2.f
+++ b/SRC/stgex2.f
@@ -103,7 +103,7 @@
 *>
 *> \param[in,out] Q
 *> \verbatim
-*>          Q is REAL array, dimension (LDZ,N)
+*>          Q is REAL array, dimension (LDQ,N)
 *>          On entry, if WANTQ = .TRUE., the orthogonal matrix Q.
 *>          On exit, the updated matrix Q.
 *>          Not referenced if WANTQ = .FALSE..
diff --git a/SRC/stgexc.f b/SRC/stgexc.f
index a1d9dbb5..0bfda7ff 100644
--- a/SRC/stgexc.f
+++ b/SRC/stgexc.f
@@ -111,7 +111,7 @@
 *>
 *> \param[in,out] Q
 *> \verbatim
-*>          Q is REAL array, dimension (LDZ,N)
+*>          Q is REAL array, dimension (LDQ,N)
 *>          On entry, if WANTQ = .TRUE., the orthogonal matrix Q.
 *>          On exit, the updated matrix Q.
 *>          If WANTQ = .FALSE., Q is not referenced.
diff --git a/SRC/zlalsd.f b/SRC/zlalsd.f
index 372b3822..01f3ac03 100644
--- a/SRC/zlalsd.f
+++ b/SRC/zlalsd.f
@@ -136,8 +136,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is COMPLEX*16 array, dimension at least
-*>         (N * NRHS).
+*>          WORK is COMPLEX*16 array, dimension (N * NRHS)
 *> \endverbatim
 *>
 *> \param[out] RWORK
diff --git a/SRC/zsytri2x.f b/SRC/zsytri2x.f
index a8e068db..38e3cb32 100644
--- a/SRC/zsytri2x.f
+++ b/SRC/zsytri2x.f
@@ -87,7 +87,7 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is COMPLEX*16 array, dimension (N+NNB+1,NNB+3)
+*>          WORK is COMPLEX*16 array, dimension (N+NB+1,NB+3)
 *> \endverbatim
 *>
 *> \param[in] NB
diff --git a/SRC/ztgex2.f b/SRC/ztgex2.f
index 3f484285..593bc247 100644
--- a/SRC/ztgex2.f
+++ b/SRC/ztgex2.f
@@ -102,7 +102,7 @@
 *>
 *> \param[in,out] Q
 *> \verbatim
-*>          Q is COMPLEX*16 array, dimension (LDZ,N)
+*>          Q is COMPLEX*16 array, dimension (LDQ,N)
 *>          If WANTQ = .TRUE, on entry, the unitary matrix Q. On exit,
 *>          the updated matrix Q.
 *>          Not referenced if WANTQ = .FALSE..
diff --git a/SRC/ztgexc.f b/SRC/ztgexc.f
index cb7b5229..d5d27a92 100644
--- a/SRC/ztgexc.f
+++ b/SRC/ztgexc.f
@@ -102,7 +102,7 @@
 *>
 *> \param[in,out] Q
 *> \verbatim
-*>          Q is COMPLEX*16 array, dimension (LDZ,N)
+*>          Q is COMPLEX*16 array, dimension (LDQ,N)
 *>          On entry, if WANTQ = .TRUE., the unitary matrix Q.
 *>          On exit, the updated matrix Q.
 *>          If WANTQ = .FALSE., Q is not referenced.
diff --git a/SRC/zuncsd.f b/SRC/zuncsd.f
index 18dce987..570b370b 100644
--- a/SRC/zuncsd.f
+++ b/SRC/zuncsd.f
@@ -188,7 +188,7 @@
 *>
 *> \param[out] U1
 *> \verbatim
-*>          U1 is COMPLEX*16 array, dimension (P,P)
+*>          U1 is COMPLEX*16 array, dimension (LDU1,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,M-P)
+*>          U2 is COMPLEX*16 array, dimension (LDU2,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,Q)
+*>          V1T is COMPLEX*16 array, dimension (LDV1T,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,M-Q)
+*>          V2T is COMPLEX*16 array, dimension (LDV2T,M-Q)
 *>          If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) unitary
 *>          matrix V2**H.
 *> \endverbatim