Skip to content

Commit c64027a

Browse files
scr2016scr2016
committed
DOWNDATED commnets in SRC/dgeqp3rk.f SRC/dlaqp2rk.f SRC/dlaqp3rk.f
1 parent 349a295 commit c64027a

3 files changed

Lines changed: 2 additions & 83 deletions

File tree

‎SRC/dgeqp3rk.f‎

Lines changed: 1 addition & 40 deletions
Original file line numberDiff line numberDiff line change
@@ -734,13 +734,6 @@ SUBROUTINE DGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
734734
WORK( J ) = DNRM2( M, A( 1, J ), 1 )
735735
WORK( N+J ) = WORK( J )
736736
END DO
737-
738-
WRITE(*,*)
739-
WRITE(*,*) "===== DGEQP3RK 2NORM ="
740-
WRITE(*,*)
741-
$ WORK( 1 ), WORK( 2 ), WORK( 3 ), WORK( 4 ),
742-
$ WORK( 5 ), WORK( 6 ), WORK( 7 ), WORK( 8 )
743-
WRITE(*,*)
744737
*
745738
* ==================================================================
746739
*
@@ -749,8 +742,6 @@ SUBROUTINE DGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
749742
*
750743
KP1 = IDAMAX( N, WORK( 1 ), 1 )
751744
MAXC2NRM = WORK( KP1 )
752-
WRITE(*,*) "======= DGEQP3RK ((( before NaN in MATRIX ))) KP1=",
753-
$ KP1
754745
*
755746
* ==================================================================.
756747
*
@@ -760,8 +751,6 @@ SUBROUTINE DGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
760751
* to the column number where the first NaN is found and return
761752
* from the routine.
762753
*
763-
WRITE(*,*) "======= DGEQP3RK ((( NaN in MATRIX ))) ===="
764-
765754
K = 0
766755
INFO = KP1
767756
*
@@ -783,8 +772,6 @@ SUBROUTINE DGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
783772
* Check is the matrix A is a zero matrix, set array TAU and
784773
* return from the routine.
785774
*
786-
WRITE(*,*) "======= DGEQP3RK ((( ZERO MATRIX ))) ===="
787-
788775
K = 0
789776
MAXC2NRMK = ZERO
790777
RELMAXC2NRMK = ZERO
@@ -947,10 +934,6 @@ SUBROUTINE DGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
947934
*
948935
* Factorize JB columns among the columns A(J:N).
949936
*
950-
WRITE(*,*)
951-
WRITE(*,*) "===== DGEQP3RK loop before block(IOFFSET, JB)=",
952-
$ J-1, JB
953-
954937
CALL DLAQP3RK( M, N_SUB, NRHS, IOFFSET, JB, ABSTOL,
955938
$ RELTOL, KP1, MAXC2NRM, A( 1, J ), LDA,
956939
$ DONE, JBF, MAXC2NRMK, RELMAXC2NRMK,
@@ -959,10 +942,6 @@ SUBROUTINE DGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
959942
$ WORK( 2*N+1 ), WORK( 2*N+JB+1 ),
960943
$ N+NRHS-J+1, IWORK, IINFO )
961944
*
962-
963-
WRITE(*,*) "======= DGEQP3RK loop after block (JBF)=",
964-
$ JBF
965-
*
966945
* Set INFO on the first occurence of Inf.
967946
*
968947
IF( IINFO.GT.N_SUB .AND. INFO.EQ.0 ) THEN
@@ -1039,9 +1018,6 @@ SUBROUTINE DGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
10391018
K = J - 1 + KF
10401019
*
10411020
* Set INFO on the first exception occurence.
1042-
*
1043-
WRITE(*,*) "======= DGEQP3RK after call to DLAQP2RK INFO=",
1044-
$ INFO
10451021
*
10461022
* Set INFO on the first exception occurence of Inf or NaN,
10471023
* (NaN takes precedence over Inf).
@@ -1078,22 +1054,7 @@ SUBROUTINE DGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
10781054
DO J = K + 1, MINMN
10791055
TAU( J ) = ZERO
10801056
END DO
1081-
1082-
1083-
WRITE(*,*)
1084-
WRITE(*,*) "===== END DGEQP3RK compute low rank ",
1085-
$ "(MAXC2NRMK, RELMAXC2NRMK)=",
1086-
$ MAXC2NRMK, RELMAXC2NRMK
1087-
1088-
ELSE
1089-
MAXC2NRMK = ZERO
1090-
RELMAXC2NRMK = ZERO
1091-
1092-
WRITE(*,*)
1093-
WRITE(*,*) "===== END DGEQP3RK compute full rank ",
1094-
$ "(MAXC2NRMK, RELMAXC2NRMK)=",
1095-
$ MAXC2NRMK, RELMAXC2NRMK
1096-
1057+
*
10971058
END IF
10981059
*
10991060
* END IF( J.LE.JMAX ) THEN

