don't subtract N_DEFLATED from istop twice #794
Merged
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Description
closes #793
During AED, eigenvalues of a subblock are calculated, some of these are deflated, others are stored to be used as shifts during the next multishift QZ iteration.
The shifts are stored in the eigenvalue array at an index corresponding to the top of the AED window. This equals
istop - nd - ns + 1
. However, after the AED call,istop
is set toistop - nd
to reflect the new active window. This means that the shifts are now stored at indexistop - ns + 1
.This issue was not detected before, because usually, when
nd
is large enough that a significant amount of wrong shifts would be used, we skip the sweep in favor of doing another AED step.Checklist