Work with eigenvalue/eigenvector solver?

[duetosymmetry]

Since this provides a variety of factorizations, I thought it might be possible to use this package to solve an infinite-dimensional eigenvalue/eigenvector problem. For example, can ArnoldiMethod.jl be made to use an infinite-dimensional factorization to find the closest eigenvalue/vector pair to some initial guess (via shift-and-invert)? I was not able to make it work out of the box, but I'm also not fluent in Julia.