Seems like GE on a 2x2 matrix should be pretty darn fast.... especially if you use scalar * vector type operations.

Refined estimates for finite element or, more generally, Galerkin approximations of the eigenvalues and eigenvectors of selfadjoint eigenvalue problems are presented. More specifically, refined ...