This work presents efficient matrix-free algorithms for evaluating the operator action in unfitted finite element discretizations, enabling high-performance computations with high-order polynomial spaces.
