In Matlab, there is pageeig function that can solve eigenvalues and eigenvectors for a huge batch of matrices efficiently. In NumPy, the linalg.eig can also accept arrays with high dimension (..., M, M) and diagonalize them in batch. I know that Parallel functions, such as ParallelMap, can accelerate the evaluation of a batch of matrices, but they are less efficient. I wonder whether there are packages for page-wise Eigensystem in Mathematica.
