Skip to content

Commit 354a16f

Browse files
jprhynejprhyne
committed
Removed mod files and extranous file changes (hopefully)
1 parent dadd80e commit 354a16f

File tree

14 files changed

+60
-59
lines changed

14 files changed

+60
-59
lines changed

‎SRC/clarft.f‎

Lines changed: 22 additions & 22 deletions
Original file line numberDiff line numberDiff line change
@@ -159,8 +159,8 @@
159159
*> \endverbatim
160160
*>
161161
* =====================================================================
162-
RECURSIVE SUBROUTINE CLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
163-
$ LDT )
162+
RECURSIVE SUBROUTINE CLARFT( DIRECT, STOREV, N, K, V, LDV,
163+
$ TAU, T, LDT )
164164
*
165165
* -- LAPACK auxiliary routine --
166166
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -319,8 +319,8 @@ RECURSIVE SUBROUTINE CLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
319319
*
320320
* T_{1,2} = T_{1,2}*V_{2,2}
321321
*
322-
CALL CTRMM('Right', 'Lower', 'No transpose', 'Unit', L, K-L,
323-
$ ONE, V(L+1, L+1), LDV, T(1, L+1), LDT)
322+
CALL CTRMM('Right', 'Lower', 'No transpose', 'Unit', L,
323+
$ K-L, ONE, V(L+1, L+1), LDV, T(1, L+1), LDT)
324324

325325
*
326326
* T_{1,2} = V_{3,1}'*V_{3,2} + T_{1,2}
@@ -410,8 +410,8 @@ RECURSIVE SUBROUTINE CLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
410410
*
411411
* T_{1,2} = T_{1,2}*V_{2,2}'
412412
*
413-
CALL CTRMM('Right', 'Upper', 'Conjugate', 'Unit', L, K-L, ONE,
414-
$ V(L+1, L+1), LDV, T(1, L+1), LDT)
413+
CALL CTRMM('Right', 'Upper', 'Conjugate', 'Unit', L, K-L,
414+
$ ONE, V(L+1, L+1), LDV, T(1, L+1), LDT)
415415

416416
*
417417
* T_{1,2} = V_{1,3}*V_{2,3}' + T_{1,2}
@@ -426,8 +426,8 @@ RECURSIVE SUBROUTINE CLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
426426
*
427427
* T_{1,2} = -T_{1,1}*T_{1,2}
428428
*
429-
CALL CTRMM('Left', 'Upper', 'No transpose', 'Non-unit', L, K-L,
430-
$ NEG_ONE, T, LDT, T(1, L+1), LDT)
429+
CALL CTRMM('Left', 'Upper', 'No transpose', 'Non-unit', L,
430+
$ K-L, NEG_ONE, T, LDT, T(1, L+1), LDT)
431431

432432
*
433433
* T_{1,2} = T_{1,2}*T_{2,2}
@@ -490,8 +490,8 @@ RECURSIVE SUBROUTINE CLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
490490
*
491491
* Compute T_{2,2} recursively
492492
*
493-
CALL CLARFT(DIRECT, STOREV, N, L, V(1, K-L+1), LDV, TAU(K-L+1),
494-
$ T(K-L+1,K-L+1), LDT)
493+
CALL CLARFT(DIRECT, STOREV, N, L, V(1, K-L+1), LDV,
494+
$ TAU(K-L+1), T(K-L+1,K-L+1), LDT)
495495
*
496496
* Compute T_{2,1}
497497
* T_{2,1} = V_{2,2}'
@@ -504,8 +504,8 @@ RECURSIVE SUBROUTINE CLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
504504
*
505505
* T_{2,1} = T_{2,1}*V_{2,1}
506506
*
507-
CALL CTRMM('Right', 'Upper', 'No transpose', 'Unit', L, K-L,
508-
$ ONE, V(N-K+1,1), LDV, T(K-L+1,1), LDT)
507+
CALL CTRMM('Right', 'Upper', 'No transpose', 'Unit', L,
508+
$ K-L, ONE, V(N-K+1,1), LDV, T(K-L+1,1), LDT)
509509

