Added SparseShiftMatrix routine