Diff for /rpl/lapack/lapack/dlasd6.f between versions 1.14 and 1.23

version 1.14, 2012/12/14 12:30:25 version 1.23, 2023/08/07 08:38:59
Line 2 Line 2
 *  *
 *  =========== DOCUMENTATION ===========  *  =========== DOCUMENTATION ===========
 *  *
 * Online html documentation available at   * Online html documentation available at
 *            http://www.netlib.org/lapack/explore-html/   *            http://www.netlib.org/lapack/explore-html/
 *  *
 *> \htmlonly  *> \htmlonly
 *> Download DLASD6 + dependencies   *> Download DLASD6 + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dlasd6.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dlasd6.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dlasd6.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dlasd6.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dlasd6.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dlasd6.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
Line 22 Line 22
 *                          IDXQ, PERM, GIVPTR, GIVCOL, LDGCOL, GIVNUM,  *                          IDXQ, PERM, GIVPTR, GIVCOL, LDGCOL, GIVNUM,
 *                          LDGNUM, POLES, DIFL, DIFR, Z, K, C, S, WORK,  *                          LDGNUM, POLES, DIFL, DIFR, Z, K, C, S, WORK,
 *                          IWORK, INFO )  *                          IWORK, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       INTEGER            GIVPTR, ICOMPQ, INFO, K, LDGCOL, LDGNUM, NL,  *       INTEGER            GIVPTR, ICOMPQ, INFO, K, LDGCOL, LDGNUM, NL,
 *      $                   NR, SQRE  *      $                   NR, SQRE
Line 35 Line 35
 *      $                   GIVNUM( LDGNUM, * ), POLES( LDGNUM, * ),  *      $                   GIVNUM( LDGNUM, * ), POLES( LDGNUM, * ),
 *      $                   VF( * ), VL( * ), WORK( * ), Z( * )  *      $                   VF( * ), VL( * ), WORK( * ), Z( * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 74 Line 74
 *>  *>
 *>       The first stage consists of deflating the size of the problem  *>       The first stage consists of deflating the size of the problem
 *>       when there are multiple singular values or if there is a zero  *>       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  *>       secular equation problem is reduced by one. This stage is
 *>       performed by the routine DLASD7.  *>       performed by the routine DLASD7.
 *>  *>
Line 164 Line 164
 *>         row.  *>         row.
 *> \endverbatim  *> \endverbatim
 *>  *>
 *> \param[out] IDXQ  *> \param[in,out] IDXQ
 *> \verbatim  *> \verbatim
 *>          IDXQ is INTEGER array, dimension ( N )  *>          IDXQ is INTEGER array, dimension ( N )
 *>         This contains the permutation which will reintegrate the  *>         This contains the permutation which will reintegrate the
Line 232 Line 232
 *> \param[out] DIFR  *> \param[out] DIFR
 *> \verbatim  *> \verbatim
 *>          DIFR is DOUBLE PRECISION array,  *>          DIFR is DOUBLE PRECISION array,
 *>                  dimension ( LDGNUM, 2 ) if ICOMPQ = 1 and  *>                   dimension ( LDDIFR, 2 ) if ICOMPQ = 1 and
 *>                  dimension ( N ) if ICOMPQ = 0.  *>                   dimension ( K ) if ICOMPQ = 0.
 *>         On exit, DIFR(I, 1) is the distance between I-th updated  *>          On exit, DIFR(I,1) = D(I) - DSIGMA(I+1), DIFR(K,1) is not
 *>         (undeflated) singular value and the I+1-th (undeflated) old  *>          defined and will not be referenced.
 *>         singular value.  
 *>  *>
 *>         If ICOMPQ = 1, DIFR(1:K,2) is an array containing the  *>          If ICOMPQ = 1, DIFR(1:K,2) is an array containing the
 *>         normalizing factors for the right singular vector matrix.  *>          normalizing factors for the right singular vector matrix.
 *>  *>
 *>         See DLASD8 for details on DIFL and DIFR.  *>         See DLASD8 for details on DIFL and DIFR.
 *> \endverbatim  *> \endverbatim
Line 293 Line 292
 *  Authors:  *  Authors:
 *  ========  *  ========
 *  *
 *> \author Univ. of Tennessee   *> \author Univ. of Tennessee
 *> \author Univ. of California Berkeley   *> \author Univ. of California Berkeley
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.   *> \author NAG Ltd.
 *  
 *> \date September 2012  
 *  *
 *> \ingroup auxOTHERauxiliary  *> \ingroup OTHERauxiliary
 *  *
 *> \par Contributors:  *> \par Contributors:
 *  ==================  *  ==================
Line 314 Line 311
      $                   LDGNUM, POLES, DIFL, DIFR, Z, K, C, S, WORK,       $                   LDGNUM, POLES, DIFL, DIFR, Z, K, C, S, WORK,
      $                   IWORK, INFO )       $                   IWORK, INFO )
 *  *
 *  -- LAPACK auxiliary routine (version 3.4.2) --  *  -- LAPACK auxiliary routine --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --  *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--  *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
 *     September 2012  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            GIVPTR, ICOMPQ, INFO, K, LDGCOL, LDGNUM, NL,        INTEGER            GIVPTR, ICOMPQ, INFO, K, LDGCOL, LDGNUM, NL,
Line 414 Line 410
       CALL DLASD8( ICOMPQ, K, D, Z, VF, VL, DIFL, DIFR, LDGNUM,        CALL DLASD8( ICOMPQ, K, D, Z, VF, VL, DIFL, DIFR, LDGNUM,
      $             WORK( ISIGMA ), WORK( IW ), INFO )       $             WORK( ISIGMA ), WORK( IW ), INFO )
 *  *
 *     Handle error returned  *     Report the possible convergence failure.
 *  *
       IF( INFO.NE.0 ) THEN        IF( INFO.NE.0 ) THEN
          CALL XERBLA( 'DLASD8', -INFO )  
          RETURN           RETURN
       END IF        END IF
 *  *

Removed from v.1.14  
changed lines
  Added in v.1.23


CVSweb interface <joel.bertrand@systella.fr>