getfem-users
[Top][All Lists]

## Re: [Getfem-users] B = curl A : how to impose a jauge?

 From: Yves Renard Subject: Re: [Getfem-users] B = curl A : how to impose a jauge? Date: Mon, 3 Mar 2008 15:40:09 +0100 User-agent: KMail/1.9.5

```On Monday 03 March 2008 07:55, Riccardo Scorretti wrote:
> Hi. I'm interested in modeling in quasi-static electromagnetism, and I
> want to represent a flux density through a vector potential (i.e. B =
> curl A). To this aim, a jauge condition has to be imposed on A. One
> convenient way to impose such a jauge is to impose A=0 over a subset of
> the edges, which is a tree (that is, all nodes of the geometry are
> "touched" without forming a closed loop).
>
> The question is: how to obtain the list of all the edges of a geometry,
> so as to build a tree?

I don't completely understand what you exactly mean by a tree in this context.
I suppose you treat a 3D problem. The edges are not stored in the mesh, so
there is no natural numbering of edges. In fact, for a mesh whose base
element is a tetrahedron, the edges are caracterised by the pairs of mesh
nodes lying on the same tetrahedron. The indices of the nodes of an element
can be obtained with the method
mesh.ind_points_of_convex(ic)
Of course, with this method, edges will be meet more than once.

Yves.

--

Yves Renard (address@hidden)       tel : (33) 04.72.43.87.08
Pole de Mathematiques, INSA de Lyon          fax : (33) 04.72.43.85.29
20, rue Albert Einstein
69621 Villeurbanne Cedex, FRANCE
http://math.univ-lyon1.fr/~renard

---------

```