[Dune] Feature request
Christian Engwer
christi at uni-hd.de
Thu May 12 17:11:37 CEST 2005
Hi,
I'm using DG myself I never got hte problems you described. The
integral over the intersection I do as follows:
1) I have intersection iterator.
2) the quadrature points are defined in local coordinates of the
intersection's reference element, i.e. a line in 2D or a triangle
or quadrilateral in 3D.
3) using intersectionSelfLocal and intersectionNeighborLocal I can
map the quadrature points onto the element's or the neighbor's
reference element.
4) I evalute the function in local coordinates of the element.
Perhaps I missed something, but where do you need the twist?
There is one thing that can be important... if you assemble one part
from the element and an other part from the neighbor then you get into
trouble if you use a non-symetric quadrature rule or you need to know
the orientation of the intersection. But this problem should not
occure as the quadrature rules in Dune are all symetric at the moment.
Christian
More information about the Dune
mailing list