Use potrf2 in Cholesky variants

angsch · angsch · commit 8e4ad311c1a5 · 2023-08-09T19:12:09.000+02:00
Replace the scalar potf2 with the recursive potrf2 implementation.
The latter is already used in the default Cholesky factorization.
diff --git a/SRC/VARIANTS/cholesky/RL/cpotrf.f b/SRC/VARIANTS/cholesky/RL/cpotrf.f
@@ -129,7 +129,7 @@ SUBROUTINE CPOTRF ( UPLO, N, A, LDA, INFO )
       EXTERNAL           LSAME, ILAENV
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CGEMM, CPOTF2, CHERK, CTRSM, XERBLA
+      EXTERNAL           CGEMM, CHERK, CPOTRF2, CTRSM, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -164,7 +164,7 @@ SUBROUTINE CPOTRF ( UPLO, N, A, LDA, INFO )
 *
 *        Use unblocked code.
 *
-         CALL CPOTF2( UPLO, N, A, LDA, INFO )
+         CALL CPOTRF2( UPLO, N, A, LDA, INFO )
       ELSE
 *
 *        Use blocked code.
@@ -180,7 +180,7 @@ SUBROUTINE CPOTRF ( UPLO, N, A, LDA, INFO )
 *
                JB = MIN( NB, N-J+1 )
 
-               CALL CPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )
+               CALL CPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )
 
                IF( INFO.NE.0 )
      $            GO TO 30
@@ -209,7 +209,7 @@ SUBROUTINE CPOTRF ( UPLO, N, A, LDA, INFO )
 *
                JB = MIN( NB, N-J+1 )
 
-               CALL CPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )
+               CALL CPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )
 
                IF( INFO.NE.0 )
      $            GO TO 30
diff --git a/SRC/VARIANTS/cholesky/RL/dpotrf.f b/SRC/VARIANTS/cholesky/RL/dpotrf.f
@@ -128,7 +128,7 @@ SUBROUTINE DPOTRF ( UPLO, N, A, LDA, INFO )
       EXTERNAL           LSAME, ILAENV
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DGEMM, DPOTF2, DSYRK, DTRSM, XERBLA
+      EXTERNAL           DGEMM, DPOTRF2, DSYRK, DTRSM, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -163,7 +163,7 @@ SUBROUTINE DPOTRF ( UPLO, N, A, LDA, INFO )
 *
 *        Use unblocked code.
 *
-         CALL DPOTF2( UPLO, N, A, LDA, INFO )
+         CALL DPOTRF2( UPLO, N, A, LDA, INFO )
       ELSE
 *
 *        Use blocked code.
@@ -179,7 +179,7 @@ SUBROUTINE DPOTRF ( UPLO, N, A, LDA, INFO )
 *
                JB = MIN( NB, N-J+1 )
 
-               CALL DPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )
+               CALL DPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )
 
                IF( INFO.NE.0 )
      $            GO TO 30
@@ -208,7 +208,7 @@ SUBROUTINE DPOTRF ( UPLO, N, A, LDA, INFO )
 *
                JB = MIN( NB, N-J+1 )
 
-               CALL DPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )
+               CALL DPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )
 
                IF( INFO.NE.0 )
      $            GO TO 30
diff --git a/SRC/VARIANTS/cholesky/RL/spotrf.f b/SRC/VARIANTS/cholesky/RL/spotrf.f
@@ -128,7 +128,7 @@ SUBROUTINE SPOTRF ( UPLO, N, A, LDA, INFO )
       EXTERNAL           LSAME, ILAENV
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SGEMM, SPOTF2, SSYRK, STRSM, XERBLA
+      EXTERNAL           SGEMM, SPOTRF2, SSYRK, STRSM, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -163,7 +163,7 @@ SUBROUTINE SPOTRF ( UPLO, N, A, LDA, INFO )
 *
 *        Use unblocked code.
 *
-         CALL SPOTF2( UPLO, N, A, LDA, INFO )
+         CALL SPOTRF2( UPLO, N, A, LDA, INFO )
       ELSE
 *
 *        Use blocked code.
@@ -179,7 +179,7 @@ SUBROUTINE SPOTRF ( UPLO, N, A, LDA, INFO )
 *
                JB = MIN( NB, N-J+1 )
 
