Hi everyone,

I am attempting to sample a few shots through the two observables

M_1 = (\sigma_x \otimes \sigma_z) and M_2 = (\sigma_z \otimes \sigma_x) simultaneously with two qubits.

The simplest circuit I can come up with is this one I believe:

```
device = qml.device("default.qubit", wires=2, shots=100)
@qml.qnode(device)
def circuit():
obs = [qml.PauliX(wires=0) @ qml.PauliZ(wires=1), qml.PauliZ(wires=0) @ qml.PauliX(wires=1)]
return list(map(qml.sample, obs))
print(qml.draw(circuit)())
```

Now on version 0.19.1 I get the following error:

```
QuantumFunctionError: Only observables that are qubit-wise commuting Pauli words can be returned on the same wire
```

I understand the commuting requirement of M_1 and M_2, and I believe this is checked since M_1 M_2 = (\sigma_x \sigma_z) \otimes (\sigma_z \sigma_x) = (-\sigma_z \sigma_x) \otimes (-\sigma_x \sigma_z) = M_2 M_1

But a similar error seems to be raised even with the following function:

```
qml.grouping.diagonalize_qwc_pauli_words(obs)
```

So I believe I am not supposed to use the framework in this way for the result I intend to obtain. Could you help me figure out what I should do?

Thank you!