version 1.18, 2018/05/29 07:18:32
|
version 1.19, 2023/08/07 08:39:34
|
Line 102
|
Line 102
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date December 2016 |
|
* |
|
*> \ingroup complex16POcomputational |
*> \ingroup complex16POcomputational |
* |
* |
* ===================================================================== |
* ===================================================================== |
SUBROUTINE ZPOTF2( UPLO, N, A, LDA, INFO ) |
SUBROUTINE ZPOTF2( UPLO, N, A, LDA, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.7.0) -- |
* -- LAPACK computational 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 |
|
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER UPLO |
CHARACTER UPLO |
Line 177
|
Line 174
|
* |
* |
* Compute U(J,J) and test for non-positive-definiteness. |
* Compute U(J,J) and test for non-positive-definiteness. |
* |
* |
AJJ = DBLE( A( J, J ) ) - ZDOTC( J-1, A( 1, J ), 1, |
AJJ = DBLE( A( J, J ) ) - DBLE( ZDOTC( J-1, A( 1, J ), 1, |
$ A( 1, J ), 1 ) |
$ A( 1, J ), 1 ) ) |
IF( AJJ.LE.ZERO.OR.DISNAN( AJJ ) ) THEN |
IF( AJJ.LE.ZERO.OR.DISNAN( AJJ ) ) THEN |
A( J, J ) = AJJ |
A( J, J ) = AJJ |
GO TO 30 |
GO TO 30 |
Line 204
|
Line 201
|
* |
* |
* Compute L(J,J) and test for non-positive-definiteness. |
* Compute L(J,J) and test for non-positive-definiteness. |
* |
* |
AJJ = DBLE( A( J, J ) ) - ZDOTC( J-1, A( J, 1 ), LDA, |
AJJ = DBLE( A( J, J ) ) - DBLE( ZDOTC( J-1, A( J, 1 ), LDA, |
$ A( J, 1 ), LDA ) |
$ A( J, 1 ), LDA ) ) |
IF( AJJ.LE.ZERO.OR.DISNAN( AJJ ) ) THEN |
IF( AJJ.LE.ZERO.OR.DISNAN( AJJ ) ) THEN |
A( J, J ) = AJJ |
A( J, J ) = AJJ |
GO TO 30 |
GO TO 30 |