xUNCSD, xORCSD; U1, U2, V1T, V2T: duplicate dimensions of square matrices

Author: thielema

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

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

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

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