Skip to content

Commit 1727eff

Browse files
iyamazakiiyamazaki
committed
add a new entry for NB of two-stage Aasen's
1 parent 5ffadec commit 1727eff

9 files changed

Lines changed: 31 additions & 21 deletions

‎SRC/chesv_aa_2stage.f‎

Lines changed: 1 addition & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -208,8 +208,7 @@ SUBROUTINE CHESV_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
208208
* ..
209209
* .. External Functions ..
210210
LOGICAL LSAME
211-
INTEGER ILAENV
212-
EXTERNAL LSAME, ILAENV
211+
EXTERNAL LSAME
213212
* ..
214213
* .. External Subroutines ..
215214
EXTERNAL XERBLA

‎SRC/chetrf_aa_2stage.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -414,7 +414,7 @@ SUBROUTINE CHETRF_AA_2STAGE( UPLO, N, A, LDA, TB, LTB, IPIV,
414414
*
415415
KB = MIN(NB, N-(J+1)*NB)
416416
CALL CLASET( 'Full', KB, NB, ZERO, ZERO,
417-
$ TB( TD+NB+1 + (J*NB)*LDTB) , LDTB-1 )
417+
$ TB( TD+NB+1 + (J*NB)*LDTB), LDTB-1 )
418418
CALL CLACPY( 'Upper', KB, NB,
419419
$ WORK, N,
420420
$ TB( TD+NB+1 + (J*NB)*LDTB ), LDTB-1 )
@@ -589,7 +589,7 @@ SUBROUTINE CHETRF_AA_2STAGE( UPLO, N, A, LDA, TB, LTB, IPIV,
589589
*
590590
KB = MIN(NB, N-(J+1)*NB)
591591
CALL CLASET( 'Full', KB, NB, ZERO, ZERO,
592-
$ TB( TD+NB+1 + (J*NB)*LDTB) , LDTB-1 )
592+
$ TB( TD+NB+1 + (J*NB)*LDTB), LDTB-1 )
593593
CALL CLACPY( 'Upper', KB, NB,
594594
$ A( (J+1)*NB+1, J*NB+1 ), LDA,
595595
$ TB( TD+NB+1 + (J*NB)*LDTB ), LDTB-1 )

‎SRC/csysv_aa_2stage.f‎

Lines changed: 1 addition & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -208,8 +208,7 @@ SUBROUTINE CSYSV_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
208208
* ..
209209
* .. External Functions ..
210210
LOGICAL LSAME
211-
INTEGER ILAENV
212-
EXTERNAL LSAME, ILAENV
211+
EXTERNAL LSAME
213212
* ..
214213
* .. External Subroutines ..
215214
EXTERNAL XERBLA

‎SRC/csytrf_aa_2stage.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -409,7 +409,7 @@ SUBROUTINE CSYTRF_AA_2STAGE( UPLO, N, A, LDA, TB, LTB, IPIV,
409409
*
410410
KB = MIN(NB, N-(J+1)*NB)
411411
CALL CLASET( 'Full', KB, NB, CZERO, CZERO,
412-
$ TB( TD+NB+1 + (J*NB)*LDTB) , LDTB-1 )
412+
$ TB( TD+NB+1 + (J*NB)*LDTB), LDTB-1 )
413413
CALL CLACPY( 'Upper', KB, NB,
414414
$ WORK, N,
415415
$ TB( TD+NB+1 + (J*NB)*LDTB ), LDTB-1 )
@@ -595,7 +595,7 @@ SUBROUTINE CSYTRF_AA_2STAGE( UPLO, N, A, LDA, TB, LTB, IPIV,
595595
*
596596
KB = MIN(NB, N-(J+1)*NB)
597597
CALL CLASET( 'Full', KB, NB, CZERO, CZERO,
598-
$ TB( TD+NB+1 + (J*NB)*LDTB) , LDTB-1 )
598+
$ TB( TD+NB+1 + (J*NB)*LDTB), LDTB-1 )
599599
CALL CLACPY( 'Upper', KB, NB,
600600
$ A( (J+1)*NB+1, J*NB+1 ), LDA,
601601
$ TB( TD+NB+1 + (J*NB)*LDTB ), LDTB-1 )

‎SRC/dsysv_aa_2stage.f‎

Lines changed: 1 addition & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -210,8 +210,7 @@ SUBROUTINE DSYSV_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
210210
* ..
211211
* .. External Functions ..
212212
LOGICAL LSAME
213-
INTEGER ILAENV
214-
EXTERNAL LSAME, ILAENV
213+
EXTERNAL LSAME
215214
* ..
216215
* .. External Subroutines ..
217216
EXTERNAL XERBLA

‎SRC/dsytrf_aa_2stage.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -402,7 +402,7 @@ SUBROUTINE DSYTRF_AA_2STAGE( UPLO, N, A, LDA, TB, LTB, IPIV,
402402
*
403403
KB = MIN(NB, N-(J+1)*NB)
404404
CALL DLASET( 'Full', KB, NB, ZERO, ZERO,
405-
$ TB( TD+NB+1 + (J*NB)*LDTB) , LDTB-1 )
405+
$ TB( TD+NB+1 + (J*NB)*LDTB), LDTB-1 )
406406
CALL DLACPY( 'Upper', KB, NB,
407407
$ WORK, N,
408408
$ TB( TD+NB+1 + (J*NB)*LDTB ), LDTB-1 )
@@ -573,7 +573,7 @@ SUBROUTINE DSYTRF_AA_2STAGE( UPLO, N, A, LDA, TB, LTB, IPIV,
573573
*
574574
KB = MIN(NB, N-(J+1)*NB)
575575
CALL DLASET( 'Full', KB, NB, ZERO, ZERO,
576-
$ TB( TD+NB+1 + (J*NB)*LDTB) , LDTB-1 )
576+
$ TB( TD+NB+1 + (J*NB)*LDTB), LDTB-1 )
577577
CALL DLACPY( 'Upper', KB, NB,
578578
$ A( (J+1)*NB+1, J*NB+1 ), LDA,
579579
$ TB( TD+NB+1 + (J*NB)*LDTB ), LDTB-1 )

