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

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

version 1.7, 2010/12/21 13:53:59 version 1.8, 2011/07/22 07:38:22
Line 1 Line 1
       SUBROUTINE ZUNMRQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC,        SUBROUTINE ZUNMRQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC,
      $                   WORK, LWORK, INFO )       $                   WORK, 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 ..
       CHARACTER          SIDE, TRANS        CHARACTER          SIDE, TRANS
Line 26 Line 26
 *  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
 *  *
 *        Q = H(1)' H(2)' . . . H(k)'  *        Q = H(1)**H H(2)**H . . . H(k)**H
 *  *
 *  as returned by ZGERQF. Q is of order M if SIDE = 'L' and of order N  *  as returned by ZGERQF. Q is of order M if SIDE = 'L' and of order N
 *  if SIDE = 'R'.  *  if SIDE = 'R'.
Line 244 Line 244
      $                   A( I, 1 ), LDA, TAU( I ), T, LDT )       $                   A( I, 1 ), LDA, TAU( I ), T, LDT )
             IF( LEFT ) THEN              IF( LEFT ) THEN
 *  *
 *              H or H' is applied to C(1:m-k+i+ib-1,1:n)  *              H or H**H is applied to C(1:m-k+i+ib-1,1:n)
 *  *
                MI = M - K + I + IB - 1                 MI = M - K + I + IB - 1
             ELSE              ELSE
 *  *
 *              H or H' is applied to C(1:m,1:n-k+i+ib-1)  *              H or H**H is applied to C(1:m,1:n-k+i+ib-1)
 *  *
                NI = N - K + I + IB - 1                 NI = N - K + I + IB - 1
             END IF              END IF
 *  *
 *           Apply H or H'  *           Apply H or H**H
 *  *
             CALL ZLARFB( SIDE, TRANST, 'Backward', 'Rowwise', MI, NI,              CALL ZLARFB( SIDE, TRANST, 'Backward', 'Rowwise', MI, NI,
      $                   IB, A( I, 1 ), LDA, T, LDT, C, LDC, WORK,       $                   IB, A( I, 1 ), LDA, T, LDT, C, LDC, WORK,
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>