Skip to content

Commit bfbd799

Browse files
scr2016scr2016
committed
DOWNDATED comments in SRC/zgeqp3rk.f SRC/zlaqp2rk.f SRC/zlaqp3rk.f
1 parent 7656583 commit bfbd799

File tree

3 files changed

+5
-152
lines changed

3 files changed

+5
-152
lines changed

‎SRC/zgeqp3rk.f

Lines changed: 2 additions & 50 deletions
Original file line numberDiff line numberDiff line change
@@ -741,23 +741,13 @@ SUBROUTINE ZGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
741741
RWORK( J ) = DZNRM2( M, A( 1, J ), 1 )
742742
RWORK( N+J ) = RWORK( J )
743743
END DO
744-
745-
WRITE(*,*)
746-
WRITE(*,*) "===== ZGEQP3RK 2NORM ="
747-
WRITE(*,*)
748-
$ RWORK( 1 ), RWORK( 2 ), RWORK( 3 ), RWORK( 4 ),
749-
$ RWORK( 5 ), RWORK( 6 ), RWORK( 7 ), RWORK( 8 )
750-
WRITE(*,*)
751744
*
752745
* ==================================================================
753746
*
754747
* Compute the pivot column index and the maximum column 2-norm
755748
* for the whole original matrix stored in A(1:M,1:N).
756749
*
757750
KP1 = IDAMAX( N, RWORK( 1 ), 1 )
758-
MAXC2NRM = RWORK( KP1 )
759-
WRITE(*,*) "======= ZGEQP3RK ((( before NaN in MATRIX ))) KP1=",
760-
$ KP1
761751
*
762752
* ==================================================================.
763753
*
@@ -767,8 +757,6 @@ SUBROUTINE ZGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
767757
* to the column number where the first NaN is found and return
768758
* from the routine.
769759
*
770-
WRITE(*,*) "======= ZGEQP3RK ((( NaN in MATRIX ))) ===="
771-
772760
K = 0
773761
INFO = KP1
774762
*
@@ -790,8 +778,6 @@ SUBROUTINE ZGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
790778
* Check is the matrix A is a zero matrix, set array TAU and
791779
* return from the routine.
792780
*
793-
WRITE(*,*) "======= ZGEQP3RK ((( ZERO MATRIX ))) ===="
794-
795781
K = 0
796782
MAXC2NRMK = ZERO
797783
RELMAXC2NRMK = ZERO
@@ -954,10 +940,6 @@ SUBROUTINE ZGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
954940
*
955941
* Factorize JB columns among the columns A(J:N).
956942
*
957-
WRITE(*,*)
958-
WRITE(*,*) "===== ZGEQP3RK loop before block(IOFFSET, JB)=",
959-
$ J-1, JB
960-
961943
CALL ZLAQP3RK( M, N_SUB, NRHS, IOFFSET, JB, ABSTOL,
962944
$ RELTOL, KP1, MAXC2NRM, A( 1, J ), LDA,
963945
$ DONE, JBF, MAXC2NRMK, RELMAXC2NRMK,
@@ -966,10 +948,6 @@ SUBROUTINE ZGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
966948
$ WORK( 1 ), WORK( JB+1 ),
967949
$ N+NRHS-J+1, IWORK, IINFO )
968950
*
969-
970-
WRITE(*,*) "======= ZGEQP3RK loop after block (JBF)=",
971-
$ JBF
972-
*
973951
* Set INFO on the first occurence of Inf.
974952
*
975953
IF( IINFO.GT.N_SUB .AND. INFO.EQ.0 ) THEN
@@ -1046,9 +1024,6 @@ SUBROUTINE ZGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
10461024
K = J - 1 + KF
10471025
*
10481026
* Set INFO on the first exception occurence.
1049-
*
1050-
WRITE(*,*) "======= ZGEQP3RK after call to ZLAQP2RK INFO=",
1051-
$ INFO
10521027
*
10531028
* Set INFO on the first exception occurence of Inf or NaN,
10541029
* (NaN takes precedence over Inf).
@@ -1085,39 +1060,16 @@ SUBROUTINE ZGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
10851060
DO J = K + 1, MINMN
10861061
TAU( J ) = CZERO
10871062
END DO
1088-
1089-
1090-
WRITE(*,*)
1091-
WRITE(*,*) "===== END ZGEQP3RK compute low rank ",
1092-
$ "(MAXC2NRMK, RELMAXC2NRMK)=",
1093-
$ MAXC2NRMK, RELMAXC2NRMK
1094-
1063+
*
10951064
ELSE
10961065
MAXC2NRMK = ZERO
10971066
RELMAXC2NRMK = ZERO
1098-
1099-
WRITE(*,*)
1100-
WRITE(*,*) "===== END ZGEQP3RK compute full rank ",
1101-
$ "(MAXC2NRMK, RELMAXC2NRMK)=",
1102-
$ MAXC2NRMK, RELMAXC2NRMK
1103-
1067+
*
11041068
END IF
11051069
*
11061070
* END IF( J.LE.JMAX ) THEN
11071071
*
11081072
END IF
1109-
1110-
1111-
DO J = 1 , MINMN
1112-
WRITE(*,*) "===== END ZGEQP3RK before return (J,TAU(J))",
1113-
$ J, TAU( J )
1114-
END DO
1115-
1116-
WRITE(*,*)
1117-
WRITE(*,*) "===== END ZGEQP3RK before return ",
1118-
$ "(K, MAXC2NRMK, RELMAXC2NRMK)=",
1119-
$ K, MAXC2NRMK, RELMAXC2NRMK
1120-
11211073
*
11221074
WORK( 1 ) = DCMPLX( LWKOPT )
11231075
*

