fix bug in mumps-assisted sparse diagonalization

When using the shift-invert method the *unshifted* matrix should
be passed to 'eigsh', but the inverse involving the *shifted* matrix
should be passed as the 'OPinv' parameter.

This did not result in a bug in the bound state solver, as we are
always trying to find eigenvalues close to 0.