in src/braket/default_simulator/noise_operations.py [0:0]
def matrices(self) -> List[np.ndarray]:
SI = np.array([[1.0, 0.0], [0.0, 1.0]], dtype=complex)
SX = np.array([[0.0, 1.0], [1.0, 0.0]], dtype=complex)
SY = np.array([[0.0, -1.0j], [1.0j, 0.0]], dtype=complex)
SZ = np.array([[1.0, 0.0], [0.0, -1.0]], dtype=complex)
K_list_single = [SI, SX, SY, SZ]
K_list = [np.kron(i, j) for i in K_list_single for j in K_list_single]
K_list[0] *= np.sqrt(1 - self._probability)
K_list[1:] = [np.sqrt(self._probability / 15) * i for i in K_list[1:]]
return K_list