How to measure pointwise error?


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?


Take a look at this answer:

With that “feature” you can at least evaluate the maximum over the set of integration points.