Skip to content

Commit d086a3d

Browse files
scr2016scr2016
committed
DOWNDATED comments in LIN/zchkqp3rk.f
1 parent bfbd799 commit d086a3d

File tree

1 file changed

+1
-140
lines changed

1 file changed

+1
-140
lines changed

‎TESTING/LIN/zchkqp3rk.f

Lines changed: 1 addition & 140 deletions
Original file line numberDiff line numberDiff line change
@@ -319,14 +319,6 @@ SUBROUTINE ZCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
319319
*
320320
IF( .NOT.DOTYPE( IMAT ) )
321321
$ CYCLE
322-
323-
WRITE(*,*) "(4) ____ ____ ____ ____ Loop for "
324-
WRITE(*,*) "IMAT, DOTYPE(IMAT), M, N, NRHS = ",
325-
$ IMAT, DOTYPE(IMAT), M, N, NRHS
326-
327-
328-
329-
330322
*
331323
* The type of distribution used to generate the random
332324
* eigen-/singular values:
@@ -363,8 +355,6 @@ SUBROUTINE ZCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
363355
DO I = 1, MINMN
364356
S( I ) = ZERO
365357
END DO
366-
367-
WRITE(*,*) "GENERATED CZERO MATRIX"
368358
*
369359
ELSE IF( (IMAT.GE.2 .AND. IMAT.LE.4 )
370360
$ .OR. (IMAT.GE.14 .AND. IMAT.LE.19 ) ) THEN
@@ -490,9 +480,6 @@ SUBROUTINE ZCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
490480
*
491481
* 1) Set the first NB_ZERO columns in COPYA(1:M,1:N)
492482
* to zero.
493-
494-
WRITE(*,*) "( M,NB_ZERO, LDA before ZLASET",
495-
$ M, NB_ZERO, LDA
496483
*
497484
CALL ZLASET( 'Full', M, NB_ZERO, CZERO, CZERO,
498485
$ COPYA, LDA )
@@ -505,21 +492,13 @@ SUBROUTINE ZCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
505492
$ ANORM, MODE, CNDNUM, DIST )
506493
*
507494
SRNAMT = 'ZLATMS'
508-
509-
WRITE(*,*) "Generate DLAMTS matrix (M,NB_GEN)=",
510-
$ M, NB_GEN
511-
495+
*
512496
IND_OFFSET_GEN = NB_ZERO * LDA
513497
*
514498
CALL ZLATMS( M, NB_GEN, DIST, ISEED, TYPE, S, MODE,
515499
$ CNDNUM, ANORM, KL, KU, 'No packing',
516500
$ COPYA( IND_OFFSET_GEN + 1 ), LDA,
517501
$ WORK, INFO )
518-
519-
WRITE(*,*) "Singular val after mat generation S="
520-
WRITE(*,*)
521-
$ S(1), S(2), S(3), S(4),
522-
$ S(5), S(6), S(7), S(8)
523502
*
524503
* Check error code from ZLATMS.
525504
*
@@ -529,18 +508,6 @@ SUBROUTINE ZCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
529508
$ NERRS, NOUT )
530509
CYCLE
531510
END IF
532-
533-
WRITE(*,*) "Matrix COPYA after",
534-
$ " generation N_ZERO and N_GEN"
535-
DO I = 1, M
536-
WRITE(*,*)
537-
$ COPYA((1-1)*LDA+I), COPYA((2-1)*LDA+I),
538-
$ COPYA((3-1)*LDA+I), COPYA((4-1)*LDA+I),
539-
$ COPYA((5-1)*LDA+I), COPYA((6-1)*LDA+I),
540-
$ COPYA((7-1)*LDA+I), COPYA((8-1)*LDA+I)
541-
END DO
542-
543-
544511
*
545512
* 3) Swap the gererated colums from the right side
546513
* NB_GEN-size block in COPYA into correct column
@@ -592,61 +559,19 @@ SUBROUTINE ZCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
592559
* The total number of singular values is MINMN.
593560
*
594561
MINMNB_GEN = MIN( M, NB_GEN )
595-
596-
WRITE(*,*) "Singular values MINB_GEN=", MINMNB_GEN
597-
WRITE(*,*) "Singular values before ordering S="
598-
WRITE(*,*)
599-
$ S(1), S(2), S(3), S(4),
600-
$ S(5), S(6), S(7), S(8)
601562
*
602563
CALL DLAORD( 'Decreasing', MINMNB_GEN, S, 1 )
603564

