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

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

version 1.7, 2010/12/21 13:53:58 version 1.8, 2011/07/22 07:38:22
Line 1 Line 1
       SUBROUTINE ZUNM2L( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC,        SUBROUTINE ZUNM2L( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC,
      $                   WORK, INFO )       $                   WORK, 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 ..
       CHARACTER          SIDE, TRANS        CHARACTER          SIDE, TRANS
Line 21 Line 21
 *  *
 *        Q * C  if SIDE = 'L' and TRANS = 'N', or  *        Q * C  if SIDE = 'L' and TRANS = 'N', or
 *  *
 *        Q'* C  if SIDE = 'L' and TRANS = 'C', or  *        Q**H* C  if SIDE = 'L' and TRANS = 'C', or
 *  *
 *        C * Q  if SIDE = 'R' and TRANS = 'N', or  *        C * Q  if SIDE = 'R' and TRANS = 'N', or
 *  *
 *        C * Q' if SIDE = 'R' and TRANS = 'C',  *        C * Q**H if SIDE = 'R' and TRANS = 'C',
 *  *
 *  where Q is a complex unitary matrix defined as the product of k  *  where Q is a complex unitary matrix defined as the product of k
 *  elementary reflectors  *  elementary reflectors
Line 39 Line 39
 *  =========  *  =========
 *  *
 *  SIDE    (input) CHARACTER*1  *  SIDE    (input) CHARACTER*1
 *          = 'L': apply Q or Q' from the Left  *          = 'L': apply Q or Q**H from the Left
 *          = 'R': apply Q or Q' from the Right  *          = 'R': apply Q or Q**H from the Right
 *  *
 *  TRANS   (input) CHARACTER*1  *  TRANS   (input) CHARACTER*1
 *          = 'N': apply Q  (No transpose)  *          = 'N': apply Q  (No transpose)
 *          = 'C': apply Q' (Conjugate transpose)  *          = 'C': apply Q**H (Conjugate transpose)
 *  *
 *  M       (input) INTEGER  *  M       (input) INTEGER
 *          The number of rows of the matrix C. M >= 0.  *          The number of rows of the matrix C. M >= 0.
Line 75 Line 75
 *  *
 *  C       (input/output) COMPLEX*16 array, dimension (LDC,N)  *  C       (input/output) COMPLEX*16 array, dimension (LDC,N)
 *          On entry, the m-by-n matrix C.  *          On entry, the m-by-n matrix C.
 *          On exit, C is overwritten by Q*C or Q'*C or C*Q' or C*Q.  *          On exit, C is overwritten by Q*C or Q**H*C or C*Q**H or C*Q.
 *  *
 *  LDC     (input) INTEGER  *  LDC     (input) INTEGER
 *          The leading dimension of the array C. LDC >= max(1,M).  *          The leading dimension of the array C. LDC >= max(1,M).
Line 168 Line 168
       DO 10 I = I1, I2, I3        DO 10 I = I1, I2, I3
          IF( LEFT ) THEN           IF( LEFT ) THEN
 *  *
 *           H(i) or H(i)' is applied to C(1:m-k+i,1:n)  *           H(i) or H(i)**H is applied to C(1:m-k+i,1:n)
 *  *
             MI = M - K + I              MI = M - K + I
          ELSE           ELSE
 *  *
 *           H(i) or H(i)' is applied to C(1:m,1:n-k+i)  *           H(i) or H(i)**H is applied to C(1:m,1:n-k+i)
 *  *
             NI = N - K + I              NI = N - K + I
          END IF           END IF
 *  *
 *        Apply H(i) or H(i)'  *        Apply H(i) or H(i)**H
 *  *
          IF( NOTRAN ) THEN           IF( NOTRAN ) THEN
             TAUI = TAU( I )              TAUI = TAU( I )
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>