version 1.4, 2010/12/21 13:53:52
|
version 1.5, 2011/07/22 07:38:18
|
Line 1
|
Line 1
|
SUBROUTINE ZLAT2C( UPLO, N, A, LDA, SA, LDSA, INFO ) |
SUBROUTINE ZLAT2C( UPLO, N, A, LDA, SA, LDSA, INFO ) |
* |
* |
* -- LAPACK PROTOTYPE auxiliary routine (version 3.1.2) -- |
* -- LAPACK PROTOTYPE auxiliary routine (version 3.3.1) -- |
* -- 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..-- |
* May 2007 |
* -- April 2011 -- |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER UPLO |
CHARACTER UPLO |
Line 56
|
Line 56
|
* PRECISION overflow threshold, in this case, the content |
* PRECISION overflow threshold, in this case, the content |
* of the UPLO part of SA in exit is unspecified. |
* of the UPLO part of SA in exit is unspecified. |
* |
* |
* ========= |
* ===================================================================== |
* |
* |
* .. Local Scalars .. |
* .. Local Scalars .. |
INTEGER I, J |
INTEGER I, J |
Line 79
|
Line 79
|
DO 20 J = 1, N |
DO 20 J = 1, N |
DO 10 I = 1, J |
DO 10 I = 1, J |
IF( ( DBLE( A( I, J ) ).LT.-RMAX ) .OR. |
IF( ( DBLE( A( I, J ) ).LT.-RMAX ) .OR. |
+ ( DBLE( A( I, J ) ).GT.RMAX ) .OR. |
$ ( DBLE( A( I, J ) ).GT.RMAX ) .OR. |
+ ( DIMAG( A( I, J ) ).LT.-RMAX ) .OR. |
$ ( DIMAG( A( I, J ) ).LT.-RMAX ) .OR. |
+ ( DIMAG( A( I, J ) ).GT.RMAX ) ) THEN |
$ ( DIMAG( A( I, J ) ).GT.RMAX ) ) THEN |
INFO = 1 |
INFO = 1 |
GO TO 50 |
GO TO 50 |
END IF |
END IF |
Line 92
|
Line 92
|
DO 40 J = 1, N |
DO 40 J = 1, N |
DO 30 I = J, N |
DO 30 I = J, N |
IF( ( DBLE( A( I, J ) ).LT.-RMAX ) .OR. |
IF( ( DBLE( A( I, J ) ).LT.-RMAX ) .OR. |
+ ( DBLE( A( I, J ) ).GT.RMAX ) .OR. |
$ ( DBLE( A( I, J ) ).GT.RMAX ) .OR. |
+ ( DIMAG( A( I, J ) ).LT.-RMAX ) .OR. |
$ ( DIMAG( A( I, J ) ).LT.-RMAX ) .OR. |
+ ( DIMAG( A( I, J ) ).GT.RMAX ) ) THEN |
$ ( DIMAG( A( I, J ) ).GT.RMAX ) ) THEN |
INFO = 1 |
INFO = 1 |
GO TO 50 |
GO TO 50 |
END IF |
END IF |