lengthy discussion/introduction to Kryolov subspaces and eigensolvers including QR and Arnoldi, implemented (apparently correctly, though error checking is still somehwat haphazard) in arnoldi_method.py for hermitian input matrices