Subroutine */
| 19467 | } /* zlarfg_ */ |
| 19468 | |
| 19469 | /* Subroutine */ int zlarft_(char *direct, char *storev, integer *n, integer * |
| 19470 | k, doublecomplex *v, integer *ldv, doublecomplex *tau, doublecomplex * |
| 19471 | t, integer *ldt) |
| 19472 | { |
| 19473 | /* System generated locals */ |
| 19474 | integer t_dim1, t_offset, v_dim1, v_offset, i__1, i__2, i__3, i__4; |
| 19475 | doublecomplex z__1; |
| 19476 | |
| 19477 | /* Local variables */ |
| 19478 | static integer i__, j, prevlastv; |
| 19479 | static doublecomplex vii; |
| 19480 | extern logical lsame_(char *, char *); |
| 19481 | extern /* Subroutine */ int zgemv_(char *, integer *, integer *, |
| 19482 | doublecomplex *, doublecomplex *, integer *, doublecomplex *, |
| 19483 | integer *, doublecomplex *, doublecomplex *, integer *); |
| 19484 | static integer lastv; |
| 19485 | extern /* Subroutine */ int ztrmv_(char *, char *, char *, integer *, |
| 19486 | doublecomplex *, integer *, doublecomplex *, integer *), zlacgv_(integer *, doublecomplex *, integer *); |
| 19487 | |
| 19488 | |
| 19489 | /* |
| 19490 | -- LAPACK auxiliary routine (version 3.2) -- |
| 19491 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 19492 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 19493 | November 2006 |
| 19494 | |
| 19495 | |
| 19496 | Purpose |
| 19497 | ======= |
| 19498 | |
| 19499 | ZLARFT forms the triangular factor T of a complex block reflector H |
| 19500 | of order n, which is defined as a product of k elementary reflectors. |
| 19501 | |
| 19502 | If DIRECT = 'F', H = H(1) H(2) . . . H(k) and T is upper triangular; |
| 19503 | |
| 19504 | If DIRECT = 'B', H = H(k) . . . H(2) H(1) and T is lower triangular. |
| 19505 | |
| 19506 | If STOREV = 'C', the vector which defines the elementary reflector |
| 19507 | H(i) is stored in the i-th column of the array V, and |
| 19508 | |
| 19509 | H = I - V * T * V' |
| 19510 | |
| 19511 | If STOREV = 'R', the vector which defines the elementary reflector |
| 19512 | H(i) is stored in the i-th row of the array V, and |
| 19513 | |
| 19514 | H = I - V' * T * V |
| 19515 | |
| 19516 | Arguments |
| 19517 | ========= |
| 19518 | |
| 19519 | DIRECT (input) CHARACTER*1 |
| 19520 | Specifies the order in which the elementary reflectors are |
| 19521 | multiplied to form the block reflector: |
| 19522 | = 'F': H = H(1) H(2) . . . H(k) (Forward) |
| 19523 | = 'B': H = H(k) . . . H(2) H(1) (Backward) |
| 19524 | |
| 19525 | STOREV (input) CHARACTER*1 |
| 19526 | Specifies how the vectors which define the elementary |
no test coverage detected