fixing lwork for Aasen's (lwork > 0).

diff --git a/SRC/chesv_aa.f b/SRC/chesv_aa.f
index 621ae1f..642c993 100644 (file)
--- a/SRC/chesv_aa.f
+++ b/SRC/chesv_aa.f
@@ -126,9 +126,9 @@
 *> \param[in] LWORK
 *> \verbatim
 *>          LWORK is INTEGER
-*>          The length of WORK.  LWORK >= 1, and for best performance
-*>          LWORK >= max(1,N*NB), where NB is the optimal blocksize for
-*>          CHETRF.
+*>          The length of WORK.  LWORK >= MAX(1,2*N,3*N-2), and for best 
+*>          performance LWORK >= MAX(1,N*NB), where NB is the optimal
+*>          blocksize for CHETRF.
 *>          for LWORK < N, TRS will be done with Level BLAS 2
 *>          for LWORK >= N, TRS will be done with Level BLAS 3
 *>

diff --git a/SRC/chetrs_aa.f b/SRC/chetrs_aa.f
index 11f19d6..ef55c8f 100644 (file)
--- a/SRC/chetrs_aa.f
+++ b/SRC/chetrs_aa.f
@@ -104,7 +104,7 @@
 *>
 *> \param[in] LWORK
 *> \verbatim
-*>          LWORK is INTEGER, LWORK >= 3*N-2.
+*>          LWORK is INTEGER, LWORK >= MAX(1,3*N-2).
 *>
 *> \param[out] INFO
 *> \verbatim
@@ -179,7 +179,7 @@
          INFO = -5
       ELSE IF( LDB.LT.MAX( 1, N ) ) THEN
          INFO = -8
-      ELSE IF( LWORK.LT.(3*N-2) .AND. .NOT.LQUERY ) THEN
+      ELSE IF( LWORK.LT.MAX( 1, 3*N-2 ) .AND. .NOT.LQUERY ) THEN
          INFO = -10
       END IF
       IF( INFO.NE.0 ) THEN

diff --git a/SRC/csytrf_aa.f b/SRC/csytrf_aa.f
index 20aa5ba..c6b7613 100644 (file)
--- a/SRC/csytrf_aa.f
+++ b/SRC/csytrf_aa.f
@@ -101,7 +101,7 @@
 *> \param[in] LWORK
 *> \verbatim
 *>          LWORK is INTEGER
-*>          The length of WORK.  LWORK >=2*N. For optimum performance
+*>          The length of WORK. LWORK >= MAX(1,2*N). For optimum performance
 *>          LWORK >= N*(1+NB), where NB is the optimal blocksize.
 *>
 *>          If LWORK = -1, then a workspace query is assumed; the routine
@@ -191,7 +191,7 @@
          INFO = -2
       ELSE IF( LDA.LT.MAX( 1, N ) ) THEN
          INFO = -4
-      ELSE IF( LWORK.LT.( 2*N ) .AND. .NOT.LQUERY ) THEN
+      ELSE IF( LWORK.LT.MAX( 1, 2*N ) .AND. .NOT.LQUERY ) THEN
          INFO = -7
       END IF
 *

diff --git a/SRC/csytrs_aa.f b/SRC/csytrs_aa.f
index 5f21078..fd75ba4 100644 (file)
--- a/SRC/csytrs_aa.f
+++ b/SRC/csytrs_aa.f
@@ -104,7 +104,7 @@
 *>
 *> \param[in] LWORK
 *> \verbatim
-*>          LWORK is INTEGER, LWORK >= 3*N-2.
+*>          LWORK is INTEGER, LWORK >= MAX(1,3*N-2).
 *>
 *> \param[out] INFO
 *> \verbatim
@@ -179,7 +179,7 @@
          INFO = -5
       ELSE IF( LDB.LT.MAX( 1, N ) ) THEN
          INFO = -8
-      ELSE IF( LWORK.LT.(3*N-2) .AND. .NOT.LQUERY ) THEN
+      ELSE IF( LWORK.LT.MAX( 1, 3*N-2 ) .AND. .NOT.LQUERY ) THEN
          INFO = -10
       END IF
       IF( INFO.NE.0 ) THEN

