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

Diff for /rpl/lapack/lapack/zlangt.f between versions 1.10 and 1.11

version 1.10, 2012/08/22 09:48:35 version 1.11, 2012/12/14 12:30:31
Line 1 Line 1
 *> \brief \b ZLANGT  *> \brief \b ZLANGT returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value of any element of a general tridiagonal matrix.
 *  *
 *  =========== DOCUMENTATION ===========  *  =========== DOCUMENTATION ===========
 *  *
Line 99 Line 99
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver 
 *> \author NAG Ltd.   *> \author NAG Ltd. 
 *  *
 *> \date November 2011  *> \date September 2012
 *  *
 *> \ingroup complex16OTHERauxiliary  *> \ingroup complex16OTHERauxiliary
 *  *
 *  =====================================================================  *  =====================================================================
       DOUBLE PRECISION FUNCTION ZLANGT( NORM, N, DL, D, DU )        DOUBLE PRECISION FUNCTION ZLANGT( NORM, N, DL, D, DU )
 *  *
 *  -- LAPACK auxiliary routine (version 3.4.0) --  *  -- LAPACK auxiliary routine (version 3.4.2) --
 *  -- 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  *     September 2012
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          NORM        CHARACTER          NORM
Line 127 Line 127
 *     ..  *     ..
 *     .. Local Scalars ..  *     .. Local Scalars ..
       INTEGER            I        INTEGER            I
       DOUBLE PRECISION   ANORM, SCALE, SUM        DOUBLE PRECISION   ANORM, SCALE, SUM, TEMP
 *     ..  *     ..
 *     .. External Functions ..  *     .. External Functions ..
       LOGICAL            LSAME        LOGICAL            LSAME, DISNAN
       EXTERNAL           LSAME        EXTERNAL           LSAME, DISNAN
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           ZLASSQ        EXTERNAL           ZLASSQ
 *     ..  *     ..
 *     .. Intrinsic Functions ..  *     .. Intrinsic Functions ..
       INTRINSIC          ABS, MAX, SQRT        INTRINSIC          ABS, SQRT
 *     ..  *     ..
 *     .. Executable Statements ..  *     .. Executable Statements ..
 *  *
Line 149 Line 149
 *  *
          ANORM = ABS( D( N ) )           ANORM = ABS( D( N ) )
          DO 10 I = 1, N - 1           DO 10 I = 1, N - 1
             ANORM = MAX( ANORM, ABS( DL( I ) ) )              IF( ANORM.LT.ABS( DL( I ) ) .OR. DISNAN( ABS( DL( I ) ) ) ) 
             ANORM = MAX( ANORM, ABS( D( I ) ) )       $           ANORM = ABS(DL(I))
             ANORM = MAX( ANORM, ABS( DU( I ) ) )              IF( ANORM.LT.ABS( D( I ) ) .OR. DISNAN( ABS( D( I ) ) ) ) 
        $           ANORM = ABS(D(I))
               IF( ANORM.LT.ABS( DU( I ) ) .OR. DISNAN (ABS( DU( I ) ) ) ) 
        $           ANORM = ABS(DU(I))
    10    CONTINUE     10    CONTINUE
       ELSE IF( LSAME( NORM, 'O' ) .OR. NORM.EQ.'1' ) THEN        ELSE IF( LSAME( NORM, 'O' ) .OR. NORM.EQ.'1' ) THEN
 *  *
Line 160 Line 163
          IF( N.EQ.1 ) THEN           IF( N.EQ.1 ) THEN
             ANORM = ABS( D( 1 ) )              ANORM = ABS( D( 1 ) )
          ELSE           ELSE
             ANORM = MAX( ABS( D( 1 ) )+ABS( DL( 1 ) ),              ANORM = ABS( D( 1 ) )+ABS( DL( 1 ) )
      $              ABS( D( N ) )+ABS( DU( N-1 ) ) )              TEMP = ABS( D( N ) )+ABS( DU( N-1 ) ) 
               IF( ANORM .LT. TEMP .OR. DISNAN( TEMP ) ) ANORM = TEMP
             DO 20 I = 2, N - 1              DO 20 I = 2, N - 1
                ANORM = MAX( ANORM, ABS( D( I ) )+ABS( DL( I ) )+                 TEMP = ABS( D( I ) )+ABS( DL( I ) )+ABS( DU( I-1 ) )
      $                 ABS( DU( I-1 ) ) )                 IF( ANORM .LT. TEMP .OR. DISNAN( TEMP ) ) ANORM = TEMP
    20       CONTINUE     20       CONTINUE
          END IF           END IF
       ELSE IF( LSAME( NORM, 'I' ) ) THEN        ELSE IF( LSAME( NORM, 'I' ) ) THEN
Line 174 Line 178
          IF( N.EQ.1 ) THEN           IF( N.EQ.1 ) THEN
             ANORM = ABS( D( 1 ) )              ANORM = ABS( D( 1 ) )
          ELSE           ELSE
             ANORM = MAX( ABS( D( 1 ) )+ABS( DU( 1 ) ),              ANORM = ABS( D( 1 ) )+ABS( DU( 1 ) )
      $              ABS( D( N ) )+ABS( DL( N-1 ) ) )              TEMP = ABS( D( N ) )+ABS( DL( N-1 ) )
               IF( ANORM .LT. TEMP .OR. DISNAN( TEMP ) ) ANORM = TEMP
             DO 30 I = 2, N - 1              DO 30 I = 2, N - 1
                ANORM = MAX( ANORM, ABS( D( I ) )+ABS( DU( I ) )+                 TEMP = ABS( D( I ) )+ABS( DU( I ) )+ABS( DL( I-1 ) )
      $                 ABS( DL( I-1 ) ) )                 IF( ANORM .LT. TEMP .OR. DISNAN( TEMP ) ) ANORM = TEMP
    30       CONTINUE     30       CONTINUE
          END IF           END IF
       ELSE IF( ( LSAME( NORM, 'F' ) ) .OR. ( LSAME( NORM, 'E' ) ) ) THEN        ELSE IF( ( LSAME( NORM, 'F' ) ) .OR. ( LSAME( NORM, 'E' ) ) ) THEN
Removed from v.1.10  
changed lines
  Added in v.1.11

CVSweb interface <joel.bertrand@systella.fr>