Re: [deal.II] FE_DGP/FE_FaceP and hp::DoFHandler

Can anyone help me with that matter? Best, Samuel On Monday, August 7, 2017 at 11:48:39 AM UTC+2, sml.i...@gmail.com wrote: > > I have tried the following implementation (for 2D computations only > hp_line_dof_identities has nonempty constraints): > > template > std::vector >

Re: [deal.II] FE_DGP/FE_FaceP and hp::DoFHandler

I have tried the following implementation (for 2D computations only hp_line_dof_identities has nonempty constraints): template std::vector FE_FaceP:: hp_line_dof_identities (const FiniteElement _other) const { // this element is continuous only for

Re: [deal.II] FE_DGP/FE_FaceP and hp::DoFHandler

Yes, that makes sense to me. But now if we have neighboring faces with differrent refinement levels (e.g. the face element on one side is only half as long as the one on the other side) the situation gets a little more complicated, right? Since then you cannot just copy the coefficients 1 to

[deal.II] FE_DGP/FE_FaceP and hp::DoFHandler

Dear all, As indicated in the documentation , the finite elements FE_DGP and FE_FaceP are not based on nodal interpolation (as eg. FE_DGQ and FE_FaceQ), but on projection. Apparently, this prevents these elements from being used

Re: [deal.II] Triangulation neighbor information when using GridIn::read_vtk

I must admit that I don't really remember ... I think it was GridOut, since I specifically wanted to export only the grid and not some solution data. What dou you mean with "GridOut should work"? On Wednesday, July 26, 2017 at 9:04:16 PM UTC+2, Wolfgang Bangerth wrote: > > On 07/26/2017 12:55

Re: [deal.II] Triangulation neighbor information when using GridIn::read_vtk

Exactly, I generated it using GridOut. I see the problem now. I replaced all duplicate vertices in the vtk file, now it works as it should. Thanks for the clarification! Best, Samuel On Wednesday, July 26, 2017 at 8:38:39 PM UTC+2, Wolfgang Bangerth wrote: > > On 07/26/2017 09:18 AM,

[deal.II] Triangulation neighbor information when using GridIn::read_vtk

Hello all I am trying to read in a triangulation from a vtk file (generated by the dealii library). It works, but it seems like no information about neighboring cells is generated: iterating through the cells and calling cell->neighbor_index(face) always gives -1, indicating that the face is

Re: [deal.II] Re: hp HDG in deal.II

Hello all Thanks a lot for your replies. I managed to submit a pull request (the video lecture was quite helpful). I hope I did it right, otherwise I'm happy for feedback. I will try to implement the error computation using Daniels suggestions. Now a follow-up question regarding development

[deal.II] Re: hp HDG in deal.II

Hi Daniel Thanks for your answer, your approach actually worked quite well. (How would I go about submitting a patch? I have absolutely no experience with these things) Now I have another question. I am trying to compute the L2 error over all the faces. For simplicity, I tested this starting

[deal.II] hp HDG in deal.II

Hello all I am attempting to implement a hp hybridized dg-method for solving ellipctic problems in deal.II. I am relying on steps 27 and 51 of the deal.II tutorial programs. However, when I try to group some FE_FaceQ elements in an hp::FeCollection and distribute the degrees of freedom, I get