Fix incorrect calls to DLASET
Reference-LAPACK issue 429
This commit is contained in:
@@ -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, 1 )
|
||||
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
|
||||
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, 1 )
|
||||
CALL SLASET( 'Full', N, 1, ZERO, ZERO, SE, 1 )
|
||||
CALL CLACPY( ' ', K+1, N, A, LDA, U, LDU )
|
||||
LH = MAX(1, 4*N)
|
||||
LW = LWORK - LH
|
||||
|
||||
Reference in New Issue
Block a user