--- /dev/null
+! { 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