510510
*
511511
* T_{2,1} = V_{2,2}'*V_{2,1} + T_{2,1}
@@ -520,8 +520,8 @@ RECURSIVE SUBROUTINE CLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
520520
*
521521
* T_{2,1} = -T_{2,2}*T_{2,1}
522522
*
523-
CALL CTRMM('Left', 'Lower', 'No transpose', 'Non-unit', L, K-L,
524-
$ NEG_ONE, T(K-L+1,K-L+1), LDT, T(K-L+1,1), LDT)
523+
CALL CTRMM('Left', 'Lower', 'No transpose', 'Non-unit', L,
524+
$ K-L, NEG_ONE, T(K-L+1,K-L+1), LDT, T(K-L+1,1), LDT)
525525
*
526526
* T_{2,1} = T_{2,1}*T_{1,1}
527527
*
@@ -585,8 +585,8 @@ RECURSIVE SUBROUTINE CLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
585585
*
586586
* Compute T_{2,2} recursively
587587
*
588-
CALL CLARFT(DIRECT, STOREV, N, L, V(K-L+1,1), LDV, TAU(K-L+1),
589-
$ T(K-L+1,K-L+1), LDT)
588+
CALL CLARFT(DIRECT, STOREV, N, L, V(K-L+1,1), LDV,
589+
$ TAU(K-L+1), T(K-L+1,K-L+1), LDT)
590590
*
591591
* Compute T_{2,1}
592592
* T_{2,1} = V_{2,2}
@@ -597,8 +597,8 @@ RECURSIVE SUBROUTINE CLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
597597
*
598598
* T_{2,1} = T_{2,1}*V_{1,2}'
599599
*
600-
CALL CTRMM('Right', 'Lower', 'Conjugate', 'Unit', L, K-L, ONE,
601-
$ V(1, N-K+1), LDV, T(K-L+1,1), LDT)
600+
CALL CTRMM('Right', 'Lower', 'Conjugate', 'Unit', L, K-L,
601+
$ ONE, V(1, N-K+1), LDV, T(K-L+1,1), LDT)
602602

603603
*
604604
* T_{2,1} = V_{2,1}*V_{1,1}' + T_{2,1}
@@ -614,13 +614,13 @@ RECURSIVE SUBROUTINE CLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
614614
*
615615
* T_{2,1} = -T_{2,2}*T_{2,1}
616616
*
617-
CALL CTRMM('Left', 'Lower', 'No tranpose', 'Non-unit', L, K-L,
618-
$ NEG_ONE, T(K-L+1,K-L+1), LDT, T(K-L+1,1), LDT)
617+
CALL CTRMM('Left', 'Lower', 'No tranpose', 'Non-unit', L,
618+
$ K-L, NEG_ONE, T(K-L+1,K-L+1), LDT, T(K-L+1,1), LDT)
619619

620620
*
621621
* T_{2,1} = T_{2,1}*T_{1,1}
622622
*
623-
CALL CTRMM('Right', 'Lower', 'No tranpose', 'Non-unit', L, K-L,
624-
$ ONE, T, LDT, T(K-L+1,1), LDT)
623+
CALL CTRMM('Right', 'Lower', 'No tranpose', 'Non-unit', L,
624+
$ K-L, ONE, T, LDT, T(K-L+1,1), LDT)
625625
END IF
626626
END SUBROUTINE

‎SRC/dgelqf.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -251,8 +251,8 @@ SUBROUTINE DGELQF( M, N, A, LDA, TAU, WORK, LWORK, INFO )
251251
* Form the triangular factor of the block reflector
252252
* H = H(i) H(i+1) . . . H(i+ib-1)
253253
*
254-
CALL DLARFT( 'Forward', 'Rowwise', N-I+1, IB, A( I, I ),
255-
$ LDA, TAU( I ), WORK, LDWORK )
254+
CALL DLARFT( 'Forward', 'Rowwise', N-I+1, IB, A( I,
255+
$ I ), LDA, TAU( I ), WORK, LDWORK )
256256
*
257257
* Apply H to A(i+ib:m,i:n) from the right
258258
*

