rpl/lapack/lapack/dlantr.f - diff

Return to dlantr.f CVS log

Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/dlantr.f between versions 1.8 and 1.19

version 1.8, 2011/11/21 20:42:56	version 1.19, 2023/08/07 08:38:55
Line 1	Line 1
*> \brief \b DLANTR	*> \brief \b DLANTR returns the value of the 1-norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a trapezoidal or triangular matrix.
*	*
* =========== DOCUMENTATION ===========	* =========== DOCUMENTATION ===========
*	*
* Online html documentation available at	* Online html documentation available at
* http://www.netlib.org/lapack/explore-html/	* http://www.netlib.org/lapack/explore-html/
*	*
*> \htmlonly	*> \htmlonly
*> Download DLANTR + dependencies	*> Download DLANTR + dependencies
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dlantr.f">	*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dlantr.f">
*> [TGZ]</a>	*> [TGZ]</a>
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dlantr.f">	*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dlantr.f">
*> [ZIP]</a>	*> [ZIP]</a>
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dlantr.f">	*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dlantr.f">
*> [TXT]</a>	*> [TXT]</a>
*> \endhtmlonly	*> \endhtmlonly
*	*
* Definition:	* Definition:
* ===========	* ===========
*	*
* DOUBLE PRECISION FUNCTION DLANTR( NORM, UPLO, DIAG, M, N, A, LDA,	* DOUBLE PRECISION FUNCTION DLANTR( NORM, UPLO, DIAG, M, N, A, LDA,
* WORK )	* WORK )
*	*
* .. Scalar Arguments ..	* .. Scalar Arguments ..
* CHARACTER DIAG, NORM, UPLO	* CHARACTER DIAG, NORM, UPLO
* INTEGER LDA, M, N	* INTEGER LDA, M, N
Line 28	Line 28
* .. Array Arguments ..	* .. Array Arguments ..
* DOUBLE PRECISION A( LDA, * ), WORK( * )	* DOUBLE PRECISION A( LDA, * ), WORK( * )
* ..	* ..
*	*
*	*
*> \par Purpose:	*> \par Purpose:
* =============	* =============
Line 128	Line 128
* Authors:	* Authors:
* ========	* ========
*	*
*> \author Univ. of Tennessee	*> \author Univ. of Tennessee
*> \author Univ. of California Berkeley	*> \author Univ. of California Berkeley
*> \author Univ. of Colorado Denver	*> \author Univ. of Colorado Denver
*> \author NAG Ltd.	*> \author NAG Ltd.
*
*> \date November 2011
*	*
*> \ingroup doubleOTHERauxiliary	*> \ingroup doubleOTHERauxiliary
*	*
Line 141	Line 139
DOUBLE PRECISION FUNCTION DLANTR( NORM, UPLO, DIAG, M, N, A, LDA,	DOUBLE PRECISION FUNCTION DLANTR( NORM, UPLO, DIAG, M, N, A, LDA,
$ WORK )	$ WORK )
*	*
* -- LAPACK auxiliary routine (version 3.4.0) --	* -- LAPACK auxiliary 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 2011
*	*
* .. Scalar Arguments ..	* .. Scalar Arguments ..
CHARACTER DIAG, NORM, UPLO	CHARACTER DIAG, NORM, UPLO
Line 169	Line 166
EXTERNAL DLASSQ	EXTERNAL DLASSQ
* ..	* ..
* .. External Functions ..	* .. External Functions ..
LOGICAL LSAME	LOGICAL LSAME, DISNAN
EXTERNAL LSAME	EXTERNAL LSAME, DISNAN
* ..	* ..
* .. Intrinsic Functions ..	* .. Intrinsic Functions ..
INTRINSIC ABS, MAX, MIN, SQRT	INTRINSIC ABS, MIN, SQRT
* ..	* ..
* .. Executable Statements ..	* .. Executable Statements ..
*	*
Line 188	Line 185
IF( LSAME( UPLO, 'U' ) ) THEN	IF( LSAME( UPLO, 'U' ) ) THEN
DO 20 J = 1, N	DO 20 J = 1, N
DO 10 I = 1, MIN( M, J-1 )	DO 10 I = 1, MIN( M, J-1 )
VALUE = MAX( VALUE, ABS( A( I, J ) ) )	SUM = ABS( A( I, J ) )
	IF( VALUE .LT. SUM .OR. DISNAN( SUM ) ) VALUE = SUM
