Kostas,
thank you for answering thу question and sorry for the delay in
communication. I am sorry for missing important point - i am in 3D space
In 2D HHO example it works fine since
# # Boundary selection
flst = m.outer_faces()
GAMMAD = 1
m.set_region(GAMMAD, flst)
managing with 1D faces that are the edges I suppose.
In 3D all_faces provides access to 2d faces. Consider 4 points simplex
that have 4 faces. I can integrate over an area of one of faces or many
using regions built with outer_faces output.
In my situation I have to make 1D integration over edges with respect to
length. The 4 point simplex has 4 faces and 6 1D edges.
At first glance outer_faces() method doesn't provide access to 1D
edges for 3D meshes and with is my problem, assuming the fact that it
returns empty array for dim=1 or dim=2
*boxm*.outer_faces()
Out[:
array([[0, 0, 0, 0, 0, 0],
[0, 1, 2, 3, 4, 5]], dtype=int32)
*boxm*.outer_faces(3)
Out:
array([[0, 0, 0, 0, 0, 0],
[0, 1, 2, 3, 4, 5]], dtype=int32)
*boxm*.outer_faces(2)
Out: array([], shape=(2, 0), dtype=int32)
*EMPTY*
*boxm*.outer_faces(1)
Out: array([], shape=(2, 0), dtype=int32)
*EMPTY*
Here *boxm *is 1 convex rectangular PARALLELEPIPED
Regards, Egor
пн, 19 июл. 2021 г. в 15:31, Konstantinos Poulios <[email protected]>:
> 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 <[email protected]>
> 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
>>
>>