-               CALL SPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )
+               CALL SPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )
 
                IF( INFO.NE.0 )
      $            GO TO 30
@@ -208,7 +208,7 @@ SUBROUTINE SPOTRF ( UPLO, N, A, LDA, INFO )
 *
                JB = MIN( NB, N-J+1 )
 
-               CALL SPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )
+               CALL SPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )
 
                IF( INFO.NE.0 )
      $            GO TO 30
diff --git a/SRC/VARIANTS/cholesky/RL/zpotrf.f b/SRC/VARIANTS/cholesky/RL/zpotrf.f
@@ -129,7 +129,7 @@ SUBROUTINE ZPOTRF ( UPLO, N, A, LDA, INFO )
       EXTERNAL           LSAME, ILAENV
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           ZGEMM, ZPOTF2, ZHERK, ZTRSM, XERBLA
+      EXTERNAL           XERBLA, ZGEMM, ZHERK, ZPOTRF2, ZTRSM
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -164,7 +164,7 @@ SUBROUTINE ZPOTRF ( UPLO, N, A, LDA, INFO )
 *
 *        Use unblocked code.
 *
-         CALL ZPOTF2( UPLO, N, A, LDA, INFO )
+         CALL ZPOTRF2( UPLO, N, A, LDA, INFO )
       ELSE
 *
 *        Use blocked code.
@@ -180,7 +180,7 @@ SUBROUTINE ZPOTRF ( UPLO, N, A, LDA, INFO )
 *
                JB = MIN( NB, N-J+1 )
 
-               CALL ZPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )
+               CALL ZPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )
 
                IF( INFO.NE.0 )
      $            GO TO 30
@@ -209,7 +209,7 @@ SUBROUTINE ZPOTRF ( UPLO, N, A, LDA, INFO )
 *
                JB = MIN( NB, N-J+1 )
 
-               CALL ZPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )
+               CALL ZPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )
 
                IF( INFO.NE.0 )
      $            GO TO 30
diff --git a/SRC/VARIANTS/cholesky/TOP/cpotrf.f b/SRC/VARIANTS/cholesky/TOP/cpotrf.f
@@ -129,7 +129,7 @@ SUBROUTINE CPOTRF ( UPLO, N, A, LDA, INFO )
       EXTERNAL           LSAME, ILAENV
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CGEMM, CPOTF2, CHERK, CTRSM, XERBLA
+      EXTERNAL           CGEMM, CHERK, CPOTRF2, CTRSM, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -164,7 +164,7 @@ SUBROUTINE CPOTRF ( UPLO, N, A, LDA, INFO )
 *
 *        Use unblocked code.
 *
-         CALL CPOTF2( UPLO, N, A, LDA, INFO )
+         CALL CPOTRF2( UPLO, N, A, LDA, INFO )
       ELSE
 *
 *        Use blocked code.
@@ -189,7 +189,7 @@ SUBROUTINE CPOTRF ( UPLO, N, A, LDA, INFO )
 *              Update and factorize the current diagonal block and test
 *              for non-positive-definiteness.
 *
-               CALL CPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )
+               CALL CPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )
                IF( INFO.NE.0 )
      $            GO TO 30
 
@@ -216,7 +216,7 @@ SUBROUTINE CPOTRF ( UPLO, N, A, LDA, INFO )
 *              Update and factorize the current diagonal block and test
 *              for non-positive-definiteness.
 *
-               CALL CPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )
+               CALL CPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )
                IF( INFO.NE.0 )
      $            GO TO 30
 
diff --git a/SRC/VARIANTS/cholesky/TOP/dpotrf.f b/SRC/VARIANTS/cholesky/TOP/dpotrf.f
@@ -128,7 +128,7 @@ SUBROUTINE DPOTRF ( UPLO, N, A, LDA, INFO )
       EXTERNAL           LSAME, ILAENV
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DGEMM, DPOTF2, DSYRK, DTRSM, XERBLA
+      EXTERNAL           DGEMM, DPOTRF2, DSYRK, DTRSM, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -163,7 +163,7 @@ SUBROUTINE DPOTRF ( UPLO, N, A, LDA, INFO )
 *
 *        Use unblocked code.
 *
-         CALL DPOTF2( UPLO, N, A, LDA, INFO )
+         CALL DPOTRF2( UPLO, N, A, LDA, INFO )
       ELSE
 *
 *        Use blocked code.
