% fun_laplacian.m — The surface Laplacian (Laplace-Beltrami operator).
%
% Higher-order operators are built by composing the basic ones. lap(f)
% is the surface Laplacian of f.

mip load --install flatironinstitute/flatironinstitute/surfacefun

p = 16;
nref = 2;
dom = surfacemesh.sphere(p + 1, nref);

f = surfacefun(@(x, y, z) cos(6*x) .* y + exp(z), dom);

figure(1);
plot(lap(f)), colorbar
title('\Delta_\Gamma f');