From 1b029a0edc2336b0d4e71bc9ed8e2b36cee331ff Mon Sep 17 00:00:00 2001
From: Henning Thielemann <git@henning-thielemann.de>
Date: Tue, 6 Jun 2017 19:26:29 +0200
Subject: xUNCSD, xORCSD: use LDx for leading dimensions, although the arrays
 contain square matrices

This fixes commit 366f9b3541705231cfbb1724b221b8857ac768a8.
---
 SRC/cuncsd.f | 8 ++++----
 SRC/dorcsd.f | 8 ++++----
 SRC/sorcsd.f | 8 ++++----
 SRC/zuncsd.f | 8 ++++----
 4 files changed, 16 insertions(+), 16 deletions(-)

(limited to 'SRC')

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/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/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/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
-- 
cgit v1.2.3