rpl/lapack/lapack/dlansf.f - diff

Return to dlansf.f CVS log

Up to [local] / rpl / lapack / lapack

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

version 1.10, 2012/08/22 09:48:18	version 1.11, 2012/12/14 12:30:23
Line 1	Line 1
*> \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 ===========	* =========== DOCUMENTATION ===========
*	*
Line 115	Line 115
*> \author Univ. of Colorado Denver	*> \author Univ. of Colorado Denver
*> \author NAG Ltd.	*> \author NAG Ltd.
*	*
*> \date November 2011	*> \date September 2012
*	*
*> \ingroup doubleOTHERcomputational	*> \ingroup doubleOTHERcomputational
*	*
Line 209	Line 209
* =====================================================================	* =====================================================================
DOUBLE PRECISION FUNCTION DLANSF( NORM, TRANSR, UPLO, N, A, WORK )	DOUBLE PRECISION FUNCTION DLANSF( NORM, TRANSR, UPLO, N, A, WORK )
*	*
* -- LAPACK computational routine (version 3.4.0) --	* -- LAPACK computational 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, TRANSR, UPLO	CHARACTER NORM, TRANSR, UPLO
Line 230	Line 230
* ..	* ..
* .. Local Scalars ..	* .. Local Scalars ..
INTEGER I, J, IFM, ILU, NOE, N1, K, L, LDA	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 ..	* .. External Functions ..
LOGICAL LSAME	LOGICAL LSAME, DISNAN
INTEGER IDAMAX	EXTERNAL LSAME, DISNAN
EXTERNAL LSAME, IDAMAX
* ..	* ..
* .. External Subroutines ..	* .. External Subroutines ..
EXTERNAL DLASSQ	EXTERNAL DLASSQ
Line 299	Line 298
* A is n by k	* A is n by k
DO J = 0, K - 1	DO J = 0, K - 1
DO I = 0, N - 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
END DO	END DO
ELSE	ELSE
* xpose case; A is k by n	* xpose case; A is k by n
DO J = 0, N - 1	DO J = 0, N - 1
DO I = 0, K - 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 DO	END DO
END IF	END IF
Line 316	Line 319
* A is n+1 by k	* A is n+1 by k
DO J = 0, K - 1	DO J = 0, K - 1
DO I = 0, N	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
END DO	END DO
ELSE	ELSE
* xpose case; A is k by n+1	* xpose case; A is k by n+1
DO J = 0, N	DO J = 0, N
DO I = 0, K - 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 DO	END DO
END IF	END IF
Line 369	Line 376
WORK( J ) = WORK( J ) + S	WORK( J ) = WORK( J ) + S
END DO	END DO
10 CONTINUE	10 CONTINUE
I = IDAMAX( N, WORK, 1 )	VALUE = WORK( 0 )
VALUE = WORK( I-1 )	DO I = 1, N-1
	TEMP = WORK( I )
	IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
	$ VALUE = TEMP
	END DO
ELSE	ELSE
* ilu = 1	* ilu = 1
K = K + 1	K = K + 1
Line 407	Line 418
END DO	END DO
WORK( J ) = WORK( J ) + S	WORK( J ) = WORK( J ) + S
END DO	END DO
I = IDAMAX( N, WORK, 1 )	VALUE = WORK( 0 )
VALUE = WORK( I-1 )	DO I = 1, N-1
	TEMP = WORK( I )
	IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
	$ VALUE = TEMP
	END DO
END IF	END IF
ELSE	ELSE
* n is even	* n is even
Line 441	Line 456
END DO	END DO
WORK( J ) = WORK( J ) + S	WORK( J ) = WORK( J ) + S
END DO	END DO
I = IDAMAX( N, WORK, 1 )	VALUE = WORK( 0 )
VALUE = WORK( I-1 )	DO I = 1, N-1
	TEMP = WORK( I )
	IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
	$ VALUE = TEMP
	END DO
ELSE	ELSE
* ilu = 1	* ilu = 1
DO I = K, N - 1	DO I = K, N - 1
Line 475	Line 494
END DO	END DO
WORK( J ) = WORK( J ) + S	WORK( J ) = WORK( J ) + S
END DO	END DO
I = IDAMAX( N, WORK, 1 )	VALUE = WORK( 0 )
VALUE = WORK( I-1 )	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	END IF
ELSE	ELSE
Line 537	Line 560
END DO	END DO
WORK( J ) = WORK( J ) + S	WORK( J ) = WORK( J ) + S
END DO	END DO
I = IDAMAX( N, WORK, 1 )	VALUE = WORK( 0 )
VALUE = WORK( I-1 )	DO I = 1, N-1
	TEMP = WORK( I )
	IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
	$ VALUE = TEMP
	END DO
ELSE	ELSE
* ilu=1	* ilu=1
K = K + 1	K = K + 1
Line 598	Line 625
END DO	END DO
WORK( J ) = WORK( J ) + S	WORK( J ) = WORK( J ) + S
END DO	END DO
I = IDAMAX( N, WORK, 1 )	VALUE = WORK( 0 )
VALUE = WORK( I-1 )	DO I = 1, N-1
	TEMP = WORK( I )
	IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
	$ VALUE = TEMP
	END DO
END IF	END IF
ELSE	ELSE
* n is even	* n is even
Line 667	Line 698
* A(k-1,k-1)	* A(k-1,k-1)
S = S + AA	S = S + AA
WORK( I ) = WORK( I ) + S	WORK( I ) = WORK( I ) + S
I = IDAMAX( N, WORK, 1 )	VALUE = WORK( 0 )
VALUE = WORK( I-1 )	DO I = 1, N-1
	TEMP = WORK( I )
	IF( VALUE .LT. TEMP .OR. DISNAN( TEMP ) )
	$ VALUE = TEMP
	END DO
ELSE	ELSE
* ilu=1	* ilu=1
DO I = K, N - 1	DO I = K, N - 1
Line 736	Line 771
END DO	END DO
WORK( J-1 ) = WORK( J-1 ) + S	WORK( J-1 ) = WORK( J-1 ) + S
END DO	END DO
I = IDAMAX( N, WORK, 1 )	VALUE = WORK( 0 )
VALUE = WORK( I-1 )	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	END IF
END IF	END IF

CVSweb interface <joel.bertrand@systella.fr>

Removed from v.1.10
changed lines
	Added in v.1.11