‎SRC/dgeqlf.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -256,8 +256,8 @@ SUBROUTINE DGEQLF( M, N, A, LDA, TAU, WORK, LWORK, INFO )
256256
* Form the triangular factor of the block reflector
257257
* H = H(i+ib-1) . . . H(i+1) H(i)
258258
*
259-
CALL DLARFT( 'Backward', 'Columnwise', M-K+I+IB-1,
260-
$ IB, A( 1, N-K+I ), LDA, TAU( I ), WORK,
259+
CALL DLARFT( 'Backward', 'Columnwise', M-K+I+IB-1, IB,
260+
$ A( 1, N-K+I ), LDA, TAU( I ), WORK,
261261
$ LDWORK )
262262
*
263263
* Apply H**T to A(1:m-k+i+ib-1,1:n-k+i-1) from the left

‎SRC/dgerqf.f‎

Lines changed: 1 addition & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -257,8 +257,7 @@ SUBROUTINE DGERQF( M, N, A, LDA, TAU, WORK, LWORK, INFO )
257257
* H = H(i+ib-1) . . . H(i+1) H(i)
258258
*
259259
CALL DLARFT( 'Backward', 'Rowwise', N-K+I+IB-1, IB,
260-
$ A( M-K+I, 1 ), LDA, TAU( I ), WORK,
261-
$ LDWORK )
260+
$ A( M-K+I, 1 ), LDA, TAU( I ), WORK, LDWORK )
262261
*
263262
* Apply H to A(1:m-k+i-1,1:n-k+i+ib-1) from the right
264263
*

‎SRC/dlarft.f‎

Lines changed: 22 additions & 22 deletions
Original file line numberDiff line numberDiff line change
@@ -159,8 +159,8 @@
159159
*> \endverbatim
160160
*>
161161
* =====================================================================
162-
RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
163-
$ LDT )
162+
RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV,
163+
$ TAU, T, LDT )
164164
*
165165
* -- LAPACK auxiliary routine --
166166
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -315,8 +315,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
315315
*
316316
* T_{1,2} = T_{1,2}*V_{2,2}
317317
*
318-
CALL DTRMM('Right', 'Lower', 'No transpose', 'Unit', L, K-L,
319-
$ ONE, V(L+1, L+1), LDV, T(1, L+1), LDT)
318+
CALL DTRMM('Right', 'Lower', 'No transpose', 'Unit', L,
319+
$ K-L, ONE, V(L+1, L+1), LDV, T(1, L+1), LDT)
320320

321321
*
322322
* T_{1,2} = V_{3,1}'*V_{3,2} + T_{1,2}
@@ -406,8 +406,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
406406
*
407407
* T_{1,2} = T_{1,2}*V_{2,2}'
408408
*
409-
CALL DTRMM('Right', 'Upper', 'Transpose', 'Unit', L, K-L, ONE,
410-
$ V(L+1, L+1), LDV, T(1, L+1), LDT)
409+
CALL DTRMM('Right', 'Upper', 'Transpose', 'Unit', L, K-L,
410+
$ ONE, V(L+1, L+1), LDV, T(1, L+1), LDT)
411411

412412
*
413413
* T_{1,2} = V_{1,3}*V_{2,3}' + T_{1,2}
@@ -422,8 +422,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
422422
*
423423
* T_{1,2} = -T_{1,1}*T_{1,2}
424424
*
425-
CALL DTRMM('Left', 'Upper', 'No transpose', 'Non-unit', L, K-L,
426-
$ NEG_ONE, T, LDT, T(1, L+1), LDT)
425+
CALL DTRMM('Left', 'Upper', 'No transpose', 'Non-unit', L,
426+
$ K-L, NEG_ONE, T, LDT, T(1, L+1), LDT)
427427

