This is a read only copy of the old FEniCS QA forum. Please visit the new QA forum to ask questions

edge moments for nedelec elements in 3D

+3 votes

Dear FENICS experts,

I'm new in FENICs and I'm interfacing dolfin with PCBDDC solver in PETSc (see here ) for solving curl-curl problems in 3D with nedelec elements. For optimal convergence of the method, I need to generate a Function defined on the nedelec elements which defines the quadrature rule for the following integrals

$$\int_E {\bf w} \cdot {\bf t_e} ds $$ and $$\int_E s {\bf w} \cdot {\bf t_e} ds $$

where E is a set of element edges (namely those located on the interface between subdomains), w is a function in the Nedelec space, t_e is the tangent unit vector on each edge and s is an arc length parameters. What is the best way of implementing them? For 3D problems, the ds you can use for integration is for facets, thus not for edges.

Thanks,
Stefano

asked May 18, 2014 by szampini FEniCS Novice (510 points)
edited May 18, 2014 by szampini

1 Answer

0 votes

Integrating over mesh entities of arbitrary dimension isn't supported yet. It would be useful, and we plan to support it (see https://bitbucket.org/fenics-project/dolfin/issue/106/generalized-assembly-over-any-entity-type), but there is not timeline as yet.

answered May 20, 2014 by Garth N. Wells FEniCS Expert (35,930 points)
...