An algorithm is presented for the simulation of reaction–diffusion systems on complex geometries, providing insight on how the interplay of cell geometry and biochemistry can control polarity in living cells.