Dear Egor Vtorushin, Have a look at the outer_faces method: [image: image.png] if you put only one element in CVIDs you will get what you want, I guess.
I am not sure how efficient it is to calculate these integrals for all elements one by one though. "mesh.all_faces" might also be relevant for your question. You can also have a look at HHO examples like https://git.savannah.nongnu.org/cgit/getfem.git/tree/interface/tests/python/demo_elasticity_HHO.py Best regards Kostas On Mon, Jul 19, 2021 at 8:59 AM Egor Vtorushin <vtorus...@gmail.com> wrote: > Dear Yves, > Do you have any hints how to manage 1D integral over a loop containing > simplex' or quad' *edges*. > Let I have an electrical field(or current ) computed in simplex or quad > mesh nodes. Then i need to calculate the field circulation over edge-loop > for each *face *according to Stokes' theorem/formula > There is an obsoleted function > edges(*CVLST=None*, **args*) > <http://getfem.org/python/cmdref_Mesh.html#getfem.Mesh.edges> > > Synopsis: [E,C] = Mesh.edges(self [, CVLST][, ‘merge’]) > > [OBSOLETE FUNCTION! will be removed in a future release] > > Return the list of edges of mesh M for the convexes.... > Maybe there is a way to utilize the Mesh.edges function or another way? > Regards, Egor Vtorushin > >