1% vec_curl.m — The surface curl of a vector field.
2%
3% curl of a surfacefunv is again a surfacefunv.
5mip load --install flatironinstitute/flatironinstitute/surfacefun
7rng(0)
8p = 16;
9nref = 2;
10dom = surfacemesh.blob(p + 1, nref);
12v = surfacefunv(@(x, y, z) cos(2*x), ...
13 @(x, y, z) sin(4*y), ...
14 @(x, y, z) sin(3*z), dom);
16figure(1);
17quiver(curl(v), 0.1, 6)
18title('Surface curl of a vector field');