‎SRC/dlaqp2rk.f‎

Lines changed: 0 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -389,10 +389,6 @@ SUBROUTINE DLAQP2RK( M, N, NRHS, IOFFSET, KMAX, ABSTOL, RELTOL,
389389
KMAX = MIN( KMAX, MINMNFACT )
390390
TOL3Z = SQRT( DLAMCH( 'Epsilon' ) )
391391
HUGEVAL = DLAMCH( 'Overflow' )
392-
393-
WRITE(*,*) "####_#### Enter DLAQP2RK "
394-
WRITE(*,*) " (M, N, NRHS, IOFFSET, KMAX, KP1, MAXC2NRM)",
395-
$ M, N, NRHS, IOFFSET, KMAX, KP1, MAXC2NRM
396392
*
397393
* Compute the factorization, KK is the lomn loop index.
398394
*
@@ -410,8 +406,6 @@ SUBROUTINE DLAQP2RK( M, N, NRHS, IOFFSET, KMAX, ABSTOL, RELTOL,
410406
*
411407

412408
KP = KP1
413-
414-
WRITE(*,*) "## ## DLAQP2RK IOFFSET=0, first column KP=",KP
415409
*
416410
* ============================================================
417411
*
@@ -577,8 +571,6 @@ SUBROUTINE DLAQP2RK( M, N, NRHS, IOFFSET, KMAX, ABSTOL, RELTOL,
577571
ELSE
578572
TAU( KK ) = ZERO
579573
END IF
580-
581-
WRITE(*,*) "## DLAQP2RK (KK, TAU(Kk) )=", KK, TAU(KK)
582574
*
583575
* Check if TAU(KK) contains NaN, set INFO parameter
584576
* to the column number where NaN is found and return from
@@ -593,7 +585,6 @@ SUBROUTINE DLAQP2RK( M, N, NRHS, IOFFSET, KMAX, ABSTOL, RELTOL,
593585
IF( DISNAN( TAU(KK) ) ) THEN
594586
K = KK - 1
595587
INFO = KK
596-
WRITE(*,*) "## ## DLAQP2RK ((TAU is NaN)) (KK, INFO)", KK, INFO
597588
*
598589
* Set MAXC2NRMK and RELMAXC2NRMK to NaN.
599590
*

‎SRC/dlaqp3rk.f‎

Lines changed: 1 addition & 34 deletions
Original file line numberDiff line numberDiff line change
@@ -444,12 +444,6 @@ SUBROUTINE DLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
444444
NB = MIN( NB, MINMNFACT )
445445
TOL3Z = SQRT( DLAMCH( 'Epsilon' ) )
446446
HUGEVAL = DLAMCH( 'Overflow' )
447-
448-
WRITE(*,*) "$$$$_$$$$ Enter DLAQP3RK "
449-
WRITE(*,*) " (M, N, NRHS, IOFFSET, NB, KP1, MAXC2NRM)",
450-
$ M, N, NRHS, IOFFSET, NB, KP1, MAXC2NRM
451-
452-
453447
*
454448
* Compute factorization in a while loop over NB columns,
455449
* K is the column index in the block A(1:M,1:N).
@@ -548,10 +542,6 @@ SUBROUTINE DLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
548542
* routine.
549543
*
550544
IF( MAXC2NRMK.EQ.ZERO ) THEN
551-
552-
553-
WRITE(*,*) "$$$$$$ DLAQP3RK zero submatrix, IOFFSET, K= ",
554-
$ IOFFSET, K
555545
*
556546
DONE = .TRUE.
557547
*
@@ -561,10 +551,6 @@ SUBROUTINE DLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
561551
* Set IF, the number of processed rows in the block, which
562552
* is the same as the number of processed rows in
563553
* the original whole matrix A_orig.
564-
*
565-
WRITE(*,*)
566-
$ "$$$$$$$$ DLAQP3RK zero submatrix (ABSTOL, K)= ",
567-
$ ABSTOL, K
568554
*
569555
KB = K - 1
570556
IF = I - 1
@@ -582,12 +568,6 @@ SUBROUTINE DLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
582568
* A(I+1:M,1:KB) * F(N+1:N+NRHS,1:KB)**T.
583569
*
584570
IF( NRHS.GT.0 .AND. KB.LT.(M-IOFFSET) ) THEN
585-
586-
587-
WRITE(*,*) "$$$$$$$$$$ DLAQP3RK zero submatrix",
588-
$ " block reflector (M-IF, NRHS, KB)",
589-
$ M-IF, NRHS, KB
590-
591571
CALL DGEMM( 'No transpose', 'Transpose',
592572
$ M-IF, NRHS, KB, -ONE, A( IF+1, 1 ), LDA,
593573
$ F( N+1, 1 ), LDF, ONE, A( IF+1, N+1 ), LDA )
@@ -652,13 +632,7 @@ SUBROUTINE DLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
652632
*
653633
KB = K - 1
654634
IF = I - 1
655-
656-
WRITE(*,*) "$$$$$$$$$$ DLAQP3RK condition for",
657-
$ " ABSTOL or RELTOL (ABSTOL, RELTOL),",
658-
$ " (MAXC2NRMK, RELMAXC2NRMK)",
659-
$ ABSTOL, RELTOL, MAXC2NRMK, RELMAXC2NRMK
660-
661-
635+
*
662636
* Apply the block reflector to the residual of the
663637
* matrix A and the residual of the right hand sides B, if
664638
* the residual matrix and and/or the residual of the right
@@ -670,11 +644,6 @@ SUBROUTINE DLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
670644
* A(IF+1:M,1:KB) * F(KB+1:N+NRHS,1:KB)**T.
671645
*
672646
IF( KB.LT.MINMNUPDT ) THEN
673-
674-
WRITE(*,*) "$$$$$$$$$$ DLAQP3RK ABSTOL or RELTOL",
675-
$ " block reflector (M-IF, N+NRHS-KB, KB)",
676-
$ M-IF, N+NRHS-KB, KB
677-
678647
CALL DGEMM( 'No transpose', 'Transpose',
679648
$ M-IF, N+NRHS-KB, KB,-ONE, A( IF+1, 1 ), LDA,
680649
$ F( KB+1, 1 ), LDF, ONE, A( IF+1, KB+1 ), LDA )
@@ -767,8 +736,6 @@ SUBROUTINE DLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
767736
KB = K - 1
768737
IF = I - 1
769738
INFO = K
770-
771-
WRITE(*,*) "$$ $$ DLAQP3RK ((TAU is NaN)) (K, INFO)", K, INFO
772739
*
773740
* Set MAXC2NRMK and RELMAXC2NRMK to NaN.
774741
*

0 commit comments

Comments
 (0)