Converge >= and > usage in docs and implementation of GESVJ

angsch · Sep 17, 2024 · 18d002a · 18d002a
1 parent 9c0ef66
commit 18d002a
Show file tree

Hide file tree

Showing 4 changed files with 14 additions and 14 deletions.
diff --git a/SRC/cgesvj.f b/SRC/cgesvj.f
@@ -101,7 +101,7 @@
 *> \param[in] M
 *> \verbatim
 *>          M is INTEGER
-*>          The number of rows of the input matrix A. 1/SLAMCH('E') > M >= 0.
+*>          The number of rows of the input matrix A. 1/SLAMCH('E') >= M >= 0.
 *> \endverbatim
 *>
 *> \param[in] N
@@ -217,7 +217,7 @@
 *>          LWORK >= 1, if MIN(M,N) = 0, and LWORK >= M+N, otherwise.
 *>
 *>          If on entry LWORK = -1, then a workspace query is assumed and
-*>          no computation is done; CWORK(1) is set to the minial (and optimal)
+*>          no computation is done; CWORK(1) is set to the minimal (and optimal)
 *>          length of CWORK.
 *> \endverbatim
 *>
@@ -258,7 +258,7 @@
 *>         LRWORK >= 1, if MIN(M,N) = 0, and LRWORK >= MAX(6,N), otherwise
 *>
 *>         If on entry LRWORK = -1, then a workspace query is assumed and
-*>         no computation is done; RWORK(1) is set to the minial (and optimal)
+*>         no computation is done; RWORK(1) is set to the minimal (and optimal)
 *>         length of RWORK.
 *> \endverbatim
 *>
