This is a direct port from the C code presented in Stam's 2003 paper, with a small convenience wrapper.

Example:

const {Fluid} = require('stam-stable-fluids')
const size = 128
const f = new Fluid(size)
for (let y = 0; y < size; y++)
  for (let x = 0; x < size/2; x++)
    f.addDensity(x, y, 10)
for (let y = 0; y < size; y++)
  f.addForce(10, y, 10, 0)
f.step(1/60)

const [vx, vy] = f.velocityAt(0, 0)
const d = f.densityAt(0, 0)

Make a new fluid simulation with width and height equal to size. If you want a goopier fluid, like honey, set viscosity to a small positive number, like 0.0001 (larger numbers will probably make the fluid so "thick" it will barely move). diffusion will make the "stuff" in the fluid spread out over time. The default for both viscosity and diffusion, if not provided, is 0.

Step the simulation by dt (e.g. 1/60).

Returns [vx, vy] representing the velocity of the fluid at (x, y).

Returns the density of "stuff" at (x, y).

Adds an external force to the fluid, e.g. from a UI or other system. The force field will be cleared after every step().

Adds an external density flow to the fluid. d is the rate of flow. The density flow field will be cleared after every step().