Subroutine */
| 14105 | } /* dlahqr_ */ |
| 14106 | |
| 14107 | /* Subroutine */ int dlahr2_(integer *n, integer *k, integer *nb, doublereal * |
| 14108 | a, integer *lda, doublereal *tau, doublereal *t, integer *ldt, |
| 14109 | doublereal *y, integer *ldy) |
| 14110 | { |
| 14111 | /* System generated locals */ |
| 14112 | integer a_dim1, a_offset, t_dim1, t_offset, y_dim1, y_offset, i__1, i__2, |
| 14113 | i__3; |
| 14114 | doublereal d__1; |
| 14115 | |
| 14116 | /* Local variables */ |
| 14117 | static integer i__; |
| 14118 | static doublereal ei; |
| 14119 | extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, |
| 14120 | integer *), dgemm_(char *, char *, integer *, integer *, integer * |
| 14121 | , doublereal *, doublereal *, integer *, doublereal *, integer *, |
| 14122 | doublereal *, doublereal *, integer *), dgemv_( |
| 14123 | char *, integer *, integer *, doublereal *, doublereal *, integer |
| 14124 | *, doublereal *, integer *, doublereal *, doublereal *, integer *), dcopy_(integer *, doublereal *, integer *, doublereal *, |
| 14125 | integer *), dtrmm_(char *, char *, char *, char *, integer *, |
| 14126 | integer *, doublereal *, doublereal *, integer *, doublereal *, |
| 14127 | integer *), daxpy_(integer *, |
| 14128 | doublereal *, doublereal *, integer *, doublereal *, integer *), |
| 14129 | dtrmv_(char *, char *, char *, integer *, doublereal *, integer *, |
| 14130 | doublereal *, integer *), dlarfg_( |
| 14131 | integer *, doublereal *, doublereal *, integer *, doublereal *), |
| 14132 | dlacpy_(char *, integer *, integer *, doublereal *, integer *, |
| 14133 | doublereal *, integer *); |
| 14134 | |
| 14135 | |
| 14136 | /* |
| 14137 | -- LAPACK auxiliary routine (version 3.2.1) -- |
| 14138 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 14139 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 14140 | -- April 2009 -- |
| 14141 | |
| 14142 | |
| 14143 | Purpose |
| 14144 | ======= |
| 14145 | |
| 14146 | DLAHR2 reduces the first NB columns of A real general n-BY-(n-k+1) |
| 14147 | matrix A so that elements below the k-th subdiagonal are zero. The |
| 14148 | reduction is performed by an orthogonal similarity transformation |
| 14149 | Q' * A * Q. The routine returns the matrices V and T which determine |
| 14150 | Q as a block reflector I - V*T*V', and also the matrix Y = A * V * T. |
| 14151 | |
| 14152 | This is an auxiliary routine called by DGEHRD. |
| 14153 | |
| 14154 | Arguments |
| 14155 | ========= |
| 14156 | |
| 14157 | N (input) INTEGER |
| 14158 | The order of the matrix A. |
| 14159 | |
| 14160 | K (input) INTEGER |
| 14161 | The offset for the reduction. Elements below the k-th |
| 14162 | subdiagonal in the first NB columns are reduced to zero. |
| 14163 | K < N. |
| 14164 |