@@ -236,13 +236,14 @@ SUBROUTINE CLAHRD( N, K, NB, A, LDA, TAU, T, LDT, Y, LDY )
236
236
*
237
237
* w := T**H *w
238
238
*
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 )
241
241
*
242
242
* b2 := b2 - V2*w
243
243
*
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 )
246
247
*
247
248
* b1 := b1 - V1*w
248
249
*
@@ -263,20 +264,20 @@ SUBROUTINE CLAHRD( N, K, NB, A, LDA, TAU, T, LDT, Y, LDY )
263
264
*
264
265
* Compute Y(1:n,i)
265
266
*
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 )
268
269
CALL CGEMV( ' Conjugate transpose' , N- K- I+1 , I-1 , ONE,
269
270
$ A( K+ I, 1 ), LDA, A( K+ I, I ), 1 , ZERO, T( 1 , I ),
270
271
$ 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 )
273
274
CALL CSCAL( N, TAU( I ), Y( 1 , I ), 1 )
274
275
*
275
276
* Compute T(1:i,i)
276
277
*
277
278
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 )
280
281
T( I, I ) = TAU( I )
281
282
*
282
283
10 CONTINUE
0 commit comments