428428
*
429429
* T_{1,2} = T_{1,2}*T_{2,2}
@@ -486,8 +486,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
486486
*
487487
* Compute T_{2,2} recursively
488488
*
489-
CALL DLARFT(DIRECT, STOREV, N, L, V(1, K-L+1), LDV, TAU(K-L+1),
490-
$ T(K-L+1,K-L+1), LDT)
489+
CALL DLARFT(DIRECT, STOREV, N, L, V(1, K-L+1), LDV,
490+
$ TAU(K-L+1), T(K-L+1,K-L+1), LDT)
491491
*
492492
* Compute T_{2,1}
493493
* T_{2,1} = V_{2,2}'
@@ -500,8 +500,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
500500
*
501501
* T_{2,1} = T_{2,1}*V_{2,1}
502502
*
503-
CALL DTRMM('Right', 'Upper', 'No transpose', 'Unit', L, K-L,
504-
$ ONE, V(N-K+1,1), LDV, T(K-L+1,1), LDT)
503+
CALL DTRMM('Right', 'Upper', 'No transpose', 'Unit', L,
504+
$ K-L, ONE, V(N-K+1,1), LDV, T(K-L+1,1), LDT)
505505

506506
*
507507
* T_{2,1} = V_{2,2}'*V_{2,1} + T_{2,1}
@@ -516,8 +516,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
516516
*
517517
* T_{2,1} = -T_{2,2}*T_{2,1}
518518
*
519-
CALL DTRMM('Left', 'Lower', 'No transpose', 'Non-unit', L, K-L,
520-
$ NEG_ONE, T(K-L+1,K-L+1), LDT, T(K-L+1,1), LDT)
519+
CALL DTRMM('Left', 'Lower', 'No transpose', 'Non-unit', L,
520+
$ K-L, NEG_ONE, T(K-L+1,K-L+1), LDT, T(K-L+1,1), LDT)
521521
*
522522
* T_{2,1} = T_{2,1}*T_{1,1}
523523
*
@@ -581,8 +581,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
581581
*
582582
* Compute T_{2,2} recursively
583583
*
584-
CALL DLARFT(DIRECT, STOREV, N, L, V(K-L+1,1), LDV, TAU(K-L+1),
585-
$ T(K-L+1,K-L+1), LDT)
584+
CALL DLARFT(DIRECT, STOREV, N, L, V(K-L+1,1), LDV,
585+
$ TAU(K-L+1), T(K-L+1,K-L+1), LDT)
586586
*
587587
* Compute T_{2,1}
588588
* T_{2,1} = V_{2,2}
@@ -593,8 +593,8 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
593593
*
594594
* T_{2,1} = T_{2,1}*V_{1,2}'
595595
*
596-
CALL DTRMM('Right', 'Lower', 'Transpose', 'Unit', L, K-L, ONE,
597-
$ V(1, N-K+1), LDV, T(K-L+1,1), LDT)
596+
CALL DTRMM('Right', 'Lower', 'Transpose', 'Unit', L, K-L,
597+
$ ONE, V(1, N-K+1), LDV, T(K-L+1,1), LDT)
598598

599599
*
600600
* T_{2,1} = V_{2,1}*V_{1,1}' + T_{2,1}
@@ -610,13 +610,13 @@ RECURSIVE SUBROUTINE DLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T,
610610
*
611611
* T_{2,1} = -T_{2,2}*T_{2,1}
612612
*
613-
CALL DTRMM('Left', 'Lower', 'No tranpose', 'Non-unit', L, K-L,
614-
$ NEG_ONE, T(K-L+1,K-L+1), LDT, T(K-L+1,1), LDT)
613+
CALL DTRMM('Left', 'Lower', 'No tranpose', 'Non-unit', L,
614+
$ K-L, NEG_ONE, T(K-L+1,K-L+1), LDT, T(K-L+1,1), LDT)
615615

