! { dg-do compile } ! { dg-options "-O3 -ffast-math" } SUBROUTINE ZUNG2L( M, N, K, A, LDA, TAU, WORK, INFO ) COMPLEX*16 A( LDA, * ), TAU( * ), WORK( * ) IF( M.LT.0 ) THEN END IF CALL ZLARF( 'LEFT', M-N+II, II-1, A( 1, II ), 1, TAU( I ), A, $ LDA, WORK ) CALL ZSCAL( M-N+II-1, -TAU( I ), A( 1, II ), 1 ) A( M-N+II, II ) = ONE - TAU( I ) END SUBROUTINE ZLARF( SIDE, M, N, V, INCV, TAU, C, LDC, WORK ) CHARACTER SIDE(*) LOGICAL LSAME COMPLEX*16 C( LDC, * ), V(*), WORK(*), TAU IF( LSAME( SIDE, 'L' ) ) THEN IF( TAU.NE.ZERO ) THEN CALL ZGEMV( 'CONJUGATE TRANSPOSE', M, N, ONE, C, LDC, V, $ INCV, ZERO, WORK, 1 ) END IF END IF END LOGICAL FUNCTION LSAME( CA, CB ) CHARACTER CA(*), CB(*) END SUBROUTINE ZGEMV ( TRANS, M, N, ALPHA, A, LDA, X, INCX, $ BETA, Y, INCY ) COMPLEX*16 A( LDA, * ), X( * ), Y( * ) CHARACTER TRANS(*) LOGICAL LSAME IF( LSAME( TRANS, 'N' ) )THEN IF( INCY.EQ.1 )THEN IF( X( JX ).NE.ZERO )THEN Y( I ) = Y( I ) + TEMP*A( I, J ) END IF END IF END IF END SUBROUTINE ZSCAL(N,ZA,ZX,INCX) COMPLEX*16 ZA,ZX(1) IF( N.LE.0 .OR. INCX.LE.0 )RETURN 20 DO 30 I = 1,N ZX(I) = ZA*ZX(I) 30 CONTINUE END