def phase(self, mu, phi): """Non-scattering phase matrix. Returns : null phase matrix """ npol = 2 return np.zeroes((npol * len(mu), npol * len(mu)))