def METHOD_NAME(self):
    if isinstance(self.exponent, sympy.Basic):
        return NotImplemented
    size = 1 << self._num_qubits
    return np.diag([1j ** (4 * i / size * self.exponent) for i in range(size)])