Consider i have a 4x4 matrix

```
H = array([[-2., 1., 0., 0.],
[ 1., -2., 1., 0.],
[ 0., 1., -2., 1.],
[ 0., 0., 1., -2.]]
```

i can decomose this matrix to unitaries with PauliX, Identity, PauliY and PauliZ with `qml.utils.decompose_hamiltonian`

. I want to decompose M with these matrices

```
C = Matrix([[0, 1],[0, 0]])
S = Matrix([[0, 0],[1, 0]])
```

where are non-unitary. These matrices are spin rising and spin lowering operators, too. How can i decompose M using Identity, C, and S with `qml.utils.decompose_hamiltonian`

?