version 1.7, 2010/12/21 13:53:29
|
version 1.8, 2011/07/22 07:38:06
|
Line 1
|
Line 1
|
SUBROUTINE DLAIC1( JOB, J, X, SEST, W, GAMMA, SESTPR, S, C ) |
SUBROUTINE DLAIC1( JOB, J, X, SEST, W, GAMMA, SESTPR, S, C ) |
* |
* |
* -- 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 .. |
INTEGER J, JOB |
INTEGER J, JOB |
Line 27
|
Line 27
|
* [ s*x ] |
* [ s*x ] |
* xhat = [ c ] |
* xhat = [ c ] |
* is an approximate singular vector of |
* is an approximate singular vector of |
* [ L 0 ] |
* [ L 0 ] |
* Lhat = [ w' gamma ] |
* Lhat = [ w**T gamma ] |
* in the sense that |
* in the sense that |
* twonorm(Lhat*xhat) = sestpr. |
* twonorm(Lhat*xhat) = sestpr. |
* |
* |
* Depending on JOB, an estimate for the largest or smallest singular |
* Depending on JOB, an estimate for the largest or smallest singular |
* value is computed. |
* value is computed. |
* |
* |
* Note that [s c]' and sestpr**2 is an eigenpair of the system |
* Note that [s c]**T and sestpr**2 is an eigenpair of the system |
* |
* |
* diag(sest*sest, 0) + [alpha gamma] * [ alpha ] |
* diag(sest*sest, 0) + [alpha gamma] * [ alpha ] |
* [ gamma ] |
* [ gamma ] |
* |
* |
* where alpha = x'*w. |
* where alpha = x**T*w. |
* |
* |
* Arguments |
* Arguments |
* ========= |
* ========= |