616616
*
617617
* T_{2,1} = T_{2,1}*T_{1,1}
618618
*
619-
CALL DTRMM('Right', 'Lower', 'No tranpose', 'Non-unit', L, K-L,
620-
$ ONE, T, LDT, T(K-L+1,1), LDT)
619+
CALL DTRMM('Right', 'Lower', 'No tranpose', 'Non-unit', L,
620+
$ K-L, ONE, T, LDT, T(K-L+1,1), LDT)
621621
END IF
622622
END SUBROUTINE

‎SRC/dorglq.f‎

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -253,7 +253,8 @@ SUBROUTINE DORGLQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
253253
* Form the triangular factor of the block reflector
254254
* H = H(i) H(i+1) . . . H(i+ib-1)
255255
*
256-
CALL DLARFT( 'Forward', 'Rowwise', N-I+1, IB, A( I, I ),
256+
CALL DLARFT( 'Forward', 'Rowwise', N-I+1, IB, A( I,
257+
$ I ),
257258
$ LDA, TAU( I ), WORK, LDWORK )
258259
*
259260
* Apply H**T to A(i+ib:m,i:n) from the right

‎SRC/dorgql.f‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -260,7 +260,7 @@ SUBROUTINE DORGQL( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
260260
* Form the triangular factor of the block reflector
261261
* H = H(i+ib-1) . . . H(i+1) H(i)
262262
*
263-
CALL DLARFT( 'Backward', 'Columnwise', M-K+I+IB-1, IB,
263+
CALL DLARFT( 'Backward', 'Columnwise', M-K+I+IB-1, IB,
264264
$ A( 1, N-K+I ), LDA, TAU( I ), WORK, LDWORK )
265265
*
266266
* Apply H to A(1:m-k+i+ib-1,1:n-k+i-1) from the left

‎SRC/dorgrq.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -261,8 +261,8 @@ SUBROUTINE DORGRQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
261261
* Form the triangular factor of the block reflector
262262
* H = H(i+ib-1) . . . H(i+1) H(i)
263263
*
264-
CALL DLARFT( 'Backward', 'Rowwise', N-K+I+IB-1, IB,
265-
$ A( II, 1 ), LDA, TAU( I ), WORK, LDWORK)
264+
CALL DLARFT( 'Backward', 'Rowwise', N-K+I+IB-1, IB,
265+
$ A( II, 1 ), LDA, TAU( I ), WORK, LDWORK )
266266
*
267267
* Apply H**T to A(1:m-k+i-1,1:n-k+i+ib-1) from the right
268268
*

‎SRC/dormlq.f‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -316,7 +316,7 @@ SUBROUTINE DORMLQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC,
316316
* Form the triangular factor of the block reflector
317317
* H = H(i) H(i+1) . . . H(i+ib-1)
318318
*
319-
CALL DLARFT( 'Forward', 'Rowwise', NQ-I+1, IB, A( I, I ),
319+
CALL DLARFT( 'Forward', 'Rowwise', NQ-I+1, IB, A( I, I ),
320320
$ LDA, TAU( I ), WORK( IWT ), LDT )
321321
IF( LEFT ) THEN
322322
*

‎SRC/dormql.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -310,8 +310,8 @@ SUBROUTINE DORMQL( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC,
310310
* Form the triangular factor of the block reflector
311311
* H = H(i+ib-1) . . . H(i+1) H(i)
312312
*
313-
CALL DLARFT( 'Backward', 'Columnwise', NQ-K+I+IB-1, IB,
314-
$ A( 1, I ), LDA, TAU( I ), WORK( IWT ), LDT)
313+
CALL DLARFT( 'Backward', 'Columnwise', NQ-K+I+IB-1, IB,
314+
$ A( 1, I ), LDA, TAU( I ), WORK( IWT ), LDT )
315315
IF( LEFT ) THEN
316316
*
317317
* H or H**T is applied to C(1:m-k+i+ib-1,1:n)

0 commit comments

Comments
 (0)