JuliaLinearAlgebra / ArnoldiMethod.jl

The Arnoldi Method with Krylov-Schur restart, natively in Julia.
https://julialinearalgebra.github.io/ArnoldiMethod.jl/dev
MIT License
96 stars 18 forks source link

Allow passing an initial partial Schur decomposition #143

Closed haampie closed 7 months ago

haampie commented 7 months ago

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; ...).

codecov[bot] commented 7 months ago

Codecov Report

All modified and coverable lines are covered by tests :white_check_mark:

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

Additional details and impacted files ```diff @@ 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 ```

:umbrella: View full report in Codecov by Sentry.
:loudspeaker: Have feedback on the report? Share it here.