decouple multiple mass (multishift) solver from IO routines and generalize

add temporary CGMMS io routine to invert_eo until spinor IO is refurbished and generalized