diff --git a/SRC/dsysv_aa.f b/SRC/dsysv_aa.f
index cb6b28e..055097f 100644 (file)
--- a/SRC/dsysv_aa.f
+++ b/SRC/dsysv_aa.f
@@ -126,8 +126,8 @@
 *> \param[in] LWORK
 *> \verbatim
 *>          LWORK is INTEGER
-*>          The length of WORK.  LWORK >= MAX(2*N, 3*N-2), and for
-*>          the best performance, LWORK >= max(1,N*NB), where NB is
+*>          The length of WORK.  LWORK >= MAX(1,2*N,3*N-2), and for
+*>          the best performance, LWORK >= MAX(1,N*NB), where NB is
 *>          the optimal blocksize for DSYTRF_AA.
 *>
 *>          If LWORK = -1, then a workspace query is assumed; the routine

diff --git a/SRC/dsytrf_aa.f b/SRC/dsytrf_aa.f
index 07919a2..0e03880 100644 (file)
--- a/SRC/dsytrf_aa.f
+++ b/SRC/dsytrf_aa.f
@@ -101,7 +101,7 @@
 *> \param[in] LWORK
 *> \verbatim
 *>          LWORK is INTEGER
-*>          The length of WORK.  LWORK >=2*N. For optimum performance
+*>          The length of WORK.  LWORK >= MAX(1,2*N). For optimum performance
 *>          LWORK >= N*(1+NB), where NB is the optimal blocksize.
 *>
 *>          If LWORK = -1, then a workspace query is assumed; the routine
@@ -191,7 +191,7 @@
          INFO = -2
       ELSE IF( LDA.LT.MAX( 1, N ) ) THEN
          INFO = -4
-      ELSE IF( LWORK.LT.( 2*N ) .AND. .NOT.LQUERY ) THEN
+      ELSE IF( LWORK.LT.MAX( 1, 2*N ) .AND. .NOT.LQUERY ) THEN
          INFO = -7
       END IF
 *

diff --git a/SRC/dsytrs_aa.f b/SRC/dsytrs_aa.f
index ddb9d3f..6c56b91 100644 (file)
--- a/SRC/dsytrs_aa.f
+++ b/SRC/dsytrs_aa.f
@@ -104,7 +104,7 @@
 *>
 *> \param[in] LWORK
 *> \verbatim
-*>          LWORK is INTEGER, LWORK >= 3*N-2.
+*>          LWORK is INTEGER, LWORK >= MAX(1,3*N-2).
 *>
 *> \param[out] INFO
 *> \verbatim
@@ -179,7 +179,7 @@
          INFO = -5
       ELSE IF( LDB.LT.MAX( 1, N ) ) THEN
          INFO = -8
-      ELSE IF( LWORK.LT.(3*N-2) .AND. .NOT.LQUERY ) THEN
+      ELSE IF( LWORK.LT.MAX( 1, 3*N-2 ) .AND. .NOT.LQUERY ) THEN
          INFO = -10
       END IF
       IF( INFO.NE.0 ) THEN

diff --git a/SRC/ssysv_aa.f b/SRC/ssysv_aa.f
index 3dcc7cb..d61a346 100644 (file)
--- a/SRC/ssysv_aa.f
+++ b/SRC/ssysv_aa.f
@@ -126,8 +126,8 @@
 *> \param[in] LWORK
 *> \verbatim
 *>          LWORK is INTEGER
-*>          The length of WORK.  LWORK >= MAX(2*N, 3*N-2), and for
-*>          the best performance, LWORK >= max(1,N*NB), where NB is
+*>          The length of WORK.  LWORK >= MAX(1,2*N,3*N-2), and for
+*>          the best performance, LWORK >= MAX(1,N*NB), where NB is
 *>          the optimal blocksize for SSYTRF_AA.
 *>
 *>          If LWORK = -1, then a workspace query is assumed; the routine