@@ -189,7 +189,7 @@ SUBROUTINE DPOTRF ( UPLO, N, A, LDA, INFO )
 *              Update and factorize the current diagonal block and test
 *              for non-positive-definiteness.
 *
-               CALL DPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )
+               CALL DPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )
                IF( INFO.NE.0 )
      $            GO TO 30
 
@@ -217,7 +217,7 @@ SUBROUTINE DPOTRF ( UPLO, N, A, LDA, INFO )
 *              Update and factorize the current diagonal block and test
 *              for non-positive-definiteness.
 *
-               CALL DPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )
+               CALL DPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )
                IF( INFO.NE.0 )
      $            GO TO 30
 
diff --git a/SRC/VARIANTS/cholesky/TOP/spotrf.f b/SRC/VARIANTS/cholesky/TOP/spotrf.f
@@ -128,7 +128,7 @@ SUBROUTINE SPOTRF ( UPLO, N, A, LDA, INFO )
       EXTERNAL           LSAME, ILAENV
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SGEMM, SPOTF2, SSYRK, STRSM, XERBLA
+      EXTERNAL           SGEMM, SPOTRF2, SSYRK, STRSM, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -163,7 +163,7 @@ SUBROUTINE SPOTRF ( UPLO, N, A, LDA, INFO )
 *
 *        Use unblocked code.
 *
-         CALL SPOTF2( UPLO, N, A, LDA, INFO )
+         CALL SPOTRF2( UPLO, N, A, LDA, INFO )
       ELSE
 *
 *        Use blocked code.
@@ -189,7 +189,7 @@ SUBROUTINE SPOTRF ( UPLO, N, A, LDA, INFO )
 *              Update and factorize the current diagonal block and test
 *              for non-positive-definiteness.
 *
-               CALL SPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )
+               CALL SPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )
                IF( INFO.NE.0 )
      $            GO TO 30
 
@@ -216,7 +216,7 @@ SUBROUTINE SPOTRF ( UPLO, N, A, LDA, INFO )
 *              Update and factorize the current diagonal block and test
 *              for non-positive-definiteness.
 *
-               CALL SPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )
+               CALL SPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )
                IF( INFO.NE.0 )
      $            GO TO 30
 
diff --git a/SRC/VARIANTS/cholesky/TOP/zpotrf.f b/SRC/VARIANTS/cholesky/TOP/zpotrf.f
@@ -129,7 +129,7 @@ SUBROUTINE ZPOTRF ( UPLO, N, A, LDA, INFO )
       EXTERNAL           LSAME, ILAENV
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           ZGEMM, ZPOTF2, ZHERK, ZTRSM, XERBLA
+      EXTERNAL           XERBLA, ZGEMM, ZHERK, ZPOTRF2, ZTRSM
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -164,7 +164,7 @@ SUBROUTINE ZPOTRF ( UPLO, N, A, LDA, INFO )
 *
 *        Use unblocked code.
 *
-         CALL ZPOTF2( UPLO, N, A, LDA, INFO )
+         CALL ZPOTRF2( UPLO, N, A, LDA, INFO )
       ELSE
 *
 *        Use blocked code.
@@ -189,7 +189,7 @@ SUBROUTINE ZPOTRF ( UPLO, N, A, LDA, INFO )
 *              Update and factorize the current diagonal block and test
 *              for non-positive-definiteness.
 *
-               CALL ZPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )
+               CALL ZPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )
                IF( INFO.NE.0 )
      $            GO TO 30
 
@@ -216,7 +216,7 @@ SUBROUTINE ZPOTRF ( UPLO, N, A, LDA, INFO )
 *              Update and factorize the current diagonal block and test
 *              for non-positive-definiteness.
 *
-               CALL ZPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )
+               CALL ZPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )
                IF( INFO.NE.0 )
      $            GO TO 30
 

Original file line number	Diff line number	Diff line change
`@@ -129,7 +129,7 @@ SUBROUTINE CPOTRF ( UPLO, N, A, LDA, INFO )`
`129`	`129`	`EXTERNAL LSAME, ILAENV`
`130`	`130`	`* ..`
`131`	`131`	`* .. External Subroutines ..`
`132`		`- EXTERNAL CGEMM, CPOTF2, CHERK, CTRSM, XERBLA`
	`132`	`+ EXTERNAL CGEMM, CHERK, CPOTRF2, CTRSM, XERBLA`
