A modified ADI method for computing the “best least-squares solution” of an incompatible system (A×I+I×B)x=g