diff --git a/SRC/ssytrf_aa.f b/SRC/ssytrf_aa.f
index 13498c9..a22ff05 100644 (file)
--- a/SRC/ssytrf_aa.f
+++ b/SRC/ssytrf_aa.f
@@ -101,7 +101,7 @@
 *> \param[in] LWORK
 *> \verbatim
 *>          LWORK is INTEGER
-*>          The length of WORK.  LWORK >=2*N. For optimum performance
+*>          The length of WORK.  LWORK >= MAX(1,2*N). For optimum performance
 *>          LWORK >= N*(1+NB), where NB is the optimal blocksize.
 *>
 *>          If LWORK = -1, then a workspace query is assumed; the routine
@@ -191,7 +191,7 @@
          INFO = -2
       ELSE IF( LDA.LT.MAX( 1, N ) ) THEN
          INFO = -4
-      ELSE IF( LWORK.LT.( 2*N ) .AND. .NOT.LQUERY ) THEN
+      ELSE IF( LWORK.LT.MAX( 1, 2*N ) .AND. .NOT.LQUERY ) THEN
          INFO = -7
       END IF
 *

diff --git a/SRC/ssytrs_aa.f b/SRC/ssytrs_aa.f
index 06c793a..911016e 100644 (file)
--- a/SRC/ssytrs_aa.f
+++ b/SRC/ssytrs_aa.f
@@ -104,7 +104,7 @@
 *>
 *> \param[in] LWORK
 *> \verbatim
-*>          LWORK is INTEGER, LWORK >= 3*N-2.
+*>          LWORK is INTEGER, LWORK >= MAX(1,3*N-2).
 *>
 *> \param[out] INFO
 *> \verbatim
@@ -179,7 +179,7 @@
          INFO = -5
       ELSE IF( LDB.LT.MAX( 1, N ) ) THEN
          INFO = -8
-      ELSE IF( LWORK.LT.(3*N-2) .AND. .NOT.LQUERY ) THEN
+      ELSE IF( LWORK.LT.MAX( 1, 3*N-2 ) .AND. .NOT.LQUERY ) THEN
          INFO = -10
       END IF
       IF( INFO.NE.0 ) THEN

diff --git a/SRC/zhesv_aa.f b/SRC/zhesv_aa.f
index 285c971..dffa475 100644 (file)
--- a/SRC/zhesv_aa.f
+++ b/SRC/zhesv_aa.f
@@ -126,9 +126,9 @@
 *> \param[in] LWORK
 *> \verbatim
 *>          LWORK is INTEGER
-*>          The length of WORK.  LWORK >= 1, and for best performance
-*>          LWORK >= max(1,N*NB), where NB is the optimal blocksize for
-*>          ZHETRF.
+*>          The length of WORK.  LWORK >= MAX(1,2*N,3*N-2), and for best 
+*>          performance LWORK >= max(1,N*NB), where NB is the optimal
+*>          blocksize for ZHETRF.
 *>          for LWORK < N, TRS will be done with Level BLAS 2
 *>          for LWORK >= N, TRS will be done with Level BLAS 3
 *>

diff --git a/SRC/zhetrf_aa.f b/SRC/zhetrf_aa.f
index adaae91..a3dd095 100644 (file)
--- a/SRC/zhetrf_aa.f
+++ b/SRC/zhetrf_aa.f
@@ -101,7 +101,7 @@
 *> \param[in] LWORK
 *> \verbatim
 *>          LWORK is INTEGER
-*>          The length of WORK.  LWORK >= 2*N. For optimum performance
+*>          The length of WORK. LWORK >= MAX(1,2*N). For optimum performance
 *>          LWORK >= N*(1+NB), where NB is the optimal blocksize.
 *>
 *>          If LWORK = -1, then a workspace query is assumed; the routine
@@ -191,7 +191,7 @@
          INFO = -2
       ELSE IF( LDA.LT.MAX( 1, N ) ) THEN
          INFO = -4
-      ELSE IF( LWORK.LT.( 2*N ) .AND. .NOT.LQUERY ) THEN
+      ELSE IF( LWORK.LT.MAX( 1, 2*N ) .AND. .NOT.LQUERY ) THEN
          INFO = -7
       END IF
 *

