getfem-users
[Top][All Lists]

## Re: (E field | current) circulation

 From: Konstantinos Poulios Subject: Re: (E field | current) circulation Date: Mon, 19 Jul 2021 10:31:22 +0200

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