Fixed comments on factorization in refinement routines. Found when coding in C++

diff --git a/SRC/chesvxx.f b/SRC/chesvxx.f
index 9d412f2..ac3486d 100644 (file)
--- a/SRC/chesvxx.f
+++ b/SRC/chesvxx.f
@@ -510,7 +510,7 @@
 *
       IF( NOFACT .OR. EQUIL ) THEN
 *
-*        Compute the LU factorization of A.
+*        Compute the LDL^T or UDU^T factorization of A.
 *
          CALL CLACPY( UPLO, N, N, A, LDA, AF, LDAF )
          CALL CHETRF( UPLO, N, AF, LDAF, IPIV, WORK, 5*MAX(1,N), INFO )

diff --git a/SRC/cposvxx.f b/SRC/cposvxx.f
index 38e8d8f..4c91e23 100644 (file)
--- a/SRC/cposvxx.f
+++ b/SRC/cposvxx.f
@@ -504,7 +504,7 @@
 *
       IF( NOFACT .OR. EQUIL ) THEN
 *
-*        Compute the LU factorization of A.
+*        Compute the Cholesky factorization of A.
 *
          CALL CLACPY( UPLO, N, N, A, LDA, AF, LDAF )
          CALL CPOTRF( UPLO, N, AF, LDAF, INFO )

diff --git a/SRC/csysvxx.f b/SRC/csysvxx.f
index 7693b88..8acce8c 100644 (file)
--- a/SRC/csysvxx.f
+++ b/SRC/csysvxx.f
@@ -515,7 +515,7 @@
 *
       IF( NOFACT .OR. EQUIL ) THEN
 *
-*        Compute the LU factorization of A.
+*        Compute the LDL^T or UDU^T factorization of A.
 *
          CALL CLACPY( UPLO, N, N, A, LDA, AF, LDAF )
          CALL CSYTRF( UPLO, N, AF, LDAF, IPIV, WORK, 5*MAX(1,N), INFO )

diff --git a/SRC/dposvxx.f b/SRC/dposvxx.f
index 7164ae3..da5dddd 100644 (file)
--- a/SRC/dposvxx.f
+++ b/SRC/dposvxx.f
@@ -503,7 +503,7 @@
 *
       IF( NOFACT .OR. EQUIL ) THEN
 *
-*        Compute the LU factorization of A.
+*        Compute the Cholesky factorization of A.
 *
          CALL DLACPY( UPLO, N, N, A, LDA, AF, LDAF )
          CALL DPOTRF( UPLO, N, AF, LDAF, INFO )

diff --git a/SRC/dsysvxx.f b/SRC/dsysvxx.f
index 943559d..5125fd5 100644 (file)
--- a/SRC/dsysvxx.f
+++ b/SRC/dsysvxx.f
@@ -510,7 +510,7 @@
 *
       IF( NOFACT .OR. EQUIL ) THEN
 *
-*        Compute the LU factorization of A.
+*        Compute the LDL^T or UDU^T factorization of A.
 *
          CALL DLACPY( UPLO, N, N, A, LDA, AF, LDAF )
          CALL DSYTRF( UPLO, N, AF, LDAF, IPIV, WORK, 5*MAX(1,N), INFO )

diff --git a/SRC/sposvxx.f b/SRC/sposvxx.f
index 21b1c20..aa6b750 100644 (file)
--- a/SRC/sposvxx.f
+++ b/SRC/sposvxx.f
@@ -506,7 +506,7 @@
 *
       IF( NOFACT .OR. EQUIL ) THEN
 *
-*        Compute the LU factorization of A.
+*        Compute the Cholesky factorization of A.
 *
          CALL SLACPY( UPLO, N, N, A, LDA, AF, LDAF )
          CALL SPOTRF( UPLO, N, AF, LDAF, INFO )

diff --git a/SRC/ssysvxx.f b/SRC/ssysvxx.f
index db43aa3..50d4850 100644 (file)
--- a/SRC/ssysvxx.f
+++ b/SRC/ssysvxx.f
@@ -513,7 +513,7 @@
 *
       IF( NOFACT .OR. EQUIL ) THEN
 *
-*        Compute the LU factorization of A.
+*        Compute the LDL^T or UDU^T factorization of A.
 *
          CALL SLACPY( UPLO, N, N, A, LDA, AF, LDAF )
          CALL SSYTRF( UPLO, N, AF, LDAF, IPIV, WORK, 5*MAX(1,N), INFO )

diff --git a/SRC/zhesvxx.f b/SRC/zhesvxx.f
index b6412ca..deb1190 100644 (file)
--- a/SRC/zhesvxx.f
+++ b/SRC/zhesvxx.f
@@ -507,7 +507,7 @@
 *
       IF( NOFACT .OR. EQUIL ) THEN
 *
-*        Compute the LU factorization of A.
+*        Compute the LDL^T or UDU^T factorization of A.
 *
          CALL ZLACPY( UPLO, N, N, A, LDA, AF, LDAF )
          CALL ZHETRF( UPLO, N, AF, LDAF, IPIV, WORK, 5*MAX(1,N), INFO )

diff --git a/SRC/zposvxx.f b/SRC/zposvxx.f
index bd202ff..7a60fe5 100644 (file)
--- a/SRC/zposvxx.f
+++ b/SRC/zposvxx.f
@@ -501,7 +501,7 @@
 *
       IF( NOFACT .OR. EQUIL ) THEN
 *
-*        Compute the LU factorization of A.
+*        Compute the Cholesky factorization of A.
 *
          CALL ZLACPY( UPLO, N, N, A, LDA, AF, LDAF )
          CALL ZPOTRF( UPLO, N, AF, LDAF, INFO )

diff --git a/SRC/zsysvxx.f b/SRC/zsysvxx.f
index f5464b2..f8e5ced 100644 (file)
--- a/SRC/zsysvxx.f
+++ b/SRC/zsysvxx.f
@@ -512,7 +512,7 @@
 *
       IF( NOFACT .OR. EQUIL ) THEN
 *
-*        Compute the LU factorization of A.
+*        Compute the LDL^T or UDU^T factorization of A.
 *
          CALL ZLACPY( UPLO, N, N, A, LDA, AF, LDAF )
          CALL ZSYTRF( UPLO, N, AF, LDAF, IPIV, WORK, 5*MAX(1,N), INFO )