I’ve been able to implement a first version of the Arnoldi method for approximation of the largest eigenvalues of a large matrix. Looks good (and fast 
).
I’ll try to make a version of it working as eigen, and another version working as eig() in the math parser.
Stay tuned!