@@ -459,7 +459,7 @@ SUBROUTINE CGESVJ( JOBA, JOBU, JOBV, M, N, A, LDA, SVA, MV, V,
       ELSE IF( ( RSVEC .AND. ( LDV.LT.N ) ) .OR.
      $          ( APPLV .AND. ( LDV.LT.MV ) ) ) THEN
          INFO = -11
-      ELSE IF( UCTOL .AND. ( RWORK( 1 ).LE.ONE ) ) THEN
+      ELSE IF( UCTOL .AND. ( RWORK( 1 ).LT.ONE ) ) THEN
          INFO = -12
       ELSE IF( LWORK.LT.LWMIN .AND. ( .NOT.LQUERY ) ) THEN
          INFO = -13

diff --git a/SRC/dgesvj.f b/SRC/dgesvj.f
@@ -103,7 +103,7 @@
 *> \param[in] M
 *> \verbatim
 *>          M is INTEGER
-*>          The number of rows of the input matrix A. 1/DLAMCH('E') > M >= 0.
+*>          The number of rows of the input matrix A. 1/DLAMCH('E') >= M >= 0.
 *> \endverbatim
 *>
 *> \param[in] N
@@ -243,7 +243,7 @@
 *>          LWORK >= 1, if MIN(M,N) = 0, and LWORK >= MAX(6,M+N), otherwise.
 *>
 *>          If on entry LWORK = -1, then a workspace query is assumed and
-*>          no computation is done; WORK(1) is set to the minial (and optimal)
+*>          no computation is done; WORK(1) is set to the minimal (and optimal)
 *>          length of WORK.
 *> \endverbatim
 *>
@@ -442,7 +442,7 @@ SUBROUTINE DGESVJ( JOBA, JOBU, JOBV, M, N, A, LDA, SVA, MV, V,
       ELSE IF( ( RSVEC .AND. ( LDV.LT.N ) ) .OR.
      $         ( APPLV .AND. ( LDV.LT.MV ) ) ) THEN
          INFO = -11
-      ELSE IF( UCTOL .AND. ( WORK( 1 ).LE.ONE ) ) THEN
+      ELSE IF( UCTOL .AND. ( WORK( 1 ).LT.ONE ) ) THEN
          INFO = -12
       ELSE IF( LWORK.LT.LWMIN .AND. ( .NOT.LQUERY ) ) THEN
          INFO = -13

diff --git a/SRC/sgesvj.f b/SRC/sgesvj.f
@@ -103,7 +103,7 @@
 *> \param[in] M
 *> \verbatim
 *>          M is INTEGER
-*>          The number of rows of the input matrix A. 1/SLAMCH('E') > M >= 0.
+*>          The number of rows of the input matrix A. 1/SLAMCH('E') >= M >= 0.
 *> \endverbatim
 *>
 *> \param[in] N
@@ -243,7 +243,7 @@
 *>          LWORK >= 1, if MIN(M,N) = 0, and LWORK >= MAX(6,M+N), otherwise.
 *>
 *>          If on entry LWORK = -1, then a workspace query is assumed and
-*>          no computation is done; WORK(1) is set to the minial (and optimal)
+*>          no computation is done; WORK(1) is set to the minimal (and optimal)
 *>          length of WORK.
 *> \endverbatim
 *>
@@ -428,7 +428,7 @@ SUBROUTINE SGESVJ( JOBA, JOBU, JOBV, M, N, A, LDA, SVA, MV, V,
       ELSE IF( ( RSVEC .AND. ( LDV.LT.N ) ) .OR.
      $         ( APPLV .AND. ( LDV.LT.MV ) ) ) THEN
          INFO = -11
-      ELSE IF( UCTOL .AND. ( WORK( 1 ).LE.ONE ) ) THEN
+      ELSE IF( UCTOL .AND. ( WORK( 1 ).LT.ONE ) ) THEN
          INFO = -12
       ELSE IF( LWORK.LT.LWMIN .AND. ( .NOT.LQUERY ) ) THEN
          INFO = -13

diff --git a/SRC/zgesvj.f b/SRC/zgesvj.f
@@ -101,7 +101,7 @@
 *> \param[in] M
 *> \verbatim
 *>          M is INTEGER
-*>          The number of rows of the input matrix A. 1/DLAMCH('E') > M >= 0.
+*>          The number of rows of the input matrix A. 1/DLAMCH('E') >= M >= 0.
 *> \endverbatim
 *>
 *> \param[in] N
@@ -217,7 +217,7 @@
 *>          LWORK >= 1, if MIN(M,N) = 0, and LWORK >= M+N, otherwise.
 *>
 *>          If on entry LWORK = -1, then a workspace query is assumed and
-*>          no computation is done; CWORK(1) is set to the minial (and optimal)
+*>          no computation is done; CWORK(1) is set to the minimal (and optimal)
 *>          length of CWORK.
 *> \endverbatim
 *>
@@ -258,7 +258,7 @@
 *>         LRWORK >= 1, if MIN(M,N) = 0, and LRWORK >= MAX(6,N), otherwise.
 *>
 *>         If on entry LRWORK = -1, then a workspace query is assumed and
-*>         no computation is done; RWORK(1) is set to the minial (and optimal)
+*>         no computation is done; RWORK(1) is set to the minimal (and optimal)
 *>         length of RWORK.
 *> \endverbatim
 *>
@@ -460,7 +460,7 @@ SUBROUTINE ZGESVJ( JOBA, JOBU, JOBV, M, N, A, LDA, SVA, MV, V,
       ELSE IF( ( RSVEC .AND. ( LDV.LT.N ) ) .OR.
      $          ( APPLV .AND. ( LDV.LT.MV ) ) ) THEN
          INFO = -11
-      ELSE IF( UCTOL .AND. ( RWORK( 1 ).LE.ONE ) ) THEN
+      ELSE IF( UCTOL .AND. ( RWORK( 1 ).LT.ONE ) ) THEN
          INFO = -12
       ELSE IF( LWORK.LT.LWMIN .AND. ( .NOT.LQUERY ) ) THEN
          INFO = -13