See: http://icl.utk.edu/lapack-forum/viewtopic.php?t=4830

The variables NDIV, NITER and NFAIL are INOUT.  These variables have to be
initialized (e.g., in LASQ2) before callin LASQ3.  LASQ3 increments these
counters.  These counters do not play a role in the algorithm; they just give
statistics.) Thanks to Tracey Brendan for the report.

diff --git a/SRC/dlasq3.f b/SRC/dlasq3.f
index 412c8ca..3ad4856 100644 (file)
--- a/SRC/dlasq3.f
+++ b/SRC/dlasq3.f
@@ -97,22 +97,22 @@
 *>         Maximum value of q.
 *> \endverbatim
 *>
-*> \param[out] NFAIL
+*> \param[in,out] NFAIL
 *> \verbatim
 *>          NFAIL is INTEGER
-*>         Number of times shift was too big.
+*>         Increment NFAIL by 1 each time the shift was too big.
 *> \endverbatim
 *>
-*> \param[out] ITER
+*> \param[in,out] ITER
 *> \verbatim
 *>          ITER is INTEGER
-*>         Number of iterations.
+*>         Increment ITER by 1 for each iteration.
 *> \endverbatim
 *>
-*> \param[out] NDIV
+*> \param[in,out] NDIV
 *> \verbatim
 *>          NDIV is INTEGER
-*>         Number of divisions.
+*>         Increment NDIV by 1 for each division.
 *> \endverbatim
 *>
 *> \param[in] IEEE

diff --git a/SRC/slasq3.f b/SRC/slasq3.f
index c7e7fdc..e310b67 100644 (file)
--- a/SRC/slasq3.f
+++ b/SRC/slasq3.f
@@ -97,22 +97,22 @@
 *>         Maximum value of q.
 *> \endverbatim
 *>
-*> \param[out] NFAIL
+*> \param[in,out] NFAIL
 *> \verbatim
 *>          NFAIL is INTEGER
-*>         Number of times shift was too big.
+*>         Increment NFAIL by 1 each time the shift was too big.
 *> \endverbatim
 *>
-*> \param[out] ITER
+*> \param[in,out] ITER
 *> \verbatim
 *>          ITER is INTEGER
-*>         Number of iterations.
+*>         Increment ITER by 1 for each iteration.
 *> \endverbatim
 *>
-*> \param[out] NDIV
+*> \param[in,out] NDIV
 *> \verbatim
 *>          NDIV is INTEGER
-*>         Number of divisions.
+*>         Increment NDIV by 1 for each division.
 *> \endverbatim
 *>
 *> \param[in] IEEE