--- rpl/lapack/lapack/dlansf.f	2012/07/31 11:06:36	1.9
+++ rpl/lapack/lapack/dlansf.f	2017/06/17 10:53:54	1.15
@@ -1,25 +1,25 @@
-*> \brief \b DLANSF
+*> \brief \b DLANSF returns the value of the 1-norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a symmetric matrix in RFP format.
 *
 *  =========== DOCUMENTATION ===========
 *
-* Online html documentation available at 
-*            http://www.netlib.org/lapack/explore-html/ 
+* Online html documentation available at
+*            http://www.netlib.org/lapack/explore-html/
 *
 *> \htmlonly
-*> Download DLANSF + dependencies 
-*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dlansf.f"> 
-*> [TGZ]</a> 
-*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dlansf.f"> 
-*> [ZIP]</a> 
-*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dlansf.f"> 
+*> Download DLANSF + dependencies
+*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dlansf.f">
+*> [TGZ]</a>
+*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dlansf.f">
+*> [ZIP]</a>
+*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dlansf.f">
 *> [TXT]</a>
-*> \endhtmlonly 
+*> \endhtmlonly
 *
 *  Definition:
 *  ===========
 *
 *       DOUBLE PRECISION FUNCTION DLANSF( NORM, TRANSR, UPLO, N, A, WORK )
-* 
+*
 *       .. Scalar Arguments ..
 *       CHARACTER          NORM, TRANSR, UPLO
 *       INTEGER            N
@@ -27,7 +27,7 @@
 *       .. Array Arguments ..
 *       DOUBLE PRECISION   A( 0: * ), WORK( 0: * )
 *       ..
-*  
+*
 *
 *> \par Purpose:
 *  =============
@@ -110,12 +110,12 @@
 *  Authors:
 *  ========
 *
-*> \author Univ. of Tennessee 
-*> \author Univ. of California Berkeley 
-*> \author Univ. of Colorado Denver 
-*> \author NAG Ltd. 
+*> \author Univ. of Tennessee
+*> \author Univ. of California Berkeley
+*> \author Univ. of Colorado Denver
+*> \author NAG Ltd.
 *
-*> \date November 2011
+*> \date December 2016
 *
 *> \ingroup doubleOTHERcomputational
 *
@@ -209,10 +209,10 @@
 *  =====================================================================
       DOUBLE PRECISION FUNCTION DLANSF( NORM, TRANSR, UPLO, N, A, WORK )
 *
-*  -- LAPACK computational routine (version 3.4.0) --
+*  -- LAPACK computational routine (version 3.7.0) --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-*     November 2011
+*     December 2016
 *
 *     .. Scalar Arguments ..
       CHARACTER          NORM, TRANSR, UPLO
@@ -230,12 +230,11 @@
 *     ..
 *     .. Local Scalars ..
       INTEGER            I, J, IFM, ILU, NOE, N1, K, L, LDA
-      DOUBLE PRECISION   SCALE, S, VALUE, AA
+      DOUBLE PRECISION   SCALE, S, VALUE, AA, TEMP
 *     ..
 *     .. External Functions ..
-      LOGICAL            LSAME
-      INTEGER            IDAMAX
-      EXTERNAL           LSAME, IDAMAX
+      LOGICAL            LSAME, DISNAN
+      EXTERNAL           LSAME, DISNAN
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           DLASSQ
@@ -299,14 +298,18 @@
 *           A is n by k
                DO J = 0, K - 1
                   DO I = 0, N - 1
-                     VALUE = MAX( VALUE, ABS( A( I+J*LDA ) ) )
+                     TEMP = ABS( A( I+J*LDA ) )
+                     IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
+     $                    VALUE = TEMP
                   END DO
                END DO
             ELSE
 *              xpose case; A is k by n
                DO J = 0, N - 1
                   DO I = 0, K - 1
-                     VALUE = MAX( VALUE, ABS( A( I+J*LDA ) ) )
+                     TEMP = ABS( A( I+J*LDA ) )
+                     IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
+     $                    VALUE = TEMP
                   END DO
                END DO
             END IF
@@ -316,14 +319,18 @@
 *              A is n+1 by k
                DO J = 0, K - 1
                   DO I = 0, N
