If the matricies are very large and sparse, you may be better off using a sparse format and using something like power method. This may also be a good job for a little inline C.