rpl/lapack/lapack/dlanst.f - diff

Return to dlanst.f CVS log

Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/dlanst.f between versions 1.1.1.1 and 1.12

-version 1.1.1.1, 2010/01/26 15:22:45
+version 1.12, 2012/12/14 14:22:33
  Line 1
+ *> \brief \b DLANST returns the value of the 1-norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a real symmetric tridiagonal matrix.
+ *
+ *  =========== DOCUMENTATION ===========
+ *
+ * Online html documentation available at
+ *            http://www.netlib.org/lapack/explore-html/
+ *
+ *> \htmlonly
+ *> Download DLANST + dependencies
+ *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dlanst.f">
+ *> [TGZ]</a>
+ *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dlanst.f">
+ *> [ZIP]</a>
+ *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dlanst.f">
+ *> [TXT]</a>
+ *> \endhtmlonly
+ *
+ *  Definition:
+ *  ===========
+ *
+ *       DOUBLE PRECISION FUNCTION DLANST( NORM, N, D, E )
+ *
+ *       .. Scalar Arguments ..
+ *       CHARACTER          NORM
+ *       INTEGER            N
+ *       ..
+ *       .. Array Arguments ..
+ *       DOUBLE PRECISION   D( * ), E( * )
+ *       ..
+ *
+ *
+ *> \par Purpose:
+ *  =============
+ *>
+ *> \verbatim
+ *>
+ *> DLANST  returns the value of the one norm,  or the Frobenius norm, or
+ *> the  infinity norm,  or the  element of  largest absolute value  of a
+ *> real symmetric tridiagonal matrix A.
+ *> \endverbatim
+ *>
+ *> \return DLANST
+ *> \verbatim
+ *>
+ *>    DLANST = ( max(abs(A(i,j))), NORM = 'M' or 'm'
+ *>             (
+ *>             ( norm1(A),         NORM = '1', 'O' or 'o'
+ *>             (
+ *>             ( normI(A),         NORM = 'I' or 'i'
+ *>             (
+ *>             ( normF(A),         NORM = 'F', 'f', 'E' or 'e'
+ *>
+ *> where  norm1  denotes the  one norm of a matrix (maximum column sum),
+ *> normI  denotes the  infinity norm  of a matrix  (maximum row sum) and
+ *> normF  denotes the  Frobenius norm of a matrix (square root of sum of
+ *> squares).  Note that  max(abs(A(i,j)))  is not a consistent matrix norm.
+ *> \endverbatim
+ *
+ *  Arguments:
+ *  ==========
+ *
+ *> \param[in] NORM
+ *> \verbatim
+ *>          NORM is CHARACTER*1
+ *>          Specifies the value to be returned in DLANST as described
+ *>          above.
+ *> \endverbatim
+ *>
+ *> \param[in] N
+ *> \verbatim
+ *>          N is INTEGER
+ *>          The order of the matrix A.  N >= 0.  When N = 0, DLANST is
+ *>          set to zero.
+ *> \endverbatim
+ *>
+ *> \param[in] D
+ *> \verbatim
+ *>          D is DOUBLE PRECISION array, dimension (N)
+ *>          The diagonal elements of A.
+ *> \endverbatim
+ *>
+ *> \param[in] E
+ *> \verbatim
+ *>          E is DOUBLE PRECISION array, dimension (N-1)
+ *>          The (n-1) sub-diagonal or super-diagonal elements of A.
+ *> \endverbatim
+ *
+ *  Authors:
+ *  ========
+ *
+ *> \author Univ. of Tennessee
+ *> \author Univ. of California Berkeley
+ *> \author Univ. of Colorado Denver
+ *> \author NAG Ltd.
+ *
+ *> \date September 2012
+ *
+ *> \ingroup auxOTHERauxiliary
+ *
+ *  =====================================================================
        DOUBLE PRECISION FUNCTION DLANST( NORM, N, D, E )
  *
- *  -- LAPACK auxiliary routine (version 3.2) --
+ *  -- LAPACK auxiliary routine (version 3.4.2) --
  *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
  *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
- *     November 2006
+ *     September 2012
  *
  *     .. Scalar Arguments ..
        CHARACTER          NORM
- Line 13
+ Line 113
        DOUBLE PRECISION   D( * ), E( * )
  *     ..
  *
- *  Purpose
- *  =======
- *
- *  DLANST  returns the value of the one norm,  or the Frobenius norm, or
- *  the  infinity norm,  or the  element of  largest absolute value  of a
- *  real symmetric tridiagonal matrix A.
- *
- *  Description
- *  ===========
- *
- *  DLANST returns the value
- *
- *     DLANST = ( max(abs(A(i,j))), NORM = 'M' or 'm'
- *              (
- *              ( norm1(A),         NORM = '1', 'O' or 'o'
- *              (
- *              ( normI(A),         NORM = 'I' or 'i'
- *              (
- *              ( normF(A),         NORM = 'F', 'f', 'E' or 'e'
- *
- *  where  norm1  denotes the  one norm of a matrix (maximum column sum),
- *  normI  denotes the  infinity norm  of a matrix  (maximum row sum) and
- *  normF  denotes the  Frobenius norm of a matrix (square root of sum of
- *  squares).  Note that  max(abs(A(i,j)))  is not a consistent matrix norm.
- *
- *  Arguments
- *  =========
- *
- *  NORM    (input) CHARACTER*1
- *          Specifies the value to be returned in DLANST as described
- *          above.
- *
- *  N       (input) INTEGER
- *          The order of the matrix A.  N >= 0.  When N = 0, DLANST is
- *          set to zero.
- *
- *  D       (input) DOUBLE PRECISION array, dimension (N)
- *          The diagonal elements of A.
- *
- *  E       (input) DOUBLE PRECISION array, dimension (N-1)
- *          The (n-1) sub-diagonal or super-diagonal elements of A.
- *
  *  =====================================================================
  *
  *     .. Parameters ..
- Line 66
+ Line 124
        DOUBLE PRECISION   ANORM, SCALE, SUM
  *     ..
  *     .. External Functions ..
-       LOGICAL            LSAME
+       LOGICAL            LSAME, DISNAN
-       EXTERNAL           LSAME
+       EXTERNAL           LSAME, DISNAN
  *     ..
  *     .. External Subroutines ..
        EXTERNAL           DLASSQ
  *     ..
  *     .. Intrinsic Functions ..
-       INTRINSIC          ABS, MAX, SQRT
+       INTRINSIC          ABS, SQRT
  *     ..
  *     .. Executable Statements ..
  *
- Line 85
+ Line 143
  *
           ANORM = ABS( D( N ) )
           DO 10 I = 1, N - 1
-             ANORM = MAX( ANORM, ABS( D( I ) ) )
+             SUM = ABS( D( I ) )
-             ANORM = MAX( ANORM, ABS( E( I ) ) )
+             IF( ANORM .LT. SUM .OR. DISNAN( SUM ) ) ANORM = SUM
+             SUM = ABS( E( I ) )
+             IF( ANORM .LT. SUM .OR. DISNAN( SUM ) ) ANORM = SUM
    CONTINUE
        ELSE IF( LSAME( NORM, 'O' ) .OR. NORM.EQ.'1' .OR.
       $         LSAME( NORM, 'I' ) ) THEN
- Line 96
+ Line 156
           IF( N.EQ.1 ) THEN
              ANORM = ABS( D( 1 ) )
           ELSE
-             ANORM = MAX( ABS( D( 1 ) )+ABS( E( 1 ) ),
+             ANORM = ABS( D( 1 ) )+ABS( E( 1 ) )
-      $              ABS( E( N-1 ) )+ABS( D( N ) ) )
+             SUM = ABS( E( N-1 ) )+ABS( D( N ) )
+             IF( ANORM .LT. SUM .OR. DISNAN( SUM ) ) ANORM = SUM
              DO 20 I = 2, N - 1
-                ANORM = MAX( ANORM, ABS( D( I ) )+ABS( E( I ) )+
+                SUM = ABS( D( I ) )+ABS( E( I ) )+ABS( E( I-1 ) )
-      $                 ABS( E( I-1 ) ) )
+                IF( ANORM .LT. SUM .OR. DISNAN( SUM ) ) ANORM = SUM
       CONTINUE
           END IF
        ELSE IF( ( LSAME( NORM, 'F' ) ) .OR. ( LSAME( NORM, 'E' ) ) ) THEN

CVSweb interface <joel.bertrand@systella.fr>

Removed from v.1.1.1.1
changed lines
	Added in v.1.12