Currently running uniform shuffling (q = 0). The q > 0 deformation will modify the hole-filling step to introduce correlations based on partition-valued stopping probabilities. For a working q-Whittaker implementation, see q-RSK Sampling of Domino Tilings.
Forward EKLP Shuffling builds random tilings of Aztec diamonds via the map An → An+1:
Enable "Granular steps" to see each phase separately. Bad blocks are highlighted in red before deletion.
The q-Whittaker deformation of domino shuffling modifies Step 3 (Fill holes) to introduce correlations between adjacent empty blocks based on partition-valued stopping probabilities.
For a working implementation using the partition-based q-Whittaker algorithm (via RSK growth diagrams), see:
The challenge for direct shuffling is translating partition indices (where islands are detected as consecutive i with μi − κi = 1) into geometric empty-block positions.
References:
Link to code
(Pure JavaScript implementation of EKLP shuffling)