`133`	`133`	`* ..`
`134`	`134`	`* .. Intrinsic Functions ..`
`135`	`135`	`INTRINSIC MAX, MIN`
`@@ -164,7 +164,7 @@ SUBROUTINE CPOTRF ( UPLO, N, A, LDA, INFO )`
`164`	`164`	`*`
`165`	`165`	`* Use unblocked code.`
`166`	`166`	`*`
`167`		`- CALL CPOTF2( UPLO, N, A, LDA, INFO )`
	`167`	`+ CALL CPOTRF2( UPLO, N, A, LDA, INFO )`
`168`	`168`	`ELSE`
`169`	`169`	`*`
`170`	`170`	`* Use blocked code.`
`@@ -180,7 +180,7 @@ SUBROUTINE CPOTRF ( UPLO, N, A, LDA, INFO )`
`180`	`180`	`*`
`181`	`181`	`JB = MIN( NB, N-J+1 )`
`182`	`182`
`183`		`- CALL CPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )`
	`183`	`+ CALL CPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )`
`184`	`184`
`185`	`185`	`IF( INFO.NE.0 )`
`186`	`186`	`$ GO TO 30`
`@@ -209,7 +209,7 @@ SUBROUTINE CPOTRF ( UPLO, N, A, LDA, INFO )`
`209`	`209`	`*`
`210`	`210`	`JB = MIN( NB, N-J+1 )`
`211`	`211`
`212`		`- CALL CPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )`
	`212`	`+ CALL CPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )`
`213`	`213`
`214`	`214`	`IF( INFO.NE.0 )`
`215`	`215`	`$ GO TO 30`
Original file line number	Diff line number	Diff line change
`@@ -128,7 +128,7 @@ SUBROUTINE DPOTRF ( UPLO, N, A, LDA, INFO )`
`128`	`128`	`EXTERNAL LSAME, ILAENV`
`129`	`129`	`* ..`
`130`	`130`	`* .. External Subroutines ..`
`131`		`- EXTERNAL DGEMM, DPOTF2, DSYRK, DTRSM, XERBLA`
	`131`	`+ EXTERNAL DGEMM, DPOTRF2, DSYRK, DTRSM, XERBLA`
`132`	`132`	`* ..`
`133`	`133`	`* .. Intrinsic Functions ..`
`134`	`134`	`INTRINSIC MAX, MIN`
`@@ -163,7 +163,7 @@ SUBROUTINE DPOTRF ( UPLO, N, A, LDA, INFO )`
`163`	`163`	`*`
`164`	`164`	`* Use unblocked code.`
`165`	`165`	`*`
`166`		`- CALL DPOTF2( UPLO, N, A, LDA, INFO )`
	`166`	`+ CALL DPOTRF2( UPLO, N, A, LDA, INFO )`
`167`	`167`	`ELSE`
`168`	`168`	`*`
`169`	`169`	`* Use blocked code.`
`@@ -179,7 +179,7 @@ SUBROUTINE DPOTRF ( UPLO, N, A, LDA, INFO )`
`179`	`179`	`*`
`180`	`180`	`JB = MIN( NB, N-J+1 )`
`181`	`181`
`182`		`- CALL DPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )`
	`182`	`+ CALL DPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )`
`183`	`183`
`184`	`184`	`IF( INFO.NE.0 )`
`185`	`185`	`$ GO TO 30`
`@@ -208,7 +208,7 @@ SUBROUTINE DPOTRF ( UPLO, N, A, LDA, INFO )`
`208`	`208`	`*`
`209`	`209`	`JB = MIN( NB, N-J+1 )`
`210`	`210`
`211`		`- CALL DPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )`
	`211`	`+ CALL DPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )`
`212`	`212`
`213`	`213`	`IF( INFO.NE.0 )`
`214`	`214`	`$ GO TO 30`
Original file line number	Diff line number	Diff line change
`@@ -128,7 +128,7 @@ SUBROUTINE SPOTRF ( UPLO, N, A, LDA, INFO )`
`128`	`128`	`EXTERNAL LSAME, ILAENV`
`129`	`129`	`* ..`
`130`	`130`	`* .. External Subroutines ..`
`131`		`- EXTERNAL SGEMM, SPOTF2, SSYRK, STRSM, XERBLA`
	`131`	`+ EXTERNAL SGEMM, SPOTRF2, SSYRK, STRSM, XERBLA`