10 CONTINUE	10 CONTINUE
20 CONTINUE	20 CONTINUE
ELSE	ELSE
DO 40 J = 1, N	DO 40 J = 1, N
DO 30 I = J + 1, M	DO 30 I = J + 1, M
VALUE = MAX( VALUE, ABS( A( I, J ) ) )	SUM = ABS( A( I, J ) )
	IF( VALUE .LT. SUM .OR. DISNAN( SUM ) ) VALUE = SUM
30 CONTINUE	30 CONTINUE
40 CONTINUE	40 CONTINUE
END IF	END IF
Line 203	Line 202
IF( LSAME( UPLO, 'U' ) ) THEN	IF( LSAME( UPLO, 'U' ) ) THEN
DO 60 J = 1, N	DO 60 J = 1, N
DO 50 I = 1, MIN( M, J )	DO 50 I = 1, MIN( M, J )
VALUE = MAX( VALUE, ABS( A( I, J ) ) )	SUM = ABS( A( I, J ) )
	IF( VALUE .LT. SUM .OR. DISNAN( SUM ) ) VALUE = SUM
50 CONTINUE	50 CONTINUE
60 CONTINUE	60 CONTINUE
ELSE	ELSE
DO 80 J = 1, N	DO 80 J = 1, N
DO 70 I = J, M	DO 70 I = J, M
VALUE = MAX( VALUE, ABS( A( I, J ) ) )	SUM = ABS( A( I, J ) )
	IF( VALUE .LT. SUM .OR. DISNAN( SUM ) ) VALUE = SUM
70 CONTINUE	70 CONTINUE
80 CONTINUE	80 CONTINUE
END IF	END IF
Line 233	Line 234
SUM = SUM + ABS( A( I, J ) )	SUM = SUM + ABS( A( I, J ) )
100 CONTINUE	100 CONTINUE
END IF	END IF
VALUE = MAX( VALUE, SUM )	IF( VALUE .LT. SUM .OR. DISNAN( SUM ) ) VALUE = SUM
110 CONTINUE	110 CONTINUE
ELSE	ELSE
DO 140 J = 1, N	DO 140 J = 1, N
Line 248	Line 249
SUM = SUM + ABS( A( I, J ) )	SUM = SUM + ABS( A( I, J ) )
130 CONTINUE	130 CONTINUE
END IF	END IF
VALUE = MAX( VALUE, SUM )	IF( VALUE .LT. SUM .OR. DISNAN( SUM ) ) VALUE = SUM
140 CONTINUE	140 CONTINUE
END IF	END IF
ELSE IF( LSAME( NORM, 'I' ) ) THEN	ELSE IF( LSAME( NORM, 'I' ) ) THEN
Line 277	Line 278
END IF	END IF
ELSE	ELSE
IF( LSAME( DIAG, 'U' ) ) THEN	IF( LSAME( DIAG, 'U' ) ) THEN
DO 210 I = 1, N	DO 210 I = 1, MIN( M, N )
WORK( I ) = ONE	WORK( I ) = ONE
210 CONTINUE	210 CONTINUE
DO 220 I = N + 1, M	DO 220 I = N + 1, M
Line 301	Line 302
END IF	END IF
VALUE = ZERO	VALUE = ZERO
DO 280 I = 1, M	DO 280 I = 1, M
VALUE = MAX( VALUE, WORK( I ) )	SUM = WORK( I )
	IF( VALUE .LT. SUM .OR. DISNAN( SUM ) ) VALUE = SUM
280 CONTINUE	280 CONTINUE
ELSE IF( ( LSAME( NORM, 'F' ) ) .OR. ( LSAME( NORM, 'E' ) ) ) THEN	ELSE IF( ( LSAME( NORM, 'F' ) ) .OR. ( LSAME( NORM, 'E' ) ) ) THEN
*	*

CVSweb interface <joel.bertrand@systella.fr>

Removed from v.1.8
changed lines
	Added in v.1.19