Generic surface integral implementation
I believe we currently don't have the functionality to perform generic surface integrals on the domain, the idea is to mark the boundaries of a macro primitive and perform surface integrals if needed on them. The marker can be the existing NeumannBoundary
markers or some more generic markers for surfaces inside the domain. The forms needs to be modified a bit to accommodate the cell coordinates and the facet coordinated where the surface integral needs to be evaluated.
Edited by Ponsuganth Ilangovan Ponkumar Ilango