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