diff --git a/SRC/zhetrs_aa.f b/SRC/zhetrs_aa.f
index e5487ff..0a02b8a 100644 (file)
--- a/SRC/zhetrs_aa.f
+++ b/SRC/zhetrs_aa.f
@@ -105,7 +105,7 @@
 *>
 *> \param[in] LWORK
 *> \verbatim
-*>          LWORK is INTEGER, LWORK >= 3*N-2.
+*>          LWORK is INTEGER, LWORK >= MAX(1,3*N-2).
 *>
 *> \param[out] INFO
 *> \verbatim
@@ -180,7 +180,7 @@
          INFO = -5
       ELSE IF( LDB.LT.MAX( 1, N ) ) THEN
          INFO = -8
-      ELSE IF( LWORK.LT.(3*N-2) .AND. .NOT.LQUERY ) THEN
+      ELSE IF( LWORK.LT.MAX( 1, 3*N-2 ) .AND. .NOT.LQUERY ) THEN
          INFO = -10
       END IF
       IF( INFO.NE.0 ) THEN

diff --git a/SRC/zsysv_aa.f b/SRC/zsysv_aa.f
index f1d1f76..6c76714 100644 (file)
--- a/SRC/zsysv_aa.f
+++ b/SRC/zsysv_aa.f
@@ -126,8 +126,8 @@
 *> \param[in] LWORK
 *> \verbatim
 *>          LWORK is INTEGER
-*>          The length of WORK.  LWORK >= MAX(2*N, 3*N-2), and for
-*>          the best performance, LWORK >= max(1,N*NB), where NB is
+*>          The length of WORK.  LWORK >= MAX(1,2*N,3*N-2), and for
+*>          the best performance, LWORK >= MAX(1,N*NB), where NB is
 *>          the optimal blocksize for ZSYTRF_AA.
 *>
 *>          If LWORK = -1, then a workspace query is assumed; the routine

diff --git a/SRC/zsytrf_aa.f b/SRC/zsytrf_aa.f
index 4e9851a..f82e513 100644 (file)
--- a/SRC/zsytrf_aa.f
+++ b/SRC/zsytrf_aa.f
@@ -101,7 +101,7 @@
 *> \param[in] LWORK
 *> \verbatim
 *>          LWORK is INTEGER
-*>          The length of WORK.  LWORK >=2*N. For optimum performance
+*>          The length of WORK. LWORK >=MAX(1,2*N). For optimum performance
 *>          LWORK >= N*(1+NB), where NB is the optimal blocksize.
 *>
 *>          If LWORK = -1, then a workspace query is assumed; the routine
@@ -191,7 +191,7 @@
          INFO = -2
       ELSE IF( LDA.LT.MAX( 1, N ) ) THEN
          INFO = -4
-      ELSE IF( LWORK.LT.( 2*N ) .AND. .NOT.LQUERY ) THEN
+      ELSE IF( LWORK.LT.MAX( 1, 2*N ) .AND. .NOT.LQUERY ) THEN
          INFO = -7
       END IF
 *

diff --git a/SRC/zsytrs_aa.f b/SRC/zsytrs_aa.f
index b9a2ab8..cae83a7 100644 (file)
--- a/SRC/zsytrs_aa.f
+++ b/SRC/zsytrs_aa.f
@@ -104,7 +104,7 @@
 *>
 *> \param[in] LWORK
 *> \verbatim
-*>          LWORK is INTEGER, LWORK >= 3*N-2.
+*>          LWORK is INTEGER, LWORK >= MAX(1,3*N-2).
 *>
 *> \param[out] INFO
 *> \verbatim
@@ -179,7 +179,7 @@
          INFO = -5
       ELSE IF( LDB.LT.MAX( 1, N ) ) THEN
          INFO = -8
-      ELSE IF( LWORK.LT.(3*N-2) .AND. .NOT.LQUERY ) THEN
+      ELSE IF( LWORK.LT.MAX( 1, 3*N-2 ) .AND. .NOT.LQUERY ) THEN
          INFO = -10
       END IF
       IF( INFO.NE.0 ) THEN