|
version 1.19, 2020/05/21 21:46:08
|
version 1.20, 2023/08/07 08:39:30
|
|
Line 116
|
Line 116
|
| *> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
| *> \author NAG Ltd. |
*> \author NAG Ltd. |
| * |
* |
| *> \date December 2016 |
|
| * |
|
| *> \ingroup complex16SYauxiliary |
*> \ingroup complex16SYauxiliary |
| * |
* |
| * ===================================================================== |
* ===================================================================== |
| DOUBLE PRECISION FUNCTION ZLANSY( NORM, UPLO, N, A, LDA, WORK ) |
DOUBLE PRECISION FUNCTION ZLANSY( NORM, UPLO, N, A, LDA, WORK ) |
| * |
* |
| * -- LAPACK auxiliary routine (version 3.7.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..-- |
| * December 2016 |
|
| * |
* |
| IMPLICIT NONE |
|
| * .. Scalar Arguments .. |
* .. Scalar Arguments .. |
| CHARACTER NORM, UPLO |
CHARACTER NORM, UPLO |
| INTEGER LDA, N |
INTEGER LDA, N |
|
Line 146
|
Line 142
|
| * .. |
* .. |
| * .. Local Scalars .. |
* .. Local Scalars .. |
| INTEGER I, J |
INTEGER I, J |
| DOUBLE PRECISION ABSA, SUM, VALUE |
DOUBLE PRECISION ABSA, SCALE, SUM, VALUE |
| * .. |
|
| * .. Local Arrays .. |
|
| DOUBLE PRECISION SSQ( 2 ), COLSSQ( 2 ) |
|
| * .. |
* .. |
| * .. External Functions .. |
* .. External Functions .. |
| LOGICAL LSAME, DISNAN |
LOGICAL LSAME, DISNAN |
| EXTERNAL LSAME, DISNAN |
EXTERNAL LSAME, DISNAN |
| * .. |
* .. |
| * .. External Subroutines .. |
* .. External Subroutines .. |
| EXTERNAL ZLASSQ, DCOMBSSQ |
EXTERNAL ZLASSQ |
| * .. |
* .. |
| * .. Intrinsic Functions .. |
* .. Intrinsic Functions .. |
| INTRINSIC ABS, SQRT |
INTRINSIC ABS, SQRT |
|
Line 222
|
Line 215
|
| ELSE IF( ( LSAME( NORM, 'F' ) ) .OR. ( LSAME( NORM, 'E' ) ) ) THEN |
ELSE IF( ( LSAME( NORM, 'F' ) ) .OR. ( LSAME( NORM, 'E' ) ) ) THEN |
| * |
* |
| * Find normF(A). |
* Find normF(A). |
| * SSQ(1) is scale |
|
| * SSQ(2) is sum-of-squares |
|
| * For better accuracy, sum each column separately. |
|
| * |
|
| SSQ( 1 ) = ZERO |
|
| SSQ( 2 ) = ONE |
|
| * |
|
| * Sum off-diagonals |
|
| * |
* |
| |
SCALE = ZERO |
| |
SUM = ONE |
| IF( LSAME( UPLO, 'U' ) ) THEN |
IF( LSAME( UPLO, 'U' ) ) THEN |
| DO 110 J = 2, N |
DO 110 J = 2, N |
| COLSSQ( 1 ) = ZERO |
CALL ZLASSQ( J-1, A( 1, J ), 1, SCALE, SUM ) |
| COLSSQ( 2 ) = ONE |
|
| CALL ZLASSQ( J-1, A( 1, J ), 1, COLSSQ(1), COLSSQ(2) ) |
|
| CALL DCOMBSSQ( SSQ, COLSSQ ) |
|
| 110 CONTINUE |
110 CONTINUE |
| ELSE |
ELSE |
| DO 120 J = 1, N - 1 |
DO 120 J = 1, N - 1 |
| COLSSQ( 1 ) = ZERO |
CALL ZLASSQ( N-J, A( J+1, J ), 1, SCALE, SUM ) |
| COLSSQ( 2 ) = ONE |
|
| CALL ZLASSQ( N-J, A( J+1, J ), 1, COLSSQ(1), COLSSQ(2) ) |
|
| CALL DCOMBSSQ( SSQ, COLSSQ ) |
|
| 120 CONTINUE |
120 CONTINUE |
| END IF |
END IF |
| SSQ( 2 ) = 2*SSQ( 2 ) |
SUM = 2*SUM |
| * |
CALL ZLASSQ( N, A, LDA+1, SCALE, SUM ) |
| * Sum diagonal |
VALUE = SCALE*SQRT( SUM ) |
| * |
|
| COLSSQ( 1 ) = ZERO |
|
| COLSSQ( 2 ) = ONE |
|
| CALL ZLASSQ( N, A, LDA+1, COLSSQ( 1 ), COLSSQ( 2 ) ) |
|
| CALL DCOMBSSQ( SSQ, COLSSQ ) |
|
| VALUE = SSQ( 1 )*SQRT( SSQ( 2 ) ) |
|
| END IF |
END IF |
| * |
* |
| ZLANSY = VALUE |
ZLANSY = VALUE |
![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to zlansy.f CVS log ![[TXT]](/cvsweb/icons/text.gif){kind=icon}
![[DIR]](/cvsweb/icons/dir.gif){kind=icon}
Up to [local] / rpl / lapack / lapack