Add "S_inv" keyword to propagator_operator

If S_inv is set to an inverse operator, it should be use instead of S.inverse_multiply/times in

D^-1 = S_inv + M

Usecase: Sometimes only a non-invertable smoothness enforcing operator S_inv \propto k^2 or k^4 should be used, or an explicitly coded pixel space operator.