I’m trying to use eigs and I’m running into arparck errors. More specifically, when eigs doesn’t converge the function itself doesn’t throw or return some kind of flag (even though that is what the documentation seemed to suggest, it returns an nconved value after all), but the arparck wrapper throws an unspecified error (1).
This is cryptic and I was wondering if this could be improved a bit? A few months earlier I ran into a similar unspecified error (-1), and that was because my particular matrix had “Inf” values.
More importantly, in my particular use-case I’m solving the eigenvalue problem a number of times, everytime increasing the desired precision.
So really, it not converging at first is not that big of a problem, but I’d still like to retrieve the partially converged eigenvectors. Is there any way to do this? Ideally I’d get the partially converged result along with a bound on the error.