Skip to content

Add scalar and block Lanczos matrix-function application - #191

Open
rileyjmurray wants to merge 1 commit into
split-138/symmetric-operatorsfrom
split-138/lanczos-fa
Open

rileyjmurray wants to merge 1 commit into
split-138/symmetric-operatorsfrom
split-138/lanczos-fa

Conversation

@rileyjmurray

@rileyjmurray rileyjmurray commented Oct 1, 2026 •

Copy link
Copy Markdown
Contributor

LanczosFA and BlockLanczosFA approximate f(A)B from matrix products with a symmetric operator. The scalar implementation processes right-hand sides independently; the block implementation builds a shared Krylov basis. Both support fixed-depth evaluation and reusable workspaces.

This extracts the matrix-function application components from #138, based on #190. It includes the shared workspace helpers and a dense eigendecomposition oracle for tests. Optional CGS2 reorthogonalization is in #195.

Validation. All 508 core Release tests passed locally with Spack Clang 19.1.3 on macOS. Component tests compare against independently evaluated dense matrix functions, including full-depth recovery. Driver integration remains in #138.

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