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 e67ea96..8fe22ab 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 6ae1a20..a33bc37 100644 (file)
--- a/SRC/dgetc2.f
+++ b/SRC/dgetc2.f
@@ -204,6 +204,11 @@
          A( N, N ) = SMIN
       END IF
 *
+*     Set last pivots to N
+*
+      IPIV( N ) = N
+      JPIV( N ) = N
+*
       RETURN
 *
 *     End of DGETC2

diff --git a/SRC/sgetc2.f b/SRC/sgetc2.f
index 7e9d485..8d688de 100644 (file)
--- a/SRC/sgetc2.f
+++ b/SRC/sgetc2.f
@@ -204,6 +204,11 @@
          A( N, N ) = SMIN
       END IF
 *
+*     Set last pivots to N
+*
+      IPIV( N ) = N
+      JPIV( N ) = N
+*
       RETURN
 *
 *     End of SGETC2

diff --git a/SRC/zgetc2.f b/SRC/zgetc2.f
index 277ac7c..c205c75 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