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

Diff for /rpl/lapack/lapack/zhetf2_rook.f between versions 1.5 and 1.6

version 1.5, 2018/05/29 07:18:20 version 1.6, 2023/08/07 08:39:24
Line 129 Line 129
 *> \author Univ. of Colorado Denver  *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.  *> \author NAG Ltd.
 *  *
 *> \date November 2013  
 *  
 *> \ingroup complex16HEcomputational  *> \ingroup complex16HEcomputational
 *  *
 *> \par Further Details:  *> \par Further Details:
Line 194 Line 192
 *  =====================================================================  *  =====================================================================
       SUBROUTINE ZHETF2_ROOK( UPLO, N, A, LDA, IPIV, INFO )        SUBROUTINE ZHETF2_ROOK( UPLO, N, A, LDA, IPIV, INFO )
 *  *
 *  -- LAPACK computational routine (version 3.5.0) --  *  -- LAPACK computational 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 2013  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          UPLO        CHARACTER          UPLO
Line 360 Line 357
 *  *
 *                 Case(2)  *                 Case(2)
 *                 Equivalent to testing for  *                 Equivalent to testing for
 *                 ABS( REAL( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX  *                 ABS( DBLE( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX
 *                 (used to handle NaN and Inf)  *                 (used to handle NaN and Inf)
 *  *
                   IF( .NOT.( ABS( DBLE( A( IMAX, IMAX ) ) )                    IF( .NOT.( ABS( DBLE( A( IMAX, IMAX ) ) )
Line 531 Line 528
 *                 D = |A12|  *                 D = |A12|
                   D = DLAPY2( DBLE( A( K-1, K ) ),                    D = DLAPY2( DBLE( A( K-1, K ) ),
      $                DIMAG( A( K-1, K ) ) )       $                DIMAG( A( K-1, K ) ) )
                   D11 = A( K, K ) / D                    D11 = DBLE( A( K, K ) / D )
                   D22 = A( K-1, K-1 ) / D                    D22 = DBLE( A( K-1, K-1 ) / D )
                   D12 = A( K-1, K ) / D                    D12 = A( K-1, K ) / D
                   TT = ONE / ( D11*D22-ONE )                    TT = ONE / ( D11*D22-ONE )
 *  *
Line 672 Line 669
 *  *
 *                 Case(2)  *                 Case(2)
 *                 Equivalent to testing for  *                 Equivalent to testing for
 *                 ABS( REAL( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX  *                 ABS( DBLE( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX
 *                 (used to handle NaN and Inf)  *                 (used to handle NaN and Inf)
 *  *
                   IF( .NOT.( ABS( DBLE( A( IMAX, IMAX ) ) )                    IF( .NOT.( ABS( DBLE( A( IMAX, IMAX ) ) )
Removed from v.1.5  
changed lines
  Added in v.1.6

CVSweb interface <joel.bertrand@systella.fr>