I think it actually solves a generalised eigenvalue problem:
though for Fourier it might reduce to a standard eigenvalue problem. Maybe the key is pruneeigs which checks the eigenvectors for convergence?
I think it actually solves a generalised eigenvalue problem:
though for Fourier it might reduce to a standard eigenvalue problem. Maybe the key is pruneeigs which checks the eigenvectors for convergence?