diff options
69 files changed, 346 insertions, 291 deletions
diff --git a/SRC/chfrk.f b/SRC/chfrk.f index caa71844..1a3d6567 100644 --- a/SRC/chfrk.f +++ b/SRC/chfrk.f @@ -39,11 +39,11 @@ * Arguments * ========== * -* TRANSR - (input) CHARACTER. +* TRANSR (input) CHARACTER * = 'N': The Normal Form of RFP A is stored; * = 'C': The Conjugate-transpose Form of RFP A is stored. * -* UPLO - (input) CHARACTER. +* UPLO (input) CHARACTER * On entry, UPLO specifies whether the upper or lower * triangular part of the array C is to be referenced as * follows: @@ -56,7 +56,7 @@ * * Unchanged on exit. * -* TRANS - (input) CHARACTER. +* TRANS (input) CHARACTER * On entry, TRANS specifies the operation to be performed as * follows: * @@ -66,41 +66,42 @@ * * Unchanged on exit. * -* N - (input) INTEGER. +* N (input) INTEGER * On entry, N specifies the order of the matrix C. N must be * at least zero. * Unchanged on exit. * -* K - (input) INTEGER. +* K (input) INTEGER * On entry with TRANS = 'N' or 'n', K specifies the number * of columns of the matrix A, and on entry with * TRANS = 'C' or 'c', K specifies the number of rows of the * matrix A. K must be at least zero. * Unchanged on exit. * -* ALPHA - (input) REAL. +* ALPHA (input) REAL * On entry, ALPHA specifies the scalar alpha. * Unchanged on exit. * -* A - (input) COMPLEX array of DIMENSION ( LDA, ka ), where KA +* A (input) COMPLEX array, dimension (LDA,ka) +* where KA * is K when TRANS = 'N' or 'n', and is N otherwise. Before * entry with TRANS = 'N' or 'n', the leading N--by--K part of * the array A must contain the matrix A, otherwise the leading * K--by--N part of the array A must contain the matrix A. * Unchanged on exit. * -* LDA - (input) INTEGER. +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. When TRANS = 'N' or 'n' * then LDA must be at least max( 1, n ), otherwise LDA must * be at least max( 1, k ). * Unchanged on exit. * -* BETA - (input) REAL. +* BETA (input) REAL * On entry, BETA specifies the scalar beta. * Unchanged on exit. * -* C - (input/output) COMPLEX array, dimension ( N*(N+1)/2 ). +* C (input/output) COMPLEX array, dimension (N*(N+1)/2) * On entry, the matrix A in RFP Format. RFP Format is * described by TRANSR, UPLO and N. Note that the imaginary * parts of the diagonal elements need not be set, they are diff --git a/SRC/chla_transtype.f b/SRC/chla_transtype.f index d1f44eed..a1eccc5c 100644 --- a/SRC/chla_transtype.f +++ b/SRC/chla_transtype.f @@ -4,6 +4,7 @@ * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- * October 2008 +* * .. Scalar Arguments .. INTEGER TRANS * .. diff --git a/SRC/cla_gbamv.f b/SRC/cla_gbamv.f index db54ccab..4615a21e 100644 --- a/SRC/cla_gbamv.f +++ b/SRC/cla_gbamv.f @@ -10,7 +10,7 @@ * -- Univ. of California Berkeley and NAG Ltd. -- * IMPLICIT NONE -* .. +* * .. Scalar Arguments .. REAL ALPHA, BETA INTEGER INCX, INCY, LDAB, M, N, KL, KU, TRANS @@ -43,7 +43,7 @@ * Arguments * ========== * -* TRANS - INTEGER +* TRANS (input) INTEGER * On entry, TRANS specifies the operation to be performed as * follows: * @@ -53,38 +53,38 @@ * * Unchanged on exit. * -* M - INTEGER +* M (input) INTEGER * On entry, M specifies the number of rows of the matrix A. * M must be at least zero. * Unchanged on exit. * -* N - INTEGER +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. * -* KL - INTEGER +* KL (input) INTEGER * The number of subdiagonals within the band of A. KL >= 0. * -* KU - INTEGER +* KU (input) INTEGER * The number of superdiagonals within the band of A. KU >= 0. * -* ALPHA - REAL +* ALPHA (input) REAL * On entry, ALPHA specifies the scalar alpha. * Unchanged on exit. * -* A - REAL array of DIMENSION ( LDA, n ) +* A (input) REAL array, dimension (LDA,n) * Before entry, the leading m by n part of the array A must * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, m ). * Unchanged on exit. * -* X - REAL array of DIMENSION at least +* X (input) REAL array, dimension at least * ( 1 + ( n - 1 )*abs( INCX ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( m - 1 )*abs( INCX ) ) otherwise. @@ -92,17 +92,17 @@ * vector x. * Unchanged on exit. * -* INCX - INTEGER +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. * -* BETA - REAL +* BETA (input) REAL * On entry, BETA specifies the scalar beta. When BETA is * supplied as zero then Y need not be set on input. * Unchanged on exit. * -* Y - REAL array of DIMENSION at least +* Y (input/output) REAL array, dimension at least * ( 1 + ( m - 1 )*abs( INCY ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( n - 1 )*abs( INCY ) ) otherwise. @@ -110,7 +110,7 @@ * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/cla_gbrpvgrw.f b/SRC/cla_gbrpvgrw.f index 20a3f15a..6f2115a7 100644 --- a/SRC/cla_gbrpvgrw.f +++ b/SRC/cla_gbrpvgrw.f @@ -98,4 +98,4 @@ END IF END DO CLA_GBRPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/cla_geamv.f b/SRC/cla_geamv.f index caee2470..bb2129f0 100644 --- a/SRC/cla_geamv.f +++ b/SRC/cla_geamv.f @@ -44,7 +44,7 @@ * Arguments * ========== * -* TRANS - INTEGER +* TRANS (input) INTEGER * On entry, TRANS specifies the operation to be performed as * follows: * @@ -54,32 +54,32 @@ * * Unchanged on exit. * -* M - INTEGER +* M (input) INTEGER * On entry, M specifies the number of rows of the matrix A. * M must be at least zero. * Unchanged on exit. * -* N - INTEGER +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. * -* ALPHA - REAL +* ALPHA (input) REAL * On entry, ALPHA specifies the scalar alpha. * Unchanged on exit. * -* A - COMPLEX array of DIMENSION ( LDA, n ) +* A (input) COMPLEX array, dimension (LDA,n) * Before entry, the leading m by n part of the array A must * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, m ). * Unchanged on exit. * -* X - COMPLEX array of DIMENSION at least +* X (input) COMPLEX array, dimension * ( 1 + ( n - 1 )*abs( INCX ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( m - 1 )*abs( INCX ) ) otherwise. @@ -87,17 +87,17 @@ * vector x. * Unchanged on exit. * -* INCX - INTEGER +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. * -* BETA - REAL +* BETA (input) REAL * On entry, BETA specifies the scalar beta. When BETA is * supplied as zero then Y need not be set on input. * Unchanged on exit. * -* Y - REAL array of DIMENSION at least +* Y (input/output) REAL array, dimension * ( 1 + ( m - 1 )*abs( INCY ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( n - 1 )*abs( INCY ) ) otherwise. @@ -105,7 +105,7 @@ * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/cla_heamv.f b/SRC/cla_heamv.f index d806ea29..eae5279e 100644 --- a/SRC/cla_heamv.f +++ b/SRC/cla_heamv.f @@ -42,7 +42,7 @@ * Arguments * ========== * -* UPLO - INTEGER +* UPLO (input) INTEGER * On entry, UPLO specifies whether the upper or lower * triangular part of the array A is to be referenced as * follows: @@ -55,12 +55,12 @@ * * Unchanged on exit. * -* N - INTEGER. +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. * -* ALPHA - REAL . +* ALPHA (input) REAL . * On entry, ALPHA specifies the scalar alpha. * Unchanged on exit. * @@ -69,35 +69,35 @@ * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER. +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, n ). * Unchanged on exit. * -* X - COMPLEX array of DIMENSION at least +* X (input) COMPLEX array, dimension * ( 1 + ( n - 1 )*abs( INCX ) ) * Before entry, the incremented array X must contain the * vector x. * Unchanged on exit. * -* INCX - INTEGER. +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. * -* BETA - REAL . +* BETA (input) REAL . * On entry, BETA specifies the scalar beta. When BETA is * supplied as zero then Y need not be set on input. * Unchanged on exit. * -* Y - REAL array of DIMENSION at least +* Y (input/output) REAL array, dimension * ( 1 + ( n - 1 )*abs( INCY ) ) * Before entry with BETA non-zero, the incremented array Y * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER. +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/cla_herpvgrw.f b/SRC/cla_herpvgrw.f index 6ef497dd..c94fac56 100644 --- a/SRC/cla_herpvgrw.f +++ b/SRC/cla_herpvgrw.f @@ -257,4 +257,4 @@ END IF CLA_HERPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/cla_lin_berr.f b/SRC/cla_lin_berr.f index 97118d1a..21daed27 100644 --- a/SRC/cla_lin_berr.f +++ b/SRC/cla_lin_berr.f @@ -48,7 +48,7 @@ * the matrix abs(op(A_s))*abs(Y) + abs(B_s). The matrices A, Y, and B * are from iterative refinement (see cla_gerfsx_extended.f). * -* RES (output) COMPLEX array, dimension (NRHS) +* BERR (output) COMPLEX array, dimension (NRHS) * The componentwise relative backward error from the formula above. * * ===================================================================== @@ -94,4 +94,4 @@ * END DO END DO - END SUBROUTINE + END diff --git a/SRC/cla_porpvgrw.f b/SRC/cla_porpvgrw.f index ded8ad6a..4aaba479 100644 --- a/SRC/cla_porpvgrw.f +++ b/SRC/cla_porpvgrw.f @@ -148,4 +148,4 @@ END IF CLA_PORPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/cla_rpvgrw.f b/SRC/cla_rpvgrw.f index 20d27ea3..d93040f8 100644 --- a/SRC/cla_rpvgrw.f +++ b/SRC/cla_rpvgrw.f @@ -84,4 +84,4 @@ END IF END DO CLA_RPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/cla_syamv.f b/SRC/cla_syamv.f index 3ea05931..61aa2a46 100644 --- a/SRC/cla_syamv.f +++ b/SRC/cla_syamv.f @@ -43,7 +43,7 @@ * Arguments * ========== * -* UPLO - INTEGER +* UPLO (input) INTEGER * On entry, UPLO specifies whether the upper or lower * triangular part of the array A is to be referenced as * follows: @@ -56,12 +56,12 @@ * * Unchanged on exit. * -* N - INTEGER. +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. * -* ALPHA - REAL . +* ALPHA (input) REAL . * On entry, ALPHA specifies the scalar alpha. * Unchanged on exit. * @@ -70,35 +70,35 @@ * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER. +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, n ). * Unchanged on exit. * -* X - COMPLEX array of DIMENSION at least +* X (input) COMPLEX array, dimension * ( 1 + ( n - 1 )*abs( INCX ) ) * Before entry, the incremented array X must contain the * vector x. * Unchanged on exit. * -* INCX - INTEGER. +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. * -* BETA - REAL . +* BETA (input) REAL . * On entry, BETA specifies the scalar beta. When BETA is * supplied as zero then Y need not be set on input. * Unchanged on exit. * -* Y - REAL array of DIMENSION at least +* Y (input/output) REAL array, dimension * ( 1 + ( n - 1 )*abs( INCY ) ) * Before entry with BETA non-zero, the incremented array Y * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER. +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/cla_syrpvgrw.f b/SRC/cla_syrpvgrw.f index d50676f9..9f7993d3 100644 --- a/SRC/cla_syrpvgrw.f +++ b/SRC/cla_syrpvgrw.f @@ -258,4 +258,4 @@ END IF CLA_SYRPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/cla_wwaddw.f b/SRC/cla_wwaddw.f index f5abec76..ed65a60c 100644 --- a/SRC/cla_wwaddw.f +++ b/SRC/cla_wwaddw.f @@ -31,10 +31,13 @@ * N (input) INTEGER * The length of vectors X, Y, and W. * -* X, Y (input/output) COMPLEX array, length N -* The doubled-single accumulation vector. +* X (input/output) COMPLEX array, dimension (N) +* The first part of the doubled-single accumulation vector. * -* W (input) COMPLEX array, length N +* Y (input/output) COMPLEX array, dimension (N) +* The second part of the doubled-single accumulation vector. +* +* W (input) COMPLEX array, dimension (N) * The vector to be added. * * ===================================================================== diff --git a/SRC/ctfsm.f b/SRC/ctfsm.f index 996d8da6..ec8db864 100644 --- a/SRC/ctfsm.f +++ b/SRC/ctfsm.f @@ -40,11 +40,11 @@ * Arguments * ========== * -* TRANSR - (input) CHARACTER +* TRANSR (input) CHARACTER * = 'N': The Normal Form of RFP A is stored; * = 'C': The Conjugate-transpose Form of RFP A is stored. * -* SIDE - (input) CHARACTER +* SIDE (input) CHARACTER * On entry, SIDE specifies whether op( A ) appears on the left * or right of X as follows: * @@ -54,7 +54,7 @@ * * Unchanged on exit. * -* UPLO - (input) CHARACTER +* UPLO (input) CHARACTER * On entry, UPLO specifies whether the RFP matrix A came from * an upper or lower triangular matrix as follows: * UPLO = 'U' or 'u' RFP A came from an upper triangular matrix @@ -62,7 +62,7 @@ * * Unchanged on exit. * -* TRANS - (input) CHARACTER +* TRANS (input) CHARACTER * On entry, TRANS specifies the form of op( A ) to be used * in the matrix multiplication as follows: * @@ -72,7 +72,7 @@ * * Unchanged on exit. * -* DIAG - (input) CHARACTER +* DIAG (input) CHARACTER * On entry, DIAG specifies whether or not RFP A is unit * triangular as follows: * @@ -83,23 +83,23 @@ * * Unchanged on exit. * -* M - (input) INTEGER. +* M (input) INTEGER * On entry, M specifies the number of rows of B. M must be at * least zero. * Unchanged on exit. * -* N - (input) INTEGER. +* N (input) INTEGER * On entry, N specifies the number of columns of B. N must be * at least zero. * Unchanged on exit. * -* ALPHA - (input) COMPLEX. +* ALPHA (input) COMPLEX * On entry, ALPHA specifies the scalar alpha. When alpha is * zero then A is not referenced and B need not be set before * entry. * Unchanged on exit. * -* A - (input) COMPLEX array, dimension ( N*(N+1)/2 ); +* A (input) COMPLEX array, dimension (N*(N+1)/2) * NT = N*(N+1)/2. On entry, the matrix A in RFP Format. * RFP Format is described by TRANSR, UPLO and N as follows: * If TRANSR='N' then RFP A is (0:N,0:K-1) when N is even; @@ -115,12 +115,12 @@ * even and is N when is odd. * See the Note below for more details. Unchanged on exit. * -* B - (input/ouptut) COMPLEX array, DIMENSION ( LDB, N ) +* B (input/ouptut) COMPLEX array, dimension (LDB,N) * Before entry, the leading m by n part of the array B must * contain the right-hand side matrix B, and on exit is * overwritten by the solution matrix X. * -* LDB - (input) INTEGER. +* LDB (input) INTEGER * On entry, LDB specifies the first dimension of B as declared * in the calling (sub) program. LDB must be at least * max( 1, m ). diff --git a/SRC/disnan.f b/SRC/disnan.f index d460d5b7..8e4e2a37 100644 --- a/SRC/disnan.f +++ b/SRC/disnan.f @@ -1,4 +1,4 @@ - LOGICAL FUNCTION DISNAN(DIN) + LOGICAL FUNCTION DISNAN( DIN) * * -- LAPACK auxiliary routine (version 3.2) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- @@ -6,7 +6,7 @@ * November 2006 * * .. Scalar Arguments .. - DOUBLE PRECISION DIN + DOUBLE PRECISION DIN * .. * * Purpose diff --git a/SRC/dla_gbamv.f b/SRC/dla_gbamv.f index 93471443..20e4644a 100644 --- a/SRC/dla_gbamv.f +++ b/SRC/dla_gbamv.f @@ -42,7 +42,7 @@ * Arguments * ========== * -* TRANS - INTEGER +* TRANS (input) INTEGER * On entry, TRANS specifies the operation to be performed as * follows: * @@ -52,20 +52,20 @@ * * Unchanged on exit. * -* M - INTEGER +* M (input) INTEGER * On entry, M specifies the number of rows of the matrix A. * M must be at least zero. * Unchanged on exit. * -* N - INTEGER +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. * -* KL - INTEGER +* KL (input) INTEGER * The number of subdiagonals within the band of A. KL >= 0. * -* KU - INTEGER +* KU (input) INTEGER * The number of superdiagonals within the band of A. KU >= 0. * * ALPHA - DOUBLE PRECISION @@ -77,13 +77,13 @@ * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, m ). * Unchanged on exit. * -* X - DOUBLE PRECISION array of DIMENSION at least +* X (input) DOUBLE PRECISION array, dimension * ( 1 + ( n - 1 )*abs( INCX ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( m - 1 )*abs( INCX ) ) otherwise. @@ -91,7 +91,7 @@ * vector x. * Unchanged on exit. * -* INCX - INTEGER +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. @@ -101,7 +101,7 @@ * supplied as zero then Y need not be set on input. * Unchanged on exit. * -* Y - DOUBLE PRECISION array of DIMENSION at least +* Y (input/output) DOUBLE PRECISION array, dimension * ( 1 + ( m - 1 )*abs( INCY ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( n - 1 )*abs( INCY ) ) otherwise. @@ -109,7 +109,7 @@ * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/dla_gbrpvgrw.f b/SRC/dla_gbrpvgrw.f index 7bac0d16..15757e1a 100644 --- a/SRC/dla_gbrpvgrw.f +++ b/SRC/dla_gbrpvgrw.f @@ -91,4 +91,4 @@ END IF END DO DLA_GBRPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/dla_geamv.f b/SRC/dla_geamv.f index a047f139..03bc06f4 100644 --- a/SRC/dla_geamv.f +++ b/SRC/dla_geamv.f @@ -42,7 +42,7 @@ * Arguments * ========== * -* TRANS - INTEGER +* TRANS (input) INTEGER * On entry, TRANS specifies the operation to be performed as * follows: * @@ -52,12 +52,12 @@ * * Unchanged on exit. * -* M - INTEGER +* M (input) INTEGER * On entry, M specifies the number of rows of the matrix A. * M must be at least zero. * Unchanged on exit. * -* N - INTEGER +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. @@ -71,13 +71,13 @@ * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, m ). * Unchanged on exit. * -* X - DOUBLE PRECISION array of DIMENSION at least +* X (input) DOUBLE PRECISION array, dimension * ( 1 + ( n - 1 )*abs( INCX ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( m - 1 )*abs( INCX ) ) otherwise. @@ -85,7 +85,7 @@ * vector x. * Unchanged on exit. * -* INCX - INTEGER +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. @@ -104,7 +104,7 @@ * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/dla_lin_berr.f b/SRC/dla_lin_berr.f index a23405a2..e204b69a 100644 --- a/SRC/dla_lin_berr.f +++ b/SRC/dla_lin_berr.f @@ -51,7 +51,7 @@ * the matrix abs(op(A_s))*abs(Y) + abs(B_s). The matrices A, Y, and B * are from iterative refinement (see dla_gerfsx_extended.f). * -* RES (output) DOUBLE PRECISION array, dimension (NRHS) +* BERR (output) DOUBLE PRECISION array, dimension (NRHS) * The component-wise relative backward error from the formula above. * * ===================================================================== @@ -90,4 +90,4 @@ * END DO END DO - END SUBROUTINE + END diff --git a/SRC/dla_porpvgrw.f b/SRC/dla_porpvgrw.f index 069e2d60..de4512a7 100644 --- a/SRC/dla_porpvgrw.f +++ b/SRC/dla_porpvgrw.f @@ -142,4 +142,4 @@ END IF DLA_PORPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/dla_rpvgrw.f b/SRC/dla_rpvgrw.f index f2a80678..e052dc0e 100644 --- a/SRC/dla_rpvgrw.f +++ b/SRC/dla_rpvgrw.f @@ -77,4 +77,4 @@ END IF END DO DLA_RPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/dla_syamv.f b/SRC/dla_syamv.f index f0ce64fd..42cb0989 100644 --- a/SRC/dla_syamv.f +++ b/SRC/dla_syamv.f @@ -41,7 +41,7 @@ * Arguments * ========== * -* UPLO - INTEGER +* UPLO (input) INTEGER * On entry, UPLO specifies whether the upper or lower * triangular part of the array A is to be referenced as * follows: @@ -54,7 +54,7 @@ * * Unchanged on exit. * -* N - INTEGER. +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. @@ -68,19 +68,19 @@ * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER. +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, n ). * Unchanged on exit. * -* X - DOUBLE PRECISION array of DIMENSION at least +* X (input) DOUBLE PRECISION array, dimension * ( 1 + ( n - 1 )*abs( INCX ) ) * Before entry, the incremented array X must contain the * vector x. * Unchanged on exit. * -* INCX - INTEGER. +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. @@ -90,13 +90,13 @@ * supplied as zero then Y need not be set on input. * Unchanged on exit. * -* Y - DOUBLE PRECISION array of DIMENSION at least +* Y (input/output) DOUBLE PRECISION array, dimension * ( 1 + ( n - 1 )*abs( INCY ) ) * Before entry with BETA non-zero, the incremented array Y * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER. +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/dla_syrpvgrw.f b/SRC/dla_syrpvgrw.f index d6ea5c22..e570c908 100644 --- a/SRC/dla_syrpvgrw.f +++ b/SRC/dla_syrpvgrw.f @@ -248,4 +248,4 @@ END IF DLA_SYRPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/dla_wwaddw.f b/SRC/dla_wwaddw.f index c80c40f3..08f85fc7 100644 --- a/SRC/dla_wwaddw.f +++ b/SRC/dla_wwaddw.f @@ -31,10 +31,13 @@ * N (input) INTEGER * The length of vectors X, Y, and W. * -* X, Y (input/output) DOUBLE PRECISION array, length N -* The doubled-single accumulation vector. +* X (input/output) DOUBLE PRECISION array, dimension (N) +* The first part of the doubled-single accumulation vector. * -* W (input) DOUBLE PRECISION array, length N +* Y (input/output) DOUBLE PRECISION array, dimension (N) +* The second part of the doubled-single accumulation vector. +* +* W (input) DOUBLE PRECISION array, dimension (N) * The vector to be added. * * ===================================================================== diff --git a/SRC/dlaisnan.f b/SRC/dlaisnan.f index d6f8bd6b..b9ecc43b 100644 --- a/SRC/dlaisnan.f +++ b/SRC/dlaisnan.f @@ -1,4 +1,4 @@ - LOGICAL FUNCTION DLAISNAN(DIN1,DIN2) + LOGICAL FUNCTION DLAISNAN( DIN1, DIN2) * * -- LAPACK auxiliary routine (version 3.2) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- @@ -6,7 +6,7 @@ * November 2006 * * .. Scalar Arguments .. - DOUBLE PRECISION DIN1,DIN2 + DOUBLE PRECISION DIN1, DIN2 * .. * * Purpose @@ -30,6 +30,7 @@ * ========= * * DIN1 (input) DOUBLE PRECISION +* * DIN2 (input) DOUBLE PRECISION * Two numbers to compare for inequality. * diff --git a/SRC/dlarrf.f b/SRC/dlarrf.f index 7d092fd6..635cf11b 100644 --- a/SRC/dlarrf.f +++ b/SRC/dlarrf.f @@ -48,15 +48,18 @@ * CLEND (input) INTEGER * The index of the last eigenvalue in the cluster. * -* W (input) DOUBLE PRECISION array, dimension >= (CLEND-CLSTRT+1) +* W (input) DOUBLE PRECISION array, dimension +* dimension is >= (CLEND-CLSTRT+1) * The eigenvalue APPROXIMATIONS of L D L^T in ascending order. * W( CLSTRT ) through W( CLEND ) form the cluster of relatively * close eigenalues. * -* WGAP (input/output) DOUBLE PRECISION array, dimension >= (CLEND-CLSTRT+1) +* WGAP (input/output) DOUBLE PRECISION array, dimension +* dimension is >= (CLEND-CLSTRT+1) * The separation from the right neighbor eigenvalue in W. * -* WERR (input) DOUBLE PRECISION array, dimension >= (CLEND-CLSTRT+1) +* WERR (input) DOUBLE PRECISION array, dimension +* dimension is >= (CLEND-CLSTRT+1) * WERR contain the semiwidth of the uncertainty * interval of the corresponding eigenvalue APPROXIMATION in W * @@ -64,6 +67,7 @@ * Gerschgorin intervals * * CLGAPL (input) DOUBLE PRECISION +* * CLGAPR (input) DOUBLE PRECISION * absolute gap on each end of the cluster. * Set by the calling routine to protect against shifts too close diff --git a/SRC/dlasdt.f b/SRC/dlasdt.f index 38aa6dad..22629432 100644 --- a/SRC/dlasdt.f +++ b/SRC/dlasdt.f @@ -40,7 +40,7 @@ * NDIMR (output) INTEGER array, dimension ( N ) * On exit, row dimensions of right children. * -* MSUB (input) INTEGER. +* MSUB (input) INTEGER * On entry, the maximum row dimension each subproblem at the * bottom of the tree can be of. * diff --git a/SRC/dlasq3.f b/SRC/dlasq3.f index e5ea244d..d93a5013 100644 --- a/SRC/dlasq3.f +++ b/SRC/dlasq3.f @@ -74,7 +74,20 @@ * TTYPE (input/output) INTEGER * Shift type. * -* DMIN1, DMIN2, DN, DN1, DN2, G, TAU (input/output) DOUBLE PRECISION +* DMIN1 (input/output) DOUBLE PRECISION +* +* DMIN2 (input/output) DOUBLE PRECISION +* +* DN (input/output) DOUBLE PRECISION +* +* DN1 (input/output) DOUBLE PRECISION +* +* DN2 (input/output) DOUBLE PRECISION +* +* G (input/output) DOUBLE PRECISION +* +* TAU (input/output) DOUBLE PRECISION +* * These are passed as arguments in order to save their values * between calls to DLASQ3. * diff --git a/SRC/dsfrk.f b/SRC/dsfrk.f index 30f82015..aa1518af 100644 --- a/SRC/dsfrk.f +++ b/SRC/dsfrk.f @@ -39,11 +39,11 @@ * Arguments * ========== * -* TRANSR (input) CHARACTER +* TRANSR (input) CHARACTER * = 'N': The Normal Form of RFP A is stored; * = 'T': The Transpose Form of RFP A is stored. * -* UPLO - (input) CHARACTER +* UPLO (input) CHARACTER * On entry, UPLO specifies whether the upper or lower * triangular part of the array C is to be referenced as * follows: @@ -56,7 +56,7 @@ * * Unchanged on exit. * -* TRANS - (input) CHARACTER +* TRANS (input) CHARACTER * On entry, TRANS specifies the operation to be performed as * follows: * @@ -66,42 +66,43 @@ * * Unchanged on exit. * -* N - (input) INTEGER. +* N (input) INTEGER * On entry, N specifies the order of the matrix C. N must be * at least zero. * Unchanged on exit. * -* K - (input) INTEGER. +* K (input) INTEGER * On entry with TRANS = 'N' or 'n', K specifies the number * of columns of the matrix A, and on entry with TRANS = 'T' * or 't', K specifies the number of rows of the matrix A. K * must be at least zero. * Unchanged on exit. * -* ALPHA - (input) DOUBLE PRECISION. +* ALPHA (input) DOUBLE PRECISION * On entry, ALPHA specifies the scalar alpha. * Unchanged on exit. * -* A - (input) DOUBLE PRECISION array of DIMENSION ( LDA, ka ), where KA +* A (input) DOUBLE PRECISION array, dimension (LDA,ka) +* where KA * is K when TRANS = 'N' or 'n', and is N otherwise. Before * entry with TRANS = 'N' or 'n', the leading N--by--K part of * the array A must contain the matrix A, otherwise the leading * K--by--N part of the array A must contain the matrix A. * Unchanged on exit. * -* LDA - (input) INTEGER. +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. When TRANS = 'N' or 'n' * then LDA must be at least max( 1, n ), otherwise LDA must * be at least max( 1, k ). * Unchanged on exit. * -* BETA - (input) DOUBLE PRECISION. +* BETA (input) DOUBLE PRECISION * On entry, BETA specifies the scalar beta. * Unchanged on exit. * * -* C - (input/output) DOUBLE PRECISION array, dimension ( NT ); +* C (input/output) DOUBLE PRECISION array, dimension (NT) * NT = N*(N+1)/2. On entry, the symmetric matrix C in RFP * Format. RFP Format is described by TRANSR, UPLO and N. * diff --git a/SRC/dtfsm.f b/SRC/dtfsm.f index 668b3e8d..c2e290f3 100644 --- a/SRC/dtfsm.f +++ b/SRC/dtfsm.f @@ -40,11 +40,11 @@ * Arguments * ========== * -* TRANSR - (input) CHARACTER +* TRANSR (input) CHARACTER * = 'N': The Normal Form of RFP A is stored; * = 'T': The Transpose Form of RFP A is stored. * -* SIDE - (input) CHARACTER +* SIDE (input) CHARACTER * On entry, SIDE specifies whether op( A ) appears on the left * or right of X as follows: * @@ -54,7 +54,7 @@ * * Unchanged on exit. * -* UPLO - (input) CHARACTER +* UPLO (input) CHARACTER * On entry, UPLO specifies whether the RFP matrix A came from * an upper or lower triangular matrix as follows: * UPLO = 'U' or 'u' RFP A came from an upper triangular matrix @@ -62,7 +62,7 @@ * * Unchanged on exit. * -* TRANS - (input) CHARACTER +* TRANS (input) CHARACTER * On entry, TRANS specifies the form of op( A ) to be used * in the matrix multiplication as follows: * @@ -72,7 +72,7 @@ * * Unchanged on exit. * -* DIAG - (input) CHARACTER +* DIAG (input) CHARACTER * On entry, DIAG specifies whether or not RFP A is unit * triangular as follows: * @@ -83,23 +83,23 @@ * * Unchanged on exit. * -* M - (input) INTEGER. +* M (input) INTEGER * On entry, M specifies the number of rows of B. M must be at * least zero. * Unchanged on exit. * -* N - (input) INTEGER. +* N (input) INTEGER * On entry, N specifies the number of columns of B. N must be * at least zero. * Unchanged on exit. * -* ALPHA - (input) DOUBLE PRECISION. +* ALPHA (input) DOUBLE PRECISION * On entry, ALPHA specifies the scalar alpha. When alpha is * zero then A is not referenced and B need not be set before * entry. * Unchanged on exit. * -* A - (input) DOUBLE PRECISION array, dimension (NT); +* A (input) DOUBLE PRECISION array, dimension (NT) * NT = N*(N+1)/2. On entry, the matrix A in RFP Format. * RFP Format is described by TRANSR, UPLO and N as follows: * If TRANSR='N' then RFP A is (0:N,0:K-1) when N is even; @@ -115,12 +115,12 @@ * even and is N when is odd. * See the Note below for more details. Unchanged on exit. * -* B - (input/ouptut) DOUBLE PRECISION array, DIMENSION (LDB,N) +* B (input/output) DOUBLE PRECISION array, dimension (LDB,N) * Before entry, the leading m by n part of the array B must * contain the right-hand side matrix B, and on exit is * overwritten by the solution matrix X. * -* LDB - (input) INTEGER. +* LDB (input) INTEGER * On entry, LDB specifies the first dimension of B as declared * in the calling (sub) program. LDB must be at least * max( 1, m ). diff --git a/SRC/ilaclc.f b/SRC/ilaclc.f index 577d5282..b2e026a7 100644 --- a/SRC/ilaclc.f +++ b/SRC/ilaclc.f @@ -1,4 +1,4 @@ - INTEGER FUNCTION ILACLC(M, N, A, LDA) + INTEGER FUNCTION ILACLC( M, N, A, LDA ) IMPLICIT NONE * * -- LAPACK auxiliary routine (version 3.2.1) -- @@ -60,4 +60,4 @@ END DO END IF RETURN - END FUNCTION + END diff --git a/SRC/ilaclr.f b/SRC/ilaclr.f index 61682191..06536b1b 100644 --- a/SRC/ilaclr.f +++ b/SRC/ilaclr.f @@ -1,4 +1,4 @@ - INTEGER FUNCTION ILACLR(M, N, A, LDA) + INTEGER FUNCTION ILACLR( M, N, A, LDA ) IMPLICIT NONE * * -- LAPACK auxiliary routine (version 3.2.1) -- @@ -62,4 +62,4 @@ END DO END IF RETURN - END FUNCTION + END diff --git a/SRC/iladlc.f b/SRC/iladlc.f index 66b141ee..0097dc60 100644 --- a/SRC/iladlc.f +++ b/SRC/iladlc.f @@ -1,4 +1,4 @@ - INTEGER FUNCTION ILADLC(M, N, A, LDA) + INTEGER FUNCTION ILADLC( M, N, A, LDA ) IMPLICIT NONE * * -- LAPACK auxiliary routine (version 3.2.1) -- @@ -60,4 +60,4 @@ END DO END IF RETURN - END FUNCTION + END diff --git a/SRC/iladlr.f b/SRC/iladlr.f index 8b12984b..36bde74d 100644 --- a/SRC/iladlr.f +++ b/SRC/iladlr.f @@ -1,4 +1,4 @@ - INTEGER FUNCTION ILADLR(M, N, A, LDA) + INTEGER FUNCTION ILADLR( M, N, A, LDA ) IMPLICIT NONE * * -- LAPACK auxiliary routine (version 3.2.1) -- @@ -62,4 +62,4 @@ END DO END IF RETURN - END FUNCTION + END diff --git a/SRC/ilaslc.f b/SRC/ilaslc.f index b7f77e8d..78487f7e 100644 --- a/SRC/ilaslc.f +++ b/SRC/ilaslc.f @@ -1,4 +1,4 @@ - INTEGER FUNCTION ILASLC(M, N, A, LDA) + INTEGER FUNCTION ILASLC( M, N, A, LDA ) IMPLICIT NONE * * -- LAPACK auxiliary routine (version 3.2.1) -- @@ -60,4 +60,4 @@ END DO END IF RETURN - END FUNCTION + END diff --git a/SRC/ilaslr.f b/SRC/ilaslr.f index afd3dba9..94052172 100644 --- a/SRC/ilaslr.f +++ b/SRC/ilaslr.f @@ -1,4 +1,4 @@ - INTEGER FUNCTION ILASLR(M, N, A, LDA) + INTEGER FUNCTION ILASLR( M, N, A, LDA ) IMPLICIT NONE * * -- LAPACK auxiliary routine (version 3.2.1) -- @@ -62,4 +62,4 @@ END DO END IF RETURN - END FUNCTION + END diff --git a/SRC/ilazlc.f b/SRC/ilazlc.f index 229e0997..ae616b10 100644 --- a/SRC/ilazlc.f +++ b/SRC/ilazlc.f @@ -1,4 +1,4 @@ - INTEGER FUNCTION ILAZLC(M, N, A, LDA) + INTEGER FUNCTION ILAZLC( M, N, A, LDA ) IMPLICIT NONE * * -- LAPACK auxiliary routine (version 3.2.1) -- @@ -60,4 +60,4 @@ END DO END IF RETURN - END FUNCTION + END diff --git a/SRC/ilazlr.f b/SRC/ilazlr.f index 1113fd2f..4615fadb 100644 --- a/SRC/ilazlr.f +++ b/SRC/ilazlr.f @@ -1,4 +1,4 @@ - INTEGER FUNCTION ILAZLR(M, N, A, LDA) + INTEGER FUNCTION ILAZLR( M, N, A, LDA ) IMPLICIT NONE * * -- LAPACK auxiliary routine (version 3.2.1) -- @@ -62,4 +62,4 @@ END DO END IF RETURN - END FUNCTION + END diff --git a/SRC/sisnan.f b/SRC/sisnan.f index 4fdf5f06..0b46abc3 100644 --- a/SRC/sisnan.f +++ b/SRC/sisnan.f @@ -1,4 +1,4 @@ - LOGICAL FUNCTION SISNAN(SIN) + LOGICAL FUNCTION SISNAN( SIN) * * -- LAPACK auxiliary routine (version 3.2) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- @@ -6,7 +6,7 @@ * November 2006 * * .. Scalar Arguments .. - REAL SIN + REAL SIN * .. * * Purpose diff --git a/SRC/sla_gbamv.f b/SRC/sla_gbamv.f index d004243f..0c0f1c50 100644 --- a/SRC/sla_gbamv.f +++ b/SRC/sla_gbamv.f @@ -42,7 +42,7 @@ * Arguments * ========== * -* TRANS - INTEGER +* TRANS (input) INTEGER * On entry, TRANS specifies the operation to be performed as * follows: * @@ -52,23 +52,23 @@ * * Unchanged on exit. * -* M - INTEGER +* M (input) INTEGER * On entry, M specifies the number of rows of the matrix A. * M must be at least zero. * Unchanged on exit. * -* N - INTEGER +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. * -* KL - INTEGER +* KL (input) INTEGER * The number of subdiagonals within the band of A. KL >= 0. * -* KU - INTEGER +* KU (input) INTEGER * The number of superdiagonals within the band of A. KU >= 0. * -* ALPHA - REAL +* ALPHA (input) REAL * On entry, ALPHA specifies the scalar alpha. * Unchanged on exit. * @@ -77,13 +77,13 @@ * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, m ). * Unchanged on exit. * -* X - REAL array of DIMENSION at least +* X (input) REAL array, dimension * ( 1 + ( n - 1 )*abs( INCX ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( m - 1 )*abs( INCX ) ) otherwise. @@ -91,17 +91,17 @@ * vector x. * Unchanged on exit. * -* INCX - INTEGER +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. * -* BETA - REAL +* BETA (input) REAL * On entry, BETA specifies the scalar beta. When BETA is * supplied as zero then Y need not be set on input. * Unchanged on exit. * -* Y - REAL array of DIMENSION at least +* Y (input/output) REAL array, dimension * ( 1 + ( m - 1 )*abs( INCY ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( n - 1 )*abs( INCY ) ) otherwise. @@ -109,7 +109,7 @@ * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/sla_gbrpvgrw.f b/SRC/sla_gbrpvgrw.f index 98b389d6..0ddbb629 100644 --- a/SRC/sla_gbrpvgrw.f +++ b/SRC/sla_gbrpvgrw.f @@ -91,4 +91,4 @@ END IF END DO SLA_GBRPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/sla_geamv.f b/SRC/sla_geamv.f index cf2ab0e1..dc319cd9 100644 --- a/SRC/sla_geamv.f +++ b/SRC/sla_geamv.f @@ -42,7 +42,7 @@ * Arguments * ========== * -* TRANS - INTEGER +* TRANS (input) INTEGER * On entry, TRANS specifies the operation to be performed as * follows: * @@ -52,17 +52,17 @@ * * Unchanged on exit. * -* M - INTEGER +* M (input) INTEGER * On entry, M specifies the number of rows of the matrix A. * M must be at least zero. * Unchanged on exit. * -* N - INTEGER +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. * -* ALPHA - REAL +* ALPHA (input) REAL * On entry, ALPHA specifies the scalar alpha. * Unchanged on exit. * @@ -71,13 +71,13 @@ * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, m ). * Unchanged on exit. * -* X - REAL array of DIMENSION at least +* X (input) REAL array, dimension * ( 1 + ( n - 1 )*abs( INCX ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( m - 1 )*abs( INCX ) ) otherwise. @@ -85,12 +85,12 @@ * vector x. * Unchanged on exit. * -* INCX - INTEGER +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. * -* BETA - REAL +* BETA (input) REAL * On entry, BETA specifies the scalar beta. When BETA is * supplied as zero then Y need not be set on input. * Unchanged on exit. @@ -104,7 +104,7 @@ * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/sla_lin_berr.f b/SRC/sla_lin_berr.f index 7a67e5f8..d228187b 100644 --- a/SRC/sla_lin_berr.f +++ b/SRC/sla_lin_berr.f @@ -51,7 +51,7 @@ * the matrix abs(op(A_s))*abs(Y) + abs(B_s). The matrices A, Y, and B * are from iterative refinement (see sla_gerfsx_extended.f). * -* RES (output) REAL array, dimension (NRHS) +* BERR (output) REAL array, dimension (NRHS) * The componentwise relative backward error from the formula above. * * ===================================================================== @@ -90,4 +90,4 @@ * END DO END DO - END SUBROUTINE + END diff --git a/SRC/sla_porpvgrw.f b/SRC/sla_porpvgrw.f index 5ea9578c..203db17f 100644 --- a/SRC/sla_porpvgrw.f +++ b/SRC/sla_porpvgrw.f @@ -141,4 +141,4 @@ END IF SLA_PORPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/sla_rpvgrw.f b/SRC/sla_rpvgrw.f index 327e3527..6c025675 100644 --- a/SRC/sla_rpvgrw.f +++ b/SRC/sla_rpvgrw.f @@ -77,4 +77,4 @@ END IF END DO SLA_RPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/sla_syamv.f b/SRC/sla_syamv.f index f6a659b1..55eee229 100644 --- a/SRC/sla_syamv.f +++ b/SRC/sla_syamv.f @@ -41,7 +41,7 @@ * Arguments * ========== * -* UPLO - INTEGER +* UPLO (input) INTEGER * On entry, UPLO specifies whether the upper or lower * triangular part of the array A is to be referenced as * follows: @@ -54,12 +54,12 @@ * * Unchanged on exit. * -* N - INTEGER. +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. * -* ALPHA - REAL . +* ALPHA (input) REAL . * On entry, ALPHA specifies the scalar alpha. * Unchanged on exit. * @@ -68,35 +68,35 @@ * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER. +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, n ). * Unchanged on exit. * -* X - REAL array of DIMENSION at least +* X (input) REAL array, dimension * ( 1 + ( n - 1 )*abs( INCX ) ) * Before entry, the incremented array X must contain the * vector x. * Unchanged on exit. * -* INCX - INTEGER. +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. * -* BETA - REAL . +* BETA (input) REAL . * On entry, BETA specifies the scalar beta. When BETA is * supplied as zero then Y need not be set on input. * Unchanged on exit. * -* Y - REAL array of DIMENSION at least +* Y (input/output) REAL array, dimension * ( 1 + ( n - 1 )*abs( INCY ) ) * Before entry with BETA non-zero, the incremented array Y * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER. +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/sla_syrpvgrw.f b/SRC/sla_syrpvgrw.f index 16c1714a..592495d7 100644 --- a/SRC/sla_syrpvgrw.f +++ b/SRC/sla_syrpvgrw.f @@ -248,4 +248,4 @@ END IF SLA_SYRPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/sla_wwaddw.f b/SRC/sla_wwaddw.f index 88152557..ede96709 100644 --- a/SRC/sla_wwaddw.f +++ b/SRC/sla_wwaddw.f @@ -31,10 +31,13 @@ * N (input) INTEGER * The length of vectors X, Y, and W. * -* X, Y (input/output) REAL array, length N -* The doubled-single accumulation vector. +* X (input/output) REAL array, dimension (N) +* The first part of the doubled-single accumulation vector. * -* W (input) REAL array, length N +* Y (input/output) REAL array, dimension (N) +* The second part of the doubled-single accumulation vector. +* +* W (input) REAL array, dimension (N) * The vector to be added. * * ===================================================================== diff --git a/SRC/slaisnan.f b/SRC/slaisnan.f index cc8cf4f3..d14cbe64 100644 --- a/SRC/slaisnan.f +++ b/SRC/slaisnan.f @@ -1,4 +1,4 @@ - LOGICAL FUNCTION SLAISNAN(SIN1,SIN2) + LOGICAL FUNCTION SLAISNAN( SIN1, SIN2) * * -- LAPACK auxiliary routine (version 3.2) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- @@ -6,7 +6,7 @@ * November 2006 * * .. Scalar Arguments .. - REAL SIN1,SIN2 + REAL SIN1, SIN2 * .. * * Purpose @@ -30,6 +30,7 @@ * ========= * * SIN1 (input) REAL +* * SIN2 (input) REAL * Two numbers to compare for inequality. * diff --git a/SRC/slaneg.f b/SRC/slaneg.f index f1e8aa88..2b6f7bd9 100644 --- a/SRC/slaneg.f +++ b/SRC/slaneg.f @@ -1,6 +1,5 @@ - FUNCTION SLANEG( N, D, LLD, SIGMA, PIVMIN, R ) + INTEGER FUNCTION SLANEG( N, D, LLD, SIGMA, PIVMIN, R ) IMPLICIT NONE - INTEGER SLANEG * * -- LAPACK auxiliary routine (version 3.2) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- diff --git a/SRC/slarrf.f b/SRC/slarrf.f index a8d8e174..1d74b15e 100644 --- a/SRC/slarrf.f +++ b/SRC/slarrf.f @@ -32,14 +32,14 @@ * N (input) INTEGER * The order of the matrix (subblock, if the matrix splitted). * -* D (input) REAL array, dimension (N) +* D (input) REAL array, dimension (N) * The N diagonal elements of the diagonal matrix D. * -* L (input) REAL array, dimension (N-1) +* L (input) REAL array, dimension (N-1) * The (N-1) subdiagonal elements of the unit bidiagonal * matrix L. * -* LD (input) REAL array, dimension (N-1) +* LD (input) REAL array, dimension (N-1) * The (N-1) elements L(i)*D(i). * * CLSTRT (input) INTEGER @@ -48,22 +48,29 @@ * CLEND (input) INTEGER * The index of the last eigenvalue in the cluster. * -* W (input) REAL array, dimension >= (CLEND-CLSTRT+1) +* W (input) REAL array, dimension +* dimension is >= (CLEND-CLSTRT+1) * The eigenvalue APPROXIMATIONS of L D L^T in ascending order. * W( CLSTRT ) through W( CLEND ) form the cluster of relatively * close eigenalues. * -* WGAP (input/output) REAL array, dimension >= (CLEND-CLSTRT+1) +* WGAP (input/output) REAL array, dimension +* dimension is >= (CLEND-CLSTRT+1) * The separation from the right neighbor eigenvalue in W. * -* WERR (input) REAL array, dimension >= (CLEND-CLSTRT+1) +* WERR (input) REAL array, dimension +* dimension is >= (CLEND-CLSTRT+1) * WERR contain the semiwidth of the uncertainty * interval of the corresponding eigenvalue APPROXIMATION in W * -* SPDIAM (input) estimate of the spectral diameter obtained from the +* SPDIAM (input) REAL +* estimate of the spectral diameter obtained from the * Gerschgorin intervals * -* CLGAPL, CLGAPR (input) absolute gap on each end of the cluster. +* CLGAPL (input) REAL +* +* CLGAPR (input) REAL +* absolute gap on each end of the cluster. * Set by the calling routine to protect against shifts too close * to eigenvalues outside the cluster. * @@ -80,7 +87,7 @@ * The first (N-1) elements of LPLUS contain the subdiagonal * elements of the unit bidiagonal matrix L(+). * -* WORK (workspace) REAL array, dimension (2*N) +* WORK (workspace) REAL array, dimension (2*N) * Workspace. * * Further Details diff --git a/SRC/slasdt.f b/SRC/slasdt.f index fbff52d9..18f37430 100644 --- a/SRC/slasdt.f +++ b/SRC/slasdt.f @@ -40,7 +40,7 @@ * NDIMR (output) INTEGER array, dimension ( N ) * On exit, row dimensions of right children. * -* MSUB (input) INTEGER. +* MSUB (input) INTEGER * On entry, the maximum row dimension each subproblem at the * bottom of the tree can be of. * diff --git a/SRC/slasq3.f b/SRC/slasq3.f index d74cdb64..d677fddb 100644 --- a/SRC/slasq3.f +++ b/SRC/slasq3.f @@ -74,7 +74,20 @@ * TTYPE (input/output) INTEGER * Shift type. * -* DMIN1, DMIN2, DN, DN1, DN2, G, TAU (input/output) REAL +* DMIN1 (input/output) REAL +* +* DMIN2 (input/output) REAL +* +* DN (input/output) REAL +* +* DN1 (input/output) REAL +* +* DN2 (input/output) REAL +* +* G (input/output) REAL +* +* TAU (input/output) REAL +* * These are passed as arguments in order to save their values * between calls to SLASQ3. * diff --git a/SRC/ssfrk.f b/SRC/ssfrk.f index a42ccbf1..75ef1d6f 100644 --- a/SRC/ssfrk.f +++ b/SRC/ssfrk.f @@ -39,11 +39,11 @@ * Arguments * ========== * -* TRANSR (input) CHARACTER +* TRANSR (input) CHARACTER * = 'N': The Normal Form of RFP A is stored; * = 'T': The Transpose Form of RFP A is stored. * -* UPLO - (input) CHARACTER +* UPLO (input) CHARACTER * On entry, UPLO specifies whether the upper or lower * triangular part of the array C is to be referenced as * follows: @@ -56,7 +56,7 @@ * * Unchanged on exit. * -* TRANS - (input) CHARACTER +* TRANS (input) CHARACTER * On entry, TRANS specifies the operation to be performed as * follows: * @@ -66,42 +66,43 @@ * * Unchanged on exit. * -* N - (input) INTEGER. +* N (input) INTEGER * On entry, N specifies the order of the matrix C. N must be * at least zero. * Unchanged on exit. * -* K - (input) INTEGER. +* K (input) INTEGER * On entry with TRANS = 'N' or 'n', K specifies the number * of columns of the matrix A, and on entry with TRANS = 'T' * or 't', K specifies the number of rows of the matrix A. K * must be at least zero. * Unchanged on exit. * -* ALPHA - (input) REAL. +* ALPHA (input) REAL * On entry, ALPHA specifies the scalar alpha. * Unchanged on exit. * -* A - (input) REAL array of DIMENSION ( LDA, ka ), where KA +* A (input) REAL array of DIMENSION (LDA,ka) +* where KA * is K when TRANS = 'N' or 'n', and is N otherwise. Before * entry with TRANS = 'N' or 'n', the leading N--by--K part of * the array A must contain the matrix A, otherwise the leading * K--by--N part of the array A must contain the matrix A. * Unchanged on exit. * -* LDA - (input) INTEGER. +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. When TRANS = 'N' or 'n' * then LDA must be at least max( 1, n ), otherwise LDA must * be at least max( 1, k ). * Unchanged on exit. * -* BETA - (input) REAL. +* BETA (input) REAL * On entry, BETA specifies the scalar beta. * Unchanged on exit. * * -* C - (input/output) REAL array, dimension ( NT ); +* C (input/output) REAL array, dimension (NT) * NT = N*(N+1)/2. On entry, the symmetric matrix C in RFP * Format. RFP Format is described by TRANSR, UPLO and N. * diff --git a/SRC/stfsm.f b/SRC/stfsm.f index d06c5e92..c9cab77e 100644 --- a/SRC/stfsm.f +++ b/SRC/stfsm.f @@ -40,11 +40,11 @@ * Arguments * ========== * -* TRANSR - (input) CHARACTER +* TRANSR (input) CHARACTER * = 'N': The Normal Form of RFP A is stored; * = 'T': The Transpose Form of RFP A is stored. * -* SIDE - (input) CHARACTER +* SIDE (input) CHARACTER * On entry, SIDE specifies whether op( A ) appears on the left * or right of X as follows: * @@ -54,7 +54,7 @@ * * Unchanged on exit. * -* UPLO - (input) CHARACTER +* UPLO (input) CHARACTER * On entry, UPLO specifies whether the RFP matrix A came from * an upper or lower triangular matrix as follows: * UPLO = 'U' or 'u' RFP A came from an upper triangular matrix @@ -62,7 +62,7 @@ * * Unchanged on exit. * -* TRANS - (input) CHARACTER +* TRANS (input) CHARACTER * On entry, TRANS specifies the form of op( A ) to be used * in the matrix multiplication as follows: * @@ -72,7 +72,7 @@ * * Unchanged on exit. * -* DIAG - (input) CHARACTER +* DIAG (input) CHARACTER * On entry, DIAG specifies whether or not RFP A is unit * triangular as follows: * @@ -83,23 +83,23 @@ * * Unchanged on exit. * -* M - (input) INTEGER. +* M (input) INTEGER * On entry, M specifies the number of rows of B. M must be at * least zero. * Unchanged on exit. * -* N - (input) INTEGER. +* N (input) INTEGER * On entry, N specifies the number of columns of B. N must be * at least zero. * Unchanged on exit. * -* ALPHA - (input) REAL. +* ALPHA (input) REAL * On entry, ALPHA specifies the scalar alpha. When alpha is * zero then A is not referenced and B need not be set before * entry. * Unchanged on exit. * -* A - (input) REAL array, dimension (NT); +* A (input) REAL array, dimension (NT) * NT = N*(N+1)/2. On entry, the matrix A in RFP Format. * RFP Format is described by TRANSR, UPLO and N as follows: * If TRANSR='N' then RFP A is (0:N,0:K-1) when N is even; @@ -115,12 +115,12 @@ * even and is N when is odd. * See the Note below for more details. Unchanged on exit. * -* B - (input/ouptut) REAL array, DIMENSION (LDB,N) +* B (input/ouptut) REAL array, DIMENSION (LDB,N) * Before entry, the leading m by n part of the array B must * contain the right-hand side matrix B, and on exit is * overwritten by the solution matrix X. * -* LDB - (input) INTEGER. +* LDB (input) INTEGER * On entry, LDB specifies the first dimension of B as declared * in the calling (sub) program. LDB must be at least * max( 1, m ). diff --git a/SRC/xerbla_array.f b/SRC/xerbla_array.f index 3e2ccd14..3743ab1e 100644 --- a/SRC/xerbla_array.f +++ b/SRC/xerbla_array.f @@ -1,4 +1,4 @@ - SUBROUTINE XERBLA_ARRAY(SRNAME_ARRAY, SRNAME_LEN, INFO) + SUBROUTINE XERBLA_ARRAY( SRNAME_ARRAY, SRNAME_LEN, INFO) * * -- LAPACK auxiliary routine (version 3.2.1) -- * diff --git a/SRC/zhfrk.f b/SRC/zhfrk.f index 44cba760..46b97239 100644 --- a/SRC/zhfrk.f +++ b/SRC/zhfrk.f @@ -39,11 +39,11 @@ * Arguments * ========== * -* TRANSR (input) CHARACTER. +* TRANSR (input) CHARACTER * = 'N': The Normal Form of RFP A is stored; * = 'C': The Conjugate-transpose Form of RFP A is stored. * -* UPLO - (input) CHARACTER. +* UPLO (input) CHARACTER * On entry, UPLO specifies whether the upper or lower * triangular part of the array C is to be referenced as * follows: @@ -56,7 +56,7 @@ * * Unchanged on exit. * -* TRANS - (input) CHARACTER. +* TRANS (input) CHARACTER * On entry, TRANS specifies the operation to be performed as * follows: * @@ -66,41 +66,42 @@ * * Unchanged on exit. * -* N - (input) INTEGER. +* N (input) INTEGER * On entry, N specifies the order of the matrix C. N must be * at least zero. * Unchanged on exit. * -* K - (input) INTEGER. +* K (input) INTEGER * On entry with TRANS = 'N' or 'n', K specifies the number * of columns of the matrix A, and on entry with * TRANS = 'C' or 'c', K specifies the number of rows of the * matrix A. K must be at least zero. * Unchanged on exit. * -* ALPHA - (input) DOUBLE PRECISION. +* ALPHA (input) DOUBLE PRECISION * On entry, ALPHA specifies the scalar alpha. * Unchanged on exit. * -* A - (input) COMPLEX*16 array of DIMENSION ( LDA, ka ), where KA +* A (input) COMPLEX*16 array of DIMENSION (LDA,ka) +* where KA * is K when TRANS = 'N' or 'n', and is N otherwise. Before * entry with TRANS = 'N' or 'n', the leading N--by--K part of * the array A must contain the matrix A, otherwise the leading * K--by--N part of the array A must contain the matrix A. * Unchanged on exit. * -* LDA - (input) INTEGER. +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. When TRANS = 'N' or 'n' * then LDA must be at least max( 1, n ), otherwise LDA must * be at least max( 1, k ). * Unchanged on exit. * -* BETA - (input) DOUBLE PRECISION. +* BETA (input) DOUBLE PRECISION * On entry, BETA specifies the scalar beta. * Unchanged on exit. * -* C - (input/output) COMPLEX*16 array, dimension ( N*(N+1)/2 ). +* C (input/output) COMPLEX*16 array, dimension (N*(N+1)/2) * On entry, the matrix A in RFP Format. RFP Format is * described by TRANSR, UPLO and N. Note that the imaginary * parts of the diagonal elements need not be set, they are diff --git a/SRC/zla_gbamv.f b/SRC/zla_gbamv.f index f465ac25..e480f4c6 100644 --- a/SRC/zla_gbamv.f +++ b/SRC/zla_gbamv.f @@ -43,7 +43,7 @@ * Arguments * ========== * -* TRANS - INTEGER +* TRANS (input) INTEGER * On entry, TRANS specifies the operation to be performed as * follows: * @@ -53,20 +53,20 @@ * * Unchanged on exit. * -* M - INTEGER +* M (input) INTEGER * On entry, M specifies the number of rows of the matrix A. * M must be at least zero. * Unchanged on exit. * -* N - INTEGER +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. * -* KL - INTEGER +* KL (input) INTEGER * The number of subdiagonals within the band of A. KL >= 0. * -* KU - INTEGER +* KU (input) INTEGER * The number of superdiagonals within the band of A. KU >= 0. * * ALPHA - DOUBLE PRECISION @@ -78,13 +78,13 @@ * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, m ). * Unchanged on exit. * -* X - DOUBLE PRECISION array of DIMENSION at least +* X (input) DOUBLE PRECISION array, dimension * ( 1 + ( n - 1 )*abs( INCX ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( m - 1 )*abs( INCX ) ) otherwise. @@ -92,7 +92,7 @@ * vector x. * Unchanged on exit. * -* INCX - INTEGER +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. @@ -102,7 +102,7 @@ * supplied as zero then Y need not be set on input. * Unchanged on exit. * -* Y - DOUBLE PRECISION array of DIMENSION at least +* Y (input/output) DOUBLE PRECISION array, dimension * ( 1 + ( m - 1 )*abs( INCY ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( n - 1 )*abs( INCY ) ) otherwise. @@ -110,7 +110,7 @@ * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/zla_gbrpvgrw.f b/SRC/zla_gbrpvgrw.f index d7cd007b..e2342bf7 100644 --- a/SRC/zla_gbrpvgrw.f +++ b/SRC/zla_gbrpvgrw.f @@ -98,4 +98,4 @@ END IF END DO ZLA_GBRPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/zla_geamv.f b/SRC/zla_geamv.f index 953d5c68..2d3d056e 100644 --- a/SRC/zla_geamv.f +++ b/SRC/zla_geamv.f @@ -44,7 +44,7 @@ * Arguments * ========== * -* TRANS - INTEGER +* TRANS (input) INTEGER * On entry, TRANS specifies the operation to be performed as * follows: * @@ -54,12 +54,12 @@ * * Unchanged on exit. * -* M - INTEGER +* M (input) INTEGER * On entry, M specifies the number of rows of the matrix A. * M must be at least zero. * Unchanged on exit. * -* N - INTEGER +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. @@ -73,7 +73,7 @@ * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, m ). @@ -87,7 +87,7 @@ * vector x. * Unchanged on exit. * -* INCX - INTEGER +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. @@ -97,7 +97,7 @@ * supplied as zero then Y need not be set on input. * Unchanged on exit. * -* Y - DOUBLE PRECISION array of DIMENSION at least +* Y (input/output) DOUBLE PRECISION array, dimension * ( 1 + ( m - 1 )*abs( INCY ) ) when TRANS = 'N' or 'n' * and at least * ( 1 + ( n - 1 )*abs( INCY ) ) otherwise. @@ -105,7 +105,7 @@ * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/zla_heamv.f b/SRC/zla_heamv.f index 9aeeb8cc..a9f96b6c 100644 --- a/SRC/zla_heamv.f +++ b/SRC/zla_heamv.f @@ -42,7 +42,7 @@ * Arguments * ========== * -* UPLO - INTEGER +* UPLO (input) INTEGER * On entry, UPLO specifies whether the upper or lower * triangular part of the array A is to be referenced as * follows: @@ -55,7 +55,7 @@ * * Unchanged on exit. * -* N - INTEGER. +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. @@ -69,7 +69,7 @@ * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER. +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, n ). @@ -81,7 +81,7 @@ * vector x. * Unchanged on exit. * -* INCX - INTEGER. +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. @@ -91,13 +91,13 @@ * supplied as zero then Y need not be set on input. * Unchanged on exit. * -* Y - DOUBLE PRECISION array of DIMENSION at least +* Y (input/output) DOUBLE PRECISION array, dimension * ( 1 + ( n - 1 )*abs( INCY ) ) * Before entry with BETA non-zero, the incremented array Y * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER. +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/zla_herpvgrw.f b/SRC/zla_herpvgrw.f index 8bd30b84..b76f8d59 100644 --- a/SRC/zla_herpvgrw.f +++ b/SRC/zla_herpvgrw.f @@ -257,4 +257,4 @@ END IF ZLA_HERPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/zla_lin_berr.f b/SRC/zla_lin_berr.f index eda99e6e..35288618 100644 --- a/SRC/zla_lin_berr.f +++ b/SRC/zla_lin_berr.f @@ -48,7 +48,7 @@ * the matrix abs(op(A_s))*abs(Y) + abs(B_s). The matrices A, Y, and B * are from iterative refinement (see zla_gerfsx_extended.f). * -* RES (output) COMPLEX*16 array, dimension (NRHS) +* BERR (output) COMPLEX*16 array, dimension (NRHS) * The componentwise relative backward error from the formula above. * * ===================================================================== @@ -94,4 +94,4 @@ * END DO END DO - END SUBROUTINE + END diff --git a/SRC/zla_porpvgrw.f b/SRC/zla_porpvgrw.f index 6c7b3647..99d6a3e8 100644 --- a/SRC/zla_porpvgrw.f +++ b/SRC/zla_porpvgrw.f @@ -149,4 +149,4 @@ END IF ZLA_PORPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/zla_rpvgrw.f b/SRC/zla_rpvgrw.f index fd9c717e..6590d3ec 100644 --- a/SRC/zla_rpvgrw.f +++ b/SRC/zla_rpvgrw.f @@ -84,4 +84,4 @@ END IF END DO ZLA_RPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/zla_syamv.f b/SRC/zla_syamv.f index 53c8b698..0d18add4 100644 --- a/SRC/zla_syamv.f +++ b/SRC/zla_syamv.f @@ -43,7 +43,7 @@ * Arguments * ========== * -* UPLO - INTEGER +* UPLO (input) INTEGER * On entry, UPLO specifies whether the upper or lower * triangular part of the array A is to be referenced as * follows: @@ -56,7 +56,7 @@ * * Unchanged on exit. * -* N - INTEGER. +* N (input) INTEGER * On entry, N specifies the number of columns of the matrix A. * N must be at least zero. * Unchanged on exit. @@ -70,7 +70,7 @@ * contain the matrix of coefficients. * Unchanged on exit. * -* LDA - INTEGER. +* LDA (input) INTEGER * On entry, LDA specifies the first dimension of A as declared * in the calling (sub) program. LDA must be at least * max( 1, n ). @@ -82,7 +82,7 @@ * vector x. * Unchanged on exit. * -* INCX - INTEGER. +* INCX (input) INTEGER * On entry, INCX specifies the increment for the elements of * X. INCX must not be zero. * Unchanged on exit. @@ -92,13 +92,13 @@ * supplied as zero then Y need not be set on input. * Unchanged on exit. * -* Y - DOUBLE PRECISION array of DIMENSION at least +* Y (input/output) DOUBLE PRECISION array, dimension * ( 1 + ( n - 1 )*abs( INCY ) ) * Before entry with BETA non-zero, the incremented array Y * must contain the vector y. On exit, Y is overwritten by the * updated vector y. * -* INCY - INTEGER. +* INCY (input) INTEGER * On entry, INCY specifies the increment for the elements of * Y. INCY must not be zero. * Unchanged on exit. diff --git a/SRC/zla_syrpvgrw.f b/SRC/zla_syrpvgrw.f index c0e02134..d155696b 100644 --- a/SRC/zla_syrpvgrw.f +++ b/SRC/zla_syrpvgrw.f @@ -258,4 +258,4 @@ END IF ZLA_SYRPVGRW = RPVGRW - END FUNCTION + END diff --git a/SRC/zla_wwaddw.f b/SRC/zla_wwaddw.f index 244afc07..6cd28e7c 100644 --- a/SRC/zla_wwaddw.f +++ b/SRC/zla_wwaddw.f @@ -31,10 +31,13 @@ * N (input) INTEGER * The length of vectors X, Y, and W. * -* X, Y (input/output) COMPLEX*16 array, length N -* The doubled-single accumulation vector. +* X (input/output) COMPLEX*16 array, dimension (N) +* The first part of the doubled-single accumulation vector. * -* W (input) COMPLEX*16 array, length N +* Y (input/output) COMPLEX*16 array, dimension (N) +* The second part of the doubled-single accumulation vector. +* +* W (input) COMPLEX*16 array, dimension (N) * The vector to be added. * * ===================================================================== diff --git a/SRC/ztfsm.f b/SRC/ztfsm.f index 9e5d3efc..295880b4 100644 --- a/SRC/ztfsm.f +++ b/SRC/ztfsm.f @@ -40,11 +40,11 @@ * Arguments * ========== * -* TRANSR - (input) CHARACTER +* TRANSR (input) CHARACTER * = 'N': The Normal Form of RFP A is stored; * = 'C': The Conjugate-transpose Form of RFP A is stored. * -* SIDE - (input) CHARACTER +* SIDE (input) CHARACTER * On entry, SIDE specifies whether op( A ) appears on the left * or right of X as follows: * @@ -54,7 +54,7 @@ * * Unchanged on exit. * -* UPLO - (input) CHARACTER +* UPLO (input) CHARACTER * On entry, UPLO specifies whether the RFP matrix A came from * an upper or lower triangular matrix as follows: * UPLO = 'U' or 'u' RFP A came from an upper triangular matrix @@ -62,7 +62,7 @@ * * Unchanged on exit. * -* TRANS - (input) CHARACTER +* TRANS (input) CHARACTER * On entry, TRANS specifies the form of op( A ) to be used * in the matrix multiplication as follows: * @@ -72,7 +72,7 @@ * * Unchanged on exit. * -* DIAG - (input) CHARACTER +* DIAG (input) CHARACTER * On entry, DIAG specifies whether or not RFP A is unit * triangular as follows: * @@ -83,23 +83,23 @@ * * Unchanged on exit. * -* M - (input) INTEGER. +* M (input) INTEGER * On entry, M specifies the number of rows of B. M must be at * least zero. * Unchanged on exit. * -* N - (input) INTEGER. +* N (input) INTEGER * On entry, N specifies the number of columns of B. N must be * at least zero. * Unchanged on exit. * -* ALPHA - (input) COMPLEX*16. +* ALPHA (input) COMPLEX*16 * On entry, ALPHA specifies the scalar alpha. When alpha is * zero then A is not referenced and B need not be set before * entry. * Unchanged on exit. * -* A - (input) COMPLEX*16 array, dimension ( N*(N+1)/2 ); +* A (input) COMPLEX*16 array, dimension (N*(N+1)/2) * NT = N*(N+1)/2. On entry, the matrix A in RFP Format. * RFP Format is described by TRANSR, UPLO and N as follows: * If TRANSR='N' then RFP A is (0:N,0:K-1) when N is even; @@ -115,12 +115,12 @@ * even and is N when is odd. * See the Note below for more details. Unchanged on exit. * -* B - (input/ouptut) COMPLEX*16 array, DIMENSION ( LDB, N) +* B (input/output) COMPLEX*16 array, dimension (LDB,N) * Before entry, the leading m by n part of the array B must * contain the right-hand side matrix B, and on exit is * overwritten by the solution matrix X. * -* LDB - (input) INTEGER. +* LDB (input) INTEGER * On entry, LDB specifies the first dimension of B as declared * in the calling (sub) program. LDB must be at least * max( 1, m ). |