version 1.7, 2010/12/21 13:53:31
|
version 1.8, 2011/07/22 07:38:07
|
Line 1
|
Line 1
|
SUBROUTINE DLAQTR( LTRAN, LREAL, N, T, LDT, B, W, SCALE, X, WORK, |
SUBROUTINE DLAQTR( LTRAN, LREAL, N, T, LDT, B, W, SCALE, X, WORK, |
$ INFO ) |
$ INFO ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.2) -- |
* -- LAPACK auxiliary routine (version 3.3.1) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* November 2006 |
* -- April 2011 -- |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
LOGICAL LREAL, LTRAN |
LOGICAL LREAL, LTRAN |
Line 36
|
Line 36
|
* [ . ] |
* [ . ] |
* [ w ] |
* [ w ] |
* |
* |
* op(A) = A or A', A' denotes the conjugate transpose of |
* op(A) = A or A**T, A**T denotes the transpose of |
* matrix A. |
* matrix A. |
* |
* |
* On input, X = [ c ]. On output, X = [ p ]. |
* On input, X = [ c ]. On output, X = [ p ]. |
Line 51
|
Line 51
|
* LTRAN (input) LOGICAL |
* LTRAN (input) LOGICAL |
* On entry, LTRAN specifies the option of conjugate transpose: |
* On entry, LTRAN specifies the option of conjugate transpose: |
* = .FALSE., op(T+i*B) = T+i*B, |
* = .FALSE., op(T+i*B) = T+i*B, |
* = .TRUE., op(T+i*B) = (T+i*B)'. |
* = .TRUE., op(T+i*B) = (T+i*B)**T. |
* |
* |
* LREAL (input) LOGICAL |
* LREAL (input) LOGICAL |
* On entry, LREAL specifies the input matrix structure: |
* On entry, LREAL specifies the input matrix structure: |
Line 290
|
Line 290
|
* |
* |
ELSE |
ELSE |
* |
* |
* Solve T'*p = scale*c |
* Solve T**T*p = scale*c |
* |
* |
JNEXT = 1 |
JNEXT = 1 |
DO 40 J = 1, N |
DO 40 J = 1, N |
Line 532
|
Line 532
|
* |
* |
ELSE |
ELSE |
* |
* |
* Solve (T + iB)'*(p+iq) = c+id |
* Solve (T + iB)**T*(p+iq) = c+id |
* |
* |
JNEXT = 1 |
JNEXT = 1 |
DO 80 J = 1, N |
DO 80 J = 1, N |