Subroutine */
| 10074 | } /* zlaed7_ */ |
| 10075 | |
| 10076 | /* Subroutine */ int zlaed8_(integer *k, integer *n, integer *qsiz, |
| 10077 | doublecomplex *q, integer *ldq, doublereal *d__, doublereal *rho, |
| 10078 | integer *cutpnt, doublereal *z__, doublereal *dlamda, doublecomplex * |
| 10079 | q2, integer *ldq2, doublereal *w, integer *indxp, integer *indx, |
| 10080 | integer *indxq, integer *perm, integer *givptr, integer *givcol, |
| 10081 | doublereal *givnum, integer *info) |
| 10082 | { |
| 10083 | /* System generated locals */ |
| 10084 | integer q_dim1, q_offset, q2_dim1, q2_offset, i__1; |
| 10085 | doublereal d__1; |
| 10086 | |
| 10087 | /* Local variables */ |
| 10088 | static doublereal c__; |
| 10089 | static integer i__, j; |
| 10090 | static doublereal s, t; |
| 10091 | static integer k2, n1, n2, jp, n1p1; |
| 10092 | static doublereal eps, tau, tol; |
| 10093 | static integer jlam, imax, jmax; |
| 10094 | extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, |
| 10095 | integer *), dcopy_(integer *, doublereal *, integer *, doublereal |
| 10096 | *, integer *), zdrot_(integer *, doublecomplex *, integer *, |
| 10097 | doublecomplex *, integer *, doublereal *, doublereal *), zcopy_( |
| 10098 | integer *, doublecomplex *, integer *, doublecomplex *, integer *) |
| 10099 | ; |
| 10100 | |
| 10101 | extern integer idamax_(integer *, doublereal *, integer *); |
| 10102 | extern /* Subroutine */ int dlamrg_(integer *, integer *, doublereal *, |
| 10103 | integer *, integer *, integer *), xerbla_(char *, integer *), zlacpy_(char *, integer *, integer *, doublecomplex *, |
| 10104 | integer *, doublecomplex *, integer *); |
| 10105 | |
| 10106 | |
| 10107 | /* |
| 10108 | -- LAPACK routine (version 3.2.2) -- |
| 10109 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 10110 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 10111 | June 2010 |
| 10112 | |
| 10113 | |
| 10114 | Purpose |
| 10115 | ======= |
| 10116 | |
| 10117 | ZLAED8 merges the two sets of eigenvalues together into a single |
| 10118 | sorted set. Then it tries to deflate the size of the problem. |
| 10119 | There are two ways in which deflation can occur: when two or more |
| 10120 | eigenvalues are close together or if there is a tiny element in the |
| 10121 | Z vector. For each such occurrence the order of the related secular |
| 10122 | equation problem is reduced by one. |
| 10123 | |
| 10124 | Arguments |
| 10125 | ========= |
| 10126 | |
| 10127 | K (output) INTEGER |
| 10128 | Contains the number of non-deflated eigenvalues. |
| 10129 | This is the order of the related secular equation. |
| 10130 | |
| 10131 | N (input) INTEGER |
| 10132 | The dimension of the symmetric tridiagonal matrix. N >= 0. |
| 10133 |
no test coverage detected