Wrappers for basis orthogonalization/normalization

nekStab / LightKrylov

Lightweight implementation of Krylov subspace techniques in Fortran.

BSD 3-Clause "New" or "Revised" License

18 stars 0 forks source link

Closed Simkern closed 3 weeks ago

Simkern commented 3 weeks ago

Added and deployed orthogonalization/normalization wrappers throughout the code.

New definitions:

orthogonalize_against_basis
orthonormalize_basis

Multiple dispatch wrappers for norml and log2. innerprod, linear_combination, axpby_basis and copy_basis exported directly in LightKrylov.
Switch from eigenvalue- to svd-based algorithm for sqrtm. Alternative kept for now.