Diff for /rpl/lapack/lapack/dgeesx.f between versions 1.9 and 1.20

version 1.9, 2011/11/21 20:42:50 version 1.20, 2023/08/07 08:38:48
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 DGEESX + dependencies   *> Download DGEESX + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dgeesx.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dgeesx.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dgeesx.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dgeesx.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dgeesx.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dgeesx.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
Line 21 Line 21
 *       SUBROUTINE DGEESX( JOBVS, SORT, SELECT, SENSE, N, A, LDA, SDIM,  *       SUBROUTINE DGEESX( JOBVS, SORT, SELECT, SENSE, N, A, LDA, SDIM,
 *                          WR, WI, VS, LDVS, RCONDE, RCONDV, WORK, LWORK,  *                          WR, WI, VS, LDVS, RCONDE, RCONDV, WORK, LWORK,
 *                          IWORK, LIWORK, BWORK, INFO )  *                          IWORK, LIWORK, BWORK, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       CHARACTER          JOBVS, SENSE, SORT  *       CHARACTER          JOBVS, SENSE, SORT
 *       INTEGER            INFO, LDA, LDVS, LIWORK, LWORK, N, SDIM  *       INTEGER            INFO, LDA, LDVS, LIWORK, LWORK, N, SDIM
Line 37 Line 37
 *       LOGICAL            SELECT  *       LOGICAL            SELECT
 *       EXTERNAL           SELECT  *       EXTERNAL           SELECT
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 90 Line 90
 *>  *>
 *> \param[in] SELECT  *> \param[in] SELECT
 *> \verbatim  *> \verbatim
 *>          SELECT is procedure) LOGICAL FUNCTION of two DOUBLE PRECISION arguments  *>          SELECT is a LOGICAL FUNCTION of two DOUBLE PRECISION arguments
 *>          SELECT must be declared EXTERNAL in the calling subroutine.  *>          SELECT must be declared EXTERNAL in the calling subroutine.
 *>          If SORT = 'S', SELECT is used to select eigenvalues to sort  *>          If SORT = 'S', SELECT is used to select eigenvalues to sort
 *>          to the top left of the Schur form.  *>          to the top left of the Schur form.
Line 267 Line 267
 *  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 November 2011  
 *  *
 *> \ingroup doubleGEeigen  *> \ingroup doubleGEeigen
 *  *
Line 281 Line 279
      $                   WR, WI, VS, LDVS, RCONDE, RCONDV, WORK, LWORK,       $                   WR, WI, VS, LDVS, RCONDE, RCONDV, WORK, LWORK,
      $                   IWORK, LIWORK, BWORK, INFO )       $                   IWORK, LIWORK, BWORK, INFO )
 *  *
 *  -- LAPACK driver routine (version 3.4.0) --  *  -- LAPACK driver 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..--
 *     November 2011  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          JOBVS, SENSE, SORT        CHARACTER          JOBVS, SENSE, SORT
Line 385 Line 382
 *  *
             CALL DHSEQR( 'S', JOBVS, N, 1, N, A, LDA, WR, WI, VS, LDVS,              CALL DHSEQR( 'S', JOBVS, N, 1, N, A, LDA, WR, WI, VS, LDVS,
      $             WORK, -1, IEVAL )       $             WORK, -1, IEVAL )
             HSWORK = WORK( 1 )              HSWORK = INT( WORK( 1 ) )
 *  *
             IF( .NOT.WANTVS ) THEN              IF( .NOT.WANTVS ) THEN
                MAXWRK = MAX( MAXWRK, N + HSWORK )                 MAXWRK = MAX( MAXWRK, N + HSWORK )
Line 583 Line 580
                      IF( N.GT.I+1 )                       IF( N.GT.I+1 )
      $                  CALL DSWAP( N-I-1, A( I, I+2 ), LDA,       $                  CALL DSWAP( N-I-1, A( I, I+2 ), LDA,
      $                              A( I+1, I+2 ), LDA )       $                              A( I+1, I+2 ), LDA )
                      CALL DSWAP( N, VS( 1, I ), 1, VS( 1, I+1 ), 1 )                       IF( WANTVS ) THEN
                          CALL DSWAP( N, VS( 1, I ), 1, VS( 1, I+1 ), 1 )
                        END IF
                      A( I, I+1 ) = A( I+1, I )                       A( I, I+1 ) = A( I+1, I )
                      A( I+1, I ) = ZERO                       A( I+1, I ) = ZERO
                   END IF                    END IF

Removed from v.1.9  
changed lines
  Added in v.1.20


CVSweb interface <joel.bertrand@systella.fr>