From 423cba722e6518d9e86363960fff36e6981039e7 Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Tue, 2 Jun 2020 10:26:57 +1000 Subject: [PATCH] code simplified: re-computation avoided --- SRC/cunmlq.f | 2 +- SRC/cunmql.f | 2 +- SRC/cunmqr.f | 2 +- SRC/cunmrq.f | 2 +- SRC/cunmrz.f | 2 +- SRC/dormlq.f | 2 +- SRC/dormql.f | 2 +- SRC/dormqr.f | 2 +- SRC/dormrq.f | 2 +- SRC/dormrz.f | 2 +- SRC/sormlq.f | 2 +- SRC/sormql.f | 2 +- SRC/sormqr.f | 2 +- SRC/sormrq.f | 2 +- SRC/sormrz.f | 2 +- SRC/zunmlq.f | 2 +- SRC/zunmql.f | 2 +- SRC/zunmqr.f | 2 +- SRC/zunmrq.f | 2 +- SRC/zunmrz.f | 2 +- 20 files changed, 20 insertions(+), 20 deletions(-) diff --git a/SRC/cunmlq.f b/SRC/cunmlq.f index 3975b000b0..25a4107708 100644 --- a/SRC/cunmlq.f +++ b/SRC/cunmlq.f @@ -271,7 +271,7 @@ SUBROUTINE CUNMLQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'CUNMLQ', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/cunmql.f b/SRC/cunmql.f index 19ae9c2061..3c71660663 100644 --- a/SRC/cunmql.f +++ b/SRC/cunmql.f @@ -270,7 +270,7 @@ SUBROUTINE CUNMQL( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.(NW*NB+TSIZE) ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'CUNMQL', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/cunmqr.f b/SRC/cunmqr.f index 23624222c1..7e59d71297 100644 --- a/SRC/cunmqr.f +++ b/SRC/cunmqr.f @@ -265,7 +265,7 @@ SUBROUTINE CUNMQR( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'CUNMQR', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/cunmrq.f b/SRC/cunmrq.f index 902e279279..5a233f6040 100644 --- a/SRC/cunmrq.f +++ b/SRC/cunmrq.f @@ -269,7 +269,7 @@ SUBROUTINE CUNMRQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'CUNMRQ', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/cunmrz.f b/SRC/cunmrz.f index 76f8ef6cba..8e06f2329a 100644 --- a/SRC/cunmrz.f +++ b/SRC/cunmrz.f @@ -294,7 +294,7 @@ SUBROUTINE CUNMRZ( SIDE, TRANS, M, N, K, L, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'CUNMRQ', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/dormlq.f b/SRC/dormlq.f index b93e6e28b2..ef039285ab 100644 --- a/SRC/dormlq.f +++ b/SRC/dormlq.f @@ -264,7 +264,7 @@ SUBROUTINE DORMLQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'DORMLQ', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/dormql.f b/SRC/dormql.f index fb8d619800..7c9f189e0d 100644 --- a/SRC/dormql.f +++ b/SRC/dormql.f @@ -266,7 +266,7 @@ SUBROUTINE DORMQL( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'DORMQL', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/dormqr.f b/SRC/dormqr.f index 4403421b6f..4d0bae3a5f 100644 --- a/SRC/dormqr.f +++ b/SRC/dormqr.f @@ -263,7 +263,7 @@ SUBROUTINE DORMQR( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'DORMQR', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/dormrq.f b/SRC/dormrq.f index c26e5af096..43422f9583 100644 --- a/SRC/dormrq.f +++ b/SRC/dormrq.f @@ -267,7 +267,7 @@ SUBROUTINE DORMRQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'DORMRQ', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/dormrz.f b/SRC/dormrz.f index 3ec12aaa18..0dea1b5ff1 100644 --- a/SRC/dormrz.f +++ b/SRC/dormrz.f @@ -291,7 +291,7 @@ SUBROUTINE DORMRZ( SIDE, TRANS, M, N, K, L, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'DORMRQ', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/sormlq.f b/SRC/sormlq.f index 679a3a28e4..ee996e5602 100644 --- a/SRC/sormlq.f +++ b/SRC/sormlq.f @@ -266,7 +266,7 @@ SUBROUTINE SORMLQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'SORMLQ', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/sormql.f b/SRC/sormql.f index 56df279fab..72a8d22eea 100644 --- a/SRC/sormql.f +++ b/SRC/sormql.f @@ -268,7 +268,7 @@ SUBROUTINE SORMQL( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'SORMQL', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/sormqr.f b/SRC/sormqr.f index 828ab7d426..5d4256f09c 100644 --- a/SRC/sormqr.f +++ b/SRC/sormqr.f @@ -265,7 +265,7 @@ SUBROUTINE SORMQR( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'SORMQR', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/sormrq.f b/SRC/sormrq.f index 602fcb9633..62fcdacdb9 100644 --- a/SRC/sormrq.f +++ b/SRC/sormrq.f @@ -269,7 +269,7 @@ SUBROUTINE SORMRQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'SORMRQ', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/sormrz.f b/SRC/sormrz.f index 522cb2fce7..cdadd62b51 100644 --- a/SRC/sormrz.f +++ b/SRC/sormrz.f @@ -290,7 +290,7 @@ SUBROUTINE SORMRZ( SIDE, TRANS, M, N, K, L, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'SORMRQ', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/zunmlq.f b/SRC/zunmlq.f index c67ee50e06..f2fb315440 100644 --- a/SRC/zunmlq.f +++ b/SRC/zunmlq.f @@ -264,7 +264,7 @@ SUBROUTINE ZUNMLQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'ZUNMLQ', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/zunmql.f b/SRC/zunmql.f index a5e12f757e..06353a0c75 100644 --- a/SRC/zunmql.f +++ b/SRC/zunmql.f @@ -266,7 +266,7 @@ SUBROUTINE ZUNMQL( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'ZUNMQL', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/zunmqr.f b/SRC/zunmqr.f index 48fc7e2ece..2ae205f4fd 100644 --- a/SRC/zunmqr.f +++ b/SRC/zunmqr.f @@ -263,7 +263,7 @@ SUBROUTINE ZUNMQR( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'ZUNMQR', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/zunmrq.f b/SRC/zunmrq.f index ff61923c89..e14e93bf2e 100644 --- a/SRC/zunmrq.f +++ b/SRC/zunmrq.f @@ -267,7 +267,7 @@ SUBROUTINE ZUNMRQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'ZUNMRQ', SIDE // TRANS, M, N, K, $ -1 ) ) diff --git a/SRC/zunmrz.f b/SRC/zunmrz.f index 0410851907..1ec5532846 100644 --- a/SRC/zunmrz.f +++ b/SRC/zunmrz.f @@ -295,7 +295,7 @@ SUBROUTINE ZUNMRZ( SIDE, TRANS, M, N, K, L, A, LDA, TAU, C, LDC, NBMIN = 2 LDWORK = NW IF( NB.GT.1 .AND. NB.LT.K ) THEN - IF( LWORK.LT.NW*NB+TSIZE ) THEN + IF( LWORK.LT.LWKOPT ) THEN NB = (LWORK-TSIZE) / LDWORK NBMIN = MAX( 2, ILAENV( 2, 'ZUNMRQ', SIDE // TRANS, M, N, K, $ -1 ) )