1 2 3 4 5 6
void streamfunction (vector u, scalar psi) { scalar omega[]; vorticity (u, omega); boundary ({psi, omega}); poisson (psi, omega); }