Dear Jan, a discontinuous RT element is actually quite simple. We to convert the current FE_RaviartThomas to a base class and derive two classes, the new FE_RaviartThomas and FE_DGRaviartThomas. Their main difference would be the dpo_vector, which locates all degrees of freedom in the interior of the cell for the DG version. We will need something like this later this spring, so I might implement it. But I don't know your time horizon. If you would be willing to test the new element, that would help a lot.
As for the manifolds: the mapping class needs to implement the Piola transform, which is currently true only for the mappings of codimension 0. Implementing the basic transformation should not be too hard, once you have decided whether the result is a 2D or a 3D vector. I was not involved in the codimension 1 project, so I am not sure what to do there. Another question in this respect concerns how to handle non-differential transitions between cells. Final question: we have started implementing FE_FaceQ, which is the multiplier space for RT. Basic tests have been made, but we have not run a hybridized code yet. I hope this helps. Best, Guido Jan Brezina wrote: > Dear all, > > I'm DEAL newbie interested in using DEAL for mixed-hybrid formulation of > a porous media problem. To this end I need discontinuous RT elements, > i.e. base functions with support only on one cell. I suppose that > FE_RaviartThomas is suitable only for the continuous case. > > Further I want to make some simulations also on a 2d manifold in 3D > space so I need something like FE_RaviartThomas<2,3>, but there is a > remark that it is not implemented. > > My question is: Is it a big deal to implement these two things and what > classes are involved. > > Second question is about Lagrange multipliers I need to enclose > hybrid/discontinuous system. Is there a way how to distribute one degree > of freedom on every edge - face common to two cells (or more with > adaptivity) ? > > Thank you very much for any suggestion. > > Jan > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
