  We extend thresholding methods for numerical realization of mean curvature
flow on obstacles to the anisotropic setting where interfacial energy depends
on the orientation of the interface. This type of schemes treats the interface
implicitly, which supports natural implementation of topology changes, such as
merging and splitting, and makes the approach attractive for applications in
material science. The main tool in the new scheme are convolution kernels
developed in previous studies, that approximate the given anisotropy in a
nonlocal way. We provide a detailed report on the numerical properties of the
proposed algorithm.
