Re: (E field | current) circulation

[Konstantinos Poulios]

Dear Egor Vtorushin,

Have a look at the outer_faces method:

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 <vtorushin@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)

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