[BACK]Return to dgeesx.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / lapack

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

version 1.16, 2017/06/17 10:53:47 version 1.20, 2023/08/07 08:38:48
Line 272 Line 272
 *> \author Univ. of Colorado Denver  *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.  *> \author NAG Ltd.
 *  *
 *> \date June 2016  
 *  
 *> \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.7.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..--
 *     June 2016  
 *  *
 *     .. 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.16  
changed lines
  Added in v.1.20

CVSweb interface <joel.bertrand@systella.fr>