Subroutine */
| 10016 | } /* claed7_ */ |
| 10017 | |
| 10018 | /* Subroutine */ int claed8_(integer *k, integer *n, integer *qsiz, complex * |
| 10019 | q, integer *ldq, real *d__, real *rho, integer *cutpnt, real *z__, |
| 10020 | real *dlamda, complex *q2, integer *ldq2, real *w, integer *indxp, |
| 10021 | integer *indx, integer *indxq, integer *perm, integer *givptr, |
| 10022 | integer *givcol, real *givnum, integer *info) |
| 10023 | { |
| 10024 | /* System generated locals */ |
| 10025 | integer q_dim1, q_offset, q2_dim1, q2_offset, i__1; |
| 10026 | real r__1; |
| 10027 | |
| 10028 | /* Local variables */ |
| 10029 | static real c__; |
| 10030 | static integer i__, j; |
| 10031 | static real s, t; |
| 10032 | static integer k2, n1, n2, jp, n1p1; |
| 10033 | static real eps, tau, tol; |
| 10034 | static integer jlam, imax, jmax; |
| 10035 | extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *), |
| 10036 | ccopy_(integer *, complex *, integer *, complex *, integer *), |
| 10037 | csrot_(integer *, complex *, integer *, complex *, integer *, |
| 10038 | real *, real *), scopy_(integer *, real *, integer *, real *, |
| 10039 | integer *); |
| 10040 | extern doublereal slapy2_(real *, real *), slamch_(char *); |
| 10041 | extern /* Subroutine */ int clacpy_(char *, integer *, integer *, complex |
| 10042 | *, integer *, complex *, integer *), xerbla_(char *, |
| 10043 | integer *); |
| 10044 | extern integer isamax_(integer *, real *, integer *); |
| 10045 | extern /* Subroutine */ int slamrg_(integer *, integer *, real *, integer |
| 10046 | *, integer *, integer *); |
| 10047 | |
| 10048 | |
| 10049 | /* |
| 10050 | -- LAPACK routine (version 3.2.2) -- |
| 10051 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 10052 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 10053 | June 2010 |
| 10054 | |
| 10055 | |
| 10056 | Purpose |
| 10057 | ======= |
| 10058 | |
| 10059 | CLAED8 merges the two sets of eigenvalues together into a single |
| 10060 | sorted set. Then it tries to deflate the size of the problem. |
| 10061 | There are two ways in which deflation can occur: when two or more |
| 10062 | eigenvalues are close together or if there is a tiny element in the |
| 10063 | Z vector. For each such occurrence the order of the related secular |
| 10064 | equation problem is reduced by one. |
| 10065 | |
| 10066 | Arguments |
| 10067 | ========= |
| 10068 | |
| 10069 | K (output) INTEGER |
| 10070 | Contains the number of non-deflated eigenvalues. |
| 10071 | This is the order of the related secular equation. |
| 10072 | |
| 10073 | N (input) INTEGER |
| 10074 | The dimension of the symmetric tridiagonal matrix. N >= 0. |
| 10075 |
no test coverage detected