To update the Krylov subspace representation of a matrix A
after adding d*A(:,k)
to the k-th column of A
, you can update the Hessenberg
matrix H
and the matrix of basis vectors V
without having to directly store or work with the matrix A
.
You can achieve this by simply updating the Hessenberg matrix H
, and the basis vectors matrix V
using the following operations:
main.m268 chars6 lines
Make sure to update the corresponding k-th column of the matrix H
by adding the appropriate value as well.
gistlibby LogSnag