pytenet.krylov

Krylov subspace algorithms.

Functions

arnoldi_iteration(afunc, vstart, numiter)

Perform a "matrix free" Arnoldi iteration.

eigh_krylov(afunc, vstart, numiter, numeig)

Compute Krylov subspace approximation of eigenvalues and vectors.

eigh_tridiag(d, e)

Solve the eigenvalue problem for a real symmetric tridiagonal matrix.

expm_krylov(afunc, vec, dt, numiter[, hermitian])

Compute Krylov subspace approximation of the matrix exponential applied to input vector: expm(dt*A)*vec.

lanczos_iteration(afunc, vstart, numiter)

Perform a "matrix free" Lanczos iteration.