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

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

version 1.7, 2010/12/21 13:53:52 version 1.8, 2011/07/22 07:38:18
Line 1 Line 1
       SUBROUTINE ZLATZM( SIDE, M, N, V, INCV, TAU, C1, C2, LDC, WORK )        SUBROUTINE ZLATZM( SIDE, M, N, V, INCV, TAU, C1, C2, LDC, WORK )
 *  *
 *  -- 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        CHARACTER          SIDE
Line 21 Line 21
 *  *
 *  ZLATZM applies a Householder matrix generated by ZTZRQF to a matrix.  *  ZLATZM applies a Householder matrix generated by ZTZRQF to a matrix.
 *  *
 *  Let P = I - tau*u*u',   u = ( 1 ),  *  Let P = I - tau*u*u**H,   u = ( 1 ),
 *                              ( v )  *                                ( v )
 *  where v is an (m-1) vector if SIDE = 'L', or a (n-1) vector if  *  where v is an (m-1) vector if SIDE = 'L', or a (n-1) vector if
 *  SIDE = 'R'.  *  SIDE = 'R'.
 *  *
Line 112 Line 112
 *  *
       IF( LSAME( SIDE, 'L' ) ) THEN        IF( LSAME( SIDE, 'L' ) ) THEN
 *  *
 *        w :=  conjg( C1 + v' * C2 )  *        w :=  ( C1 + v**H * C2 )**H
 *  *
          CALL ZCOPY( N, C1, LDC, WORK, 1 )           CALL ZCOPY( N, C1, LDC, WORK, 1 )
          CALL ZLACGV( N, WORK, 1 )           CALL ZLACGV( N, WORK, 1 )
          CALL ZGEMV( 'Conjugate transpose', M-1, N, ONE, C2, LDC, V,           CALL ZGEMV( 'Conjugate transpose', M-1, N, ONE, C2, LDC, V,
      $               INCV, ONE, WORK, 1 )       $               INCV, ONE, WORK, 1 )
 *  *
 *        [ C1 ] := [ C1 ] - tau* [ 1 ] * w'  *        [ C1 ] := [ C1 ] - tau* [ 1 ] * w**H
 *        [ C2 ]    [ C2 ]        [ v ]  *        [ C2 ]    [ C2 ]        [ v ]
 *  *
          CALL ZLACGV( N, WORK, 1 )           CALL ZLACGV( N, WORK, 1 )
Line 134 Line 134
          CALL ZGEMV( 'No transpose', M, N-1, ONE, C2, LDC, V, INCV, ONE,           CALL ZGEMV( 'No transpose', M, N-1, ONE, C2, LDC, V, INCV, ONE,
      $               WORK, 1 )       $               WORK, 1 )
 *  *
 *        [ C1, C2 ] := [ C1, C2 ] - tau* w * [ 1 , v']  *        [ C1, C2 ] := [ C1, C2 ] - tau* w * [ 1 , v**H]
 *  *
          CALL ZAXPY( M, -TAU, WORK, 1, C1, 1 )           CALL ZAXPY( M, -TAU, WORK, 1, C1, 1 )
          CALL ZGERC( M, N-1, -TAU, WORK, 1, V, INCV, C2, LDC )           CALL ZGERC( M, N-1, -TAU, WORK, 1, V, INCV, C2, LDC )
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>