Skip to content

stedc scales only if norm is tiny or huge (and not unconditionally)#1053

Draft
angsch wants to merge 1 commit intoReference-LAPACK:masterfrom
angsch:stedc
Draft

stedc scales only if norm is tiny or huge (and not unconditionally)#1053
angsch wants to merge 1 commit intoReference-LAPACK:masterfrom
angsch:stedc

Conversation

@angsch
Copy link
Collaborator

@angsch angsch commented Sep 18, 2024

Description

Replace unconditional scaling

lapack/SRC/dstedc.f

Lines 380 to 387 in 2cbabee

ORGNRM = DLANST( 'M', M, D( START ), E( START ) )
CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, M, 1, D( START ),
$ M,
$ INFO )
CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, M-1, 1,
$ E( START ),
$ M-1, INFO )
*

with conditional scaling as done in, for example, steqr.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

2 participants