We numerically solve the Schrödinger-Helmholtz equation

in two dimensions, via a standard pseudo-spectral method, at resolution 128².

Freely-evolving simulation

Forced-dissipated simulation