version 1.5, 2010/08/07 13:22:22
|
version 1.8, 2011/07/22 07:38:08
|
Line 201
|
Line 201
|
END IF |
END IF |
ELSE |
ELSE |
* |
* |
* Form P', determined by a call to DGEBRD to reduce a k-by-n |
* Form P**T, determined by a call to DGEBRD to reduce a k-by-n |
* matrix |
* matrix |
* |
* |
IF( K.LT.N ) THEN |
IF( K.LT.N ) THEN |
Line 215
|
Line 215
|
* If k >= n, assume m = n |
* If k >= n, assume m = n |
* |
* |
* Shift the vectors which define the elementary reflectors one |
* Shift the vectors which define the elementary reflectors one |
* row downward, and set the first row and column of P' to |
* row downward, and set the first row and column of P**T to |
* those of the unit matrix |
* those of the unit matrix |
* |
* |
A( 1, 1 ) = ONE |
A( 1, 1 ) = ONE |
Line 230
|
Line 230
|
60 CONTINUE |
60 CONTINUE |
IF( N.GT.1 ) THEN |
IF( N.GT.1 ) THEN |
* |
* |
* Form P'(2:n,2:n) |
* Form P**T(2:n,2:n) |
* |
* |
CALL DORGLQ( N-1, N-1, N-1, A( 2, 2 ), LDA, TAU, WORK, |
CALL DORGLQ( N-1, N-1, N-1, A( 2, 2 ), LDA, TAU, WORK, |
$ LWORK, IINFO ) |
$ LWORK, IINFO ) |