@@ -236,13 +236,14 @@ SUBROUTINE CLAHRD( N, K, NB, A, LDA, TAU, T, LDT, Y, LDY )
236236*
237237* w := T**H *w
238238*
239- CALL CTRMV( ' Upper' , ' Conjugate transpose' , ' Non-unit' , I -1 ,
240- $ T, LDT, T( 1 , NB ), 1 )
239+ CALL CTRMV( ' Upper' , ' Conjugate transpose' , ' Non-unit' ,
240+ $ I -1 , T, LDT, T( 1 , NB ), 1 )
241241*
242242* b2 := b2 - V2*w
243243*
244- CALL CGEMV( ' No transpose' , N- K- I+1 , I-1 , - ONE, A( K+ I, 1 ),
245- $ LDA, T( 1 , NB ), 1 , ONE, A( K+ I, I ), 1 )
244+ CALL CGEMV( ' No transpose' , N- K- I+1 , I-1 , - ONE,
245+ $ A( K+ I, 1 ), LDA, T( 1 , NB ), 1 , ONE,
246+ $ A( K+ I, I ), 1 )
246247*
247248* b1 := b1 - V1*w
248249*
@@ -263,20 +264,20 @@ SUBROUTINE CLAHRD( N, K, NB, A, LDA, TAU, T, LDT, Y, LDY )
263264*
264265* Compute Y(1:n,i)
265266*
266- CALL CGEMV( ' No transpose' , N, N- K- I+1 , ONE, A( 1 , I+1 ), LDA,
267- $ A( K+ I, I ), 1 , ZERO, Y( 1 , I ), 1 )
267+ CALL CGEMV( ' No transpose' , N, N- K- I+1 , ONE, A( 1 , I+1 ),
268+ $ LDA, A( K+ I, I ), 1 , ZERO, Y( 1 , I ), 1 )
268269 CALL CGEMV( ' Conjugate transpose' , N- K- I+1 , I-1 , ONE,
269270 $ A( K+ I, 1 ), LDA, A( K+ I, I ), 1 , ZERO, T( 1 , I ),
270271 $ 1 )
271- CALL CGEMV( ' No transpose' , N, I-1 , - ONE, Y, LDY, T( 1 , I ), 1 ,
272- $ ONE, Y( 1 , I ), 1 )
272+ CALL CGEMV( ' No transpose' , N, I-1 , - ONE, Y, LDY,
273+ $ T( 1 , I ), 1 , ONE, Y( 1 , I ), 1 )
273274 CALL CSCAL( N, TAU( I ), Y( 1 , I ), 1 )
274275*
275276* Compute T(1:i,i)
276277*
277278 CALL CSCAL( I-1 , - TAU( I ), T( 1 , I ), 1 )
278- CALL CTRMV( ' Upper' , ' No transpose' , ' Non-unit' , I-1 , T, LDT,
279- $ T( 1 , I ), 1 )
279+ CALL CTRMV( ' Upper' , ' No transpose' , ' Non-unit' , I-1 , T,
280+ $ LDT, T( 1 , I ), 1 )
280281 T( I, I ) = TAU( I )
281282*
282283 10 CONTINUE
0 commit comments