@@ -324,8 +324,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV,
324324*
325325* T_{1,2} = T_{1,2}*V_{2,2}
326326*
327- CALL DTRMM(' Right' , ' Lower' , ' No transpose' , ' Unit' , L, K - L,
328- $ ONE, V(L+1 , L+1 ), LDV, T(1 , L+1 ), LDT)
327+ CALL DTRMM(' Right' , ' Lower' , ' No transpose' , ' Unit' , L,
328+ $ K - L, ONE, V(L+1 , L+1 ), LDV, T(1 , L+1 ), LDT)
329329
330330*
331331* T_{1,2} = V_{3,1}'*V_{3,2} + T_{1,2}
@@ -415,8 +415,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV,
415415*
416416* T_{1,2} = T_{1,2}*V_{2,2}'
417417*
418- CALL DTRMM(' Right' , ' Upper' , ' Transpose' , ' Unit' , L, K- L, ONE,
419- $ V(L+1 , L+1 ), LDV, T(1 , L+1 ), LDT)
418+ CALL DTRMM(' Right' , ' Upper' , ' Transpose' , ' Unit' , L, K- L,
419+ $ ONE, V(L+1 , L+1 ), LDV, T(1 , L+1 ), LDT)
420420
421421*
422422* T_{1,2} = V_{1,3}*V_{2,3}' + T_{1,2}
@@ -431,8 +431,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV,
431431*
432432* T_{1,2} = -T_{1,1}*T_{1,2}
433433*
434- CALL DTRMM(' Left' , ' Upper' , ' No transpose' , ' Non-unit' , L, K - L,
435- $ NEG_ONE, T, LDT, T(1 , L+1 ), LDT)
434+ CALL DTRMM(' Left' , ' Upper' , ' No transpose' , ' Non-unit' , L,
435+ $ K - L, NEG_ONE, T, LDT, T(1 , L+1 ), LDT)
436436
437437*
438438* T_{1,2} = T_{1,2}*T_{2,2}
@@ -603,8 +603,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV,
603603*
604604* T_{2,1} = T_{2,1}*V_{2,1}
605605*
606- CALL DTRMM(' Right' , ' Upper' , ' No transpose' , ' Unit' , L, K - L,
607- $ ONE, V(N- K+1 ,1 ), LDV, T(K- L+1 ,1 ), LDT)
606+ CALL DTRMM(' Right' , ' Upper' , ' No transpose' , ' Unit' , L,
607+ $ K - L, ONE, V(N- K+1 ,1 ), LDV, T(K- L+1 ,1 ), LDT)
608608
609609*
610610* T_{2,1} = V_{2,2}'*V_{2,1} + T_{2,1}
@@ -619,8 +619,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV,
619619*
620620* T_{2,1} = -T_{2,2}*T_{2,1}
621621*
622- CALL DTRMM(' Left' , ' Lower' , ' No transpose' , ' Non-unit' , L, K - L,
623- $ NEG_ONE, T(K- L+1 ,K- L+1 ), LDT, T(K- L+1 ,1 ), LDT)
622+ CALL DTRMM(' Left' , ' Lower' , ' No transpose' , ' Non-unit' , L,
623+ $ K - L, NEG_ONE, T(K- L+1 ,K- L+1 ), LDT, T(K- L+1 ,1 ), LDT)
624624*
625625* T_{2,1} = T_{2,1}*T_{1,1}
626626*
@@ -694,8 +694,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV,
694694*
695695* T_{2,1} = T_{2,1}*V_{1,2}'
696696*
697- CALL DTRMM(' Right' , ' Lower' , ' Transpose' , ' Unit' , L, K- L, ONE,
698- $ V(1 , N- K+1 ), LDV, T(K- L+1 ,1 ), LDT)
697+ CALL DTRMM(' Right' , ' Lower' , ' Transpose' , ' Unit' , L, K- L,
698+ $ ONE, V(1 , N- K+1 ), LDV, T(K- L+1 ,1 ), LDT)
699699
700700*
701701* T_{2,1} = V_{2,1}*V_{1,1}' + T_{2,1}
@@ -711,14 +711,14 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV,
711711*
712712* T_{2,1} = -T_{2,2}*T_{2,1}
713713*
714- CALL DTRMM(' Left' , ' Lower' , ' No tranpose' , ' Non-unit' , L, K - L,
715- $ NEG_ONE, T(K- L+1 ,K- L+1 ), LDT, T(K- L+1 ,1 ), LDT)
714+ CALL DTRMM(' Left' , ' Lower' , ' No tranpose' , ' Non-unit' , L,
715+ $ K - L, NEG_ONE, T(K- L+1 ,K- L+1 ), LDT, T(K- L+1 ,1 ), LDT)
716716
717717*
718718* T_{2,1} = T_{2,1}*T_{1,1}
719719*
720- CALL DTRMM(' Right' , ' Lower' , ' No tranpose' , ' Non-unit' , L, K - L,
721- $ ONE, T, LDT, T(K- L+1 ,1 ), LDT)
720+ CALL DTRMM(' Right' , ' Lower' , ' No tranpose' , ' Non-unit' , L,
721+ $ K - L, ONE, T, LDT, T(K- L+1 ,1 ), LDT)
722722 ELSE IF (RQT) THEN
723723* Break V apart into 6 components
724724*
0 commit comments