--- rpl/lapack/lapack/dlasd6.f 2011/11/21 20:42:58 1.11 +++ rpl/lapack/lapack/dlasd6.f 2016/08/27 15:34:31 1.19 @@ -1,4 +1,4 @@ -*> \brief \b DLASD6 +*> \brief \b DLASD6 computes the SVD of an updated upper bidiagonal matrix obtained by merging two smaller ones by appending a row. Used by sbdsdc. * * =========== DOCUMENTATION =========== * @@ -74,7 +74,7 @@ *> *> The first stage consists of deflating the size of the problem *> when there are multiple singular values or if there is a zero -*> in the Z vector. For each such occurence the dimension of the +*> in the Z vector. For each such occurrence the dimension of the *> secular equation problem is reduced by one. This stage is *> performed by the routine DLASD7. *> @@ -164,7 +164,7 @@ *> row. *> \endverbatim *> -*> \param[out] IDXQ +*> \param[in,out] IDXQ *> \verbatim *> IDXQ is INTEGER array, dimension ( N ) *> This contains the permutation which will reintegrate the @@ -232,14 +232,13 @@ *> \param[out] DIFR *> \verbatim *> DIFR is DOUBLE PRECISION array, -*> dimension ( LDGNUM, 2 ) if ICOMPQ = 1 and -*> dimension ( N ) if ICOMPQ = 0. -*> On exit, DIFR(I, 1) is the distance between I-th updated -*> (undeflated) singular value and the I+1-th (undeflated) old -*> singular value. +*> dimension ( LDDIFR, 2 ) if ICOMPQ = 1 and +*> dimension ( K ) if ICOMPQ = 0. +*> On exit, DIFR(I,1) = D(I) - DSIGMA(I+1), DIFR(K,1) is not +*> defined and will not be referenced. *> -*> If ICOMPQ = 1, DIFR(1:K,2) is an array containing the -*> normalizing factors for the right singular vector matrix. +*> If ICOMPQ = 1, DIFR(1:K,2) is an array containing the +*> normalizing factors for the right singular vector matrix. *> *> See DLASD8 for details on DIFL and DIFR. *> \endverbatim @@ -298,7 +297,7 @@ *> \author Univ. of Colorado Denver *> \author NAG Ltd. * -*> \date November 2011 +*> \date June 2016 * *> \ingroup auxOTHERauxiliary * @@ -314,10 +313,10 @@ $ LDGNUM, POLES, DIFL, DIFR, Z, K, C, S, WORK, $ IWORK, INFO ) * -* -- LAPACK auxiliary routine (version 3.4.0) -- +* -- LAPACK auxiliary routine (version 3.6.1) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* November 2011 +* June 2016 * * .. Scalar Arguments .. INTEGER GIVPTR, ICOMPQ, INFO, K, LDGCOL, LDGNUM, NL, @@ -414,10 +413,9 @@ CALL DLASD8( ICOMPQ, K, D, Z, VF, VL, DIFL, DIFR, LDGNUM, $ WORK( ISIGMA ), WORK( IW ), INFO ) * -* Handle error returned +* Report the possible convergence failure. * IF( INFO.NE.0 ) THEN - CALL XERBLA( 'DLASD8', -INFO ) RETURN END IF *