[BACK]Return to zggqrf.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/zggqrf.f between versions 1.7 and 1.8

version 1.7, 2010/12/21 13:53:45 version 1.8, 2011/07/22 07:38:14
Line 1 Line 1
       SUBROUTINE ZGGQRF( N, M, P, A, LDA, TAUA, B, LDB, TAUB, WORK,        SUBROUTINE ZGGQRF( N, M, P, A, LDA, TAUA, B, LDB, TAUB, WORK,
      $                   LWORK, INFO )       $                   LWORK, INFO )
 *  *
 *  -- LAPACK routine (version 3.2) --  *  -- LAPACK 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            INFO, LDA, LDB, LWORK, M, N, P        INTEGER            INFO, LDA, LDB, LWORK, M, N, P
Line 40 Line 40
 *  In particular, if B is square and nonsingular, the GQR factorization  *  In particular, if B is square and nonsingular, the GQR factorization
 *  of A and B implicitly gives the QR factorization of inv(B)*A:  *  of A and B implicitly gives the QR factorization of inv(B)*A:
 *  *
 *               inv(B)*A = Z'*(inv(T)*R)  *               inv(B)*A = Z**H * (inv(T)*R)
 *  *
 *  where inv(B) denotes the inverse of the matrix B, and Z' denotes the  *  where inv(B) denotes the inverse of the matrix B, and Z**H denotes the
 *  conjugate transpose of matrix Z.  *  conjugate transpose of matrix Z.
 *  *
 *  Arguments  *  Arguments
Line 119 Line 119
 *  *
 *  Each H(i) has the form  *  Each H(i) has the form
 *  *
 *     H(i) = I - taua * v * v'  *     H(i) = I - taua * v * v**H
 *  *
 *  where taua is a complex scalar, and v is a complex vector with  *  where taua is a complex scalar, and v is a complex vector with
 *  v(1:i-1) = 0 and v(i) = 1; v(i+1:n) is stored on exit in A(i+1:n,i),  *  v(1:i-1) = 0 and v(i) = 1; v(i+1:n) is stored on exit in A(i+1:n,i),
Line 133 Line 133
 *  *
 *  Each H(i) has the form  *  Each H(i) has the form
 *  *
 *     H(i) = I - taub * v * v'  *     H(i) = I - taub * v * v**H
 *  *
 *  where taub is a complex scalar, and v is a complex vector with  *  where taub is a complex scalar, and v is a complex vector with
 *  v(p-k+i+1:p) = 0 and v(p-k+i) = 1; v(1:p-k+i-1) is stored on exit in  *  v(p-k+i+1:p) = 0 and v(p-k+i) = 1; v(1:p-k+i-1) is stored on exit in
Line 194 Line 194
       CALL ZGEQRF( N, M, A, LDA, TAUA, WORK, LWORK, INFO )        CALL ZGEQRF( N, M, A, LDA, TAUA, WORK, LWORK, INFO )
       LOPT = WORK( 1 )        LOPT = WORK( 1 )
 *  *
 *     Update B := Q'*B.  *     Update B := Q**H*B.
 *  *
       CALL ZUNMQR( 'Left', 'Conjugate Transpose', N, P, MIN( N, M ), A,        CALL ZUNMQR( 'Left', 'Conjugate Transpose', N, P, MIN( N, M ), A,
      $             LDA, TAUA, B, LDB, WORK, LWORK, INFO )       $             LDA, TAUA, B, LDB, WORK, LWORK, INFO )
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>