I’m using a cut finite element method for solving Stokes equations. I want to investigate the pointwise divergence of the velocity solution. In another library, I do this by iterating over all quadrature points (both on cut elements, noncut elements, and on element faces) and taking the maximum value of the divergence of the solution. Can something like this be done in ngsxfem?