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

Allow passing an initial partial Schur decomposition #143

Merged

Conversation

haampie
Copy link
Member

@haampie haampie commented Feb 22, 2024

This commit allows users to

  1. avoid allocations by reusing existing arrays/buffers
  2. specify custom matrix types for V, H and temporaries
  3. starting the algorithm from an existing partial Schur decomposition

by exporting ArnoldiWorkspace, which holds the major arrays, and
partialschur!(A, ::ArnoldiWorkspace; ...).

Copy link

codecov bot commented Feb 22, 2024

Codecov Report

All modified and coverable lines are covered by tests ✅

Comparison is base (3cb957d) 95.91% compared to head (0be7991) 96.75%.

Additional details and impacted files
@@            Coverage Diff             @@
##           master     #143      +/-   ##
==========================================
+ Coverage   95.91%   96.75%   +0.83%     
==========================================
  Files          10       10              
  Lines         856      800      -56     
==========================================
- Hits          821      774      -47     
+ Misses         35       26       -9     

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

@haampie haampie merged commit e7141e6 into JuliaLinearAlgebra:master Feb 22, 2024
13 checks passed
@haampie haampie deleted the fix/pre-allocate-workspace branch February 22, 2024 09:59
haampie added a commit that referenced this pull request Feb 22, 2024
This commit allows users to

1. avoid allocations by reusing existing arrays/buffers
2. specify custom matrix types for V, H and temporaries
3. starting the algorithm from an existing partial Schur decomposition

by exporting `ArnoldiWorkspace`, which holds the major arrays, and
`partialschur!(A, ::ArnoldiWorkspace; ...)`.

Co-authored-by: Samuel Powell <[email protected]>
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