‎SRC/ilaenv.f‎

Lines changed: 19 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -176,8 +176,8 @@ INTEGER FUNCTION ILAENV( ISPEC, NAME, OPTS, N1, N2, N3, N4 )
176176
*
177177
* .. Local Scalars ..
178178
INTEGER I, IC, IZ, NB, NBMIN, NX
179-
LOGICAL CNAME, SNAME
180-
CHARACTER C1*1, C2*2, C4*2, C3*3, SUBNAM*6
179+
LOGICAL CNAME, SNAME, TWOSTAGE
180+
CHARACTER C1*1, C2*2, C4*2, C3*3, SUBNAM*16
181181
* ..
182182
* .. Intrinsic Functions ..
183183
INTRINSIC CHAR, ICHAR, INT, MIN, REAL
@@ -256,6 +256,8 @@ INTEGER FUNCTION ILAENV( ISPEC, NAME, OPTS, N1, N2, N3, N4 )
256256
C2 = SUBNAM( 2: 3 )
257257
C3 = SUBNAM( 4: 6 )
258258
C4 = C3( 2: 3 )
259+
TWOSTAGE = LEN( SUBNAM ).GE.11
260+
$ .AND. SUBNAM( 11: 11 ).EQ.'2'
259261
*
260262
GO TO ( 50, 60, 70 )ISPEC
261263
*
@@ -359,9 +361,17 @@ INTEGER FUNCTION ILAENV( ISPEC, NAME, OPTS, N1, N2, N3, N4 )
359361
ELSE IF( C2.EQ.'SY' ) THEN
360362
IF( C3.EQ.'TRF' ) THEN
361363
IF( SNAME ) THEN
362-
NB = 64
364+
IF( TWOSTAGE ) THEN
365+
NB = 192
366+
ELSE
367+
NB = 64
368+
END IF
363369
ELSE
364-
NB = 64
370+
IF( TWOSTAGE ) THEN
371+
NB = 192
372+
ELSE
373+
NB = 64
374+
END IF
365375
END IF
366376
ELSE IF( SNAME .AND. C3.EQ.'TRD' ) THEN
367377
NB = 32
@@ -370,7 +380,11 @@ INTEGER FUNCTION ILAENV( ISPEC, NAME, OPTS, N1, N2, N3, N4 )
370380
END IF
371381
ELSE IF( CNAME .AND. C2.EQ.'HE' ) THEN
372382
IF( C3.EQ.'TRF' ) THEN
373-
NB = 64
383+
IF( TWOSTAGE ) THEN
384+
NB = 192
385+
ELSE
386+
NB = 64
387+
END IF
374388
ELSE IF( C3.EQ.'TRD' ) THEN
375389
NB = 32
376390
ELSE IF( C3.EQ.'GST' ) THEN

‎SRC/ssysv_aa_2stage.f‎

Lines changed: 1 addition & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -208,8 +208,7 @@ SUBROUTINE SSYSV_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
208208
* ..
209209
* .. External Functions ..
210210
LOGICAL LSAME
211-
INTEGER ILAENV
212-
EXTERNAL LSAME, ILAENV
211+
EXTERNAL LSAME
213212
* ..
214213
* .. External Subroutines ..
215214
EXTERNAL XERBLA

‎SRC/ssytrf_aa_2stage.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -402,7 +402,7 @@ SUBROUTINE SSYTRF_AA_2STAGE( UPLO, N, A, LDA, TB, LTB, IPIV,
402402
*
403403
KB = MIN(NB, N-(J+1)*NB)
404404
CALL SLASET( 'Full', KB, NB, ZERO, ZERO,
405-
$ TB( TD+NB+1 + (J*NB)*LDTB) , LDTB-1 )
405+
$ TB( TD+NB+1 + (J*NB)*LDTB), LDTB-1 )
406406
CALL SLACPY( 'Upper', KB, NB,
407407
$ WORK, N,
408408
$ TB( TD+NB+1 + (J*NB)*LDTB ), LDTB-1 )
@@ -573,7 +573,7 @@ SUBROUTINE SSYTRF_AA_2STAGE( UPLO, N, A, LDA, TB, LTB, IPIV,
573573
*
574574
KB = MIN(NB, N-(J+1)*NB)
575575
CALL SLASET( 'Full', KB, NB, ZERO, ZERO,
576-
$ TB( TD+NB+1 + (J*NB)*LDTB) , LDTB-1 )
576+
$ TB( TD+NB+1 + (J*NB)*LDTB), LDTB-1 )
577577
CALL SLACPY( 'Upper', KB, NB,
578578
$ A( (J+1)*NB+1, J*NB+1 ), LDA,
579579
$ TB( TD+NB+1 + (J*NB)*LDTB ), LDTB-1 )

0 commit comments

Comments
 (0)