phynexis-deps/lapack
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Merge pull request #438 from h-vetinari/backport_openblas_2778

Browse Source 
Fix some incorrect SLASET/DLASET calls (includes backport of openblas#2778)
This commit is contained in:
julielangou
2020-09-07 18:21:04 -07:00
committed by GitHub
parent 9b136e860a bffac66e45
commit cc6e21e3c6
8 changed files with 32 additions and 32 deletions
Download Patch File Download Diff File Expand all files Collapse all files
TESTING/EIG/cchkhb2stg.f
+4 -4
View File
@@ -680,8 +680,8 @@
* the one from above. Compare it with D1 computed
* using the DSBTRD.
*
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL CLACPY( ' ', K+1, N, A, LDA, U, LDU )
LH = MAX(1, 4*N)
LW = LWORK - LH
@@ -753,8 +753,8 @@
* the one from above. Compare it with D1 computed
* using the DSBTRD.
*
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL CLACPY( ' ', K+1, N, A, LDA, U, LDU )
LH = MAX(1, 4*N)
LW = LWORK - LH
TESTING/EIG/cchkst2stg.f
+4 -4
View File
@@ -1014,8 +1014,8 @@
* the one from above. Compare it with D1 computed
* using the 1-stage.
*
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL CLACPY( 'U', N, N, A, LDA, V, LDU )
LH = MAX(1, 4*N)
LW = LWORK - LH
@@ -1048,8 +1048,8 @@
* the one from above. Compare it with D1 computed
* using the 1-stage.
*
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL CLACPY( 'L', N, N, A, LDA, V, LDU )
CALL CHETRD_2STAGE( 'N', "L", N, V, LDU, SD, SE, TAU,
$ WORK, LH, WORK( LH+1 ), LW, IINFO )
TESTING/EIG/dchksb2stg.f
+4 -4
View File
@@ -670,8 +670,8 @@
* the one from above. Compare it with D1 computed
* using the DSBTRD.
*
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL DLACPY( ' ', K+1, N, A, LDA, U, LDU )
LH = MAX(1, 4*N)
LW = LWORK - LH
@@ -743,8 +743,8 @@
* the one from above. Compare it with D1 computed
* using the DSBTRD.
*
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL DLACPY( ' ', K+1, N, A, LDA, U, LDU )
LH = MAX(1, 4*N)
LW = LWORK - LH
TESTING/EIG/dchkst2stg.f
+4 -4
View File
@@ -999,8 +999,8 @@
* the one from above. Compare it with D1 computed
* using the 1-stage.
*
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL DLACPY( "U", N, N, A, LDA, V, LDU )
LH = MAX(1, 4*N)
LW = LWORK - LH
@@ -1032,8 +1032,8 @@
* the one from above. Compare it with D1 computed
* using the 1-stage.
*
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL DLACPY( "L", N, N, A, LDA, V, LDU )
CALL DSYTRD_2STAGE( 'N', "L", N, V, LDU, SD, SE, TAU,
$ WORK, LH, WORK( LH+1 ), LW, IINFO )
TESTING/EIG/schksb2stg.f
+4 -4
View File
@@ -670,8 +670,8 @@
* the one from above. Compare it with D1 computed
* using the SSBTRD.
*
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL SLACPY( ' ', K+1, N, A, LDA, U, LDU )
LH = MAX(1, 4*N)
LW = LWORK - LH
@@ -743,8 +743,8 @@
* the one from above. Compare it with D1 computed
* using the SSBTRD.
*
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL SLACPY( ' ', K+1, N, A, LDA, U, LDU )
LH = MAX(1, 4*N)
LW = LWORK - LH
TESTING/EIG/schkst2stg.f
+4 -4
View File
@@ -999,8 +999,8 @@
* the one from above. Compare it with D1 computed
* using the 1-stage.
*
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL SLACPY( "U", N, N, A, LDA, V, LDU )
LH = MAX(1, 4*N)
LW = LWORK - LH
@@ -1032,8 +1032,8 @@
* the one from above. Compare it with D1 computed
* using the 1-stage.
*
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL SLACPY( "L", N, N, A, LDA, V, LDU )
CALL SSYTRD_2STAGE( 'N', "L", N, V, LDU, SD, SE, TAU,
$ WORK, LH, WORK( LH+1 ), LW, IINFO )
TESTING/EIG/zchkhb2stg.f
+4 -4
View File
@@ -680,8 +680,8 @@
* the one from above. Compare it with D1 computed
* using the DSBTRD.
*
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL ZLACPY( ' ', K+1, N, A, LDA, U, LDU )
LH = MAX(1, 4*N)
LW = LWORK - LH
@@ -753,8 +753,8 @@
* the one from above. Compare it with D1 computed
* using the DSBTRD.
*
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL ZLACPY( ' ', K+1, N, A, LDA, U, LDU )
LH = MAX(1, 4*N)
LW = LWORK - LH
TESTING/EIG/zchkst2stg.f
+4 -4
View File
@@ -1014,8 +1014,8 @@
* the one from above. Compare it with D1 computed
* using the 1-stage.
*
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL ZLACPY( 'U', N, N, A, LDA, V, LDU )
LH = MAX(1, 4*N)
LW = LWORK - LH
@@ -1048,8 +1048,8 @@
* the one from above. Compare it with D1 computed
* using the 1-stage.
*
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SD, N )
CALL DLASET( 'Full', N, 1, ZERO, ZERO, SE, N )
CALL ZLACPY( 'L', N, N, A, LDA, V, LDU )
CALL ZHETRD_2STAGE( 'N', "L", N, V, LDU, SD, SE, TAU,
$ WORK, LH, WORK( LH+1 ), LW, IINFO )