I'm looking to implement spatially correlated large-area errors similar to this paper, where a local patch of qubits experience random $I$,$X$, $Y$, $Z$ errors (simultaneously) with some probability $p(r)$, where $r$ is related to the size of the patch.
I can employ the CORRELATED_ERROR
and ELSE_CORRELATED_ERROR
instructions to implement this, but because I have to specify each Pauli product, this seemingly becomes intractable as the number of Pauli products grows as $4^N$ for an $N$-qubit error.
Is there some functionality of Stim or a workaround that I'm overlooking that would make this implementable? Or would I have to resort stim.FlipSimulator
or stim.TableauSimulator
?