‎SRC/zlaqp2rk.f

Lines changed: 0 additions & 35 deletions
Original file line numberDiff line numberDiff line change
@@ -394,14 +394,6 @@ SUBROUTINE ZLAQP2RK( M, N, NRHS, IOFFSET, KMAX, ABSTOL, RELTOL,
394394
KMAX = MIN( KMAX, MINMNFACT )
395395
TOL3Z = SQRT( DLAMCH( 'Epsilon' ) )
396396
HUGEVAL = DLAMCH( 'Overflow' )
397-
398-
WRITE(*,*)
399-
WRITE(*,*) "####_#### Enter ZLAQP2RK "
400-
WRITE(*,*) " (M, N, NRHS, IOFFSET, KMAX, KP1,",
401-
$ " MAXC2NRM)",
402-
$ M, N, NRHS, IOFFSET, KMAX, KP1, MAXC2NRM
403-
WRITE(*,*)
404-
405397
*
406398
* Compute the factorization, KK is the lomn loop index.
407399
*
@@ -417,11 +409,7 @@ SUBROUTINE ZLAQP2RK( M, N, NRHS, IOFFSET, KMAX, ABSTOL, RELTOL,
417409
* therefore we use the computed KP1 and MAXC2NRM from the
418410
* main routine.
419411
*
420-
421412
KP = KP1
422-
423-
WRITE(*,*) " _#### ZLAQP2RK IOFFSET=0,
424-
$ first column KP=",KP
425413
*
426414
* ============================================================
427415
*
@@ -587,8 +575,6 @@ SUBROUTINE ZLAQP2RK( M, N, NRHS, IOFFSET, KMAX, ABSTOL, RELTOL,
587575
ELSE
588576
TAU( KK ) = CZERO
589577
END IF
590-
591-
WRITE(*,*) " _#### ZLAQP2RK (K, TAU(Kk) )=", KK, TAU(KK)
592578
*
593579
* Check if TAU(KK) contains NaN, set INFO parameter
594580
* to the column number where NaN is found and return from
@@ -611,9 +597,6 @@ SUBROUTINE ZLAQP2RK( M, N, NRHS, IOFFSET, KMAX, ABSTOL, RELTOL,
611597
IF( DISNAN( TAUNAN ) ) THEN
612598
K = KK - 1
613599
INFO = KK
614-
WRITE(*,*) " _#### ZLAQP2RK ((TAU is NaN))",
615-
$ " (KK, INFO),TAU",
616-
$ KK, INFO, TAU(KK)
617600
*
618601
* Set MAXC2NRMK and RELMAXC2NRMK to NaN.
619602
*
@@ -701,10 +684,6 @@ SUBROUTINE ZLAQP2RK( M, N, NRHS, IOFFSET, KMAX, ABSTOL, RELTOL,
701684
* factorized, we need to set MAXC2NRMK and RELMAXC2NRMK before
702685
* we return.
703686
*
704-
705-
WRITE(*,*) " _#### ZLAQP2RK after loop (K,MINMNFACT)", K,
706-
$ MINMNFACT
707-
708687
IF( K.LT.MINMNFACT ) THEN
709688
*
710689
JMAXC2NRM = K + IDAMAX( N-K, VN1( K+1 ), 1 )
@@ -728,20 +707,6 @@ SUBROUTINE ZLAQP2RK( M, N, NRHS, IOFFSET, KMAX, ABSTOL, RELTOL,
728707
DO J = K + 1, MINMNFACT
729708
TAU( J ) = CZERO
730709
END DO
731-
732-
DO J = 1 , MINMNFACT
733-
WRITE(*,*) " _#### END ZLAQP2RK before return",
734-
$ " (J,TAU(J))",
735-
$ J, TAU( J )
736-
END DO
737-
738-
739-
WRITE(*,*)
740-
WRITE(*,*) " _#### Exit ZLAQP2RK "
741-
WRITE(*,*) " (K, MAXC2NRMK, RELMAXC2NRMK)",
742-
$ K, MAXC2NRMK, RELMAXC2NRMK
743-
WRITE(*,*)
744-
745710
*
746711
RETURN
747712
*

‎SRC/zlaqp3rk.f

Lines changed: 3 additions & 67 deletions
Original file line numberDiff line numberDiff line change
@@ -448,12 +448,6 @@ SUBROUTINE ZLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
448448
NB = MIN( NB, MINMNFACT )
449449
TOL3Z = SQRT( DLAMCH( 'Epsilon' ) )
450450
HUGEVAL = DLAMCH( 'Overflow' )
451-
452-
WRITE(*,*) "$$$$_$$$$ Enter ZLAQP3RK "
453-
WRITE(*,*) " (M, N, NRHS, IOFFSET, NB, KP1, MAXC2NRM)",
454-
$ M, N, NRHS, IOFFSET, NB, KP1, MAXC2NRM
455-
456-
457451
*
458452
* Compute factorization in a while loop over NB columns,
459453
* K is the column index in the block A(1:M,1:N).
@@ -552,10 +546,6 @@ SUBROUTINE ZLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
552546
* routine.
553547
*
554548
IF( MAXC2NRMK.EQ.ZERO ) THEN
555-
556-
557-
WRITE(*,*) "$$$$$$ ZLAQP3RK zero submatrix, IOFFSET, K= ",
558-
$ IOFFSET, K
559549
*
560550
DONE = .TRUE.
561551
*
@@ -565,10 +555,6 @@ SUBROUTINE ZLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
565555
* Set IF, the number of processed rows in the block, which
566556
* is the same as the number of processed rows in
567557
* the original whole matrix A_orig.
568-
*
569-
WRITE(*,*)
570-
$ "$$$$$$$$ ZLAQP3RK zero submatrix (ABSTOL, K)= ",
571-
$ ABSTOL, K
572558
*
573559
KB = K - 1
574560
IF = I - 1
@@ -586,12 +572,7 @@ SUBROUTINE ZLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
586572
* A(I+1:M,1:KB) * F(N+1:N+NRHS,1:KB)**H.
587573
*
588574
IF( NRHS.GT.0 .AND. KB.LT.(M-IOFFSET) ) THEN
589-
590-
591-
WRITE(*,*) "$$$$$$$$$$ ZLAQP3RK zero submatrix",
592-
$ " block reflector (M-IF, NRHS, KB)",
593-
$ M-IF, NRHS, KB
594-
575+
*
595576
CALL ZGEMM( 'No transpose', 'Conjugate transpose',
596577
$ M-IF, NRHS, KB, -CONE, A( IF+1, 1 ), LDA,
597578
$ F( N+1, 1 ), LDF, CONE, A( IF+1, N+1 ), LDA )
@@ -656,13 +637,7 @@ SUBROUTINE ZLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
656637
*
657638
KB = K - 1
658639
IF = I - 1
659-
660-
WRITE(*,*) "$$$$$$$$$$ ZLAQP3RK condition for",
661-
$ " ABSTOL or RELTOL (ABSTOL, RELTOL),",
662-
$ " (MAXC2NRMK, RELMAXC2NRMK)",
663-
$ ABSTOL, RELTOL, MAXC2NRMK, RELMAXC2NRMK
664-
665-
640+
*
666641
* Apply the block reflector to the residual of the
667642
* matrix A and the residual of the right hand sides B, if
668643
* the residual matrix and and/or the residual of the right
@@ -674,11 +649,7 @@ SUBROUTINE ZLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
674649
* A(IF+1:M,1:KB) * F(KB+1:N+NRHS,1:KB)**H.
675650
*
676651
IF( KB.LT.MINMNUPDT ) THEN
677-
678-
WRITE(*,*) "$$$$$$$$$$ ZLAQP3RK ABSTOL or RELTOL",
679-
$ " block reflector (M-IF, N+NRHS-KB, KB)",
680-
$ M-IF, N+NRHS-KB, KB
681-
652+
*
682653
CALL ZGEMM( 'No transpose', 'Conjugate transpose',
683654
$ M-IF, N+NRHS-KB, KB,-CONE, A( IF+1, 1 ), LDA,
684655
$ F( KB+1, 1 ), LDF, CONE, A( IF+1, KB+1 ), LDA )
@@ -785,9 +756,6 @@ SUBROUTINE ZLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
785756
KB = K - 1
786757
IF = I - 1
787758
INFO = K
788-
789-
WRITE(*,*) "## ZLAQP2RK ((TAU is NaN)) (K, INFO),TAU",
790-
$ K, INFO, TAU(K)
791759
*
792760
* Set MAXC2NRMK and RELMAXC2NRMK to NaN.
793761
*
@@ -934,11 +902,6 @@ SUBROUTINE ZLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
934902
* A(IF+1:M,K+1:N+NRHS) := A(IF+1:M,KB+1:N+NRHS) -
935903
* A(IF+1:M,1:KB) * F(KB+1:N+NRHS,1:KB)**H.
936904
*
937-
WRITE(*,*) " _$$$$ ZLAQP3RK A before ZGEMM"
938-
WRITE(*,*) " ( KB, IF, MINMNUPDT, NRHS )",
939-
$ KB, IF, MINMNUPDT, NRHS
940-
941-
942905
IF( KB.LT.MINMNUPDT ) THEN
943906
CALL ZGEMM( 'No transpose', 'Conjugate transpose',
944907
$ M-IF, N+NRHS-KB, KB, -CONE, A( IF+1, 1 ), LDA,
@@ -949,20 +912,7 @@ SUBROUTINE ZLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
949912
* Loop over the index of the difficult columns from the largest
950913
* to the smallest index.
951914
*
952-
953-
WRITE(*,*) " _$$$$ ZLAQP3RK A after ZGEMM"
954-
DO J = 1, M
955-
WRITE(*,*)
956-
$ A(J,1), A(J,2),
957-
$ A(J,3), A(J,4),
958-
$ A(J, 5), A(J,6)
959-
END DO
960-
WRITE(*,*)
961-
962-
963915
DO WHILE( LSTICC.GT.0 )
964-
965-
WRITE(*,*) " _$$$$ ZLAQP3RK inside bad norm updating"
966916
*
967917
* LSTICC is the index of the last difficult column is greater
968918
* than 1.
@@ -986,20 +936,6 @@ SUBROUTINE ZLAQP3RK( M, N, NRHS, IOFFSET, NB, ABSTOL,
986936
LSTICC = ITEMP
987937
*
988938
END DO
989-
990-
DO J = 1 , MINMNFACT
991-
WRITE(*,*) " _$$$$ END ZLAQP3RK before return",
992-
$ " (J,TAU(J))",
993-
$ J, TAU( J )
994-
END DO
995-
996-
997-
WRITE(*,*)
998-
WRITE(*,*) " _$$$$ Exit ZLAQP2RK "
999-
WRITE(*,*) " (KB, MAXC2NRMK, RELMAXC2NRMK)",
1000-
$ KB, MAXC2NRMK, RELMAXC2NRMK
1001-
WRITE(*,*)
1002-
1003939
*
1004940
RETURN
1005941
*

0 commit comments

Comments
 (0)