--- rpl/lapack/lapack/zbbcsd.f 2016/08/27 15:34:43 1.12
+++ rpl/lapack/lapack/zbbcsd.f 2023/08/07 08:39:15 1.17
@@ -2,18 +2,18 @@
*
* =========== DOCUMENTATION ===========
*
-* Online html documentation available at
-* http://www.netlib.org/lapack/explore-html/
+* Online html documentation available at
+* http://www.netlib.org/lapack/explore-html/
*
*> \htmlonly
-*> Download ZBBCSD + dependencies
-*>
-*> [TGZ]
-*>
-*> [ZIP]
-*>
+*> Download ZBBCSD + dependencies
+*>
+*> [TGZ]
+*>
+*> [ZIP]
+*>
*> [TXT]
-*> \endhtmlonly
+*> \endhtmlonly
*
* Definition:
* ===========
@@ -22,7 +22,7 @@
* THETA, PHI, U1, LDU1, U2, LDU2, V1T, LDV1T,
* V2T, LDV2T, B11D, B11E, B12D, B12E, B21D, B21E,
* B22D, B22E, RWORK, LRWORK, INFO )
-*
+*
* .. Scalar Arguments ..
* CHARACTER JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS
* INTEGER INFO, LDU1, LDU2, LDV1T, LDV2T, LRWORK, M, P, Q
@@ -34,7 +34,7 @@
* COMPLEX*16 U1( LDU1, * ), U2( LDU2, * ), V1T( LDV1T, * ),
* $ V2T( LDV2T, * )
* ..
-*
+*
*
*> \par Purpose:
* =============
@@ -190,7 +190,7 @@
*>
*> \param[in,out] V2T
*> \verbatim
-*> V2T is COMPLEX*16 array, dimenison (LDV2T,M-Q)
+*> V2T is COMPLEX*16 array, dimension (LDV2T,M-Q)
*> On entry, an (M-Q)-by-(M-Q) matrix. On exit, V2T is
*> premultiplied by the conjugate transpose of the right
*> singular vector matrix common to [ B12 0 0 ; 0 -I 0 ] and
@@ -273,8 +273,8 @@
*>
*> \param[out] RWORK
*> \verbatim
-*> RWORK is DOUBLE PRECISION array, dimension (MAX(1,LWORK))
-*> On exit, if INFO = 0, WORK(1) returns the optimal LWORK.
+*> RWORK is DOUBLE PRECISION array, dimension (MAX(1,LRWORK))
+*> On exit, if INFO = 0, RWORK(1) returns the optimal LRWORK.
*> \endverbatim
*>
*> \param[in] LRWORK
@@ -317,12 +317,10 @@
* Authors:
* ========
*
-*> \author Univ. of Tennessee
-*> \author Univ. of California Berkeley
-*> \author Univ. of Colorado Denver
-*> \author NAG Ltd.
-*
-*> \date June 2016
+*> \author Univ. of Tennessee
+*> \author Univ. of California Berkeley
+*> \author Univ. of Colorado Denver
+*> \author NAG Ltd.
*
*> \ingroup complex16OTHERcomputational
*
@@ -332,10 +330,9 @@
$ V2T, LDV2T, B11D, B11E, B12D, B12E, B21D, B21E,
$ B22D, B22E, RWORK, LRWORK, INFO )
*
-* -- LAPACK computational routine (version 3.6.1) --
+* -- LAPACK computational routine --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* June 2016
*
* .. Scalar Arguments ..
CHARACTER JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS
@@ -354,12 +351,13 @@
* .. Parameters ..
INTEGER MAXITR
PARAMETER ( MAXITR = 6 )
- DOUBLE PRECISION HUNDRED, MEIGHTH, ONE, PIOVER2, TEN, ZERO
+ DOUBLE PRECISION HUNDRED, MEIGHTH, ONE, TEN, ZERO
PARAMETER ( HUNDRED = 100.0D0, MEIGHTH = -0.125D0,
- $ ONE = 1.0D0, PIOVER2 = 1.57079632679489662D0,
- $ TEN = 10.0D0, ZERO = 0.0D0 )
+ $ ONE = 1.0D0, TEN = 10.0D0, ZERO = 0.0D0 )
COMPLEX*16 NEGONECOMPLEX
PARAMETER ( NEGONECOMPLEX = (-1.0D0,0.0D0) )
+ DOUBLE PRECISION PIOVER2
+ PARAMETER ( PIOVER2 = 1.57079632679489661923132169163975144210D0 )
* ..
* .. Local Scalars ..
LOGICAL COLMAJOR, LQUERY, RESTART11, RESTART12,