Open gthyagi opened 5 months ago
I think that would work - we probably need to wrap it as we did with the volume integrals.
We could also fake it with an approximate delta function at the surface (we did this before for boundary conditions). That would be a workaround until the wrapper is implemented.
I will try the suggestions
@jcgraciosa - can we assign this to you ?
@lmoresi yes, please do.
knepley's comment for this feature Can you use https://petsc.org/main/manualpages/DMPlex/DMPlexComputeBdIntegral/ ?