1 files changed, 56 insertions, 24 deletions
diff --git a/SRC/dgehrd.f b/SRC/dgehrd.f
index b24b2e3f..653a5fb8 100644
--- a/SRC/dgehrd.f
+++ b/SRC/dgehrd.f
@@ -45,47 +45,57 @@
 *>          The order of the matrix A.  N >= 0.
 *> \endverbatim
 *>
-*
-*  Authors
-*  =======
-*
-*> \author Univ. of Tennessee 
-*> \author Univ. of California Berkeley 
-*> \author Univ. of Colorado Denver 
-*> \author NAG Ltd. 
-*
-*> \date November 2011
-*
-*> \ingroup doubleGEcomputational
-*
-*
-*  Further Details
-*  ===============
-*>\details \b Further \b Details
+*> \param[in] ILO
 *> \verbatim
-*          set to 1 and N respectively. See Further Details.
+*>          ILO is INTEGER
+*> \endverbatim
+*>
+*> \param[in] IHI
+*> \verbatim
+*>          IHI is INTEGER
+*> \endverbatim
+*> \verbatim
+*>          It is assumed that A is already upper triangular in rows
+*>          and columns 1:ILO-1 and IHI+1:N. ILO and IHI are normally
+*>          set by a previous call to DGEBAL; otherwise they should be
+*>          set to 1 and N respectively. See Further Details.
 *>          1 <= ILO <= IHI <= N, if N > 0; ILO=1 and IHI=0, if N=0.
+*> \endverbatim
 *>
-*>  A       (input/output) DOUBLE PRECISION array, dimension (LDA,N)
+*> \param[in,out] A
+*> \verbatim
+*>          A is DOUBLE PRECISION array, dimension (LDA,N)
 *>          On entry, the N-by-N general matrix to be reduced.
 *>          On exit, the upper triangle and the first subdiagonal of A
 *>          are overwritten with the upper Hessenberg matrix H, and the
 *>          elements below the first subdiagonal, with the array TAU,
 *>          represent the orthogonal matrix Q as a product of elementary
 *>          reflectors. See Further Details.
+*> \endverbatim
 *>
-*>  LDA     (input) INTEGER
+*> \param[in] LDA
+*> \verbatim
+*>          LDA is INTEGER
 *>          The leading dimension of the array A.  LDA >= max(1,N).
+*> \endverbatim
 *>
-*>  TAU     (output) DOUBLE PRECISION array, dimension (N-1)
+*> \param[out] TAU
+*> \verbatim
+*>          TAU is DOUBLE PRECISION array, dimension (N-1)
 *>          The scalar factors of the elementary reflectors (see Further
 *>          Details). Elements 1:ILO-1 and IHI:N-1 of TAU are set to
 *>          zero.
+*> \endverbatim
 *>
-*>  WORK    (workspace/output) DOUBLE PRECISION array, dimension (LWORK)
+*> \param[out] WORK
+*> \verbatim
+*>          WORK is DOUBLE PRECISION array, dimension (LWORK)
 *>          On exit, if INFO = 0, WORK(1) returns the optimal LWORK.
+*> \endverbatim
 *>
-*>  LWORK   (input) INTEGER
+*> \param[in] LWORK
+*> \verbatim
+*>          LWORK is INTEGER
 *>          The length of the array WORK.  LWORK >= max(1,N).
 *>          For optimum performance LWORK >= N*NB, where NB is the
 *>          optimal blocksize.
@@ -94,11 +104,33 @@
 *>          only calculates the optimal size of the WORK array, returns
 *>          this value as the first entry of the WORK array, and no error
 *>          message related to LWORK is issued by XERBLA.
+*> \endverbatim
 *>
-*>  INFO    (output) INTEGER
+*> \param[out] INFO
+*> \verbatim
+*>          INFO is INTEGER
 *>          = 0:  successful exit
 *>          < 0:  if INFO = -i, the i-th argument had an illegal value.
+*> \endverbatim
 *>
+*
+*  Authors
+*  =======
+*
+*> \author Univ. of Tennessee 
+*> \author Univ. of California Berkeley 
+*> \author Univ. of Colorado Denver 
+*> \author NAG Ltd. 
+*
+*> \date November 2011
+*
+*> \ingroup doubleGEcomputational
+*
+*
+*  Further Details
+*  ===============
+*>\details \b Further \b Details
+*> \verbatim
 *>
 *>  The matrix Q is represented as a product of (ihi-ilo) elementary
 *>  reflectors