added code to set the last pivots to n: ipiv(n)=jpiv(n)=n

diff --git a/SRC/cgetc2.f b/SRC/cgetc2.f
index e67ea96946e6bc30d4c193f79f395e7f0d21b21b..8fe22abce2ddab1f466742ea84cc1ef4b2e71881 100644 (file)
--- a/SRC/cgetc2.f
+++ b/SRC/cgetc2.f
@@ -203,6 +203,12 @@
          INFO = N
          A( N, N ) = CMPLX( SMIN, ZERO )
       END IF
+*
+*     Set last pivots to N
+*
+      IPIV( N ) = N
+      JPIV( N ) = N
+*
       RETURN
 *
 *     End of CGETC2

diff --git a/SRC/dgetc2.f b/SRC/dgetc2.f
index 6ae1a204b3f89ac76860c69367fd6be8f2b131f5..a33bc374839936462e71a588bcccd43387a6f72a 100644 (file)
--- a/SRC/dgetc2.f
+++ b/SRC/dgetc2.f
@@ -203,6 +203,11 @@
          INFO = N
          A( N, N ) = SMIN
       END IF
+*
+*     Set last pivots to N
+*
+      IPIV( N ) = N
+      JPIV( N ) = N
 *
       RETURN
 *

diff --git a/SRC/sgetc2.f b/SRC/sgetc2.f
index 7e9d485072c6b1bbf6f73f100f6931d238242241..8d688decb74f8345d0675797a149758a4e7d1c05 100644 (file)
--- a/SRC/sgetc2.f
+++ b/SRC/sgetc2.f
@@ -203,6 +203,11 @@
          INFO = N
          A( N, N ) = SMIN
       END IF
+*
+*     Set last pivots to N
+*
+      IPIV( N ) = N
+      JPIV( N ) = N
 *
       RETURN
 *

diff --git a/SRC/zgetc2.f b/SRC/zgetc2.f
index 277ac7ce225cbcb8659a0b3d4390959d4e96ca8f..c205c7561562f42c5d8627dfd533153524214d51 100644 (file)
--- a/SRC/zgetc2.f
+++ b/SRC/zgetc2.f
@@ -203,6 +203,12 @@
          INFO = N
          A( N, N ) = DCMPLX( SMIN, ZERO )
       END IF
+*
+*     Set last pivots to N
+*
+      IPIV( N ) = N
+      JPIV( N ) = N
+*
       RETURN
 *
 *     End of ZGETC2