-                     VALUE = MAX( VALUE, ABS( A( I+J*LDA ) ) )
+                     TEMP = ABS( A( I+J*LDA ) )
+                     IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
+     $                    VALUE = TEMP
                   END DO
                END DO
             ELSE
 *              xpose case; A is k by n+1
                DO J = 0, N
                   DO I = 0, K - 1
-                     VALUE = MAX( VALUE, ABS( A( I+J*LDA ) ) )
+                     TEMP = ABS( A( I+J*LDA ) )
+                     IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
+     $                    VALUE = TEMP
                   END DO
                END DO
             END IF
@@ -369,8 +376,12 @@
                      WORK( J ) = WORK( J ) + S
                   END DO
    10             CONTINUE
-                  I = IDAMAX( N, WORK, 1 )
-                  VALUE = WORK( I-1 )
+                  VALUE = WORK( 0 )
+                  DO I = 1, N-1
+                     TEMP = WORK( I )
+                     IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
+     $                    VALUE = TEMP
+                  END DO
                ELSE
 *                 ilu = 1
                   K = K + 1
@@ -407,8 +418,12 @@
                      END DO
                      WORK( J ) = WORK( J ) + S
                   END DO
-                  I = IDAMAX( N, WORK, 1 )
-                  VALUE = WORK( I-1 )
+                  VALUE = WORK( 0 )
+                  DO I = 1, N-1
+                     TEMP = WORK( I )
+                     IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
+     $                    VALUE = TEMP
+                  END DO
                END IF
             ELSE
 *              n is even
@@ -441,8 +456,12 @@
                      END DO
                      WORK( J ) = WORK( J ) + S
                   END DO
-                  I = IDAMAX( N, WORK, 1 )
-                  VALUE = WORK( I-1 )
+                  VALUE = WORK( 0 )
+                  DO I = 1, N-1
+                     TEMP = WORK( I )
+                     IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
+     $                    VALUE = TEMP
+                  END DO
                ELSE
 *                 ilu = 1
                   DO I = K, N - 1
@@ -475,8 +494,12 @@
                      END DO
                      WORK( J ) = WORK( J ) + S
                   END DO
-                  I = IDAMAX( N, WORK, 1 )
-                  VALUE = WORK( I-1 )
+                  VALUE = WORK( 0 )
+                  DO I = 1, N-1
+                     TEMP = WORK( I )
+                     IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
+     $                    VALUE = TEMP
+                  END DO
                END IF
             END IF
          ELSE
@@ -537,8 +560,12 @@
                      END DO
                      WORK( J ) = WORK( J ) + S
                   END DO
-                  I = IDAMAX( N, WORK, 1 )
-                  VALUE = WORK( I-1 )
+                  VALUE = WORK( 0 )
+                  DO I = 1, N-1
+                     TEMP = WORK( I )
+                     IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
+     $                    VALUE = TEMP
+                  END DO
                ELSE
 *                 ilu=1
                   K = K + 1
@@ -598,8 +625,12 @@
                      END DO
                      WORK( J ) = WORK( J ) + S
                   END DO
-                  I = IDAMAX( N, WORK, 1 )
-                  VALUE = WORK( I-1 )
+                  VALUE = WORK( 0 )
+                  DO I = 1, N-1
+                     TEMP = WORK( I )
+                     IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
+     $                    VALUE = TEMP
+                  END DO
                END IF
             ELSE
 *              n is even
@@ -667,8 +698,12 @@
 *                 A(k-1,k-1)
                   S = S + AA
                   WORK( I ) = WORK( I ) + S
-                  I = IDAMAX( N, WORK, 1 )
-                  VALUE = WORK( I-1 )
+                  VALUE = WORK( 0 )
+                  DO I = 1, N-1
+                     TEMP = WORK( I )
+                     IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
+     $                    VALUE = TEMP
+                  END DO
                ELSE
 *                 ilu=1
                   DO I = K, N - 1
@@ -736,8 +771,12 @@
                      END DO
                      WORK( J-1 ) = WORK( J-1 ) + S
                   END DO
-                  I = IDAMAX( N, WORK, 1 )
-                  VALUE = WORK( I-1 )
+                  VALUE = WORK( 0 )
+                  DO I = 1, N-1
+                     TEMP = WORK( I )
+                     IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
+     $                    VALUE = TEMP
+                  END DO
                END IF
             END IF
          END IF