summaryrefslogtreecommitdiff
path: root/SRC
diff options
context:
space:
mode:
Diffstat (limited to 'SRC')
-rw-r--r--SRC/cgbrfsx.f4
-rw-r--r--SRC/cgbsvx.f10
-rw-r--r--SRC/cgbsvxx.f10
-rw-r--r--SRC/cgesvx.f10
-rw-r--r--SRC/cgesvxx.f10
-rw-r--r--SRC/cgtsvx.f10
-rw-r--r--SRC/cherfsx.f2
-rw-r--r--SRC/chesvx.f4
-rw-r--r--SRC/chesvxx.f8
-rw-r--r--SRC/chpsvx.f4
-rw-r--r--SRC/clatbs.f2
-rw-r--r--SRC/clatps.f2
-rw-r--r--SRC/clatrs.f2
-rw-r--r--SRC/cpbsvx.f6
-rw-r--r--SRC/cporfsx.f2
-rw-r--r--SRC/cposvx.f6
-rw-r--r--SRC/cposvxx.f6
-rw-r--r--SRC/cppsvx.f6
-rw-r--r--SRC/cptsvx.f4
-rw-r--r--SRC/cspsvx.f4
-rw-r--r--SRC/csyrfsx.f2
-rw-r--r--SRC/csysvx.f4
-rw-r--r--SRC/csysvxx.f8
-rw-r--r--SRC/dgbrfsx.f4
-rw-r--r--SRC/dgbsvx.f10
-rw-r--r--SRC/dgbsvxx.f10
-rw-r--r--SRC/dgesvx.f10
-rw-r--r--SRC/dgesvxx.f10
-rw-r--r--SRC/dgtsvx.f10
-rw-r--r--SRC/dlatbs.f2
-rw-r--r--SRC/dlatps.f2
-rw-r--r--SRC/dlatrs.f2
-rw-r--r--SRC/dpbsvx.f6
-rw-r--r--SRC/dporfsx.f2
-rw-r--r--SRC/dposvx.f6
-rw-r--r--SRC/dposvxx.f6
-rw-r--r--SRC/dppsvx.f6
-rw-r--r--SRC/dptsvx.f4
-rw-r--r--SRC/dspsvx.f4
-rw-r--r--SRC/dsyrfsx.f2
-rw-r--r--SRC/dsysvx.f4
-rw-r--r--SRC/dsysvxx.f8
-rw-r--r--SRC/sgbrfsx.f4
-rw-r--r--SRC/sgbsvx.f10
-rw-r--r--SRC/sgbsvxx.f10
-rw-r--r--SRC/sgesvx.f10
-rw-r--r--SRC/sgesvxx.f10
-rw-r--r--SRC/sgtsvx.f10
-rw-r--r--SRC/slatbs.f2
-rw-r--r--SRC/slatps.f2
-rw-r--r--SRC/slatrs.f2
-rw-r--r--SRC/spbsvx.f6
-rw-r--r--SRC/sporfsx.f2
-rw-r--r--SRC/sposvx.f6
-rw-r--r--SRC/sposvxx.f6
-rw-r--r--SRC/sppsvx.f6
-rw-r--r--SRC/sptsvx.f4
-rw-r--r--SRC/sspsvx.f4
-rw-r--r--SRC/ssyrfsx.f2
-rw-r--r--SRC/ssysvx.f4
-rw-r--r--SRC/ssysvxx.f8
-rw-r--r--SRC/zgbrfsx.f4
-rw-r--r--SRC/zgbsvx.f10
-rw-r--r--SRC/zgbsvxx.f10
-rw-r--r--SRC/zgesvx.f10
-rw-r--r--SRC/zgesvxx.f10
-rw-r--r--SRC/zgtsvx.f10
-rw-r--r--SRC/zherfsx.f2
-rw-r--r--SRC/zhesvx.f4
-rw-r--r--SRC/zhesvxx.f8
-rw-r--r--SRC/zhpsvx.f4
-rw-r--r--SRC/zlatbs.f2
-rw-r--r--SRC/zlatps.f2
-rw-r--r--SRC/zlatrs.f2
-rw-r--r--SRC/zpbsvx.f6
-rw-r--r--SRC/zporfsx.f2
-rw-r--r--SRC/zposvx.f6
-rw-r--r--SRC/zposvxx.f6
-rw-r--r--SRC/zppsvx.f6
-rw-r--r--SRC/zptsvx.f4
-rw-r--r--SRC/zspsvx.f4
-rw-r--r--SRC/zsyrfsx.f2
-rw-r--r--SRC/zsysvx.f4
-rw-r--r--SRC/zsysvxx.f8
84 files changed, 234 insertions, 234 deletions
diff --git a/SRC/cgbrfsx.f b/SRC/cgbrfsx.f
index a28ddd4f..252da0d9 100644
--- a/SRC/cgbrfsx.f
+++ b/SRC/cgbrfsx.f
@@ -159,7 +159,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) REAL array, dimension (N)
+*> R is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -177,7 +177,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) REAL array, dimension (N)
+*> C is REAL array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/cgbsvx.f b/SRC/cgbsvx.f
index abf8aa74..0086f221 100644
--- a/SRC/cgbsvx.f
+++ b/SRC/cgbsvx.f
@@ -173,7 +173,7 @@
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is or output) COMPLEX array, dimension (LDAFB,N)
+*> AFB is COMPLEX array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains details of the LU factorization of the band matrix
*> A, as computed by CGBTRF. U is stored as an upper triangular
@@ -199,7 +199,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = L*U
*> as computed by CGBTRF; row i of the matrix was interchanged
@@ -216,7 +216,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -231,7 +231,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) REAL array, dimension (N)
+*> R is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -241,7 +241,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) REAL array, dimension (N)
+*> C is REAL array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/cgbsvxx.f b/SRC/cgbsvxx.f
index 3bd7db2d..18bfb913 100644
--- a/SRC/cgbsvxx.f
+++ b/SRC/cgbsvxx.f
@@ -200,7 +200,7 @@
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is or output) REAL array, dimension (LDAFB,N)
+*> AFB is REAL array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains details of the LU factorization of the band matrix
*> A, as computed by CGBTRF. U is stored as an upper triangular
@@ -227,7 +227,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by SGETRF; row i of the matrix was interchanged
@@ -244,7 +244,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -259,7 +259,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) REAL array, dimension (N)
+*> R is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -277,7 +277,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) REAL array, dimension (N)
+*> C is REAL array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/cgesvx.f b/SRC/cgesvx.f
index ca26fcca..4d41d6e8 100644
--- a/SRC/cgesvx.f
+++ b/SRC/cgesvx.f
@@ -155,7 +155,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX array, dimension (LDAF,N)
+*> AF is COMPLEX array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the factors L and U from the factorization
*> A = P*L*U as computed by CGETRF. If EQUED .ne. 'N', then
@@ -179,7 +179,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by CGETRF; row i of the matrix was interchanged
@@ -196,7 +196,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -211,7 +211,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) REAL array, dimension (N)
+*> R is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -221,7 +221,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) REAL array, dimension (N)
+*> C is REAL array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/cgesvxx.f b/SRC/cgesvxx.f
index 2051557a..3ecd635a 100644
--- a/SRC/cgesvxx.f
+++ b/SRC/cgesvxx.f
@@ -183,7 +183,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX array, dimension (LDAF,N)
+*> AF is COMPLEX array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the factors L and U from the factorization
*> A = P*L*U as computed by CGETRF. If EQUED .ne. 'N', then
@@ -207,7 +207,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by CGETRF; row i of the matrix was interchanged
@@ -224,7 +224,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -239,7 +239,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) REAL array, dimension (N)
+*> R is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -257,7 +257,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) REAL array, dimension (N)
+*> C is REAL array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/cgtsvx.f b/SRC/cgtsvx.f
index e16bc2b9..a839381b 100644
--- a/SRC/cgtsvx.f
+++ b/SRC/cgtsvx.f
@@ -134,7 +134,7 @@
*>
*> \param[in,out] DLF
*> \verbatim
-*> DLF is or output) COMPLEX array, dimension (N-1)
+*> DLF is COMPLEX array, dimension (N-1)
*> If FACT = 'F', then DLF is an input argument and on entry
*> contains the (n-1) multipliers that define the matrix L from
*> the LU factorization of A as computed by CGTTRF.
@@ -146,7 +146,7 @@
*>
*> \param[in,out] DF
*> \verbatim
-*> DF is or output) COMPLEX array, dimension (N)
+*> DF is COMPLEX array, dimension (N)
*> If FACT = 'F', then DF is an input argument and on entry
*> contains the n diagonal elements of the upper triangular
*> matrix U from the LU factorization of A.
@@ -158,7 +158,7 @@
*>
*> \param[in,out] DUF
*> \verbatim
-*> DUF is or output) COMPLEX array, dimension (N-1)
+*> DUF is COMPLEX array, dimension (N-1)
*> If FACT = 'F', then DUF is an input argument and on entry
*> contains the (n-1) elements of the first superdiagonal of U.
*>
@@ -168,7 +168,7 @@
*>
*> \param[in,out] DU2
*> \verbatim
-*> DU2 is or output) COMPLEX array, dimension (N-2)
+*> DU2 is COMPLEX array, dimension (N-2)
*> If FACT = 'F', then DU2 is an input argument and on entry
*> contains the (n-2) elements of the second superdiagonal of
*> U.
@@ -180,7 +180,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the LU factorization of A as
*> computed by CGTTRF.
diff --git a/SRC/cherfsx.f b/SRC/cherfsx.f
index 91132d35..a04f0e87 100644
--- a/SRC/cherfsx.f
+++ b/SRC/cherfsx.f
@@ -141,7 +141,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/chesvx.f b/SRC/chesvx.f
index 8b7a4bd4..9503aa10 100644
--- a/SRC/chesvx.f
+++ b/SRC/chesvx.f
@@ -133,7 +133,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX array, dimension (LDAF,N)
+*> AF is COMPLEX array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
*> to obtain the factor U or L from the factorization
@@ -153,7 +153,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block structure
*> of D, as determined by CHETRF.
diff --git a/SRC/chesvxx.f b/SRC/chesvxx.f
index afad83e2..82bf3265 100644
--- a/SRC/chesvxx.f
+++ b/SRC/chesvxx.f
@@ -181,7 +181,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX array, dimension (LDAF,N)
+*> AF is COMPLEX array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers
*> used to obtain the factor U or L from the factorization A =
@@ -201,7 +201,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block
*> structure of D, as determined by CHETRF. If IPIV(k) > 0,
@@ -220,7 +220,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Both row and column equilibration, i.e., A has been
@@ -231,7 +231,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/chpsvx.f b/SRC/chpsvx.f
index 8c553022..817856f8 100644
--- a/SRC/chpsvx.f
+++ b/SRC/chpsvx.f
@@ -124,7 +124,7 @@
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is or output) COMPLEX array, dimension (N*(N+1)/2)
+*> AFP is COMPLEX array, dimension (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
*> to obtain the factor U or L from the factorization
@@ -140,7 +140,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block structure
*> of D, as determined by CHPTRF.
diff --git a/SRC/clatbs.f b/SRC/clatbs.f
index e07f62f7..62bdb4df 100644
--- a/SRC/clatbs.f
+++ b/SRC/clatbs.f
@@ -135,7 +135,7 @@
*>
*> \param[in,out] CNORM
*> \verbatim
-*> CNORM is or output) REAL array, dimension (N)
+*> CNORM is REAL array, dimension (N)
*>
*> If NORMIN = 'Y', CNORM is an input argument and CNORM(j)
*> contains the norm of the off-diagonal part of the j-th column
diff --git a/SRC/clatps.f b/SRC/clatps.f
index 7c96a7cb..e3c8814f 100644
--- a/SRC/clatps.f
+++ b/SRC/clatps.f
@@ -123,7 +123,7 @@
*>
*> \param[in,out] CNORM
*> \verbatim
-*> CNORM is or output) REAL array, dimension (N)
+*> CNORM is REAL array, dimension (N)
*>
*> If NORMIN = 'Y', CNORM is an input argument and CNORM(j)
*> contains the norm of the off-diagonal part of the j-th column
diff --git a/SRC/clatrs.f b/SRC/clatrs.f
index 9c6048c3..5ce8cb1b 100644
--- a/SRC/clatrs.f
+++ b/SRC/clatrs.f
@@ -131,7 +131,7 @@
*>
*> \param[in,out] CNORM
*> \verbatim
-*> CNORM is or output) REAL array, dimension (N)
+*> CNORM is REAL array, dimension (N)
*>
*> If NORMIN = 'Y', CNORM is an input argument and CNORM(j)
*> contains the norm of the off-diagonal part of the j-th column
diff --git a/SRC/cpbsvx.f b/SRC/cpbsvx.f
index 409abcc5..8dd5be3d 100644
--- a/SRC/cpbsvx.f
+++ b/SRC/cpbsvx.f
@@ -160,7 +160,7 @@
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is or output) COMPLEX array, dimension (LDAFB,N)
+*> AFB is COMPLEX array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**H*U or A = L*L**H of the band matrix
@@ -186,7 +186,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Equilibration was done, i.e., A has been replaced by
@@ -197,7 +197,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The scale factors for A; not accessed if EQUED = 'N'. S is
*> an input argument if FACT = 'F'; otherwise, S is an output
*> argument. If FACT = 'F' and EQUED = 'Y', each element of S
diff --git a/SRC/cporfsx.f b/SRC/cporfsx.f
index 3889dadf..7dfb4a93 100644
--- a/SRC/cporfsx.f
+++ b/SRC/cporfsx.f
@@ -133,7 +133,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/cposvx.f b/SRC/cposvx.f
index 53bd377e..9e5b773a 100644
--- a/SRC/cposvx.f
+++ b/SRC/cposvx.f
@@ -155,7 +155,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX array, dimension (LDAF,N)
+*> AF is COMPLEX array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**H*U or A = L*L**H, in the same storage
@@ -182,7 +182,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Equilibration was done, i.e., A has been replaced by
@@ -193,7 +193,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The scale factors for A; not accessed if EQUED = 'N'. S is
*> an input argument if FACT = 'F'; otherwise, S is an output
*> argument. If FACT = 'F' and EQUED = 'Y', each element of S
diff --git a/SRC/cposvxx.f b/SRC/cposvxx.f
index 7aeee08d..def825c6 100644
--- a/SRC/cposvxx.f
+++ b/SRC/cposvxx.f
@@ -180,7 +180,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX array, dimension (LDAF,N)
+*> AF is COMPLEX array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**T*U or A = L*L**T, in the same storage
@@ -207,7 +207,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Both row and column equilibration, i.e., A has been
@@ -218,7 +218,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/cppsvx.f b/SRC/cppsvx.f
index e722e428..6a1c38c1 100644
--- a/SRC/cppsvx.f
+++ b/SRC/cppsvx.f
@@ -147,7 +147,7 @@
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is or output) COMPLEX array, dimension (N*(N+1)/2)
+*> AFP is COMPLEX array, dimension (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**H*U or A = L*L**H, in the same storage
@@ -168,7 +168,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Equilibration was done, i.e., A has been replaced by
@@ -179,7 +179,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The scale factors for A; not accessed if EQUED = 'N'. S is
*> an input argument if FACT = 'F'; otherwise, S is an output
*> argument. If FACT = 'F' and EQUED = 'Y', each element of S
diff --git a/SRC/cptsvx.f b/SRC/cptsvx.f
index 85c0e33a..9bb12c08 100644
--- a/SRC/cptsvx.f
+++ b/SRC/cptsvx.f
@@ -117,7 +117,7 @@
*>
*> \param[in,out] DF
*> \verbatim
-*> DF is or output) REAL array, dimension (N)
+*> DF is REAL array, dimension (N)
*> If FACT = 'F', then DF is an input argument and on entry
*> contains the n diagonal elements of the diagonal matrix D
*> from the L*D*L**H factorization of A.
@@ -128,7 +128,7 @@
*>
*> \param[in,out] EF
*> \verbatim
-*> EF is or output) COMPLEX array, dimension (N-1)
+*> EF is COMPLEX array, dimension (N-1)
*> If FACT = 'F', then EF is an input argument and on entry
*> contains the (n-1) subdiagonal elements of the unit
*> bidiagonal factor L from the L*D*L**H factorization of A.
diff --git a/SRC/cspsvx.f b/SRC/cspsvx.f
index df4bd39b..49729171 100644
--- a/SRC/cspsvx.f
+++ b/SRC/cspsvx.f
@@ -124,7 +124,7 @@
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is or output) COMPLEX array, dimension (N*(N+1)/2)
+*> AFP is COMPLEX array, dimension (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
*> to obtain the factor U or L from the factorization
@@ -140,7 +140,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block structure
*> of D, as determined by CSPTRF.
diff --git a/SRC/csyrfsx.f b/SRC/csyrfsx.f
index a7fb995f..0ff61b12 100644
--- a/SRC/csyrfsx.f
+++ b/SRC/csyrfsx.f
@@ -142,7 +142,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/csysvx.f b/SRC/csysvx.f
index 21397b78..05362172 100644
--- a/SRC/csysvx.f
+++ b/SRC/csysvx.f
@@ -133,7 +133,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX array, dimension (LDAF,N)
+*> AF is COMPLEX array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
*> to obtain the factor U or L from the factorization
@@ -153,7 +153,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block structure
*> of D, as determined by CSYTRF.
diff --git a/SRC/csysvxx.f b/SRC/csysvxx.f
index 79de2a79..c5654ff5 100644
--- a/SRC/csysvxx.f
+++ b/SRC/csysvxx.f
@@ -181,7 +181,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX array, dimension (LDAF,N)
+*> AF is COMPLEX array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers
*> used to obtain the factor U or L from the factorization A =
@@ -201,7 +201,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block
*> structure of D, as determined by SSYTRF. If IPIV(k) > 0,
@@ -220,7 +220,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Both row and column equilibration, i.e., A has been
@@ -231,7 +231,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/dgbrfsx.f b/SRC/dgbrfsx.f
index bdad0119..cf0bf6e8 100644
--- a/SRC/dgbrfsx.f
+++ b/SRC/dgbrfsx.f
@@ -159,7 +159,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) DOUBLE PRECISION array, dimension (N)
+*> R is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -177,7 +177,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) DOUBLE PRECISION array, dimension (N)
+*> C is DOUBLE PRECISION array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/dgbsvx.f b/SRC/dgbsvx.f
index 32c34795..b31fe12b 100644
--- a/SRC/dgbsvx.f
+++ b/SRC/dgbsvx.f
@@ -172,7 +172,7 @@
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is or output) DOUBLE PRECISION array, dimension (LDAFB,N)
+*> AFB is DOUBLE PRECISION array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains details of the LU factorization of the band matrix
*> A, as computed by DGBTRF. U is stored as an upper triangular
@@ -198,7 +198,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = L*U
*> as computed by DGBTRF; row i of the matrix was interchanged
@@ -215,7 +215,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -230,7 +230,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) DOUBLE PRECISION array, dimension (N)
+*> R is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -240,7 +240,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) DOUBLE PRECISION array, dimension (N)
+*> C is DOUBLE PRECISION array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/dgbsvxx.f b/SRC/dgbsvxx.f
index 9990638d..f4259c76 100644
--- a/SRC/dgbsvxx.f
+++ b/SRC/dgbsvxx.f
@@ -200,7 +200,7 @@
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is or output) DOUBLE PRECISION array, dimension (LDAFB,N)
+*> AFB is DOUBLE PRECISION array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains details of the LU factorization of the band matrix
*> A, as computed by DGBTRF. U is stored as an upper triangular
@@ -227,7 +227,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by DGETRF; row i of the matrix was interchanged
@@ -244,7 +244,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -259,7 +259,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) DOUBLE PRECISION array, dimension (N)
+*> R is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -277,7 +277,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) DOUBLE PRECISION array, dimension (N)
+*> C is DOUBLE PRECISION array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/dgesvx.f b/SRC/dgesvx.f
index c9e6c704..1ebd4b13 100644
--- a/SRC/dgesvx.f
+++ b/SRC/dgesvx.f
@@ -154,7 +154,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) DOUBLE PRECISION array, dimension (LDAF,N)
+*> AF is DOUBLE PRECISION array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the factors L and U from the factorization
*> A = P*L*U as computed by DGETRF. If EQUED .ne. 'N', then
@@ -178,7 +178,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by DGETRF; row i of the matrix was interchanged
@@ -195,7 +195,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -210,7 +210,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) DOUBLE PRECISION array, dimension (N)
+*> R is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -220,7 +220,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) DOUBLE PRECISION array, dimension (N)
+*> C is DOUBLE PRECISION array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/dgesvxx.f b/SRC/dgesvxx.f
index 33b83b8f..265e1e0f 100644
--- a/SRC/dgesvxx.f
+++ b/SRC/dgesvxx.f
@@ -183,7 +183,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) DOUBLE PRECISION array, dimension (LDAF,N)
+*> AF is DOUBLE PRECISION array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the factors L and U from the factorization
*> A = P*L*U as computed by DGETRF. If EQUED .ne. 'N', then
@@ -207,7 +207,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by DGETRF; row i of the matrix was interchanged
@@ -224,7 +224,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -239,7 +239,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) DOUBLE PRECISION array, dimension (N)
+*> R is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -257,7 +257,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) DOUBLE PRECISION array, dimension (N)
+*> C is DOUBLE PRECISION array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/dgtsvx.f b/SRC/dgtsvx.f
index b0ce0fbe..5ab690c5 100644
--- a/SRC/dgtsvx.f
+++ b/SRC/dgtsvx.f
@@ -133,7 +133,7 @@
*>
*> \param[in,out] DLF
*> \verbatim
-*> DLF is or output) DOUBLE PRECISION array, dimension (N-1)
+*> DLF is DOUBLE PRECISION array, dimension (N-1)
*> If FACT = 'F', then DLF is an input argument and on entry
*> contains the (n-1) multipliers that define the matrix L from
*> the LU factorization of A as computed by DGTTRF.
@@ -145,7 +145,7 @@
*>
*> \param[in,out] DF
*> \verbatim
-*> DF is or output) DOUBLE PRECISION array, dimension (N)
+*> DF is DOUBLE PRECISION array, dimension (N)
*> If FACT = 'F', then DF is an input argument and on entry
*> contains the n diagonal elements of the upper triangular
*> matrix U from the LU factorization of A.
@@ -157,7 +157,7 @@
*>
*> \param[in,out] DUF
*> \verbatim
-*> DUF is or output) DOUBLE PRECISION array, dimension (N-1)
+*> DUF is DOUBLE PRECISION array, dimension (N-1)
*> If FACT = 'F', then DUF is an input argument and on entry
*> contains the (n-1) elements of the first superdiagonal of U.
*>
@@ -167,7 +167,7 @@
*>
*> \param[in,out] DU2
*> \verbatim
-*> DU2 is or output) DOUBLE PRECISION array, dimension (N-2)
+*> DU2 is DOUBLE PRECISION array, dimension (N-2)
*> If FACT = 'F', then DU2 is an input argument and on entry
*> contains the (n-2) elements of the second superdiagonal of
*> U.
@@ -179,7 +179,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the LU factorization of A as
*> computed by DGTTRF.
diff --git a/SRC/dlatbs.f b/SRC/dlatbs.f
index f594033e..9579e8ba 100644
--- a/SRC/dlatbs.f
+++ b/SRC/dlatbs.f
@@ -134,7 +134,7 @@
*>
*> \param[in,out] CNORM
*> \verbatim
-*> CNORM is or output) DOUBLE PRECISION array, dimension (N)
+*> CNORM is DOUBLE PRECISION array, dimension (N)
*>
*> If NORMIN = 'Y', CNORM is an input argument and CNORM(j)
*> contains the norm of the off-diagonal part of the j-th column
diff --git a/SRC/dlatps.f b/SRC/dlatps.f
index ef3df08e..d3c77d13 100644
--- a/SRC/dlatps.f
+++ b/SRC/dlatps.f
@@ -121,7 +121,7 @@
*>
*> \param[in,out] CNORM
*> \verbatim
-*> CNORM is or output) DOUBLE PRECISION array, dimension (N)
+*> CNORM is DOUBLE PRECISION array, dimension (N)
*>
*> If NORMIN = 'Y', CNORM is an input argument and CNORM(j)
*> contains the norm of the off-diagonal part of the j-th column
diff --git a/SRC/dlatrs.f b/SRC/dlatrs.f
index f5f01b46..fddff337 100644
--- a/SRC/dlatrs.f
+++ b/SRC/dlatrs.f
@@ -130,7 +130,7 @@
*>
*> \param[in,out] CNORM
*> \verbatim
-*> CNORM is or output) DOUBLE PRECISION array, dimension (N)
+*> CNORM is DOUBLE PRECISION array, dimension (N)
*>
*> If NORMIN = 'Y', CNORM is an input argument and CNORM(j)
*> contains the norm of the off-diagonal part of the j-th column
diff --git a/SRC/dpbsvx.f b/SRC/dpbsvx.f
index dc57bcd3..f5ab009a 100644
--- a/SRC/dpbsvx.f
+++ b/SRC/dpbsvx.f
@@ -161,7 +161,7 @@
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is or output) DOUBLE PRECISION array, dimension (LDAFB,N)
+*> AFB is DOUBLE PRECISION array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**T*U or A = L*L**T of the band matrix
@@ -187,7 +187,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Equilibration was done, i.e., A has been replaced by
@@ -198,7 +198,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The scale factors for A; not accessed if EQUED = 'N'. S is
*> an input argument if FACT = 'F'; otherwise, S is an output
*> argument. If FACT = 'F' and EQUED = 'Y', each element of S
diff --git a/SRC/dporfsx.f b/SRC/dporfsx.f
index 27555071..bc35fc0e 100644
--- a/SRC/dporfsx.f
+++ b/SRC/dporfsx.f
@@ -134,7 +134,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/dposvx.f b/SRC/dposvx.f
index e09c623c..8bdadd5a 100644
--- a/SRC/dposvx.f
+++ b/SRC/dposvx.f
@@ -156,7 +156,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) DOUBLE PRECISION array, dimension (LDAF,N)
+*> AF is DOUBLE PRECISION array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**T*U or A = L*L**T, in the same storage
@@ -183,7 +183,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Equilibration was done, i.e., A has been replaced by
@@ -194,7 +194,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The scale factors for A; not accessed if EQUED = 'N'. S is
*> an input argument if FACT = 'F'; otherwise, S is an output
*> argument. If FACT = 'F' and EQUED = 'Y', each element of S
diff --git a/SRC/dposvxx.f b/SRC/dposvxx.f
index d6114a23..1544ec98 100644
--- a/SRC/dposvxx.f
+++ b/SRC/dposvxx.f
@@ -181,7 +181,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) DOUBLE PRECISION array, dimension (LDAF,N)
+*> AF is DOUBLE PRECISION array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**T*U or A = L*L**T, in the same storage
@@ -208,7 +208,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Both row and column equilibration, i.e., A has been
@@ -219,7 +219,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/dppsvx.f b/SRC/dppsvx.f
index ae82afd3..9ec06b31 100644
--- a/SRC/dppsvx.f
+++ b/SRC/dppsvx.f
@@ -147,7 +147,7 @@
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is or output) DOUBLE PRECISION array, dimension
+*> AFP is DOUBLE PRECISION array, dimension
*> (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
@@ -169,7 +169,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Equilibration was done, i.e., A has been replaced by
@@ -180,7 +180,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The scale factors for A; not accessed if EQUED = 'N'. S is
*> an input argument if FACT = 'F'; otherwise, S is an output
*> argument. If FACT = 'F' and EQUED = 'Y', each element of S
diff --git a/SRC/dptsvx.f b/SRC/dptsvx.f
index 3d51d35e..9643607b 100644
--- a/SRC/dptsvx.f
+++ b/SRC/dptsvx.f
@@ -116,7 +116,7 @@
*>
*> \param[in,out] DF
*> \verbatim
-*> DF is or output) DOUBLE PRECISION array, dimension (N)
+*> DF is DOUBLE PRECISION array, dimension (N)
*> If FACT = 'F', then DF is an input argument and on entry
*> contains the n diagonal elements of the diagonal matrix D
*> from the L*D*L**T factorization of A.
@@ -127,7 +127,7 @@
*>
*> \param[in,out] EF
*> \verbatim
-*> EF is or output) DOUBLE PRECISION array, dimension (N-1)
+*> EF is DOUBLE PRECISION array, dimension (N-1)
*> If FACT = 'F', then EF is an input argument and on entry
*> contains the (n-1) subdiagonal elements of the unit
*> bidiagonal factor L from the L*D*L**T factorization of A.
diff --git a/SRC/dspsvx.f b/SRC/dspsvx.f
index 402b6817..e25d5dff 100644
--- a/SRC/dspsvx.f
+++ b/SRC/dspsvx.f
@@ -123,7 +123,7 @@
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is or output) DOUBLE PRECISION array, dimension
+*> AFP is DOUBLE PRECISION array, dimension
*> (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
@@ -140,7 +140,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block structure
*> of D, as determined by DSPTRF.
diff --git a/SRC/dsyrfsx.f b/SRC/dsyrfsx.f
index 6959d039..fde12474 100644
--- a/SRC/dsyrfsx.f
+++ b/SRC/dsyrfsx.f
@@ -142,7 +142,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/dsysvx.f b/SRC/dsysvx.f
index 0855f65e..2f57cf2f 100644
--- a/SRC/dsysvx.f
+++ b/SRC/dsysvx.f
@@ -132,7 +132,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) DOUBLE PRECISION array, dimension (LDAF,N)
+*> AF is DOUBLE PRECISION array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
*> to obtain the factor U or L from the factorization
@@ -152,7 +152,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block structure
*> of D, as determined by DSYTRF.
diff --git a/SRC/dsysvxx.f b/SRC/dsysvxx.f
index 10aee2d9..f8cf7ea2 100644
--- a/SRC/dsysvxx.f
+++ b/SRC/dsysvxx.f
@@ -180,7 +180,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) DOUBLE PRECISION array, dimension (LDAF,N)
+*> AF is DOUBLE PRECISION array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers
*> used to obtain the factor U or L from the factorization A =
@@ -200,7 +200,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block
*> structure of D, as determined by DSYTRF. If IPIV(k) > 0,
@@ -219,7 +219,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Both row and column equilibration, i.e., A has been
@@ -230,7 +230,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/sgbrfsx.f b/SRC/sgbrfsx.f
index 6df4f271..08575899 100644
--- a/SRC/sgbrfsx.f
+++ b/SRC/sgbrfsx.f
@@ -159,7 +159,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) REAL array, dimension (N)
+*> R is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -177,7 +177,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) REAL array, dimension (N)
+*> C is REAL array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/sgbsvx.f b/SRC/sgbsvx.f
index 92c95d57..97a6be86 100644
--- a/SRC/sgbsvx.f
+++ b/SRC/sgbsvx.f
@@ -172,7 +172,7 @@
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is or output) REAL array, dimension (LDAFB,N)
+*> AFB is REAL array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains details of the LU factorization of the band matrix
*> A, as computed by SGBTRF. U is stored as an upper triangular
@@ -198,7 +198,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = L*U
*> as computed by SGBTRF; row i of the matrix was interchanged
@@ -215,7 +215,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -230,7 +230,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) REAL array, dimension (N)
+*> R is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -240,7 +240,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) REAL array, dimension (N)
+*> C is REAL array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/sgbsvxx.f b/SRC/sgbsvxx.f
index 7f880394..82466320 100644
--- a/SRC/sgbsvxx.f
+++ b/SRC/sgbsvxx.f
@@ -200,7 +200,7 @@
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is or output) REAL array, dimension (LDAFB,N)
+*> AFB is REAL array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains details of the LU factorization of the band matrix
*> A, as computed by SGBTRF. U is stored as an upper triangular
@@ -227,7 +227,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by SGETRF; row i of the matrix was interchanged
@@ -244,7 +244,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -259,7 +259,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) REAL array, dimension (N)
+*> R is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -277,7 +277,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) REAL array, dimension (N)
+*> C is REAL array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/sgesvx.f b/SRC/sgesvx.f
index 3e2f1b12..6def2508 100644
--- a/SRC/sgesvx.f
+++ b/SRC/sgesvx.f
@@ -154,7 +154,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) REAL array, dimension (LDAF,N)
+*> AF is REAL array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the factors L and U from the factorization
*> A = P*L*U as computed by SGETRF. If EQUED .ne. 'N', then
@@ -178,7 +178,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by SGETRF; row i of the matrix was interchanged
@@ -195,7 +195,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -210,7 +210,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) REAL array, dimension (N)
+*> R is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -220,7 +220,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) REAL array, dimension (N)
+*> C is REAL array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/sgesvxx.f b/SRC/sgesvxx.f
index bbff217e..8ed2e599 100644
--- a/SRC/sgesvxx.f
+++ b/SRC/sgesvxx.f
@@ -183,7 +183,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) REAL array, dimension (LDAF,N)
+*> AF is REAL array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the factors L and U from the factorization
*> A = P*L*U as computed by SGETRF. If EQUED .ne. 'N', then
@@ -207,7 +207,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by SGETRF; row i of the matrix was interchanged
@@ -224,7 +224,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -239,7 +239,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) REAL array, dimension (N)
+*> R is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -257,7 +257,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) REAL array, dimension (N)
+*> C is REAL array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/sgtsvx.f b/SRC/sgtsvx.f
index 535f17a8..20e8c9c4 100644
--- a/SRC/sgtsvx.f
+++ b/SRC/sgtsvx.f
@@ -133,7 +133,7 @@
*>
*> \param[in,out] DLF
*> \verbatim
-*> DLF is or output) REAL array, dimension (N-1)
+*> DLF is REAL array, dimension (N-1)
*> If FACT = 'F', then DLF is an input argument and on entry
*> contains the (n-1) multipliers that define the matrix L from
*> the LU factorization of A as computed by SGTTRF.
@@ -145,7 +145,7 @@
*>
*> \param[in,out] DF
*> \verbatim
-*> DF is or output) REAL array, dimension (N)
+*> DF is REAL array, dimension (N)
*> If FACT = 'F', then DF is an input argument and on entry
*> contains the n diagonal elements of the upper triangular
*> matrix U from the LU factorization of A.
@@ -157,7 +157,7 @@
*>
*> \param[in,out] DUF
*> \verbatim
-*> DUF is or output) REAL array, dimension (N-1)
+*> DUF is REAL array, dimension (N-1)
*> If FACT = 'F', then DUF is an input argument and on entry
*> contains the (n-1) elements of the first superdiagonal of U.
*>
@@ -167,7 +167,7 @@
*>
*> \param[in,out] DU2
*> \verbatim
-*> DU2 is or output) REAL array, dimension (N-2)
+*> DU2 is REAL array, dimension (N-2)
*> If FACT = 'F', then DU2 is an input argument and on entry
*> contains the (n-2) elements of the second superdiagonal of
*> U.
@@ -179,7 +179,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the LU factorization of A as
*> computed by SGTTRF.
diff --git a/SRC/slatbs.f b/SRC/slatbs.f
index 44a52de8..89c9f88b 100644
--- a/SRC/slatbs.f
+++ b/SRC/slatbs.f
@@ -134,7 +134,7 @@
*>
*> \param[in,out] CNORM
*> \verbatim
-*> CNORM is or output) REAL array, dimension (N)
+*> CNORM is REAL array, dimension (N)
*>
*> If NORMIN = 'Y', CNORM is an input argument and CNORM(j)
*> contains the norm of the off-diagonal part of the j-th column
diff --git a/SRC/slatps.f b/SRC/slatps.f
index 463e5b2f..95061953 100644
--- a/SRC/slatps.f
+++ b/SRC/slatps.f
@@ -121,7 +121,7 @@
*>
*> \param[in,out] CNORM
*> \verbatim
-*> CNORM is or output) REAL array, dimension (N)
+*> CNORM is REAL array, dimension (N)
*>
*> If NORMIN = 'Y', CNORM is an input argument and CNORM(j)
*> contains the norm of the off-diagonal part of the j-th column
diff --git a/SRC/slatrs.f b/SRC/slatrs.f
index a753472b..df1b4f11 100644
--- a/SRC/slatrs.f
+++ b/SRC/slatrs.f
@@ -130,7 +130,7 @@
*>
*> \param[in,out] CNORM
*> \verbatim
-*> CNORM is or output) REAL array, dimension (N)
+*> CNORM is REAL array, dimension (N)
*>
*> If NORMIN = 'Y', CNORM is an input argument and CNORM(j)
*> contains the norm of the off-diagonal part of the j-th column
diff --git a/SRC/spbsvx.f b/SRC/spbsvx.f
index eb0b9f6d..84f27443 100644
--- a/SRC/spbsvx.f
+++ b/SRC/spbsvx.f
@@ -161,7 +161,7 @@
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is or output) REAL array, dimension (LDAFB,N)
+*> AFB is REAL array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**T*U or A = L*L**T of the band matrix
@@ -187,7 +187,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Equilibration was done, i.e., A has been replaced by
@@ -198,7 +198,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The scale factors for A; not accessed if EQUED = 'N'. S is
*> an input argument if FACT = 'F'; otherwise, S is an output
*> argument. If FACT = 'F' and EQUED = 'Y', each element of S
diff --git a/SRC/sporfsx.f b/SRC/sporfsx.f
index 2ef392ff..e2e87551 100644
--- a/SRC/sporfsx.f
+++ b/SRC/sporfsx.f
@@ -134,7 +134,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/sposvx.f b/SRC/sposvx.f
index aff52230..948249e0 100644
--- a/SRC/sposvx.f
+++ b/SRC/sposvx.f
@@ -156,7 +156,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) REAL array, dimension (LDAF,N)
+*> AF is REAL array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**T*U or A = L*L**T, in the same storage
@@ -183,7 +183,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Equilibration was done, i.e., A has been replaced by
@@ -194,7 +194,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The scale factors for A; not accessed if EQUED = 'N'. S is
*> an input argument if FACT = 'F'; otherwise, S is an output
*> argument. If FACT = 'F' and EQUED = 'Y', each element of S
diff --git a/SRC/sposvxx.f b/SRC/sposvxx.f
index 906c0910..9372076e 100644
--- a/SRC/sposvxx.f
+++ b/SRC/sposvxx.f
@@ -181,7 +181,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) REAL array, dimension (LDAF,N)
+*> AF is REAL array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**T*U or A = L*L**T, in the same storage
@@ -208,7 +208,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Both row and column equilibration, i.e., A has been
@@ -219,7 +219,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The row scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/sppsvx.f b/SRC/sppsvx.f
index 10b89d11..11b5586a 100644
--- a/SRC/sppsvx.f
+++ b/SRC/sppsvx.f
@@ -147,7 +147,7 @@
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is or output) REAL array, dimension
+*> AFP is REAL array, dimension
*> (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
@@ -169,7 +169,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Equilibration was done, i.e., A has been replaced by
@@ -180,7 +180,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The scale factors for A; not accessed if EQUED = 'N'. S is
*> an input argument if FACT = 'F'; otherwise, S is an output
*> argument. If FACT = 'F' and EQUED = 'Y', each element of S
diff --git a/SRC/sptsvx.f b/SRC/sptsvx.f
index 6b8b566d..77bdf380 100644
--- a/SRC/sptsvx.f
+++ b/SRC/sptsvx.f
@@ -116,7 +116,7 @@
*>
*> \param[in,out] DF
*> \verbatim
-*> DF is or output) REAL array, dimension (N)
+*> DF is REAL array, dimension (N)
*> If FACT = 'F', then DF is an input argument and on entry
*> contains the n diagonal elements of the diagonal matrix D
*> from the L*D*L**T factorization of A.
@@ -127,7 +127,7 @@
*>
*> \param[in,out] EF
*> \verbatim
-*> EF is or output) REAL array, dimension (N-1)
+*> EF is REAL array, dimension (N-1)
*> If FACT = 'F', then EF is an input argument and on entry
*> contains the (n-1) subdiagonal elements of the unit
*> bidiagonal factor L from the L*D*L**T factorization of A.
diff --git a/SRC/sspsvx.f b/SRC/sspsvx.f
index cbfeb410..b243764d 100644
--- a/SRC/sspsvx.f
+++ b/SRC/sspsvx.f
@@ -123,7 +123,7 @@
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is or output) REAL array, dimension
+*> AFP is REAL array, dimension
*> (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
@@ -140,7 +140,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block structure
*> of D, as determined by SSPTRF.
diff --git a/SRC/ssyrfsx.f b/SRC/ssyrfsx.f
index 0e152926..611ce3a1 100644
--- a/SRC/ssyrfsx.f
+++ b/SRC/ssyrfsx.f
@@ -142,7 +142,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/ssysvx.f b/SRC/ssysvx.f
index c2cf1dc8..5b76f184 100644
--- a/SRC/ssysvx.f
+++ b/SRC/ssysvx.f
@@ -132,7 +132,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) REAL array, dimension (LDAF,N)
+*> AF is REAL array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
*> to obtain the factor U or L from the factorization
@@ -152,7 +152,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block structure
*> of D, as determined by SSYTRF.
diff --git a/SRC/ssysvxx.f b/SRC/ssysvxx.f
index 7c6261d6..89f28434 100644
--- a/SRC/ssysvxx.f
+++ b/SRC/ssysvxx.f
@@ -180,7 +180,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) REAL array, dimension (LDAF,N)
+*> AF is REAL array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers
*> used to obtain the factor U or L from the factorization A =
@@ -200,7 +200,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block
*> structure of D, as determined by SSYTRF. If IPIV(k) > 0,
@@ -219,7 +219,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Both row and column equilibration, i.e., A has been
@@ -230,7 +230,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) REAL array, dimension (N)
+*> S is REAL array, dimension (N)
*> The scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/zgbrfsx.f b/SRC/zgbrfsx.f
index 9943d824..d34950c1 100644
--- a/SRC/zgbrfsx.f
+++ b/SRC/zgbrfsx.f
@@ -159,7 +159,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) DOUBLE PRECISION array, dimension (N)
+*> R is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -177,7 +177,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) DOUBLE PRECISION array, dimension (N)
+*> C is DOUBLE PRECISION array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/zgbsvx.f b/SRC/zgbsvx.f
index cd4e381b..266be625 100644
--- a/SRC/zgbsvx.f
+++ b/SRC/zgbsvx.f
@@ -173,7 +173,7 @@
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is or output) COMPLEX*16 array, dimension (LDAFB,N)
+*> AFB is COMPLEX*16 array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains details of the LU factorization of the band matrix
*> A, as computed by ZGBTRF. U is stored as an upper triangular
@@ -199,7 +199,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = L*U
*> as computed by ZGBTRF; row i of the matrix was interchanged
@@ -216,7 +216,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -231,7 +231,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) DOUBLE PRECISION array, dimension (N)
+*> R is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -241,7 +241,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) DOUBLE PRECISION array, dimension (N)
+*> C is DOUBLE PRECISION array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/zgbsvxx.f b/SRC/zgbsvxx.f
index 22b0a4d6..70639948 100644
--- a/SRC/zgbsvxx.f
+++ b/SRC/zgbsvxx.f
@@ -200,7 +200,7 @@
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is or output) COMPLEX*16 array, dimension (LDAFB,N)
+*> AFB is COMPLEX*16 array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains details of the LU factorization of the band matrix
*> A, as computed by ZGBTRF. U is stored as an upper triangular
@@ -227,7 +227,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by DGETRF; row i of the matrix was interchanged
@@ -244,7 +244,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -259,7 +259,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) DOUBLE PRECISION array, dimension (N)
+*> R is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -277,7 +277,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) DOUBLE PRECISION array, dimension (N)
+*> C is DOUBLE PRECISION array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/zgesvx.f b/SRC/zgesvx.f
index 77bcc0e0..4fac956c 100644
--- a/SRC/zgesvx.f
+++ b/SRC/zgesvx.f
@@ -155,7 +155,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX*16 array, dimension (LDAF,N)
+*> AF is COMPLEX*16 array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the factors L and U from the factorization
*> A = P*L*U as computed by ZGETRF. If EQUED .ne. 'N', then
@@ -179,7 +179,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by ZGETRF; row i of the matrix was interchanged
@@ -196,7 +196,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -211,7 +211,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) DOUBLE PRECISION array, dimension (N)
+*> R is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -221,7 +221,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) DOUBLE PRECISION array, dimension (N)
+*> C is DOUBLE PRECISION array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/zgesvxx.f b/SRC/zgesvxx.f
index 01cd900e..06a0a1a4 100644
--- a/SRC/zgesvxx.f
+++ b/SRC/zgesvxx.f
@@ -183,7 +183,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX*16 array, dimension (LDAF,N)
+*> AF is COMPLEX*16 array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the factors L and U from the factorization
*> A = P*L*U as computed by ZGETRF. If EQUED .ne. 'N', then
@@ -207,7 +207,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by ZGETRF; row i of the matrix was interchanged
@@ -224,7 +224,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -239,7 +239,7 @@
*>
*> \param[in,out] R
*> \verbatim
-*> R is or output) DOUBLE PRECISION array, dimension (N)
+*> R is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'R' or 'B', A is
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -257,7 +257,7 @@
*>
*> \param[in,out] C
*> \verbatim
-*> C is or output) DOUBLE PRECISION array, dimension (N)
+*> C is DOUBLE PRECISION array, dimension (N)
*> The column scale factors for A. If EQUED = 'C' or 'B', A is
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
diff --git a/SRC/zgtsvx.f b/SRC/zgtsvx.f
index dc74f169..1aba5d0f 100644
--- a/SRC/zgtsvx.f
+++ b/SRC/zgtsvx.f
@@ -134,7 +134,7 @@
*>
*> \param[in,out] DLF
*> \verbatim
-*> DLF is or output) COMPLEX*16 array, dimension (N-1)
+*> DLF is COMPLEX*16 array, dimension (N-1)
*> If FACT = 'F', then DLF is an input argument and on entry
*> contains the (n-1) multipliers that define the matrix L from
*> the LU factorization of A as computed by ZGTTRF.
@@ -146,7 +146,7 @@
*>
*> \param[in,out] DF
*> \verbatim
-*> DF is or output) COMPLEX*16 array, dimension (N)
+*> DF is COMPLEX*16 array, dimension (N)
*> If FACT = 'F', then DF is an input argument and on entry
*> contains the n diagonal elements of the upper triangular
*> matrix U from the LU factorization of A.
@@ -158,7 +158,7 @@
*>
*> \param[in,out] DUF
*> \verbatim
-*> DUF is or output) COMPLEX*16 array, dimension (N-1)
+*> DUF is COMPLEX*16 array, dimension (N-1)
*> If FACT = 'F', then DUF is an input argument and on entry
*> contains the (n-1) elements of the first superdiagonal of U.
*>
@@ -168,7 +168,7 @@
*>
*> \param[in,out] DU2
*> \verbatim
-*> DU2 is or output) COMPLEX*16 array, dimension (N-2)
+*> DU2 is COMPLEX*16 array, dimension (N-2)
*> If FACT = 'F', then DU2 is an input argument and on entry
*> contains the (n-2) elements of the second superdiagonal of
*> U.
@@ -180,7 +180,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains the pivot indices from the LU factorization of A as
*> computed by ZGTTRF.
diff --git a/SRC/zherfsx.f b/SRC/zherfsx.f
index 62fdf989..6ecf5f21 100644
--- a/SRC/zherfsx.f
+++ b/SRC/zherfsx.f
@@ -141,7 +141,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/zhesvx.f b/SRC/zhesvx.f
index ce349187..cd2ecd98 100644
--- a/SRC/zhesvx.f
+++ b/SRC/zhesvx.f
@@ -133,7 +133,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX*16 array, dimension (LDAF,N)
+*> AF is COMPLEX*16 array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
*> to obtain the factor U or L from the factorization
@@ -153,7 +153,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block structure
*> of D, as determined by ZHETRF.
diff --git a/SRC/zhesvxx.f b/SRC/zhesvxx.f
index c02f925a..9e0b693d 100644
--- a/SRC/zhesvxx.f
+++ b/SRC/zhesvxx.f
@@ -181,7 +181,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX*16 array, dimension (LDAF,N)
+*> AF is COMPLEX*16 array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers
*> used to obtain the factor U or L from the factorization A =
@@ -201,7 +201,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block
*> structure of D, as determined by ZHETRF. If IPIV(k) > 0,
@@ -220,7 +220,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Both row and column equilibration, i.e., A has been
@@ -231,7 +231,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/zhpsvx.f b/SRC/zhpsvx.f
index 7f827686..d9ffb17f 100644
--- a/SRC/zhpsvx.f
+++ b/SRC/zhpsvx.f
@@ -124,7 +124,7 @@
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is or output) COMPLEX*16 array, dimension (N*(N+1)/2)
+*> AFP is COMPLEX*16 array, dimension (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
*> to obtain the factor U or L from the factorization
@@ -140,7 +140,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block structure
*> of D, as determined by ZHPTRF.
diff --git a/SRC/zlatbs.f b/SRC/zlatbs.f
index 8705d280..79448be8 100644
--- a/SRC/zlatbs.f
+++ b/SRC/zlatbs.f
@@ -135,7 +135,7 @@
*>
*> \param[in,out] CNORM
*> \verbatim
-*> CNORM is or output) DOUBLE PRECISION array, dimension (N)
+*> CNORM is DOUBLE PRECISION array, dimension (N)
*>
*> If NORMIN = 'Y', CNORM is an input argument and CNORM(j)
*> contains the norm of the off-diagonal part of the j-th column
diff --git a/SRC/zlatps.f b/SRC/zlatps.f
index 7e571c75..438b1283 100644
--- a/SRC/zlatps.f
+++ b/SRC/zlatps.f
@@ -123,7 +123,7 @@
*>
*> \param[in,out] CNORM
*> \verbatim
-*> CNORM is or output) DOUBLE PRECISION array, dimension (N)
+*> CNORM is DOUBLE PRECISION array, dimension (N)
*>
*> If NORMIN = 'Y', CNORM is an input argument and CNORM(j)
*> contains the norm of the off-diagonal part of the j-th column
diff --git a/SRC/zlatrs.f b/SRC/zlatrs.f
index 4fce05e5..77d694c8 100644
--- a/SRC/zlatrs.f
+++ b/SRC/zlatrs.f
@@ -131,7 +131,7 @@
*>
*> \param[in,out] CNORM
*> \verbatim
-*> CNORM is or output) DOUBLE PRECISION array, dimension (N)
+*> CNORM is DOUBLE PRECISION array, dimension (N)
*>
*> If NORMIN = 'Y', CNORM is an input argument and CNORM(j)
*> contains the norm of the off-diagonal part of the j-th column
diff --git a/SRC/zpbsvx.f b/SRC/zpbsvx.f
index 64a1aaba..0b920f6b 100644
--- a/SRC/zpbsvx.f
+++ b/SRC/zpbsvx.f
@@ -160,7 +160,7 @@
*>
*> \param[in,out] AFB
*> \verbatim
-*> AFB is or output) COMPLEX*16 array, dimension (LDAFB,N)
+*> AFB is COMPLEX*16 array, dimension (LDAFB,N)
*> If FACT = 'F', then AFB is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**H *U or A = L*L**H of the band matrix
@@ -186,7 +186,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Equilibration was done, i.e., A has been replaced by
@@ -197,7 +197,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The scale factors for A; not accessed if EQUED = 'N'. S is
*> an input argument if FACT = 'F'; otherwise, S is an output
*> argument. If FACT = 'F' and EQUED = 'Y', each element of S
diff --git a/SRC/zporfsx.f b/SRC/zporfsx.f
index 0687fa47..4fbb09f5 100644
--- a/SRC/zporfsx.f
+++ b/SRC/zporfsx.f
@@ -133,7 +133,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/zposvx.f b/SRC/zposvx.f
index 2dec9a09..fcf4a7df 100644
--- a/SRC/zposvx.f
+++ b/SRC/zposvx.f
@@ -155,7 +155,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX*16 array, dimension (LDAF,N)
+*> AF is COMPLEX*16 array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**H *U or A = L*L**H, in the same storage
@@ -182,7 +182,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Equilibration was done, i.e., A has been replaced by
@@ -193,7 +193,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The scale factors for A; not accessed if EQUED = 'N'. S is
*> an input argument if FACT = 'F'; otherwise, S is an output
*> argument. If FACT = 'F' and EQUED = 'Y', each element of S
diff --git a/SRC/zposvxx.f b/SRC/zposvxx.f
index 263805f9..9e0f4cb5 100644
--- a/SRC/zposvxx.f
+++ b/SRC/zposvxx.f
@@ -180,7 +180,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX*16 array, dimension (LDAF,N)
+*> AF is COMPLEX*16 array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**T*U or A = L*L**T, in the same storage
@@ -207,7 +207,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Both row and column equilibration, i.e., A has been
@@ -218,7 +218,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The row scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/zppsvx.f b/SRC/zppsvx.f
index 54ecb4c4..a6914ce2 100644
--- a/SRC/zppsvx.f
+++ b/SRC/zppsvx.f
@@ -147,7 +147,7 @@
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is or output) COMPLEX*16 array, dimension (N*(N+1)/2)
+*> AFP is COMPLEX*16 array, dimension (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the triangular factor U or L from the Cholesky
*> factorization A = U**H*U or A = L*L**H, in the same storage
@@ -168,7 +168,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Equilibration was done, i.e., A has been replaced by
@@ -179,7 +179,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The scale factors for A; not accessed if EQUED = 'N'. S is
*> an input argument if FACT = 'F'; otherwise, S is an output
*> argument. If FACT = 'F' and EQUED = 'Y', each element of S
diff --git a/SRC/zptsvx.f b/SRC/zptsvx.f
index 1dcca724..dc0fb5a9 100644
--- a/SRC/zptsvx.f
+++ b/SRC/zptsvx.f
@@ -117,7 +117,7 @@
*>
*> \param[in,out] DF
*> \verbatim
-*> DF is or output) DOUBLE PRECISION array, dimension (N)
+*> DF is DOUBLE PRECISION array, dimension (N)
*> If FACT = 'F', then DF is an input argument and on entry
*> contains the n diagonal elements of the diagonal matrix D
*> from the L*D*L**H factorization of A.
@@ -128,7 +128,7 @@
*>
*> \param[in,out] EF
*> \verbatim
-*> EF is or output) COMPLEX*16 array, dimension (N-1)
+*> EF is COMPLEX*16 array, dimension (N-1)
*> If FACT = 'F', then EF is an input argument and on entry
*> contains the (n-1) subdiagonal elements of the unit
*> bidiagonal factor L from the L*D*L**H factorization of A.
diff --git a/SRC/zspsvx.f b/SRC/zspsvx.f
index 6fb6a3ed..3679a3c8 100644
--- a/SRC/zspsvx.f
+++ b/SRC/zspsvx.f
@@ -124,7 +124,7 @@
*>
*> \param[in,out] AFP
*> \verbatim
-*> AFP is or output) COMPLEX*16 array, dimension (N*(N+1)/2)
+*> AFP is COMPLEX*16 array, dimension (N*(N+1)/2)
*> If FACT = 'F', then AFP is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
*> to obtain the factor U or L from the factorization
@@ -140,7 +140,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block structure
*> of D, as determined by ZSPTRF.
diff --git a/SRC/zsyrfsx.f b/SRC/zsyrfsx.f
index 994e2553..3ffa6774 100644
--- a/SRC/zsyrfsx.f
+++ b/SRC/zsyrfsx.f
@@ -142,7 +142,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED
diff --git a/SRC/zsysvx.f b/SRC/zsysvx.f
index 236e6d4d..3e6d6782 100644
--- a/SRC/zsysvx.f
+++ b/SRC/zsysvx.f
@@ -133,7 +133,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX*16 array, dimension (LDAF,N)
+*> AF is COMPLEX*16 array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers used
*> to obtain the factor U or L from the factorization
@@ -153,7 +153,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block structure
*> of D, as determined by ZSYTRF.
diff --git a/SRC/zsysvxx.f b/SRC/zsysvxx.f
index 81da5780..56b2fff8 100644
--- a/SRC/zsysvxx.f
+++ b/SRC/zsysvxx.f
@@ -181,7 +181,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) COMPLEX*16 array, dimension (LDAF,N)
+*> AF is COMPLEX*16 array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the block diagonal matrix D and the multipliers
*> used to obtain the factor U or L from the factorization A =
@@ -201,7 +201,7 @@
*>
*> \param[in,out] IPIV
*> \verbatim
-*> IPIV is or output) INTEGER array, dimension (N)
+*> IPIV is INTEGER array, dimension (N)
*> If FACT = 'F', then IPIV is an input argument and on entry
*> contains details of the interchanges and the block
*> structure of D, as determined by DSYTRF. If IPIV(k) > 0,
@@ -220,7 +220,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'Y': Both row and column equilibration, i.e., A has been
@@ -231,7 +231,7 @@
*>
*> \param[in,out] S
*> \verbatim
-*> S is or output) DOUBLE PRECISION array, dimension (N)
+*> S is DOUBLE PRECISION array, dimension (N)
*> The scale factors for A. If EQUED = 'Y', A is multiplied on
*> the left and right by diag(S). S is an input argument if FACT =
*> 'F'; otherwise, S is an output argument. If FACT = 'F' and EQUED