8 files changed, 47 insertions, 19 deletions

diff --git a/SRC/cgehrd.f b/SRC/cgehrd.f
index 055db08d..cbc50440 100644
--- a/SRC/cgehrd.f
+++ b/SRC/cgehrd.f
@@ -97,9 +97,9 @@
 *  modified element of the upper Hessenberg matrix H, and vi denotes an
 *  element of the vector defining H(i).
 *
-*  This file is a slight modification of LAPACK-3.0's CGEHRD
+*  This file is a slight modification of LAPACK-3.0's DGEHRD
 *  subroutine incorporating improvements proposed by Quintana-Orti and
-*  Van de Geijn (2005). 
+*  Van de Geijn (2006). (See DLAHR2.)
 *
 *  =====================================================================
 *
diff --git a/SRC/clahr2.f b/SRC/clahr2.f
index 3e97d3d5..6caaf832 100644
--- a/SRC/clahr2.f
+++ b/SRC/clahr2.f
@@ -100,11 +100,18 @@
 *  modified element of the upper Hessenberg matrix H, and vi denotes an
 *  element of the vector defining H(i).
 *
-*  This file is a slight modification of LAPACK-3.0's CLAHRD
+*  This subroutine is a slight modification of LAPACK-3.0's DLAHRD
 *  incorporating improvements proposed by Quintana-Orti and Van de
 *  Gejin. Note that the entries of A(1:K,2:NB) differ from those
-*  returned by the original LAPACK routine. This function is
-*  not backward compatible with LAPACK3.0.
+*  returned by the original LAPACK-3.0's DLAHRD routine. (This
+*  subroutine is not backward compatible with LAPACK-3.0's DLAHRD.)
+*
+*  References
+*  ==========
+*
+*  Gregorio Quintana-Orti and Robert van de Geijn, "Improving the
+*  performance of reduction to Hessenberg form," ACM Transactions on
+*  Mathematical Software, 32(2):180-194, June 2006.
 *
 *  =====================================================================
 *
diff --git a/SRC/dgehrd.f b/SRC/dgehrd.f
index 4ad82246..136f90b6 100644
--- a/SRC/dgehrd.f
+++ b/SRC/dgehrd.f
@@ -99,7 +99,7 @@
 *
 *  This file is a slight modification of LAPACK-3.0's DGEHRD
 *  subroutine incorporating improvements proposed by Quintana-Orti and
-*  Van de Geijn (2005). 
+*  Van de Geijn (2006). (See DLAHR2.)
 *
 *  =====================================================================
 *
diff --git a/SRC/dlahr2.f b/SRC/dlahr2.f
index d204394c..dd687076 100644
--- a/SRC/dlahr2.f
+++ b/SRC/dlahr2.f
@@ -100,11 +100,18 @@
 *  modified element of the upper Hessenberg matrix H, and vi denotes an
 *  element of the vector defining H(i).
 *
-*  This file is a slight modification of LAPACK-3.0's DLAHRD
+*  This subroutine is a slight modification of LAPACK-3.0's DLAHRD
 *  incorporating improvements proposed by Quintana-Orti and Van de
 *  Gejin. Note that the entries of A(1:K,2:NB) differ from those
-*  returned by the original LAPACK routine. This function is
-*  not backward compatible with LAPACK3.0.
+*  returned by the original LAPACK-3.0's DLAHRD routine. (This
+*  subroutine is not backward compatible with LAPACK-3.0's DLAHRD.)
+*
+*  References
+*  ==========
+*
+*  Gregorio Quintana-Orti and Robert van de Geijn, "Improving the
+*  performance of reduction to Hessenberg form," ACM Transactions on
+*  Mathematical Software, 32(2):180-194, June 2006.
 *
 *  =====================================================================
 *
diff --git a/SRC/sgehrd.f b/SRC/sgehrd.f
index fbe794b4..152bfb40 100644
--- a/SRC/sgehrd.f
+++ b/SRC/sgehrd.f
@@ -97,9 +97,9 @@
 *  modified element of the upper Hessenberg matrix H, and vi denotes an
 *  element of the vector defining H(i).
 *
-*  This file is a slight modification of LAPACK-3.0's SGEHRD
+*  This file is a slight modification of LAPACK-3.0's DGEHRD
 *  subroutine incorporating improvements proposed by Quintana-Orti and
-*  Van de Geijn (2005). 
+*  Van de Geijn (2006). (See DLAHR2.)
 *
 *  =====================================================================
 *
diff --git a/SRC/slahr2.f b/SRC/slahr2.f
index 0bf73d12..fa83c760 100644
--- a/SRC/slahr2.f
+++ b/SRC/slahr2.f
@@ -100,11 +100,18 @@
 *  modified element of the upper Hessenberg matrix H, and vi denotes an
 *  element of the vector defining H(i).
 *
-*  This file is a slight modification of LAPACK-3.0's SLAHRD
+*  This subroutine is a slight modification of LAPACK-3.0's DLAHRD
 *  incorporating improvements proposed by Quintana-Orti and Van de
 *  Gejin. Note that the entries of A(1:K,2:NB) differ from those
-*  returned by the original LAPACK routine. This function is
-*  not backward compatible with LAPACK3.0.
+*  returned by the original LAPACK-3.0's DLAHRD routine. (This
+*  subroutine is not backward compatible with LAPACK-3.0's DLAHRD.)
+*
+*  References
+*  ==========
+*
+*  Gregorio Quintana-Orti and Robert van de Geijn, "Improving the
+*  performance of reduction to Hessenberg form," ACM Transactions on
+*  Mathematical Software, 32(2):180-194, June 2006.
 *
 *  =====================================================================
 *
diff --git a/SRC/zgehrd.f b/SRC/zgehrd.f
index 045feaee..51f1b94e 100644
--- a/SRC/zgehrd.f
+++ b/SRC/zgehrd.f
@@ -97,9 +97,9 @@
 *  modified element of the upper Hessenberg matrix H, and vi denotes an
 *  element of the vector defining H(i).
 *
-*  This file is a slight modification of LAPACK-3.0's ZGEHRD
+*  This file is a slight modification of LAPACK-3.0's DGEHRD
 *  subroutine incorporating improvements proposed by Quintana-Orti and
-*  Van de Geijn (2005). 
+*  Van de Geijn (2006). (See DLAHR2.)
 *
 *  =====================================================================
 *
diff --git a/SRC/zlahr2.f b/SRC/zlahr2.f
index 8822ed69..ba56ed77 100644
--- a/SRC/zlahr2.f
+++ b/SRC/zlahr2.f
@@ -100,11 +100,18 @@
 *  modified element of the upper Hessenberg matrix H, and vi denotes an
 *  element of the vector defining H(i).
 *
-*  This file is a slight modification of LAPACK-3.0's ZLAHRD
+*  This subroutine is a slight modification of LAPACK-3.0's DLAHRD
 *  incorporating improvements proposed by Quintana-Orti and Van de
 *  Gejin. Note that the entries of A(1:K,2:NB) differ from those
-*  returned by the original LAPACK routine. This function is
-*  not backward compatible with LAPACK3.0.
+*  returned by the original LAPACK-3.0's DLAHRD routine. (This
+*  subroutine is not backward compatible with LAPACK-3.0's DLAHRD.)
+*
+*  References
+*  ==========
+*
+*  Gregorio Quintana-Orti and Robert van de Geijn, "Improving the
+*  performance of reduction to Hessenberg form," ACM Transactions on
+*  Mathematical Software, 32(2):180-194, June 2006.
 *
 *  =====================================================================
 *