Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Purging: remove converged, unwanted vecs #140

Merged

Conversation

haampie
Copy link
Member

@haampie haampie commented Feb 18, 2024

It looks like converged, unwanted vectors can actually be problematic for the
Krylov-Schur restart, since they have tiny residuals, and we need to construct
reflectors out of those values.

By reordering the converged, unwanted to the tail of the decomposition, they're
removed, and we don't have to compute reflectors involving their residuals.

Copy link

codecov bot commented Feb 18, 2024

Codecov Report

All modified and coverable lines are covered by tests ✅

Comparison is base (2a264db) 95.47% compared to head (9fd7e67) 95.52%.

Additional details and impacted files
@@            Coverage Diff             @@
##           master     #140      +/-   ##
==========================================
+ Coverage   95.47%   95.52%   +0.05%     
==========================================
  Files          10       10              
  Lines         839      849      +10     
==========================================
+ Hits          801      811      +10     
  Misses         38       38              

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

@haampie haampie mentioned this pull request Feb 18, 2024
@haampie haampie merged commit 55a09aa into JuliaLinearAlgebra:master Feb 18, 2024
13 checks passed
@haampie haampie deleted the fix/purge-converged-unwanted branch February 18, 2024 23:05
haampie added a commit that referenced this pull request Feb 18, 2024
It looks like converged, unwanted vectors can actually be problematic for the
Krylov-Schur restart, since they have tiny residuals, and we need to construct
reflectors out of those values.

By reordering the converged, unwanted to the tail of the decomposition, they're
removed, and we don't have to compute reflectors involving their residuals.

Should also restore the original pre-0.3 behavior of purging.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant