SORCSD: fix documentation on matrix dimensions

christoph-conrads · christoph-conrads · commit e86cbeb823fb · 2022-05-29T18:30:09.000Z
diff --git a/SRC/sorbdb2.f b/SRC/sorbdb2.f
@@ -122,14 +122,14 @@
 *>
 *> \param[out] TAUP1
 *> \verbatim
-*>          TAUP1 is REAL array, dimension (P)
+*>          TAUP1 is REAL array, dimension (P-1)
 *>           The scalar factors of the elementary reflectors that define
 *>           P1.
 *> \endverbatim
 *>
 *> \param[out] TAUP2
 *> \verbatim
-*>          TAUP2 is REAL array, dimension (M-P)
+*>          TAUP2 is REAL array, dimension (Q)
 *>           The scalar factors of the elementary reflectors that define
 *>           P2.
 *> \endverbatim
diff --git a/SRC/sorbdb4.f b/SRC/sorbdb4.f
@@ -131,7 +131,7 @@
 *>
 *> \param[out] TAUP2
 *> \verbatim
-*>          TAUP2 is REAL array, dimension (M-P)
+*>          TAUP2 is REAL array, dimension (M-Q)
 *>           The scalar factors of the elementary reflectors that define
 *>           P2.
 *> \endverbatim
diff --git a/SRC/sorcsd2by1.f b/SRC/sorcsd2by1.f
@@ -672,7 +672,7 @@ SUBROUTINE SORCSD2BY1( JOBU1, JOBU2, JOBV1T, M, P, Q, X11, LDX11,
 *        Accumulate Householder reflectors
 *
          IF( WANTU2 .AND. M-P .GT. 0 ) THEN
-            CALL SCOPY( M-P, WORK(IORBDB+P), 1, U2, 1 )
+            CALL SCOPY( M-Q, WORK(IORBDB+P), 1, U2, 1 )
          END IF
          IF( WANTU1 .AND. P .GT. 0 ) THEN
             CALL SCOPY( P, WORK(IORBDB), 1, U1, 1 )

Original file line number	Diff line number	Diff line change
`@@ -672,7 +672,7 @@ SUBROUTINE SORCSD2BY1( JOBU1, JOBU2, JOBV1T, M, P, Q, X11, LDX11,`
`672`	`672`	`* Accumulate Householder reflectors`
`673`	`673`	`*`
`674`	`674`	`IF( WANTU2 .AND. M-P .GT. 0 ) THEN`
`675`		`- CALL SCOPY( M-P, WORK(IORBDB+P), 1, U2, 1 )`
	`675`	`+ CALL SCOPY( M-Q, WORK(IORBDB+P), 1, U2, 1 )`
`676`	`676`	`END IF`
`677`	`677`	`IF( WANTU1 .AND. P .GT. 0 ) THEN`
`678`	`678`	`CALL SCOPY( P, WORK(IORBDB), 1, U1, 1 )`