Subroutine */
| 12614 | } /* dlaed9_ */ |
| 12615 | |
| 12616 | /* Subroutine */ int dlaeda_(integer *n, integer *tlvls, integer *curlvl, |
| 12617 | integer *curpbm, integer *prmptr, integer *perm, integer *givptr, |
| 12618 | integer *givcol, doublereal *givnum, doublereal *q, integer *qptr, |
| 12619 | doublereal *z__, doublereal *ztemp, integer *info) |
| 12620 | { |
| 12621 | /* System generated locals */ |
| 12622 | integer i__1, i__2, i__3; |
| 12623 | |
| 12624 | /* Local variables */ |
| 12625 | static integer i__, k, mid, ptr; |
| 12626 | extern /* Subroutine */ int drot_(integer *, doublereal *, integer *, |
| 12627 | doublereal *, integer *, doublereal *, doublereal *); |
| 12628 | static integer curr, bsiz1, bsiz2, psiz1, psiz2, zptr1; |
| 12629 | extern /* Subroutine */ int dgemv_(char *, integer *, integer *, |
| 12630 | doublereal *, doublereal *, integer *, doublereal *, integer *, |
| 12631 | doublereal *, doublereal *, integer *), dcopy_(integer *, |
| 12632 | doublereal *, integer *, doublereal *, integer *), xerbla_(char *, |
| 12633 | integer *); |
| 12634 | |
| 12635 | |
| 12636 | /* |
| 12637 | -- LAPACK routine (version 3.2.2) -- |
| 12638 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 12639 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 12640 | June 2010 |
| 12641 | |
| 12642 | |
| 12643 | Purpose |
| 12644 | ======= |
| 12645 | |
| 12646 | DLAEDA computes the Z vector corresponding to the merge step in the |
| 12647 | CURLVLth step of the merge process with TLVLS steps for the CURPBMth |
| 12648 | problem. |
| 12649 | |
| 12650 | Arguments |
| 12651 | ========= |
| 12652 | |
| 12653 | N (input) INTEGER |
| 12654 | The dimension of the symmetric tridiagonal matrix. N >= 0. |
| 12655 | |
| 12656 | TLVLS (input) INTEGER |
| 12657 | The total number of merging levels in the overall divide and |
| 12658 | conquer tree. |
| 12659 | |
| 12660 | CURLVL (input) INTEGER |
| 12661 | The current level in the overall merge routine, |
| 12662 | 0 <= curlvl <= tlvls. |
| 12663 | |
| 12664 | CURPBM (input) INTEGER |
| 12665 | The current problem in the current level in the overall |
| 12666 | merge routine (counting from upper left to lower right). |
| 12667 | |
| 12668 | PRMPTR (input) INTEGER array, dimension (N lg N) |
| 12669 | Contains a list of pointers which indicate where in PERM a |
| 12670 | level's permutation is stored. PRMPTR(i+1) - PRMPTR(i) |
| 12671 | indicates the size of the permutation and incidentally the |
| 12672 | size of the full, non-deflated problem. |
| 12673 |