def getSbus(self, buses=None):
        """ Returns the net complex bus power injection vector in p.u.
        """
        bs = self.buses if buses is None else buses
        s = array([self.s_surplus(v) / self.base_mva for v in bs])
        return s
Returns the net complex bus power injection vector in p.u.