The method used is Cellular Automata, or rather, a variation of it:

https://slackermanz.com/understanding-multiple-neighborhood-...

The simulation was done on a RX 6800XT, and ran for ~4.5 million GPU-frames, which took just over two hours at 1024x1024 resolution (~560 FPS)