ToolPythonAdvanced 2D turbulence with a pseudo-spectral solver: eddies merge into vortices K2inv = np.divide(1, K2, out=np.zeros_like(K2), where=K2 > 0)
ToolPythonIntermediate Cahn-Hilliard with a semi-implicit spectral step: how fast domains grow kinv = np.divide(1.0, kabs, out=np.zeros_like(kabs), where=kabs > 0)