604-
WRITE(*,*) "Singular values after ordering S="
605-
WRITE(*,*)
606-
$ S(1), S(2), S(3), S(4),
607-
$ S(5), S(6), S(7), S(8)
608-
609565
DO I = MINMNB_GEN+1, MINMN
610-
WRITE(*,*) "Singular values to ZERO out I=", I
611566
S( I ) = ZERO
612567
END DO
613-
614-
WRITE(*,*) "Matrix with ZERO columnms COPYA"
615-
DO I = 1, M
616-
WRITE(*,*)
617-
$ COPYA((1-1)*LDA+I), COPYA((2-1)*LDA+I),
618-
$ COPYA((3-1)*LDA+I), COPYA((4-1)*LDA+I),
619-
$ COPYA((5-1)*LDA+I), COPYA((6-1)*LDA+I),
620-
$ COPYA((7-1)*LDA+I), COPYA((8-1)*LDA+I)
621-
END DO
622-
623-
624-
WRITE(*,*) "Matrix with ZERO columns ordering of S="
625-
WRITE(*,*)
626-
$ S(1), S(2), S(3), S(4),
627-
$ S(5), S(6), S(7), S(8)
628568
*
629569
ELSE
630570
*
631571
* IF(MINMN.LT.2) skip this size for this matrix type.
632572
*
633573
CYCLE
634574
END IF
635-
*
636-
WRITE(*,*) "AFTER GENERATING COPYA"
637-
DO I = 1, M
638-
WRITE(*,*)
639-
$ COPYA((1-1)*LDA+I), COPYA((2-1)*LDA+I),
640-
$ COPYA((3-1)*LDA+I), COPYA((4-1)*LDA+I),
641-
$ COPYA((5-1)*LDA+I), COPYA((6-1)*LDA+I)
642-
** $ C,OPYA((7-1)*LDA+I), COPYA((8-1)*LDA+I)
643-
END DO
644-
WRITE(*,*) "AFTER GENERATING COPYA, S="
645-
WRITE(*,*)
646-
$ S(1), S(2), S(3), S(4),
647-
$ S(5), S(6)
648-
** $ , S(7), S(8)
649-
650575
*
651576
* Initialize a copy array for a pivot array for DGEQP3RK.
652577
*
@@ -662,22 +587,12 @@ SUBROUTINE ZCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
662587
CALL XLAENV( 1, NB )
663588
NX = NXVAL( INB )
664589
CALL XLAENV( 3, NX )
665-
*
666-
667-
WRITE(*,*) "(5) ____ ____ ____ ____ ____ Loop for NB,NX=",
668-
$ NB, NX
669-
670-
671590
*
672591
* We do MIN(M,N)+1 because we need a test for KMAX > N,
673592
* when KMAX is larger than MIN(M,N), KMAX should be
674593
* KMAX = MIN(M,N)
675594
*
676595
DO KMAX = 0, MIN(M,N)+1
677-
678-
679-
WRITE(*,*) "(6) ____ ____ ____ ____ ____" ,
680-
$ " ____ Loop for KMAX=", KMAX
681596
*
682597
* Get a working copy of COPYA into A( 1:M,1:N ).
683598
* Get a working copy of COPYB into A( 1:M, (N+1):NRHS ).
@@ -696,29 +611,6 @@ SUBROUTINE ZCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
696611
*
697612
ABSTOL = -1.0
698613
RELTOl = -1.0
699-
700-
701-
WRITE(*,*) "(BEGIN) B VALUE OF B after COPYB into B"
702-
703-
DO I = 1, LDA
704-
WRITE(*,*)
705-
$ B((1-1)*LDA+I), B((2-1)*LDA+I),
706-
$ B((3-1)*LDA+I), B((4-1)*LDA+I),
707-
$ B((5-1)*LDA+I), B((6-1)*LDA+I)
708-
* $ ,B((7-1)*LDA+I), B((8-1)*LDA+I)
709-
* $ ,B((9-1)*LDA+I), B((10-1)*LDA+I)
710-
* $ ,B((11-1)*LDA+I), B((12-1)*LDA+I)
711-
END DO
712-
713-
WRITE(*,*)
714-
WRITE(*,*) "(BEGIN) A after COPYA and COPYB into A"
715-
DO I = 1, LDA
716-
WRITE(*,*)
717-
$ A((1-1)*LDA+I), A((2-1)*LDA+I),
718-
$ A((3-1)*LDA+I), A((4-1)*LDA+I),
719-
$ A((5-1)*LDA+I), A((6-1)*LDA+I)
720-
END DO
721-
WRITE(*,*)
722614
*
723615
* Compute the QR factorization with pivoting of A
724616
*
@@ -733,31 +625,6 @@ SUBROUTINE ZCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
733625
$ RELMAXC2NRMK, IWORK( N+1 ), TAU,
734626
$ WORK, LW, RWORK, IWORK( 2*N+1 ),
735627
$ INFO )
736-
737-
WRITE(*,*)
738-
WRITE(*,*) "(END) A after ZGEQP3RK"
739-
DO I = 1, LDA
740-
WRITE(*,*)
741-
$ A((1-1)*LDA+I), A((2-1)*LDA+I),
742-
$ A((3-1)*LDA+I), A((4-1)*LDA+I),
743-
$ A((5-1)*LDA+I), A((6-1)*LDA+I)
744-
END DO
745-
WRITE(*,*)
746-
747-
WRITE(*,*) "TAU after ZGEQP3RK"
748-
WRITE(*,*)
749-
$ TAU(1), TAU(2), TAU(3), TAU(4),
750-
$ TAU(5), TAU(6)
751-
WRITE(*,*)
752-
753-
WRITE(*,*) "JPIV after ZGEQP3RK",
754-
$ IWORK(N+1), IWORK(N+2), IWORK(N+3), IWORK(N+4),
755-
$ IWORK(N+5), IWORK(N+6)
756-
WRITE(*,*)
757-
*
758-
WRITE(*,*) "INFO after ZGEQP3RK"
759-
WRITE(*,*) INFO
760-
WRITE(*,*)
761628
*
762629
* Check error code from ZGEQP3RK.
763630
*
@@ -766,12 +633,6 @@ SUBROUTINE ZCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
766633
$ M, N, NX, -1, NB, IMAT,
767634
$ NFAIL, NERRS, NOUT )
768635
*
769-
770-
771-
WRITE(*,*) "After ZGEQP3RK, KFACT, MAX2N, REL2N",
772-
$ KFACT, MAXC2NRMK, RELMAXC2NRMK
773-
774-
775636
IF( KFACT.EQ.MINMN ) THEN
776637
*
777638
* Compute test 1:

0 commit comments

Comments
 (0)