@@ -202,7 +202,7 @@ SUBROUTINE SLARF1F( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
202202* w(1:lastc,1) := C(2:lastv,1:lastc)**T * v(2:lastv,1)
203203*
204204 CALL SGEMV( ' Transpose' , LASTV - 1 , LASTC, ONE, C( 2 , 1 ),
205- $ LDC, V( 1 + INCV ), INCV, ZERO, WORK, 1 )
205+ $ LDC, V( 1 + INCV ), INCV, ZERO, WORK, 1 )
206206*
207207* w(1:lastc,1) += v(1,1) * C(1,1:lastc)**T
208208*
@@ -214,8 +214,8 @@ SUBROUTINE SLARF1F( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
214214*
215215* C(2:lastv,1:lastc) += - tau * v(2:lastv,1) * w(1:lastc,1)**T
216216*
217- CALL SGER( LASTV - 1 , LASTC, - TAU, V( 1 + INCV ), INCV, WORK,
218- $ 1 , C( 2 , 1 ), LDC )
217+ CALL SGER( LASTV - 1 , LASTC, - TAU, V( 1 + INCV ), INCV,
218+ $ WORK, 1 , C( 2 , 1 ), LDC )
219219 END IF
220220 ELSE
221221*
@@ -230,8 +230,9 @@ SUBROUTINE SLARF1F( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
230230*
231231* w(1:lastc,1) := C(1:lastc,2:lastv) * v(2:lastv,1)
232232*
233- CALL SGEMV( ' No transpose' , LASTC, LASTV - 1 , ONE,
234- $ C( 1 , 2 ), LDC, V( 1 + INCV ), INCV, ZERO, WORK, 1 )
233+ CALL SGEMV( ' No transpose' , LASTC, LASTV - 1 , ONE,
234+ $ C( 1 , 2 ), LDC, V( 1 + INCV ), INCV, ZERO,
235+ $ WORK, 1 )
235236*
236237* w(1:lastc,1) += v(1,1) * C(1:lastc,1)
237238*
@@ -244,7 +245,7 @@ SUBROUTINE SLARF1F( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
244245* C(1:lastc,2:lastv) += - tau * w(1:lastc,1) * v(2:lastv)**T
245246*
246247 CALL SGER( LASTC, LASTV - 1 , - TAU, WORK, 1 ,
247- $ V( 1 + INCV ), INCV, C( 1 , 2 ), LDC )
248+ $ V( 1 + INCV ), INCV, C( 1 , 2 ), LDC )
248249 END IF
249250 END IF
250251 RETURN
0 commit comments