Subroutine */
| 53 | static doublereal c_b2435 = .5; |
| 54 | |
| 55 | /* Subroutine */ int zgebak_(char *job, char *side, integer *n, integer *ilo, |
| 56 | integer *ihi, doublereal *scale, integer *m, doublecomplex *v, |
| 57 | integer *ldv, integer *info) |
| 58 | { |
| 59 | /* System generated locals */ |
| 60 | integer v_dim1, v_offset, i__1; |
| 61 | |
| 62 | /* Local variables */ |
| 63 | static integer i__, k; |
| 64 | static doublereal s; |
| 65 | static integer ii; |
| 66 | extern logical lsame_(char *, char *); |
| 67 | static logical leftv; |
| 68 | extern /* Subroutine */ int zswap_(integer *, doublecomplex *, integer *, |
| 69 | doublecomplex *, integer *), xerbla_(char *, integer *), |
| 70 | zdscal_(integer *, doublereal *, doublecomplex *, integer *); |
| 71 | static logical rightv; |
| 72 | |
| 73 | |
| 74 | /* |
| 75 | -- LAPACK routine (version 3.2) -- |
| 76 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 77 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 78 | November 2006 |
| 79 | |
| 80 | |
| 81 | Purpose |
| 82 | ======= |
| 83 | |
| 84 | ZGEBAK forms the right or left eigenvectors of a complex general |
| 85 | matrix by backward transformation on the computed eigenvectors of the |
| 86 | balanced matrix output by ZGEBAL. |
| 87 | |
| 88 | Arguments |
| 89 | ========= |
| 90 | |
| 91 | JOB (input) CHARACTER*1 |
| 92 | Specifies the type of backward transformation required: |
| 93 | = 'N', do nothing, return immediately; |
| 94 | = 'P', do backward transformation for permutation only; |
| 95 | = 'S', do backward transformation for scaling only; |
| 96 | = 'B', do backward transformations for both permutation and |
| 97 | scaling. |
| 98 | JOB must be the same as the argument JOB supplied to ZGEBAL. |
| 99 | |
| 100 | SIDE (input) CHARACTER*1 |
| 101 | = 'R': V contains right eigenvectors; |
| 102 | = 'L': V contains left eigenvectors. |
| 103 | |
| 104 | N (input) INTEGER |
| 105 | The number of rows of the matrix V. N >= 0. |
| 106 | |
| 107 | ILO (input) INTEGER |
| 108 | IHI (input) INTEGER |
| 109 | The integers ILO and IHI determined by ZGEBAL. |
| 110 | 1 <= ILO <= IHI <= N, if N > 0; ILO=1 and IHI=0, if N=0. |
| 111 | |
| 112 | SCALE (input) DOUBLE PRECISION array, dimension (N) |