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