`132`	`132`	`* ..`
`133`	`133`	`* .. Intrinsic Functions ..`
`134`	`134`	`INTRINSIC MAX, MIN`
`@@ -163,7 +163,7 @@ SUBROUTINE SPOTRF ( UPLO, N, A, LDA, INFO )`
`163`	`163`	`*`
`164`	`164`	`* Use unblocked code.`
`165`	`165`	`*`
`166`		`- CALL SPOTF2( UPLO, N, A, LDA, INFO )`
	`166`	`+ CALL SPOTRF2( UPLO, N, A, LDA, INFO )`
`167`	`167`	`ELSE`
`168`	`168`	`*`
`169`	`169`	`* Use blocked code.`
`@@ -179,7 +179,7 @@ SUBROUTINE SPOTRF ( UPLO, N, A, LDA, INFO )`
`179`	`179`	`*`
`180`	`180`	`JB = MIN( NB, N-J+1 )`
`181`	`181`
`182`		`- CALL SPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )`
	`182`	`+ CALL SPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )`
`183`	`183`
`184`	`184`	`IF( INFO.NE.0 )`
`185`	`185`	`$ GO TO 30`
`@@ -208,7 +208,7 @@ SUBROUTINE SPOTRF ( UPLO, N, A, LDA, INFO )`
`208`	`208`	`*`
`209`	`209`	`JB = MIN( NB, N-J+1 )`
`210`	`210`
`211`		`- CALL SPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )`
	`211`	`+ CALL SPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )`
`212`	`212`
`213`	`213`	`IF( INFO.NE.0 )`
`214`	`214`	`$ GO TO 30`
Original file line number	Diff line number	Diff line change
`@@ -129,7 +129,7 @@ SUBROUTINE ZPOTRF ( UPLO, N, A, LDA, INFO )`
`129`	`129`	`EXTERNAL LSAME, ILAENV`
`130`	`130`	`* ..`
`131`	`131`	`* .. External Subroutines ..`
`132`		`- EXTERNAL ZGEMM, ZPOTF2, ZHERK, ZTRSM, XERBLA`
	`132`	`+ EXTERNAL XERBLA, ZGEMM, ZHERK, ZPOTRF2, ZTRSM`
`133`	`133`	`* ..`
`134`	`134`	`* .. Intrinsic Functions ..`
`135`	`135`	`INTRINSIC MAX, MIN`
`@@ -164,7 +164,7 @@ SUBROUTINE ZPOTRF ( UPLO, N, A, LDA, INFO )`
`164`	`164`	`*`
`165`	`165`	`* Use unblocked code.`
`166`	`166`	`*`
`167`		`- CALL ZPOTF2( UPLO, N, A, LDA, INFO )`
	`167`	`+ CALL ZPOTRF2( UPLO, N, A, LDA, INFO )`
`168`	`168`	`ELSE`
`169`	`169`	`*`
`170`	`170`	`* Use blocked code.`
`@@ -180,7 +180,7 @@ SUBROUTINE ZPOTRF ( UPLO, N, A, LDA, INFO )`
`180`	`180`	`*`
`181`	`181`	`JB = MIN( NB, N-J+1 )`
`182`	`182`
`183`		`- CALL ZPOTF2( 'Upper', JB, A( J, J ), LDA, INFO )`
	`183`	`+ CALL ZPOTRF2( 'Upper', JB, A( J, J ), LDA, INFO )`
`184`	`184`
`185`	`185`	`IF( INFO.NE.0 )`
`186`	`186`	`$ GO TO 30`
`@@ -209,7 +209,7 @@ SUBROUTINE ZPOTRF ( UPLO, N, A, LDA, INFO )`
`209`	`209`	`*`
`210`	`210`	`JB = MIN( NB, N-J+1 )`
`211`	`211`
`212`		`- CALL ZPOTF2( 'Lower', JB, A( J, J ), LDA, INFO )`
	`212`	`+ CALL ZPOTRF2( 'Lower', JB, A( J, J ), LDA, INFO )`
`213`	`213`
`214`	`214`	`IF( INFO.NE.0 )`
`215`	`215`	`$ GO TO 30`