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

Diff for /rpl/lapack/lapack/zggev3.f between versions 1.1 and 1.2

version 1.1, 2015/11/26 11:44:22 version 1.2, 2016/08/27 15:27:12
Line 216 Line 216
       SUBROUTINE ZGGEV3( JOBVL, JOBVR, N, A, LDA, B, LDB, ALPHA, BETA,        SUBROUTINE ZGGEV3( JOBVL, JOBVR, N, A, LDA, B, LDB, ALPHA, BETA,
      $                   VL, LDVL, VR, LDVR, WORK, LWORK, RWORK, INFO )       $                   VL, LDVL, VR, LDVR, WORK, LWORK, RWORK, INFO )
 *  *
 *  -- LAPACK driver routine (version 3.6.0) --  *  -- LAPACK driver routine (version 3.6.1) --
 *  -- 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..--
 *     January 2015  *     January 2015
Line 340 Line 340
             LWKOPT = MAX( LWKOPT, N+INT( WORK( 1 ) ) )              LWKOPT = MAX( LWKOPT, N+INT( WORK( 1 ) ) )
             CALL ZHGEQZ( 'S', JOBVL, JOBVR, N, 1, N, A, LDA, B, LDB,              CALL ZHGEQZ( 'S', JOBVL, JOBVR, N, 1, N, A, LDA, B, LDB,
      $                   ALPHA, BETA, VL, LDVL, VR, LDVR, WORK, -1,       $                   ALPHA, BETA, VL, LDVL, VR, LDVR, WORK, -1,
      $                   WORK, IERR )       $                   RWORK, IERR )
             LWKOPT = MAX( LWKOPT, N+INT( WORK( 1 ) ) )              LWKOPT = MAX( LWKOPT, N+INT( WORK( 1 ) ) )
          ELSE           ELSE
             CALL ZGGHD3( JOBVL, JOBVR, N, 1, N, A, LDA, B, LDB, VL,              CALL ZGGHD3( JOBVL, JOBVR, N, 1, N, A, LDA, B, LDB, VL,
Line 348 Line 348
             LWKOPT = MAX( LWKOPT, N+INT( WORK( 1 ) ) )              LWKOPT = MAX( LWKOPT, N+INT( WORK( 1 ) ) )
             CALL ZHGEQZ( 'E', JOBVL, JOBVR, N, 1, N, A, LDA, B, LDB,              CALL ZHGEQZ( 'E', JOBVL, JOBVR, N, 1, N, A, LDA, B, LDB,
      $                   ALPHA, BETA, VL, LDVL, VR, LDVR, WORK, -1,       $                   ALPHA, BETA, VL, LDVL, VR, LDVR, WORK, -1,
      $                   WORK, IERR )       $                   RWORK, IERR )
             LWKOPT = MAX( LWKOPT, N+INT( WORK( 1 ) ) )              LWKOPT = MAX( LWKOPT, N+INT( WORK( 1 ) ) )
          END IF           END IF
          WORK( 1 ) = DCMPLX( LWKOPT )           WORK( 1 ) = DCMPLX( LWKOPT )
Removed from v.1.1  
changed lines
  Added in v.1.2

CVSweb interface <joel.bertrand@systella.fr>