On the sparse matrix computation of critical frequencies